# ../bin/fastx36 -q -k 1000 ../seq/mgstm1.esq ../seq/xurt8c.aa
FASTX compares a DNA sequence to a protein sequence data bank
 version 36.3.8i May, 2023
Please cite:
 Pearson et al, Genomics (1997) 46:24-36

Query: ../seq/mgstm1.esq
  1>>>mgstm1e  - 1087 aa
Library: ../seq/xurt8c.aa
      222 residues in     1 sequences

Statistics: (shuffled [115]) MLE statistics: Lambda= 0.1381;  K=0.0009471
 statistics sampled from 1 (1) to 115 sequences
Algorithm: FASTX (3.8 June 2014) [optimized]
Parameters: BL50 matrix (15:-5), open/ext: -12/-2, shift: -20
 ktup: 2, E-join: 0.5 (0.5), E-opt: 0.1 (0.5), width:  16
 Scan time:  0.160

The best scores are:                                      opt bits E(1)
XURT8C | 40001 | glutathione transferase (EC 2 ( 222) [f]  153 45.5 4.8e-09
XURT8C | 40001 | glutathione transferase (EC 2 ( 222) [r]   29 20.8    0.12

>>XURT8C | 40001 | glutathione transferase (EC 2.5.1.18)  (222 aa)
 initn: 135 init1: 111 opt: 153  Z-score: 194.9  bits: 45.5 E(1): 4.8e-09
Smith-Waterman score: 153; 27.0% identity (56.0% similar) in 159 aa overlap (161-620:50-207)

        170        200       230       260       290       320     
mgstm1 LDFPNLPYL-IDGSHKITQSNAILRYLARKHHLDGETEEERIRADIVENQVMDTRMQLIM
       : : ..: . :::   .::. ::: ::: :..: :.  .::.: :.  . ..:  :..: 
XURT8C LLFGQVPLVEIDG-MLLTQTRAILSYLAAKYNLYGKDLKERVRIDMYADGTQDLMMMIIG
      50        60         70        80        90       100        

         350       380       410          440          470         
mgstm1 LCYNPDFEKQKPEFLKTIPEKMKLYSEF--LGK\GHGFA---GDKVTYVDFLAYDILDQY
         ..   ::..   : .   : . .  :  . :  :: :   :......:.   . . . 
XURT8C APFKAPQEKEESLALAVKRAKNRYFPVFEKILK-DHGEAFLVGNQLSWADIQLLEAILMV
      110       120       130       140        150       160       

     500       530       560       590         
mgstm1 RMFEPKCLDAFPNLRDFLARFEGLKKISAYMKSSRYIATP
       .      :. :: :. : .:. ..  :. ... .     :
XURT8C EEVSAPVLSDFPLLQAFKTRISNIPTIKKFLQPGSQRKPP
       170       180       190       200       

>>XURT8C | 40001 | glutathione transferase (EC 2.5.1.18)  (222 aa)
rev-comp initn:  29 init1:  29 opt:  29  Z-score: 61.3  bits: 20.8 E(1): 0.12
Smith-Waterman score: 45; 28.6% identity (71.4% similar) in 21 aa overlap (267-205:3-23)

          250       220     
mgstm- VSPSRWCFRAR*RRMALLWVI
       :.:. . :..: :  .. :..
XURT8C VKPKLYYFQGRGRMEVIRWLL
             10        20   

>--
rev-comp initn:  45 init1:  26 opt:  33  Z-score: 65.7  bits: 21.6 E(1): 0.073
Smith-Waterman score: 33; 25.0% identity (75.0% similar) in 20 aa overlap (688-629:63-82)

       680       650       
mgstm- VLLVSARVARAYLLLQWAIF
       .::...:.  .::  .. ..
XURT8C MLLTQTRAILSYLAAKYNLY
             70        80  

>--
rev-comp initn:  40 init1:  26 opt:  31  Z-score: 63.5  bits: 21.2 E(1): 0.095
Smith-Waterman score: 31; 53.3% identity (60.0% similar) in 15 aa overlap (712-677:44-58)

     710          680 
mgstm- LQDPQC---GQVLLV
       ::  .:   ::: ::
XURT8C LQKDDCLLFGQVPLV
            50        

>--
rev-comp initn:  49 init1:  29 opt:  31  Z-score: 63.5  bits: 21.2 E(1): 0.095
Smith-Waterman score: 48; 31.6% identity (73.7% similar) in 19 aa overlap (415-359:89-107)

      410       380       
mgstm- RVELHFLRDGLQELWLLLL
       ::.. .  :: :.: ....
XURT8C RVRIDMYADGTQDLMMMII
       90       100       

>--
rev-comp initn:  46 init1:  29 opt:  31  Z-score: 63.5  bits: 21.2 E(1): 0.095
Smith-Waterman score: 37; 40.0% identity (53.3% similar) in 15 aa overlap (740-696:199-213)

           720        
mgstm- RPGCCPHPGPAGSPV
       .::   .: : :  :
XURT8C QPGSQRKPPPDGHYV
      200       210   



1087 residues in 1 query   sequences
222 residues in 1 library sequences
 Tcomplib [36.3.8i May, 2023] (64 proc in memory [0G])
 start: Tue May 23 09:59:24 2023 done: Tue May 23 09:59:24 2023
 Total Scan time:  0.160 Total Display time:  0.020

Function used was FASTX [36.3.8i May, 2023]