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: pF1KE6113, 134 aa
  1>>>pF1KE6113 134 - 134 aa - 134 aa
Library: human.CCDS.faa
  18511270 residues in 32554 sequences

Statistics:  Expectation_n fit: rho(ln(x))= 4.8215+/-0.000631; mu= 13.9962+/- 0.038
 mean_var=60.0462+/-11.966, 0's: 0 Z-trim(110.1): 19  B-trim: 0 in 0/53
 Lambda= 0.165513
 statistics sampled from 11358 (11372) to 11358 sequences
Algorithm: FASTA (3.7 Nov 2010) [optimized]
Parameters: BL50 matrix (15:-5), open/ext: -10/-2
 ktup: 2, E-join: 1 (0.745), E-opt: 0.2 (0.349), width:  16
 Scan time:  1.550

The best scores are:                                      opt bits E(32554)
CCDS12004.1 CYB5A gene_id:1528|Hs108|chr18         ( 134)  897 221.7 1.1e-58
CCDS12005.1 CYB5A gene_id:1528|Hs108|chr18         (  98)  656 164.1 1.8e-41
CCDS54188.1 CYB5A gene_id:1528|Hs108|chr18         ( 124)  575 144.8 1.4e-35
CCDS10880.2 CYB5B gene_id:80777|Hs108|chr16        ( 150)  457 116.7   5e-27


>>CCDS12004.1 CYB5A gene_id:1528|Hs108|chr18              (134 aa)
 initn: 897 init1: 897 opt: 897  Z-score: 1167.2  bits: 221.7 E(32554): 1.1e-58
Smith-Waterman score: 897; 100.0% identity (100.0% similar) in 134 aa overlap (1-134:1-134)

               10        20        30        40        50        60
pF1KE6 MAEQSDEAVKYYTLEEIQKHNHSKSTWLILHHKVYDLTKFLEEHPGGEEVLREQAGGDAT
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS12 MAEQSDEAVKYYTLEEIQKHNHSKSTWLILHHKVYDLTKFLEEHPGGEEVLREQAGGDAT
               10        20        30        40        50        60

               70        80        90       100       110       120
pF1KE6 ENFEDVGHSTDAREMSKTFIIGELHPDDRPKLNKPPETLITTIDSSSSWWTNWVIPAISA
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS12 ENFEDVGHSTDAREMSKTFIIGELHPDDRPKLNKPPETLITTIDSSSSWWTNWVIPAISA
               70        80        90       100       110       120

              130    
pF1KE6 VAVALMYRLYMAED
       ::::::::::::::
CCDS12 VAVALMYRLYMAED
              130    

>>CCDS12005.1 CYB5A gene_id:1528|Hs108|chr18              (98 aa)
 initn: 656 init1: 656 opt: 656  Z-score: 858.2  bits: 164.1 E(32554): 1.8e-41
Smith-Waterman score: 656; 100.0% identity (100.0% similar) in 97 aa overlap (1-97:1-97)

               10        20        30        40        50        60
pF1KE6 MAEQSDEAVKYYTLEEIQKHNHSKSTWLILHHKVYDLTKFLEEHPGGEEVLREQAGGDAT
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS12 MAEQSDEAVKYYTLEEIQKHNHSKSTWLILHHKVYDLTKFLEEHPGGEEVLREQAGGDAT
               10        20        30        40        50        60

               70        80        90       100       110       120
pF1KE6 ENFEDVGHSTDAREMSKTFIIGELHPDDRPKLNKPPETLITTIDSSSSWWTNWVIPAISA
       :::::::::::::::::::::::::::::::::::::                       
CCDS12 ENFEDVGHSTDAREMSKTFIIGELHPDDRPKLNKPPEP                      
               70        80        90                              

              130    
pF1KE6 VAVALMYRLYMAED

>>CCDS54188.1 CYB5A gene_id:1528|Hs108|chr18              (124 aa)
 initn: 575 init1: 575 opt: 575  Z-score: 752.2  bits: 144.8 E(32554): 1.4e-35
Smith-Waterman score: 790; 92.5% identity (92.5% similar) in 134 aa overlap (1-134:1-124)

               10        20        30        40        50        60
pF1KE6 MAEQSDEAVKYYTLEEIQKHNHSKSTWLILHHKVYDLTKFLEEHPGGEEVLREQAGGDAT
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS54 MAEQSDEAVKYYTLEEIQKHNHSKSTWLILHHKVYDLTKFLEEHPGGEEVLREQAGGDAT
               10        20        30        40        50        60

               70        80        90       100       110       120
pF1KE6 ENFEDVGHSTDAREMSKTFIIGELHPDDRPKLNKPPETLITTIDSSSSWWTNWVIPAISA
       ::::::::::::::::::::::::::          ::::::::::::::::::::::::
CCDS54 ENFEDVGHSTDAREMSKTFIIGELHP----------ETLITTIDSSSSWWTNWVIPAISA
               70        80                  90       100       110

              130    
pF1KE6 VAVALMYRLYMAED
       ::::::::::::::
CCDS54 VAVALMYRLYMAED
              120    

>>CCDS10880.2 CYB5B gene_id:80777|Hs108|chr16             (150 aa)
 initn: 395 init1: 368 opt: 457  Z-score: 598.7  bits: 116.7 E(32554): 5e-27
Smith-Waterman score: 457; 50.4% identity (75.9% similar) in 133 aa overlap (3-133:18-146)

                              10        20        30        40     
pF1KE6                MAEQSDEAVKYYTLEEIQKHNHSKSTWLILHHKVYDLTKFLEEHP
                        .. . .: :: :::. :.:  :  ::..: .:::.:.::.:::
CCDS10 MSGSMATAEASGSDGKGQEVETSVTYYRLEEVAKRNSLKELWLVIHGRVYDVTRFLNEHP
               10        20        30        40        50        60

          50        60        70        80         90        100   
pF1KE6 GGEEVLREQAGGDATENFEDVGHSTDAREMSKTFIIGELHPDD-RPKL-NKPPETLITTI
       :::::: :::: ::.:.:::::::.::::: : . ::..::.: .:.  .: :    .  
CCDS10 GGEEVLLEQAGVDASESFEDVGHSSDAREMLKQYYIGDIHPSDLKPESGSKDP----SKN
               70        80        90       100       110          

           110       120       130       
pF1KE6 DSSSSWWTNWVIPAISAVAVALMYRLYMAED   
       :. .: :. :..: :.:: ....:: : .:    
CCDS10 DTCKSCWAYWILPIIGAVLLGFLYRYYTSESKSS
        120       130       140       150




134 residues in 1 query   sequences
18511270 residues in 32554 library sequences
 Tcomplib [36.3.4 Apr, 2011] (8 proc)
 start: Tue Nov  8 09:33:44 2016 done: Tue Nov  8 09:33:45 2016
 Total Scan time:  1.550 Total Display time: -0.030

Function used was FASTA [36.3.4 Apr, 2011]