Практикум 9

1. Поиск гомологов в неаннотированном геноме

Для поиска выбрала β-тубулин, актин (белки цитоскелета) и гистон H3 (белок нуклеосомы). Все они впервые появляются у эукариот и несут базовые функции, поэтому должны быть консервативны и встречаться у подавляющего большинства эукариот. Использовала белки Saccharomyces cerevisiae, так как группа Aphelida родственна грибам.

β-тубулин

Представлена только часть выдачи, поскольку по ней очевидно наличие белка (e-value - машинный ноль, выравнялась почти вся последовательность)

Query= sp|P02557|TBB_YEAST Tubulin beta chain OS=Saccharomyces cerevisiae
(strain ATCC 204508 / S288c) OX=559292 GN=TUB2 PE=1 SV=2

Length=457
                                                                      Score     E
Sequences producing significant alignments:                          (Bits)  Value

  unplaced-665                                                        680     0.0   
  scaffold-26                                                         631     0.0   
  scaffold-57                                                         343     3e-105
  unplaced-5                                                          343     3e-105
  scaffold-423                                                        143     5e-42 
  scaffold-693                                                        28.5    7.8   
  scaffold-17                                                         28.5    8.0   


> unplaced-665
Length=20473

 Score = 680 bits (1754),  Expect = 0.0, Method: Compositional matrix adjust.
 Identities = 339/451 (75%), Positives = 386/451 (86%), Gaps = 22/451 (5%)
 Frame = -2

Query  1     MREIIHISTGQCGNQIGAAFWETICGEHGLDFNGTYHGHDDIQKERLN------------  48
             MREI+H+ TG CGN IGA FWE I  EHG+D NG Y G  D+Q ER+N            
Sbjct  7236  MREIVHL*TG*CGN*IGAKFWEVISDEHGIDPNGRYEGDSDLQLERINGEFLNVLFCA**  7057

Query  49    ----------VYFNEASSGKWVPRSINVDLEPGTIDAVRNSAIGNLFRPDNYIFGQSSAG  98
                       VYFNEAS GK+VPR++ VDLEPGT+D+VR    GNLFRPDN+IFGQS AG
Sbjct  7056  AFCITLLLIVVYFNEASGGKYVPRAVLVDLEPGTMDSVRAGPYGNLFRPDNFIFGQSGAG  6877

Query  99    NVWAKGHYTEGAELVDSVMDVIRREAEGCDSLQGFQIthslgggtgsgmgtllisKIREE  158
             N WAKGHYTEGAELVDSV+DV+R+EAEGCD LQGFQITHSLGGGTG+GMGTLLISKIREE
Sbjct  6876  NNWAKGHYTEGAELVDSVLDVVRKEAEGCDCLQGFQITHSLGGGTGAGMGTLLISKIREE  6697

Query  159   FPDRMMATFSVLPSPKTSDTVVEPYNATLSVHQLVEHSDETFCIDNEALYDICQRTLKLN  218
             +PDRMM TFSV+PSPK SDTVVEPYNATLSVHQLVE+SDETFCIDNEALYDIC RTLKL 
Sbjct  6696  YPDRMMCTFSVVPSPKVSDTVVEPYNATLSVHQLVENSDETFCIDNEALYDICFRTLKLT  6517

Query  219   QPSYGDLNNLVSSVMSGVTTSLRYPGQLNSDLRKLAVNLVPFPRLHFFMVGYAPLTAIGS  278
              P+YGDLN+LVS+VMSGVTTS+R+PGQLN+DLRKLAVN+VPFPRLHFFMVG+APLT+ GS
Sbjct  6516  TPTYGDLNHLVSAVMSGVTTSIRFPGQLNADLRKLAVNMVPFPRLHFFMVGFAPLTSRGS  6337

Query  279   QSFRSLTVPELTQQMFDAKNMMAAADPRNGRYLTVAAFFRGKVSVKEVEDEMHKVQSKNS  338
             Q +R+L+V ELT QMFDAKNMMAA+DPR+GRYL VAA FRGK+S+KEV+++M  VQ+KNS
Sbjct  6336  QQYRALSVAELTTQMFDAKNMMAASDPRHGRYLAVAAIFRGKMSMKEVDEQMLSVQTKNS  6157

Query  339   DYFVEWIPNNVQTAVCSVAPQGLDMAATFIANSTSIQELFKRVGDQFSAMFKRKAFLHWY  398
              YFVEWIPNNV+TAVC + P+GL M+ATFI NST+IQELFKR+ DQFS MFKRKAFLHWY
Sbjct  6156  SYFVEWIPNNVKTAVCDIPPKGLKMSATFIGNSTAIQELFKRISDQFSVMFKRKAFLHWY  5977

Query  399   TSEGMDELEFSEAESNMNDLVSEYQQYQEAT  429
             T EGMDE+EF+EAESNMNDLVSEYQQYQ+AT
Sbjct  5976  TGEGMDEMEFTEAESNMNDLVSEYQQYQDAT  5884


> scaffold-26
Length=483268

 Score = 631 bits (1628),  Expect = 0.0, Method: Compositional matrix adjust.
 Identities = 308/381 (81%), Positives = 351/381 (92%), Gaps = 0/381 (0%)
 Frame = -3

Query  49      VYFNEASSGKWVPRSINVDLEPGTIDAVRNSAIGNLFRPDNYIFGQSSAGNVWAKGHYTE  108
               VY+NEAS GK+VPR++ VDLEPGT+D+VR    GNLFRPDN+IFGQS AGN WAKGHYTE
Sbjct  109793  VYYNEASGGKYVPRAVLVDLEPGTMDSVRAGPYGNLFRPDNFIFGQSGAGNNWAKGHYTE  109614

Query  109     GAELVDSVMDVIRREAEGCDSLQGFQIthslgggtgsgmgtllisKIREEFPDRMMATFS  168
               GAELVDSV+DV+R+EAEGCD LQGFQITHSLGGGTG+GMGTLLISKIREE+PDRMM TFS
Sbjct  109613  GAELVDSVLDVVRKEAEGCDCLQGFQITHSLGGGTGAGMGTLLISKIREEYPDRMMCTFS  109434

Query  169     VLPSPKTSDTVVEPYNATLSVHQLVEHSDETFCIDNEALYDICQRTLKLNQPSYGDLNNL  228
               V+PSPK SDTVVEPYNATLSVHQLVE+SDETFCIDNEALYDIC RTLKL  P+YGDLN+L
Sbjct  109433  VVPSPKVSDTVVEPYNATLSVHQLVENSDETFCIDNEALYDICFRTLKLTTPTYGDLNHL  109254

Query  229     VSSVMSGVTTSLRYPGQLNSDLRKLAVNLVPFPRLHFFMVGYAPLTAIGSQSFRSLTVPE  288
               VS+VMSGVTTS+R+PGQLN+DLRKLAVN+VPFPRLHFFMVG+APLT+ GSQ +R+L+V E
Sbjct  109253  VSAVMSGVTTSIRFPGQLNADLRKLAVNMVPFPRLHFFMVGFAPLTSRGSQQYRALSVAE  109074

Query  289     LTQQMFDAKNMMAAADPRNGRYLTVAAFFRGKVSVKEVEDEMHKVQSKNSDYFVEWIPNN  348
               LT QMFDAKNMMAA+DPR+GRYL VAA FRGK+S+KEV+++M  VQ+KNS YFVEWIPNN
Sbjct  109073  LTTQMFDAKNMMAASDPRHGRYLAVAAIFRGKMSMKEVDEQMLSVQTKNSSYFVEWIPNN  108894

Query  349     VQTAVCSVAPQGLDMAATFIANSTSIQELFKRVGDQFSAMFKRKAFLHWYTSEGMDELEF  408
               V+TAVC + P+GL M+ATFI NST+IQELFKR+ DQFS MFKRKAFLHWYT EGMDE+EF
Sbjct  108893  VKTAVCDIPPKGLKMSATFIGNSTAIQELFKRISDQFSVMFKRKAFLHWYTGEGMDEMEF  108714

Query  409     SEAESNMNDLVSEYQQYQEAT  429
               +EAESNMNDLVSEYQQYQ+AT
Sbjct  108713  TEAESNMNDLVSEYQQYQDAT  108651

Актин

Ситуация аналогична предыдущей. Ген актина в геноме есть

Query= sp|P60010|ACT_YEAST Actin OS=Saccharomyces cerevisiae (strain ATCC
204508 / S288c) OX=559292 GN=ACT1 PE=1 SV=1

Length=375
                                                                      Score     E
Sequences producing significant alignments:                          (Bits)  Value

  scaffold-17                                                         721     0.0   
  scaffold-444                                                        721     0.0   
  scaffold-514                                                        717     0.0   
  scaffold-693                                                        376     6e-118
  scaffold-170                                                        374     2e-117
  scaffold-96                                                         356     4e-111
  scaffold-423                                                        356     4e-111
  unplaced-1003                                                       215     2e-70 
  unplaced-948                                                        215     2e-70 
  unplaced-987                                                        234     1e-68 
  unplaced-1039                                                       191     4e-61 
  unplaced-925                                                        191     4e-61 
  scaffold-277                                                        188     6e-53 
  scaffold-199                                                        180     4e-50 
  scaffold-105                                                        62.8    1e-10 
  scaffold-115                                                        45.4    4e-05 
  scaffold-361                                                        45.1    5e-05 
  scaffold-157                                                        30.4    1.7   
  scaffold-287                                                        29.3    4.1   


> scaffold-17
Length=2125590

 Score = 721 bits (1861),  Expect = 0.0, Method: Compositional matrix adjust.
 Identities = 333/375 (89%), Positives = 363/375 (97%), Gaps = 0/375 (0%)
 Frame = -3

Query  1       MDSEVAALVIDNGSGMCKAGFAGDDAPRAVFPSIVGRPRHQGIMVGMGQKDSYVGDEAQS  60
               M+ EVAALVIDNGSGMCKAGFAGDDAPRAVFPSIVGRPRHQG+MVGMGQKDSYVGDEAQS
Sbjct  893374  MEDEVAALVIDNGSGMCKAGFAGDDAPRAVFPSIVGRPRHQGVMVGMGQKDSYVGDEAQS  893195

Query  61      KRGILTLRYPIEHGIVTNWDDMEKIWHHTFYNELRVAPEEHPVLLTEAPMNPKSNREKMT  120
               KRGILTL+YPIEHGIVTNWDDMEKIWHHTFYNELRVAPEEHPVL TEAP+NPK+NREKMT
Sbjct  893194  KRGILTLKYPIEHGIVTNWDDMEKIWHHTFYNELRVAPEEHPVLCTEAPLNPKANREKMT  893015

Query  121     QIMFETFNVPAFYVSIQAVLSLYSSGRTTGIVLDSGDGVTHVVPIYAGFSLPHAILRIDL  180
               QIMFETFNVPAFYV+IQAVLSLY+SGRTTGIVLDSGDGV+H VPIY G++LPHAILR+DL
Sbjct  893014  QIMFETFNVPAFYVAIQAVLSLYASGRTTGIVLDSGDGVSHTVPIYEGYALPHAILRLDL  892835

Query  181     AGRDLTDYLMKILSERGYSFSTTAEREIVRDIKEKLCYVALDFEQEMQTAAQSSSIEKSY  240
               AGRDLTDYLMKIL+ERGYSF+TTAEREIVRDIKEKLCYVALDFE EMQTA+ SS++EKSY
Sbjct  892834  AGRDLTDYLMKILTERGYSFTTTAEREIVRDIKEKLCYVALDFE*EMQTASTSSALEKSY  892655

Query  241     ELPDGQVITIGNERFRAPEALFHPSVLGLESAGIDQTTYNSIMKCDVDVRKELYGNIVMS  300
               ELPDGQVITIGNERFR PEALF PS LG+E+AG+ +TTYNSIMKCDVD+RK+LY NIV+S
Sbjct  892654  ELPDGQVITIGNERFRCPEALFQPSFLGMEAAGVHETTYNSIMKCDVDIRKDLYANIVLS  892475

Query  301     GGTTMFPGIAERMQKEITALAPSSMKVKIIAPPERKYSVWIGGSILASLTTFQQMWISKQ  360
               GGTTM+PGIA+R+QKEITALAPSSMK+K++APPERKYSVWIGGSILASL+TFQQMWISKQ
Sbjct  892474  GGTTMYPGIADRVQKEITALAPSSMKIKVVAPPERKYSVWIGGSILASLSTFQQMWISKQ  892295

Query  361     EYDESGPSIVHHKCF  375
               EYDESGPSIVH KCF
Sbjct  892294  EYDESGPSIVHRKCF  892250


 Score = 67.0 bits (162),  Expect = 6e-12, Method: Compositional matrix adjust.
 Identities = 56/239 (23%), Positives = 106/239 (44%), Gaps = 24/239 (10%)
 Frame = +1

Query  3        SEVAALVIDNGSGMCKAGFAGDDAPRAVFPSIVGRPRHQGIMVGMGQKDSYVGDEAQSKR  62
                S+ + +VID GS   +AG+A  D P+ VF ++V +      + G    D Y G  +    
Sbjct  1426510  SQESPIVIDFGSWQLRAGYADTD-PQLVFDNLVSKTS*GIHIGGGD*CD-YSGARSA---  1426674

Query  63       GILTLRYPIEHGIVTNWDDMEKIWHHTFYNELRVAPEE----HPVLLTEAPMNPKSNREK  118
                      +     ++ +  ++E    H F N L+   E+     P++ TE   N    R  
Sbjct  1426675  ------FDFS*SLIVHPPNLELTLDHVFANLLKAQREKTLA*RPIVFTETLCNLN*CRRI  1426836

Query  119      MTQIMFETFNVPAFYVSIQAVLSLYSSGRTTG---------IVLDSGDGVTHVVPIYAGF  169
                ++++MFE ++V      + A+ S Y +G  +G         ++L SG   +H++P+  G 
Sbjct  1426837  VSELMFECYSVGRLAYGVDALFSFYHNGGFSGDDKKKKRNALILSSGHNASHIIPVIDGG  1427016

Query  170      SLPHAILRIDLAGRDLTDYLMKILSERGYSFSTTAEREIVRDIKEKLCYVALDFEQEMQ  228
                       R+D+ GR  T+++ K L  +    +          ++ + C VA D++ +M 
Sbjct  1427017  VDYSHCKRLDIGGRSSTEFMQKHLLVKYPELNQKLSLMTTD*LRVEFCRVASDYKGDMN  1427193


 Score = 58.2 bits (139),  Expect = 4e-09, Method: Compositional matrix adjust.
 Identities = 35/133 (26%), Positives = 71/133 (53%), Gaps = 7/133 (5%)
 Frame = +1

Query  248      ITIGNERFRAPEALFHPSVLGLESAGIDQ--TTYNSIMKCDVDVRKELYGNIVMSGGTTM  305
                I + ++R++ PE LF P ++G ESAG+ +  +T   + + D+  + ++  N+ ++GG TM
Sbjct  1428313  IKLSSDRYKVPEILFQPHLMGSESAGLSEMCSTLAKVYESDLS-KLQMMSNVFVTGGNTM  1428489

Query  306      FPGIAERMQKEITALAPSSMKVKIIAPPERKYSVWIGGS--ILASLTTFQ-QMWISKQEY  362
                +PG   R+  ++       + VK+    +     W G S   L  L+    +++I+++ Y
Sbjct  1428490  YPGFD*RLFNDLR*NVLCDVPVKVYRAQDSILDAWRGASKWALDHLSNKN*KVFITRK*Y  1428669

Query  363      DESGPSIVH-HKC  374
                +E G   +  H+C
Sbjct  1428670  EECGADYLQEHQC  1428708


> scaffold-444
Length=851481

 Score = 721 bits (1861),  Expect = 0.0, Method: Compositional matrix adjust.
 Identities = 333/375 (89%), Positives = 363/375 (97%), Gaps = 0/375 (0%)
 Frame = -2

Query  1       MDSEVAALVIDNGSGMCKAGFAGDDAPRAVFPSIVGRPRHQGIMVGMGQKDSYVGDEAQS  60
               M+ EVAALVIDNGSGMCKAGFAGDDAPRAVFPSIVGRPRHQG+MVGMGQKDSYVGDEAQS
Sbjct  165923  MEDEVAALVIDNGSGMCKAGFAGDDAPRAVFPSIVGRPRHQGVMVGMGQKDSYVGDEAQS  165744

Query  61      KRGILTLRYPIEHGIVTNWDDMEKIWHHTFYNELRVAPEEHPVLLTEAPMNPKSNREKMT  120
               KRGILTL+YPIEHGIVTNWDDMEKIWHHTFYNELRVAPEEHPVL TEAP+NPK+NREKMT
Sbjct  165743  KRGILTLKYPIEHGIVTNWDDMEKIWHHTFYNELRVAPEEHPVLCTEAPLNPKANREKMT  165564

Query  121     QIMFETFNVPAFYVSIQAVLSLYSSGRTTGIVLDSGDGVTHVVPIYAGFSLPHAILRIDL  180
               QIMFETFNVPAFYV+IQAVLSLY+SGRTTGIVLDSGDGV+H VPIY G++LPHAILR+DL
Sbjct  165563  QIMFETFNVPAFYVAIQAVLSLYASGRTTGIVLDSGDGVSHTVPIYEGYALPHAILRLDL  165384

Query  181     AGRDLTDYLMKILSERGYSFSTTAEREIVRDIKEKLCYVALDFEQEMQTAAQSSSIEKSY  240
               AGRDLTDYLMKIL+ERGYSF+TTAEREIVRDIKEKLCYVALDFE EMQTA+ SS++EKSY
Sbjct  165383  AGRDLTDYLMKILTERGYSFTTTAEREIVRDIKEKLCYVALDFE*EMQTASTSSALEKSY  165204

Query  241     ELPDGQVITIGNERFRAPEALFHPSVLGLESAGIDQTTYNSIMKCDVDVRKELYGNIVMS  300
               ELPDGQVITIGNERFR PEALF PS LG+E+AG+ +TTYNSIMKCDVD+RK+LY NIV+S
Sbjct  165203  ELPDGQVITIGNERFRCPEALFQPSFLGMEAAGVHETTYNSIMKCDVDIRKDLYANIVLS  165024

Query  301     GGTTMFPGIAERMQKEITALAPSSMKVKIIAPPERKYSVWIGGSILASLTTFQQMWISKQ  360
               GGTTM+PGIA+R+QKEITALAPSSMK+K++APPERKYSVWIGGSILASL+TFQQMWISKQ
Sbjct  165023  GGTTMYPGIADRVQKEITALAPSSMKIKVVAPPERKYSVWIGGSILASLSTFQQMWISKQ  164844

Query  361     EYDESGPSIVHHKCF  375
               EYDESGPSIVH KCF
Sbjct  164843  EYDESGPSIVHRKCF  164799


> scaffold-514
Length=59281

 Score = 717 bits (1851),  Expect = 0.0, Method: Compositional matrix adjust.
 Identities = 332/375 (89%), Positives = 362/375 (97%), Gaps = 0/375 (0%)
 Frame = -2

Query  1      MDSEVAALVIDNGSGMCKAGFAGDDAPRAVFPSIVGRPRHQGIMVGMGQKDSYVGDEAQS  60
              M+ EVAALVIDNGSGMCKAGFAGDDAPRAVFPSIVGRPRHQG+MVGMGQKDSYVGDEAQS
Sbjct  10368  MEDEVAALVIDNGSGMCKAGFAGDDAPRAVFPSIVGRPRHQGVMVGMGQKDSYVGDEAQS  10189

Query  61     KRGILTLRYPIEHGIVTNWDDMEKIWHHTFYNELRVAPEEHPVLLTEAPMNPKSNREKMT  120
              KRGILTL+YPIEHGIVTNWDDMEKIWHHTFYNELRVAPEEHPVL TEAP+NPK+NREKMT
Sbjct  10188  KRGILTLKYPIEHGIVTNWDDMEKIWHHTFYNELRVAPEEHPVLCTEAPLNPKANREKMT  10009

Query  121    QIMFETFNVPAFYVSIQAVLSLYSSGRTTGIVLDSGDGVTHVVPIYAGFSLPHAILRIDL  180
              QIMFETFNVPAFYV+I AVLSLY+SGRTTGIVLDSGDGV+H VPIY G++LPHAILR+DL
Sbjct  10008  QIMFETFNVPAFYVAI*AVLSLYASGRTTGIVLDSGDGVSHTVPIYEGYALPHAILRLDL  9829

Query  181    AGRDLTDYLMKILSERGYSFSTTAEREIVRDIKEKLCYVALDFEQEMQTAAQSSSIEKSY  240
              AGRDLTDYLMKIL+ERGYSF+TTAEREIVRDIKEKLCYVALDFE EMQTA+ SS++EKSY
Sbjct  9828   AGRDLTDYLMKILTERGYSFTTTAEREIVRDIKEKLCYVALDFE*EMQTASTSSALEKSY  9649

Query  241    ELPDGQVITIGNERFRAPEALFHPSVLGLESAGIDQTTYNSIMKCDVDVRKELYGNIVMS  300
              ELPDGQVITIGNERFR PEALF PS LG+E+AG+ +TTYNSIMKCDVD+RK+LY NIV+S
Sbjct  9648   ELPDGQVITIGNERFRCPEALFQPSFLGMEAAGVHETTYNSIMKCDVDIRKDLYANIVLS  9469

Query  301    GGTTMFPGIAERMQKEITALAPSSMKVKIIAPPERKYSVWIGGSILASLTTFQQMWISKQ  360
              GGTTM+PGIA+R+QKEITALAPSSMK+K++APPERKYSVWIGGSILASL+TFQQMWISKQ
Sbjct  9468   GGTTMYPGIADRVQKEITALAPSSMKIKVVAPPERKYSVWIGGSILASLSTFQQMWISKQ  9289

Query  361    EYDESGPSIVHHKCF  375
              EYDESGPSIVH KCF
Sbjct  9288   EYDESGPSIVHRKCF  9244

Гистон Н3

В этом случае выравниваний с нулевыми e-value нет, но есть несколько убедительных находок с выравниванием полной последовательности белка и низкими e-value

Query= sp|P61830|H3_YEAST Histone H3 OS=Saccharomyces cerevisiae (strain
ATCC 204508 / S288c) OX=559292 GN=HHT1 PE=1 SV=2

Length=136
                                                                      Score     E
Sequences producing significant alignments:                          (Bits)  Value

  scaffold-126                                                        232     9e-72
  scaffold-104                                                        232     9e-72
  unplaced-984                                                        174     1e-57
  scaffold-153                                                        153     1e-48
  scaffold-70                                                         58.5    2e-11
  unplaced-368                                                        57.0    1e-10
  unplaced-13                                                         29.3    0.36 
  scaffold-168                                                        27.7    1.8  


> scaffold-126
Length=92370

 Score = 232 bits (591),  Expect = 9e-72, Method: Compositional matrix adjust.
 Identities = 119/136 (88%), Positives = 127/136 (93%), Gaps = 1/136 (1%)
 Frame = -3

Query  1      MARTKQTARKSTGGKAPRKQLASKAARKSAPSTGGVKKPHRYKPGTVALREIRRFQKSTE  60
              MARTK TARKSTGGKAPRKQLA+KA RKS P+TGGVKKPHRY+PGTVALREIRR+QKSTE
Sbjct  85096  MARTK*TARKSTGGKAPRKQLATKA-RKSTPATGGVKKPHRYRPGTVALREIRRYQKSTE  84920

Query  61     LLIRKLPFQRLVREIAQDFKTDLRFQSSAIGALQESVEAYLVSLFEDTNLAAIHAKRVTI  120
              LLIRKLPFQRLVREIAQDFKTDLRF SSAIGAL E+ EAYLVSLFEDTNL AIHAKRVTI
Sbjct  84919  LLIRKLPFQRLVREIAQDFKTDLRF*SSAIGAL*EAAEAYLVSLFEDTNLCAIHAKRVTI  84740

Query  121    QKKDIKLARRLRGERS  136
                KDI+LARR+RGER+
Sbjct  84739  FPKDIQLARRIRGERA  84692


 Score = 57.0 bits (136),  Expect = 2e-10, Method: Compositional matrix adjust.
 Identities = 32/37 (86%), Positives = 34/37 (92%), Gaps = 1/37 (3%)
 Frame = +2

Query  1      MARTKQTARKSTGGKAPRKQLASKAARKSAPSTGGVK  37
              MARTK TARKSTGGKAPRKQLA+K ARKS P+TGGVK
Sbjct  92261  MARTK*TARKSTGGKAPRKQLATK-ARKSTPATGGVK  92368


> scaffold-104
Length=736462

 Score = 232 bits (591),  Expect = 9e-72, Method: Compositional matrix adjust.
 Identities = 119/136 (88%), Positives = 127/136 (93%), Gaps = 1/136 (1%)
 Frame = +3

Query  1       MARTKQTARKSTGGKAPRKQLASKAARKSAPSTGGVKKPHRYKPGTVALREIRRFQKSTE  60
               MARTK TARKSTGGKAPRKQLA+KA RKS P+TGGVKKPHRY+PGTVALREIRR+QKSTE
Sbjct  185496  MARTK*TARKSTGGKAPRKQLATKA-RKSTPATGGVKKPHRYRPGTVALREIRRYQKSTE  185672

Query  61      LLIRKLPFQRLVREIAQDFKTDLRFQSSAIGALQESVEAYLVSLFEDTNLAAIHAKRVTI  120
               LLIRKLPFQRLVREIAQDFKTDLRF SSAIGAL E+ EAYLVSLFEDTNL AIHAKRVTI
Sbjct  185673  LLIRKLPFQRLVREIAQDFKTDLRF*SSAIGAL*EAAEAYLVSLFEDTNLCAIHAKRVTI  185852

Query  121     QKKDIKLARRLRGERS  136
                 KDI+LARR+RGER+
Sbjct  185853  FPKDIQLARRIRGERA  185900


 Score = 232 bits (591),  Expect = 9e-72, Method: Compositional matrix adjust.
 Identities = 119/136 (88%), Positives = 127/136 (93%), Gaps = 1/136 (1%)
 Frame = -2

Query  1       MARTKQTARKSTGGKAPRKQLASKAARKSAPSTGGVKKPHRYKPGTVALREIRRFQKSTE  60
               MARTK TARKSTGGKAPRKQLA+KA RKS P+TGGVKKPHRY+PGTVALREIRR+QKSTE
Sbjct  178362  MARTK*TARKSTGGKAPRKQLATKA-RKSTPATGGVKKPHRYRPGTVALREIRRYQKSTE  178186

Query  61      LLIRKLPFQRLVREIAQDFKTDLRFQSSAIGALQESVEAYLVSLFEDTNLAAIHAKRVTI  120
               LLIRKLPFQRLVREIAQDFKTDLRF SSAIGAL E+ EAYLVSLFEDTNL AIHAKRVTI
Sbjct  178185  LLIRKLPFQRLVREIAQDFKTDLRF*SSAIGAL*EAAEAYLVSLFEDTNLCAIHAKRVTI  178006

Query  121     QKKDIKLARRLRGERS  136
                 KDI+LARR+RGER+
Sbjct  178005  FPKDIQLARRIRGERA  177958


> unplaced-984
Length=389

 Score = 174 bits (442),  Expect = 1e-57, Method: Compositional matrix adjust.
 Identities = 86/98 (88%), Positives = 92/98 (94%), Gaps = 0/98 (0%)
 Frame = +2

Query  39   PHRYKPGTVALREIRRFQKSTELLIRKLPFQRLVREIAQDFKTDLRFQSSAIGALQESVE  98
            PHRY+PGTVALREIRR+QKSTELLIRKLPFQRLVREIAQDFKTDLRF SSAIGAL E+ E
Sbjct  2    PHRYRPGTVALREIRRYQKSTELLIRKLPFQRLVREIAQDFKTDLRF*SSAIGAL*EAAE  181

Query  99   AYLVSLFEDTNLAAIHAKRVTIQKKDIKLARRLRGERS  136
            AYLVSLFEDTNL AIHAKRVTI  KDI+LARR+RGER+
Sbjct  182  AYLVSLFEDTNLCAIHAKRVTIFPKDIQLARRIRGERA  295


> scaffold-153
Length=574

 Score = 153 bits (387),  Expect = 1e-48, Method: Compositional matrix adjust.
 Identities = 78/89 (88%), Positives = 83/89 (93%), Gaps = 0/89 (0%)
 Frame = -2

Query  48   ALREIRRFQKSTELLIRKLPFQRLVREIAQDFKTDLRFQSSAIGALQESVEAYLVSLFED  107
            ALREIRR+QKSTELLIRKLPFQRLVREIAQDFKTDLRF SSAIGAL E+ EAYLVSLFED
Sbjct  573  ALREIRRYQKSTELLIRKLPFQRLVREIAQDFKTDLRF*SSAIGAL*EAAEAYLVSLFED  394

Query  108  TNLAAIHAKRVTIQKKDIKLARRLRGERS  136
            TNL AIHAKRVTI  KDI+LARR+RGER+
Sbjct  393  TNLCAIHAKRVTIFPKDIQLARRIRGERA  307