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

Statistics:  Expectation_n fit: rho(ln(x))= 4.9593+/-0.000881; mu= 10.3384+/- 0.053
 mean_var=50.9615+/-10.369, 0's: 0 Z-trim(103.8): 15  B-trim: 0 in 0/50
 Lambda= 0.179661
 statistics sampled from 7570 (7575) to 7570 sequences
Algorithm: FASTA (3.7 Nov 2010) [optimized]
Parameters: BL50 matrix (15:-5), open/ext: -10/-2
 ktup: 2, E-join: 1 (0.613), E-opt: 0.2 (0.233), width:  16
 Scan time:  0.970

The best scores are:                                      opt bits E(32554)
CCDS54809.1 GYPB gene_id:2994|Hs108|chr4           (  91)  537 146.7 1.9e-36
CCDS47138.1 GYPE gene_id:2996|Hs108|chr4           (  78)  276 79.0 3.9e-16
CCDS77965.1 GYPA gene_id:2993|Hs108|chr4           ( 137)  254 73.4 3.3e-14
CCDS34069.1 GYPA gene_id:2993|Hs108|chr4           ( 150)  254 73.4 3.6e-14
CCDS82959.1 GYPA gene_id:2993|Hs108|chr4           ( 117)  237 68.9 6.1e-13


>>CCDS54809.1 GYPB gene_id:2994|Hs108|chr4                (91 aa)
 initn: 537 init1: 537 opt: 537  Z-score: 767.8  bits: 146.7 E(32554): 1.9e-36
Smith-Waterman score: 537; 100.0% identity (100.0% similar) in 91 aa overlap (1-91:1-91)

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

               70        80        90 
pF1KE6 VVIILIILCVMAGIIGTILLISYSIRRLIKA
       :::::::::::::::::::::::::::::::
CCDS54 VVIILIILCVMAGIIGTILLISYSIRRLIKA
               70        80        90 

>>CCDS47138.1 GYPE gene_id:2996|Hs108|chr4                (78 aa)
 initn: 298 init1: 237 opt: 276  Z-score: 403.2  bits: 79.0 E(32554): 3.9e-16
Smith-Waterman score: 276; 65.1% identity (76.7% similar) in 86 aa overlap (1-86:1-78)

               10        20        30        40        50        60
pF1KE6 MYGKIIFVLLLSEIVSISALSTTEVAMHTSTSSSVTKSYISSQTNGETGQLVHRFTVPAP
       :::::::::::: :::::: ::: :::::::::::::::::::::: :  :.. ..    
CCDS47 MYGKIIFVLLLSGIVSISASSTTGVAMHTSTSSSVTKSYISSQTNGIT--LINWWA----
               10        20        30        40          50        

               70        80        90 
pF1KE6 VVIILIILCVMAGIIGTILLISYSIRRLIKA
         .  .:. ::  ..: :.:::: ::     
CCDS47 --MARVIFEVMLVVVGMIILISYCIR     
             60        70             

>>CCDS77965.1 GYPA gene_id:2993|Hs108|chr4                (137 aa)
 initn: 372 init1: 254 opt: 254  Z-score: 368.5  bits: 73.4 E(32554): 3.3e-14
Smith-Waterman score: 346; 67.0% identity (75.5% similar) in 106 aa overlap (1-90:1-106)

               10        20        30        40                    
pF1KE6 MYGKIIFVLLLSEIVSISALSTTEVAMHTSTSSSVTKSYISSQTNG--------------
       :::::::::::::::::::::::::::::::::::::::::::::               
CCDS77 MYGKIIFVLLLSEIVSISALSTTEVAMHTSTSSSVTKSYISSQTNDTHKRDTYAATPRAH
               10        20        30        40        50        60

         50        60          70        80        90              
pF1KE6 ETGQLVHRFTVPAP--VVIILIILCVMAGIIGTILLISYSIRRLIKA             
       :....  : . :    . : :::. ::::.:::::::::.::::::              
CCDS77 EVSEISVRTVYPPEEETEITLIIFGVMAGVIGTILLISYGIRRLIKKSPSDVKPLPSPDT
               70        80        90       100       110       120

CCDS77 DVPLSSVEIENPETSDQ
              130       

>>CCDS34069.1 GYPA gene_id:2993|Hs108|chr4                (150 aa)
 initn: 374 init1: 254 opt: 254  Z-score: 367.9  bits: 73.4 E(32554): 3.6e-14
Smith-Waterman score: 343; 61.5% identity (66.4% similar) in 122 aa overlap (1-90:1-119)

               10        20        30        40                    
pF1KE6 MYGKIIFVLLLSEIVSISALSTTEVAMHTSTSSSVTKSYISSQTN---------------
       :::::::::::::::::::::::::::::::::::::::::::::               
CCDS34 MYGKIIFVLLLSEIVSISALSTTEVAMHTSTSSSVTKSYISSQTNDTHKRDTYAATPRAH
               10        20        30        40        50        60

                           50        60        70        80        
pF1KE6 -----------------GETGQLVHRFTVPAPVVIILIILCVMAGIIGTILLISYSIRRL
                        ::  ::.:.:. :    : :::. ::::.:::::::::.::::
CCDS34 EVSEISVRTVYPPEEETGERVQLAHHFSEPE---ITLIIFGVMAGVIGTILLISYGIRRL
               70        80        90          100       110       

       90                               
pF1KE6 IKA                              
       ::                               
CCDS34 IKKSPSDVKPLPSPDTDVPLSSVEIENPETSDQ
       120       130       140       150

>>CCDS82959.1 GYPA gene_id:2993|Hs108|chr4                (117 aa)
 initn: 194 init1: 134 opt: 237  Z-score: 345.8  bits: 68.9 E(32554): 6.1e-13
Smith-Waterman score: 237; 53.3% identity (70.0% similar) in 90 aa overlap (1-90:1-86)

               10        20        30        40        50        60
pF1KE6 MYGKIIFVLLLSEIVSISALSTTEVAMHTSTSSSVTKSYISSQTNGETGQLVHRFTVPAP
       ::::::::::::.  . .. ..:  : :  .  ::   :   . .::  ::.:.:. :  
CCDS82 MYGKIIFVLLLSDTHKRDTYAATPRA-HEVSEISVRTVYPPEEETGERVQLAHHFSEPE-
               10        20         30        40        50         

               70        80        90                              
pF1KE6 VVIILIILCVMAGIIGTILLISYSIRRLIKA                             
         : :::. ::::.:::::::::.::::::                              
CCDS82 --ITLIIFGVMAGVIGTILLISYGIRRLIKKSPSDVKPLPSPDTDVPLSSVEIENPETSD
         60        70        80        90       100       110      




91 residues in 1 query   sequences
18511270 residues in 32554 library sequences
 Tcomplib [36.3.4 Apr, 2011] (8 proc)
 start: Tue Nov  8 14:14:56 2016 done: Tue Nov  8 14:14:56 2016
 Total Scan time:  0.970 Total Display time: -0.020

Function used was FASTA [36.3.4 Apr, 2011]