FASTA searches a protein or DNA sequence data bank 36.3.4 Apr, 2011
Please cite:
W.R. Pearson & D.J. Lipman PNAS (1988) 85:2444-2448

Query: pF1KB7631, 320 aa
  1>>>pF1KB7631 320 - 320 aa - 320 aa
Library: human.CCDS.faa
  18511270 residues in 32554 sequences

Statistics:  Expectation_n fit: rho(ln(x))= 7.4774+/-0.000773; mu= 8.6106+/- 0.047
 mean_var=177.3390+/-35.361, 0's: 0 Z-trim(115.9): 28  B-trim: 723 in 1/52
 Lambda= 0.096310
 statistics sampled from 16418 (16446) to 16418 sequences
Algorithm: FASTA (3.7 Nov 2010) [optimized]
Parameters: BL50 matrix (15:-5), open/ext: -10/-2
 ktup: 2, E-join: 1 (0.818), E-opt: 0.2 (0.505), width:  16
 Scan time:  2.990

The best scores are:                                      opt bits E(32554)
CCDS7826.1 MYOD1 gene_id:4654|Hs108|chr11          ( 320) 2225 320.4 1.2e-87
CCDS9020.1 MYF5 gene_id:4617|Hs108|chr12           ( 255)  564 89.5 2.9e-18
CCDS1433.1 MYOG gene_id:4656|Hs108|chr1            ( 224)  422 69.8 2.3e-12
CCDS9019.1 MYF6 gene_id:4618|Hs108|chr12           ( 242)  420 69.5   3e-12


>>CCDS7826.1 MYOD1 gene_id:4654|Hs108|chr11               (320 aa)
 initn: 2225 init1: 2225 opt: 2225  Z-score: 1687.3  bits: 320.4 E(32554): 1.2e-87
Smith-Waterman score: 2225; 100.0% identity (100.0% similar) in 320 aa overlap (1-320:1-320)

               10        20        30        40        50        60
pF1KB7 MELLSPPLRDVDLTAPDGSLCSFATTDDFYDDPCFDSPDLRFFEDLDPRLMHVGALLKPE
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS78 MELLSPPLRDVDLTAPDGSLCSFATTDDFYDDPCFDSPDLRFFEDLDPRLMHVGALLKPE
               10        20        30        40        50        60

               70        80        90       100       110       120
pF1KB7 EHSHFPAAVHPAPGAREDEHVRAPSGHHQAGRCLLWACKACKRKTTNADRRKAATMRERR
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS78 EHSHFPAAVHPAPGAREDEHVRAPSGHHQAGRCLLWACKACKRKTTNADRRKAATMRERR
               70        80        90       100       110       120

              130       140       150       160       170       180
pF1KB7 RLSKVNEAFETLKRCTSSNPNQRLPKVEILRNAIRYIEGLQALLRDQDAAPPGAAAAFYA
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS78 RLSKVNEAFETLKRCTSSNPNQRLPKVEILRNAIRYIEGLQALLRDQDAAPPGAAAAFYA
              130       140       150       160       170       180

              190       200       210       220       230       240
pF1KB7 PGPLPPGRGGEHYSGDSDASSPRSNCSDGMMDYSGPPSGARRRNCYEGAYYNEAPSEPRP
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS78 PGPLPPGRGGEHYSGDSDASSPRSNCSDGMMDYSGPPSGARRRNCYEGAYYNEAPSEPRP
              190       200       210       220       230       240

              250       260       270       280       290       300
pF1KB7 GKSAAVSSLDCLSSIVERISTESPAAPALLLADVPSESPPRRQEAAAPSEGESSGDPTQS
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS78 GKSAAVSSLDCLSSIVERISTESPAAPALLLADVPSESPPRRQEAAAPSEGESSGDPTQS
              250       260       270       280       290       300

              310       320
pF1KB7 PDAAPQCPAGANPNPIYQVL
       ::::::::::::::::::::
CCDS78 PDAAPQCPAGANPNPIYQVL
              310       320

>>CCDS9020.1 MYF5 gene_id:4617|Hs108|chr12                (255 aa)
 initn: 758 init1: 550 opt: 564  Z-score: 441.3  bits: 89.5 E(32554): 2.9e-18
Smith-Waterman score: 767; 48.7% identity (69.2% similar) in 279 aa overlap (17-294:5-247)

               10        20        30        40        50        60
pF1KB7 MELLSPPLRDVDLTAPDGSLCSFATTDDFYDDPCFDSPDLRFFEDLDPRLMHVGALLKPE
                       ::  :.:. .. :::  :. ::. .: ... ::.   ::     
CCDS90             MDVMDG--CQFSPSEYFYDGSCIPSPEGEFGDEFVPRVAAFGA-----
                             10        20        30        40      

               70        80        90       100       110       120
pF1KB7 EHSHFPAAVHPAPGAREDEHVRAPSGHHQAGRCLLWACKACKRKTTNADRRKAATMRERR
        :.   : ..   :. ::::::::.::::::.::.:::::::::.:. ::::::::::::
CCDS90 -HK---AELQ---GSDEDEHVRAPTGHHQAGHCLMWACKACKRKSTTMDRRKAATMRERR
                     50        60        70        80        90    

              130       140       150       160       170       180
pF1KB7 RLSKVNEAFETLKRCTSSNPNQRLPKVEILRNAIRYIEGLQALLRDQDAAPPGAAAAFYA
       ::.:::.:::::::::..::::::::::::::::::::.:: :::.:       .  .:.
CCDS90 RLKKVNQAFETLKRCTTTNPNQRLPKVEILRNAIRYIESLQELLREQ-------VENYYS
          100       110       120       130       140              

              190       200       210       220       230       240
pF1KB7 PGPLPPGRGGEHYSGDSDASSPRSNCSDGMMDYSGPPSGARRRNCYEGAYYNEAPSEPRP
            ::..       :. .:: ::::::: . ..:   .:. . ... :  .. .    
CCDS90 L----PGQSC------SEPTSPTSNCSDGMPECNSP-VWSRKSSTFDSIYCPDVSNVYAT
           150             160       170        180       190      

              250        260       270       280       290         
pF1KB7 GKSAAVSSLDCLSSIVERI-STESPAAPALLLADVPSESPPRRQEAAAPSEGESSGDPTQ
        :.. .:::::::.::.:: :.:.:. :   : :. : ::    ..   . : ::     
CCDS90 DKNS-LSSLDCLSNIVDRITSSEQPGLP---LQDLASLSPVASTDSQPATPGASSSRLIY
        200        210       220          230       240       250  

     300       310       320
pF1KB7 SPDAAPQCPAGANPNPIYQVL
                            
CCDS90 HVL                  
                            

>>CCDS1433.1 MYOG gene_id:4656|Hs108|chr1                 (224 aa)
 initn: 444 init1: 382 opt: 422  Z-score: 335.4  bits: 69.8 E(32554): 2.3e-12
Smith-Waterman score: 434; 41.7% identity (64.3% similar) in 199 aa overlap (23-217:4-183)

               10        20        30        40        50        60
pF1KB7 MELLSPPLRDVDLTAPDGSLCSFATTDDFYDDPCFDSPDLRFFEDLDPRLMHVGALLKPE
                             . :.  ::..:       ::..  .   .:. .. .: 
CCDS14                    MELYETSPYFYQEP-------RFYDGENYLPVHLQGF-EPP
                                  10               20        30    

               70         80        90       100       110         
pF1KB7 EHSHFPAAVHP-APGAREDEHVRAPSGHHQAGRCLLWACKACKRKTTNADRRKAATMRER
        . .   .. : :::  ::. . .:  .:  :.:: ::::.::::....:::.:::.::.
CCDS14 GYERTELTLSPEAPGPLEDKGLGTP--EHCPGQCLPWACKVCKRKSVSVDRRRAATLREK
            40        50          60        70        80        90 

     120       130       140       150       160       170         
pF1KB7 RRLSKVNEAFETLKRCTSSNPNQRLPKVEILRNAIRYIEGLQALLRDQDAAPPGAAAAFY
       :::.:::::::.::: :  :::::::::::::.::.::: ::::: . .          :
CCDS14 RRLKKVNEAFEALKRSTLLNPNQRLPKVEILRSAIQYIERLQALLSSLNQEERDLR---Y
             100       110       120       130       140           

     180       190       200          210       220       230      
pF1KB7 APGPLPPGRGGEHYSGDSDASSPRSNCSD---GMMDYSGPPSGARRRNCYEGAYYNEAPS
         :      :: . .  :. ::  ..::    . ...:. :                   
CCDS14 RGG------GGPQPGVPSECSSHSASCSPEWGSALEFSANPGDHLLTADPTDAHNLHSLT
      150             160       170       180       190       200  

        240       250       260       270       280       290      
pF1KB7 EPRPGKSAAVSSLDCLSSIVERISTESPAAPALLLADVPSESPPRRQEAAAPSEGESSGD
                                                                   
CCDS14 SIVDSITVEDVSVAFPDETMPN                                      
            210       220                                          

>>CCDS9019.1 MYF6 gene_id:4618|Hs108|chr12                (242 aa)
 initn: 472 init1: 388 opt: 420  Z-score: 333.4  bits: 69.5 E(32554): 3e-12
Smith-Waterman score: 446; 44.9% identity (58.4% similar) in 214 aa overlap (57-267:41-234)

         30        40        50        60        70        80      
pF1KB7 DDFYDDPCFDSPDLRFFEDLDPRLMHVGALLKPEEHSHFPAAVHPAPGAREDEHVRAPSG
                                     :.: . .  : :   . :   .::: :: :
CCDS90 YFFYLDGENVTLQPLEVAEGSPLYPGSDGTLSPCQDQMPPEAGSDSSG---EEHVLAPPG
               20        30        40        50           60       

            90       100       110       120       130       140   
pF1KB7 ---HHQAGRCLLWACKACKRKTTNADRRKAATMRERRRLSKVNEAFETLKRCTSSNPNQR
           :  :.::.::::.::::.. .:::::::.::::::.:.:::::.::: : .:::::
CCDS90 LQPPHCPGQCLIWACKTCKRKSAPTDRRKAATLRERRRLKKINEAFEALKRRTVANPNQR
        70        80        90       100       110       120       

           150       160       170       180       190       200   
pF1KB7 LPKVEILRNAIRYIEGLQALLRDQDAAPPGAAAAFYAPGPLPPGRGGEHYSGDSDASSPR
       ::::::::.:: ::: :: ::.  :        .   :    : .  :.  :   :.  :
CCDS90 LPKVEILRSAISYIERLQDLLHRLDQQEKMQELGV-DPFSYRPKQ--ENLEG---ADFLR
       130       140       150       160        170            180 

           210       220       230       240       250       260   
pF1KB7 SNCSDGMMDYSGPPSGARRRNCYEGAYYNEAPSEPRPGKSAAVSSLDCLSSIVERISTES
       . ::.   . :    :        ::  .          :.: ::: ::::::. ::.: 
CCDS90 T-CSSQWPSVSDHSRGLVITAKEGGASID----------SSASSSLRCLSSIVDSISSEE
              190       200                 210       220       230

           270       280       290       300       310       320
pF1KB7 PAAPALLLADVPSESPPRRQEAAAPSEGESSGDPTQSPDAAPQCPAGANPNPIYQVL
          :                                                     
CCDS90 RKLPCVEEVVEK                                             
              240                                               




320 residues in 1 query   sequences
18511270 residues in 32554 library sequences
 Tcomplib [36.3.4 Apr, 2011] (8 proc)
 start: Fri Nov  4 09:04:18 2016 done: Fri Nov  4 09:04:19 2016
 Total Scan time:  2.990 Total Display time: -0.030

Function used was FASTA [36.3.4 Apr, 2011]