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

Statistics:  Expectation_n fit: rho(ln(x))= 7.0117+/-0.000789; mu= 9.7228+/- 0.048
 mean_var=153.2191+/-32.316, 0's: 0 Z-trim(113.2): 176  B-trim: 858 in 3/51
 Lambda= 0.103614
 statistics sampled from 13679 (13868) to 13679 sequences
Algorithm: FASTA (3.7 Nov 2010) [optimized]
Parameters: BL50 matrix (15:-5), open/ext: -10/-2
 ktup: 2, E-join: 1 (0.758), E-opt: 0.2 (0.426), width:  16
 Scan time:  2.970

The best scores are:                                      opt bits E(32554)
CCDS5401.1 HOXA1 gene_id:3198|Hs108|chr7           ( 335) 2307 356.2 2.2e-98
CCDS5402.2 HOXA1 gene_id:3198|Hs108|chr7           ( 137)  836 135.9 1.8e-32
CCDS32675.1 HOXB1 gene_id:3211|Hs108|chr17         ( 301)  625 104.7   1e-22
CCDS2271.1 HOXD1 gene_id:3231|Hs108|chr2           ( 328)  525 89.8 3.4e-18


>>CCDS5401.1 HOXA1 gene_id:3198|Hs108|chr7                (335 aa)
 initn: 2307 init1: 2307 opt: 2307  Z-score: 1879.6  bits: 356.2 E(32554): 2.2e-98
Smith-Waterman score: 2307; 99.7% identity (100.0% similar) in 335 aa overlap (1-335:1-335)

               10        20        30        40        50        60
pF1KB8 MDNARMNSFLEYPILSSGDSGTCSARAYPSDHRITTFQSCAVSANSCGGDDRFLVGRGVQ
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS54 MDNARMNSFLEYPILSSGDSGTCSARAYPSDHRITTFQSCAVSANSCGGDDRFLVGRGVQ
               10        20        30        40        50        60

               70        80        90       100       110       120
pF1KB8 IGSPHHHHHHHHHHPQPATYQTSGNLGVSYSHSSCGPSYGSQNFSAPYSPYALNQEADVS
       ::::::::::::.:::::::::::::::::::::::::::::::::::::::::::::::
CCDS54 IGSPHHHHHHHHRHPQPATYQTSGNLGVSYSHSSCGPSYGSQNFSAPYSPYALNQEADVS
               70        80        90       100       110       120

              130       140       150       160       170       180
pF1KB8 GGYPQCAPAVYSGNLSSPMVQHHHHHQGYAGGAVGSPQYIHHSYGQEHQSLALATYNNSL
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS54 GGYPQCAPAVYSGNLSSPMVQHHHHHQGYAGGAVGSPQYIHHSYGQEHQSLALATYNNSL
              130       140       150       160       170       180

              190       200       210       220       230       240
pF1KB8 SPLHASHQEACRSPASETSSPAQTFDWMKVKRNPPKTGKVGEYGYLGQPNAVRTNFTTKQ
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS54 SPLHASHQEACRSPASETSSPAQTFDWMKVKRNPPKTGKVGEYGYLGQPNAVRTNFTTKQ
              190       200       210       220       230       240

              250       260       270       280       290       300
pF1KB8 LTELEKEFHFNKYLTRARRVEIAASLQLNETQVKIWFQNRRMKQKKREKEGLLPISPATP
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS54 LTELEKEFHFNKYLTRARRVEIAASLQLNETQVKIWFQNRRMKQKKREKEGLLPISPATP
              250       260       270       280       290       300

              310       320       330     
pF1KB8 PGNDEKAEESSEKSSSSPCVPSPGSSTSDTLTTSH
       :::::::::::::::::::::::::::::::::::
CCDS54 PGNDEKAEESSEKSSSSPCVPSPGSSTSDTLTTSH
              310       320       330     

>>CCDS5402.2 HOXA1 gene_id:3198|Hs108|chr7                (137 aa)
 initn: 836 init1: 836 opt: 836  Z-score: 696.3  bits: 135.9 E(32554): 1.8e-32
Smith-Waterman score: 836; 99.2% identity (100.0% similar) in 118 aa overlap (1-118:1-118)

               10        20        30        40        50        60
pF1KB8 MDNARMNSFLEYPILSSGDSGTCSARAYPSDHRITTFQSCAVSANSCGGDDRFLVGRGVQ
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS54 MDNARMNSFLEYPILSSGDSGTCSARAYPSDHRITTFQSCAVSANSCGGDDRFLVGRGVQ
               10        20        30        40        50        60

               70        80        90       100       110       120
pF1KB8 IGSPHHHHHHHHHHPQPATYQTSGNLGVSYSHSSCGPSYGSQNFSAPYSPYALNQEADVS
       ::::::::::::.:::::::::::::::::::::::::::::::::::::::::::::  
CCDS54 IGSPHHHHHHHHRHPQPATYQTSGNLGVSYSHSSCGPSYGSQNFSAPYSPYALNQEADPP
               70        80        90       100       110       120

              130       140       150       160       170       180
pF1KB8 GGYPQCAPAVYSGNLSSPMVQHHHHHQGYAGGAVGSPQYIHHSYGQEHQSLALATYNNSL
                                                                   
CCDS54 RSLSLPRIGDIFSSADF                                           
              130                                                  

>>CCDS32675.1 HOXB1 gene_id:3211|Hs108|chr17              (301 aa)
 initn: 614 init1: 377 opt: 625  Z-score: 521.4  bits: 104.7 E(32554): 1e-22
Smith-Waterman score: 731; 45.9% identity (67.1% similar) in 340 aa overlap (1-328:1-301)

               10        20        30         40        50         
pF1KB8 MDNARMNSFLEYPILSSGDSGTCSARAYPSDHRITTFQ-SCAVSANSCGGDDRFLVGRGV
       ::  :::::::::. . : :.  ::.. :     :.:  : : ...: ... :.  : : 
CCDS32 MDYNRMNSFLEYPLCNRGPSAY-SAHSAP-----TSFPPSSAQAVDSYASEGRY--GGG-
               10        20              30        40          50  

      60        70             80        90       100       110    
pF1KB8 QIGSPHHHHHHHHHHPQPAT-----YQTSGNLGVSYSHSSCGPSYGSQNFSAPYSPYALN
        ..::  ...  .   :: .     . .:.  :  :. ..:.::::      : . : :.
CCDS32 -LSSPAFQQNSGYPAQQPPSTLGVPFPSSAPSG--YAPAACSPSYG------PSQYYPLG
               60        70        80          90             100  

           120       130       140       150       160         170 
pF1KB8 Q-EADVSGGYPQCAPAVYSGNLSSPMVQHHHHHQGYAGGAVGSPQYI--HHSYGQEHQSL
       : :.:  ::: .  :. :...:..         .::..:..:   :   :  ::.:. . 
CCDS32 QSEGD--GGYFH--PSSYGAQLGG-------LSDGYGAGGAGPGPYPPQHPPYGNEQTAS
              110         120              130       140       150 

             180       190       200        210       220       230
pF1KB8 ALATYNNSLSPLHASHQEACRSPASETSSP-AQTFDWMKVKRNPPKTGKVGEYGYLGQPN
          .: . ::    ...  :  : :: ..: :.::::::::::::::.::.: : ::.:.
CCDS32 FAPAYADLLS---EDKETPC--P-SEPNTPTARTFDWMKVKRNPPKTAKVSEPG-LGSPS
             160             170       180       190        200    

              240       250       260       270       280       290
pF1KB8 AVRTNFTTKQLTELEKEFHFNKYLTRARRVEIAASLQLNETQVKIWFQNRRMKQKKREKE
       ..::::::.:::::::::::::::.:::::::::.:.:::::::::::::::::::::.:
CCDS32 GLRTNFTTRQLTELEKEFHFNKYLSRARRVEIAATLELNETQVKIWFQNRRMKQKKRERE
          210       220       230       240       250       260    

               300        310       320       330     
pF1KB8 -GLLPISPATPPGND-EKAEESSEKSSSSPCVPSPGSSTSDTLTTSH
        : .:  :: :::   : : ..:..:. .    ::.: ::       
CCDS32 EGRVP--PA-PPGCPKEAAGDASDQSTCTSPEASPSSVTS       
            270        280       290       300        

>>CCDS2271.1 HOXD1 gene_id:3231|Hs108|chr2                (328 aa)
 initn: 535 init1: 499 opt: 525  Z-score: 440.1  bits: 89.8 E(32554): 3.4e-18
Smith-Waterman score: 535; 37.0% identity (58.7% similar) in 341 aa overlap (6-328:1-328)

               10             20        30        40        50     
pF1KB8 MDNARMNSFLEYPILSS-----GDSGTCSARAYPSDHRITTFQSCAVSANSCGGDDRFLV
            :.:.:::   ::     ::  . . .   :: : ...:     .:   ::  :. 
CCDS22      MSSYLEYVSCSSSGGVGGDVLSLAPKFCRSDARPVALQPAFPLGN---GDGAFVS
                    10        20        30        40           50  

          60        70           80               90       100     
pF1KB8 GRGVQIGSPHHHHHHHHHHPQ---PATYQTS-----GNL--GVSYSHSSCGPSYGSQNFS
          .  . :         .:.   ::. : .     :    :.. . .. : .::  . .
CCDS22 CLPLAAARPSPSPPAAPARPSVPPPAAPQYAQCTLEGAYEPGAAPAAAAGGADYGFLGSG
             60        70        80        90       100       110  

          110       120        130       140       150        160  
pF1KB8 APYS-PYALNQEADVSGGYPQCAP-AVYSGNLSSPMVQHHHHHQGYAGGAVGS-PQYIHH
         :. : .:.. :: .:.. . :  ::.::. .: ... .  . ..  :  :  :  .. 
CCDS22 PAYDFPGVLGRAADDGGSHVHYATSAVFSGG-GSFLLSGQVDYAAF--GEPGPFPACLKA
            120       130       140        150         160         

            170       180       190       200       210       220  
pF1KB8 SYGQEHQSLALATYNNSLSPLHASHQEACRSPASETSSPAQTFDWMKVKRNPPKTGKVGE
       :   .  ..  :      ::  ... ..  ::::   .  .::.:::::::  : ::..:
CCDS22 SADGHPGAFQTA------SPAPGTYPKSV-SPASGLPAAFSTFEWMKVKRNASKKGKLAE
     170       180             190        200       210       220  

            230       240       250       260       270       280  
pF1KB8 YGYLGQPNAVRTNFTTKQLTELEKEFHFNKYLTRARRVEIAASLQLNETQVKIWFQNRRM
       ::  .  .:.::::.::::::::::::::::::::::.:::  :.::.::::::::::::
CCDS22 YGAASPSSAIRTNFSTKQLTELEKEFHFNKYLTRARRIEIANCLHLNDTQVKIWFQNRRM
            230       240       250       260       270       280  

            290       300       310       320       330     
pF1KB8 KQKKREKEGLLPISPATPPGNDEKAEESSEKSSSSPCVPSPGSSTSDTLTTSH
       ::::::.::::  .  . : .   .  .  :  ..:  :: ..  :       
CCDS22 KQKKREREGLLATAIPVAPLQLPLSGTTPTKFIKNPGSPSQSQEPS       
            290       300       310       320               




335 residues in 1 query   sequences
18511270 residues in 32554 library sequences
 Tcomplib [36.3.4 Apr, 2011] (8 proc)
 start: Fri Nov  4 16:40:54 2016 done: Fri Nov  4 16:40:54 2016
 Total Scan time:  2.970 Total Display time:  0.000

Function used was FASTA [36.3.4 Apr, 2011]