Bombyx BLAST Result against fruitfly

BLASTX 2.2.12 [Aug-07-2005]


Reference: Altschul, Stephen F., Thomas L. Madden, Alejandro A. Schaffer, 
Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (1997), 
"Gapped BLAST and PSI-BLAST: a new generation of protein database search
programs",  Nucleic Acids Res. 25:3389-3402.

Query= wdS20831
         (707 letters)

Database: fruitfly 
           53,049 sequences; 24,988,368 total letters

Searching..................................................done

                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

AY047538-1|AAK77270.1|  311|Drosophila melanogaster GH04846p pro...   155   7e-38
AE014297-844|AAF54312.1|  311|Drosophila melanogaster CG11968-PA...   155   7e-38
AY051510-1|AAK92934.1|  385|Drosophila melanogaster GH16429p pro...    38   0.010
AE013599-609|AAF59120.2|  385|Drosophila melanogaster CG8707-PA ...    38   0.010
AY122200-1|AAM52712.1|  261|Drosophila melanogaster LD47384p pro...    31   1.2  
AF209706-1|AAF21440.1|  261|Drosophila melanogaster Rab-related ...    31   1.2  
AE013599-3617|AAF47018.1|  261|Drosophila melanogaster CG3870-PA...    31   1.2  
D84348-1|BAA21712.1|  219|Drosophila melanogaster rab-related pr...    31   2.0  
AY071400-1|AAL49022.1|  219|Drosophila melanogaster RE48347p pro...    31   2.0  
AY071340-1|AAL48962.1|  218|Drosophila melanogaster RE37651p pro...    31   2.0  
AE014296-1412|AAF50452.1|  219|Drosophila melanogaster CG7062-PA...    31   2.0  
AE014298-1056|AAF46271.1|  218|Drosophila melanogaster CG12156-P...    30   2.7  
BT023374-1|AAY55790.1|  358|Drosophila melanogaster IP10363p pro...    30   3.5  
AE014296-330|AAF47546.1|  305|Drosophila melanogaster CG7980-PA ...    30   3.5  
D84313-1|BAA21706.1|  213|Drosophila melanogaster rab2 protein.        29   6.2  
BT003315-1|AAO25075.1|  213|Drosophila melanogaster GH01619p pro...    29   6.2  
AE013599-209|AAM70817.1|  213|Drosophila melanogaster CG3269-PA ...    29   6.2  
AB035352-1|BAA87878.1|  213|Drosophila melanogaster Drab2 protein.     29   6.2  
AE014296-740|AAF47825.3|  353|Drosophila melanogaster CG14987-PA...    29   8.2  
AE014296-739|AAS64968.1|  429|Drosophila melanogaster CG14987-PB...    29   8.2  

>AY047538-1|AAK77270.1|  311|Drosophila melanogaster GH04846p
           protein.
          Length = 311

 Score =  155 bits (375), Expect = 7e-38
 Identities = 68/85 (80%), Positives = 76/85 (89%)
 Frame = +3

Query: 255 LWDCGGQEAFMENYFASQRDNIFRNVEVLIYVFDVESREMEKDMHYYQSCLEAILQNSPE 434
           LWDCGGQE FM+ YFA QRDNIFRNVEVLIYVFDVES+E+E+D+HYYQSCLEA+LQNSPE
Sbjct: 55  LWDCGGQEGFMKQYFAQQRDNIFRNVEVLIYVFDVESQEIERDIHYYQSCLEALLQNSPE 114

Query: 435 ARIFCLVHKMDLVAEEQREEIFRAR 509
           A+IFCLVHKMDLV E  RE +F  R
Sbjct: 115 AKIFCLVHKMDLVPEGLRESVFTER 139



 Score =  109 bits (262), Expect = 4e-24
 Identities = 53/54 (98%), Positives = 53/54 (98%)
 Frame = +1

Query: 94  MKKKVLLMGKSGSGKTSMRSIIFANYIARDTRRLGATIDVEHSHVRFLGNLVLN 255
           MKKKVLLMGKSGSGKTSMRSIIFANYIARDT RLGATIDVEHSHVRFLGNLVLN
Sbjct: 1   MKKKVLLMGKSGSGKTSMRSIIFANYIARDTTRLGATIDVEHSHVRFLGNLVLN 54



 Score = 99.5 bits (237), Expect = 4e-21
 Identities = 47/65 (72%), Positives = 53/65 (81%)
 Frame = +2

Query: 512 EDLKGLSKPLECTCFRTSIWDETLYKAWSSIVYKLIPNVKALESSLKQFAAIVDADEVLL 691
           EDL  LSKP   TCFRTSIWDETLYKAWSSIV  LIPNV ALE+S+  F  +++ADEVLL
Sbjct: 141 EDLIKLSKPGNVTCFRTSIWDETLYKAWSSIVTMLIPNVAALENSVTHFGNVIEADEVLL 200

Query: 692 FERAT 706
           FE+AT
Sbjct: 201 FEKAT 205


>AE014297-844|AAF54312.1|  311|Drosophila melanogaster CG11968-PA
           protein.
          Length = 311

 Score =  155 bits (375), Expect = 7e-38
 Identities = 68/85 (80%), Positives = 76/85 (89%)
 Frame = +3

Query: 255 LWDCGGQEAFMENYFASQRDNIFRNVEVLIYVFDVESREMEKDMHYYQSCLEAILQNSPE 434
           LWDCGGQE FM+ YFA QRDNIFRNVEVLIYVFDVES+E+E+D+HYYQSCLEA+LQNSPE
Sbjct: 55  LWDCGGQEGFMKQYFAQQRDNIFRNVEVLIYVFDVESQEIERDIHYYQSCLEALLQNSPE 114

Query: 435 ARIFCLVHKMDLVAEEQREEIFRAR 509
           A+IFCLVHKMDLV E  RE +F  R
Sbjct: 115 AKIFCLVHKMDLVPEGLRESVFTER 139



 Score =  109 bits (262), Expect = 4e-24
 Identities = 53/54 (98%), Positives = 53/54 (98%)
 Frame = +1

Query: 94  MKKKVLLMGKSGSGKTSMRSIIFANYIARDTRRLGATIDVEHSHVRFLGNLVLN 255
           MKKKVLLMGKSGSGKTSMRSIIFANYIARDT RLGATIDVEHSHVRFLGNLVLN
Sbjct: 1   MKKKVLLMGKSGSGKTSMRSIIFANYIARDTTRLGATIDVEHSHVRFLGNLVLN 54



 Score = 99.5 bits (237), Expect = 4e-21
 Identities = 47/65 (72%), Positives = 53/65 (81%)
 Frame = +2

Query: 512 EDLKGLSKPLECTCFRTSIWDETLYKAWSSIVYKLIPNVKALESSLKQFAAIVDADEVLL 691
           EDL  LSKP   TCFRTSIWDETLYKAWSSIV  LIPNV ALE+S+  F  +++ADEVLL
Sbjct: 141 EDLIKLSKPGNVTCFRTSIWDETLYKAWSSIVTMLIPNVAALENSVTHFGNVIEADEVLL 200

Query: 692 FERAT 706
           FE+AT
Sbjct: 201 FEKAT 205


>AY051510-1|AAK92934.1|  385|Drosophila melanogaster GH16429p
           protein.
          Length = 385

 Score = 38.3 bits (85), Expect = 0.010
 Identities = 26/91 (28%), Positives = 49/91 (53%), Gaps = 5/91 (5%)
 Frame = +2

Query: 440 DILFSTQNGLS--GRRTKRRNI-QSKEEDLK--GLSKPLECTCFRTSIWDETLYKAWSSI 604
           ++     +G+S   +   +R+I Q   +DL   GL + +  +   TSI+D ++++A+S +
Sbjct: 153 EVFIHKVDGISDDSKMESQRDIHQRSSDDLNEAGLDQ-IHLSFHLTSIYDHSIFEAFSKV 211

Query: 605 VYKLIPNVKALESSLKQFAAIVDADEVLLFE 697
           V KLIP +  LE+ L  F      ++  LF+
Sbjct: 212 VQKLIPQLPTLENLLNIFIPNSGIEKAFLFD 242



 Score = 30.7 bits (66), Expect = 2.0
 Identities = 14/39 (35%), Positives = 24/39 (61%)
 Frame = +1

Query: 88  SKMKKKVLLMGKSGSGKTSMRSIIFANYIARDTRRLGAT 204
           S+ K ++LLMG   SGK+S++ ++F      +T  L +T
Sbjct: 37  SETKPRILLMGMRRSGKSSIQKVVFHKMSPNETLFLEST 75


>AE013599-609|AAF59120.2|  385|Drosophila melanogaster CG8707-PA
           protein.
          Length = 385

 Score = 38.3 bits (85), Expect = 0.010
 Identities = 26/91 (28%), Positives = 49/91 (53%), Gaps = 5/91 (5%)
 Frame = +2

Query: 440 DILFSTQNGLS--GRRTKRRNI-QSKEEDLK--GLSKPLECTCFRTSIWDETLYKAWSSI 604
           ++     +G+S   +   +R+I Q   +DL   GL + +  +   TSI+D ++++A+S +
Sbjct: 153 EVFIHKVDGISDDSKMESQRDIHQRSSDDLNEAGLDQ-IHLSFHLTSIYDHSIFEAFSKV 211

Query: 605 VYKLIPNVKALESSLKQFAAIVDADEVLLFE 697
           V KLIP +  LE+ L  F      ++  LF+
Sbjct: 212 VQKLIPQLPTLENLLNIFIPNSGIEKAFLFD 242



 Score = 30.7 bits (66), Expect = 2.0
 Identities = 14/39 (35%), Positives = 24/39 (61%)
 Frame = +1

Query: 88  SKMKKKVLLMGKSGSGKTSMRSIIFANYIARDTRRLGAT 204
           S+ K ++LLMG   SGK+S++ ++F      +T  L +T
Sbjct: 37  SETKPRILLMGMRRSGKSSIQKVVFHKMSPNETLFLEST 75


>AY122200-1|AAM52712.1|  261|Drosophila melanogaster LD47384p
           protein.
          Length = 261

 Score = 31.5 bits (68), Expect = 1.2
 Identities = 26/86 (30%), Positives = 40/86 (46%), Gaps = 3/86 (3%)
 Frame = +3

Query: 252 QLWDCGGQEAFMENYFASQRDNIFRNVEVLIYVFDVESREMEKDMHYYQSCLEAILQNSP 431
           Q+WD  GQE      + +     +RN    I VFD+      K     +S ++ + +N  
Sbjct: 59  QIWDTAGQER-----YRAVAPMYYRNANAAILVFDL---TQYKTFTEIKSWIQELHRNVQ 110

Query: 432 EARIFCLV-HKMDLVAEE--QREEIF 500
           +  I  LV +KMD+ A+    REE F
Sbjct: 111 DPMILTLVGNKMDMQAQRAVSREEAF 136


>AF209706-1|AAF21440.1|  261|Drosophila melanogaster Rab-related
           GTP-binding proteinChrowded protein.
          Length = 261

 Score = 31.5 bits (68), Expect = 1.2
 Identities = 26/86 (30%), Positives = 40/86 (46%), Gaps = 3/86 (3%)
 Frame = +3

Query: 252 QLWDCGGQEAFMENYFASQRDNIFRNVEVLIYVFDVESREMEKDMHYYQSCLEAILQNSP 431
           Q+WD  GQE      + +     +RN    I VFD+      K     +S ++ + +N  
Sbjct: 59  QIWDTAGQER-----YRAVAPMYYRNANAAILVFDL---TQYKTFTEIKSWIQELHRNVQ 110

Query: 432 EARIFCLV-HKMDLVAEE--QREEIF 500
           +  I  LV +KMD+ A+    REE F
Sbjct: 111 DPMILTLVGNKMDMQAQRAVSREEAF 136


>AE013599-3617|AAF47018.1|  261|Drosophila melanogaster CG3870-PA
           protein.
          Length = 261

 Score = 31.5 bits (68), Expect = 1.2
 Identities = 26/86 (30%), Positives = 40/86 (46%), Gaps = 3/86 (3%)
 Frame = +3

Query: 252 QLWDCGGQEAFMENYFASQRDNIFRNVEVLIYVFDVESREMEKDMHYYQSCLEAILQNSP 431
           Q+WD  GQE      + +     +RN    I VFD+      K     +S ++ + +N  
Sbjct: 59  QIWDTAGQER-----YRAVAPMYYRNANAAILVFDL---TQYKTFTEIKSWIQELHRNVQ 110

Query: 432 EARIFCLV-HKMDLVAEE--QREEIF 500
           +  I  LV +KMD+ A+    REE F
Sbjct: 111 DPMILTLVGNKMDMQAQRAVSREEAF 136


>D84348-1|BAA21712.1|  219|Drosophila melanogaster rab-related
           protein 3 protein.
          Length = 219

 Score = 30.7 bits (66), Expect = 2.0
 Identities = 25/95 (26%), Positives = 46/95 (48%), Gaps = 1/95 (1%)
 Frame = +3

Query: 252 QLWDCGGQEAFMENYFASQRDNIFRNVEVLIYVFDVESREMEKDMHYYQSCLEAILQNSP 431
           Q+WD  GQE      F +   + +R+   ++ V+D+  R    ++   Q  +E + + + 
Sbjct: 73  QIWDTAGQER-----FRTITQSYYRSANGVLIVYDITKRSSFSNL---QKWIEEVRRYTA 124

Query: 432 EARIFCLV-HKMDLVAEEQREEIFRARKKI*KGFP 533
              +  LV +K DL  EEQRE  F   +++ +  P
Sbjct: 125 SNVLIILVGNKCDL--EEQREVDFEEARQMCQYIP 157


>AY071400-1|AAL49022.1|  219|Drosophila melanogaster RE48347p
           protein.
          Length = 219

 Score = 30.7 bits (66), Expect = 2.0
 Identities = 25/95 (26%), Positives = 46/95 (48%), Gaps = 1/95 (1%)
 Frame = +3

Query: 252 QLWDCGGQEAFMENYFASQRDNIFRNVEVLIYVFDVESREMEKDMHYYQSCLEAILQNSP 431
           Q+WD  GQE      F +   + +R+   ++ V+D+  R    ++   Q  +E + + + 
Sbjct: 73  QIWDTAGQER-----FRTITQSYYRSANGVLIVYDITKRSSFSNL---QKWIEEVRRYTA 124

Query: 432 EARIFCLV-HKMDLVAEEQREEIFRARKKI*KGFP 533
              +  LV +K DL  EEQRE  F   +++ +  P
Sbjct: 125 SNVLIILVGNKCDL--EEQREVDFEEARQMCQYIP 157


>AY071340-1|AAL48962.1|  218|Drosophila melanogaster RE37651p
           protein.
          Length = 218

 Score = 30.7 bits (66), Expect = 2.0
 Identities = 23/83 (27%), Positives = 38/83 (45%), Gaps = 1/83 (1%)
 Frame = +3

Query: 252 QLWDCGGQEAFMENYFASQRDNIFRNVEVLIYVFDVESREMEKDMHYYQSCLEAILQNSP 431
           QLWD  GQE      F S   + +RN   ++ V+D+ +    + +  +   +EA     P
Sbjct: 62  QLWDTAGQER-----FRSITKSYYRNSVGVLLVYDISNHASFEHIPLWM--MEAQRHIEP 114

Query: 432 EARIFCLVH-KMDLVAEEQREEI 497
              +F LV  K+DL+      E+
Sbjct: 115 HRSVFALVGCKLDLINAGGHREV 137


>AE014296-1412|AAF50452.1|  219|Drosophila melanogaster CG7062-PA
           protein.
          Length = 219

 Score = 30.7 bits (66), Expect = 2.0
 Identities = 25/95 (26%), Positives = 46/95 (48%), Gaps = 1/95 (1%)
 Frame = +3

Query: 252 QLWDCGGQEAFMENYFASQRDNIFRNVEVLIYVFDVESREMEKDMHYYQSCLEAILQNSP 431
           Q+WD  GQE      F +   + +R+   ++ V+D+  R    ++   Q  +E + + + 
Sbjct: 73  QIWDTAGQER-----FRTITQSYYRSANGVLIVYDITKRSSFSNL---QKWIEEVRRYTA 124

Query: 432 EARIFCLV-HKMDLVAEEQREEIFRARKKI*KGFP 533
              +  LV +K DL  EEQRE  F   +++ +  P
Sbjct: 125 SNVLIILVGNKCDL--EEQREVDFEEARQMCQYIP 157


>AE014298-1056|AAF46271.1|  218|Drosophila melanogaster CG12156-PA
           protein.
          Length = 218

 Score = 30.3 bits (65), Expect = 2.7
 Identities = 23/83 (27%), Positives = 38/83 (45%), Gaps = 1/83 (1%)
 Frame = +3

Query: 252 QLWDCGGQEAFMENYFASQRDNIFRNVEVLIYVFDVESREMEKDMHYYQSCLEAILQNSP 431
           QLWD  GQE      F S   + +RN   ++ V+D+ +    + +  +   +EA     P
Sbjct: 62  QLWDTAGQER-----FRSITKSYYRNSVGVLLVYDISNHASFEHIPLWM--MEAQRHIEP 114

Query: 432 EARIFCLVH-KMDLVAEEQREEI 497
              +F LV  K+DL+      E+
Sbjct: 115 HRPVFALVGCKLDLINAGGHREV 137


>BT023374-1|AAY55790.1|  358|Drosophila melanogaster IP10363p
           protein.
          Length = 358

 Score = 29.9 bits (64), Expect = 3.5
 Identities = 20/82 (24%), Positives = 38/82 (46%), Gaps = 1/82 (1%)
 Frame = +3

Query: 243 FGSQLWDCGGQEAFMENYFASQRDNIFRNVEVLIYVFDVESRE-MEKDMHYYQSCLEAIL 419
           +  ++WD  GQE F     A      +RN  V++  +D+  ++ +E    +  S L    
Sbjct: 194 YSLEMWDTAGQERF--RCIAGA---YYRNASVIVVTYDMSKKDSLESAKKWLNSALNYNA 248

Query: 420 QNSPEARIFCLVHKMDLVAEEQ 485
              P   +F +  K DL+++E+
Sbjct: 249 SKRP--LVFLVGTKADLLSKEE 268


>AE014296-330|AAF47546.1|  305|Drosophila melanogaster CG7980-PA
           protein.
          Length = 305

 Score = 29.9 bits (64), Expect = 3.5
 Identities = 20/82 (24%), Positives = 38/82 (46%), Gaps = 1/82 (1%)
 Frame = +3

Query: 243 FGSQLWDCGGQEAFMENYFASQRDNIFRNVEVLIYVFDVESRE-MEKDMHYYQSCLEAIL 419
           +  ++WD  GQE F     A      +RN  V++  +D+  ++ +E    +  S L    
Sbjct: 113 YSLEMWDTAGQERF--RCIAGA---YYRNASVIVVTYDMSKKDSLESAKKWLNSALNYNA 167

Query: 420 QNSPEARIFCLVHKMDLVAEEQ 485
              P   +F +  K DL+++E+
Sbjct: 168 SKRP--LVFLVGTKADLLSKEE 187


>D84313-1|BAA21706.1|  213|Drosophila melanogaster rab2 protein.
          Length = 213

 Score = 29.1 bits (62), Expect = 6.2
 Identities = 25/84 (29%), Positives = 41/84 (48%), Gaps = 3/84 (3%)
 Frame = +3

Query: 252 QLWDCGGQEAFMENYFASQRDNIFRNVEVLIYVFDVESREMEKDMHYYQSCLEAILQNSP 431
           Q+WD  GQEA     F S   + +R     + V+D+  RE     ++  + LE   Q+S 
Sbjct: 58  QIWDTAGQEA-----FRSITRSYYRGAAGALLVYDITRRE---TFNHLTTWLEDARQHSN 109

Query: 432 EARIFCLV-HKMDLVA--EEQREE 494
              +  L+ +K DL +  E ++EE
Sbjct: 110 SNMVIMLIGNKSDLDSRREVKKEE 133


>BT003315-1|AAO25075.1|  213|Drosophila melanogaster GH01619p
           protein.
          Length = 213

 Score = 29.1 bits (62), Expect = 6.2
 Identities = 25/84 (29%), Positives = 41/84 (48%), Gaps = 3/84 (3%)
 Frame = +3

Query: 252 QLWDCGGQEAFMENYFASQRDNIFRNVEVLIYVFDVESREMEKDMHYYQSCLEAILQNSP 431
           Q+WD  GQEA     F S   + +R     + V+D+  RE     ++  + LE   Q+S 
Sbjct: 58  QIWDTAGQEA-----FRSITRSYYRGAAGALLVYDITRRE---TFNHLTTWLEDARQHSN 109

Query: 432 EARIFCLV-HKMDLVA--EEQREE 494
              +  L+ +K DL +  E ++EE
Sbjct: 110 SNMVIMLIGNKSDLDSRREVKKEE 133


>AE013599-209|AAM70817.1|  213|Drosophila melanogaster CG3269-PA
           protein.
          Length = 213

 Score = 29.1 bits (62), Expect = 6.2
 Identities = 25/84 (29%), Positives = 41/84 (48%), Gaps = 3/84 (3%)
 Frame = +3

Query: 252 QLWDCGGQEAFMENYFASQRDNIFRNVEVLIYVFDVESREMEKDMHYYQSCLEAILQNSP 431
           Q+WD  GQEA     F S   + +R     + V+D+  RE     ++  + LE   Q+S 
Sbjct: 58  QIWDTAGQEA-----FRSITRSYYRGAAGALLVYDITRRE---TFNHLTTWLEDARQHSN 109

Query: 432 EARIFCLV-HKMDLVA--EEQREE 494
              +  L+ +K DL +  E ++EE
Sbjct: 110 SNMVIMLIGNKSDLDSRREVKKEE 133


>AB035352-1|BAA87878.1|  213|Drosophila melanogaster Drab2 protein.
          Length = 213

 Score = 29.1 bits (62), Expect = 6.2
 Identities = 25/84 (29%), Positives = 41/84 (48%), Gaps = 3/84 (3%)
 Frame = +3

Query: 252 QLWDCGGQEAFMENYFASQRDNIFRNVEVLIYVFDVESREMEKDMHYYQSCLEAILQNSP 431
           Q+WD  GQEA     F S   + +R     + V+D+  RE     ++  + LE   Q+S 
Sbjct: 58  QIWDTAGQEA-----FRSITRSYYRGAAGALLVYDITRRE---TFNHLTTWLEDARQHSN 109

Query: 432 EARIFCLV-HKMDLVA--EEQREE 494
              +  L+ +K DL +  E ++EE
Sbjct: 110 SNMVIMLIGNKSDLDSRREVKKEE 133


>AE014296-740|AAF47825.3|  353|Drosophila melanogaster CG14987-PA,
           isoform A protein.
          Length = 353

 Score = 28.7 bits (61), Expect = 8.2
 Identities = 15/42 (35%), Positives = 22/42 (52%), Gaps = 1/42 (2%)
 Frame = +3

Query: 252 QLWDCGGQEAFMENYFASQRDNIFRNVEVLIY-VFDVESREM 374
           Q+ DC  +  F+ENYF  +R  +F   E   + VF +E   M
Sbjct: 100 QIRDCNLKVGFVENYFLRERQELFSVFEYRAWMVFFIEWNTM 141


>AE014296-739|AAS64968.1|  429|Drosophila melanogaster CG14987-PB
           protein.
          Length = 429

 Score = 28.7 bits (61), Expect = 8.2
 Identities = 15/42 (35%), Positives = 22/42 (52%), Gaps = 1/42 (2%)
 Frame = +3

Query: 252 QLWDCGGQEAFMENYFASQRDNIFRNVEVLIY-VFDVESREM 374
           Q+ DC  +  F+ENYF  +R  +F   E   + VF +E   M
Sbjct: 176 QIRDCNLKVGFVENYFLRERQELFSVFEYRAWMVFFIEWNTM 217


  Database: fruitfly
    Posted date:  Oct 23, 2007  1:17 PM
  Number of letters in database: 24,988,368
  Number of sequences in database:  53,049
  
Lambda     K      H
   0.318    0.134    0.401 

Gapped
Lambda     K      H
   0.279   0.0580    0.190 


Matrix: BLOSUM62
Gap Penalties: Existence: 9, Extension: 2
Number of Hits to DB: 31,283,640
Number of Sequences: 53049
Number of extensions: 661312
Number of successful extensions: 1778
Number of sequences better than 10.0: 20
Number of HSP's better than 10.0 without gapping: 1682
Number of HSP's successfully gapped in prelim test: 0
Number of HSP's that attempted gapping in prelim test: 0
Number of HSP's gapped (non-prelim): 1778
length of database: 24,988,368
effective HSP length: 83
effective length of database: 20,585,301
effective search space used: 3128965752
frameshift window, decay const: 40,  0.1
T: 12
A: 40
X1: 16 ( 7.3 bits)
X2: 37 (14.9 bits)
X3: 62 (25.0 bits)
S1: 41 (21.7 bits)