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

Statistics:  Expectation_n fit: rho(ln(x))= 4.9098+/-0.000701; mu= 13.5924+/- 0.042
 mean_var=51.9932+/-10.555, 0's: 0 Z-trim(108.2): 18  B-trim: 0 in 0/51
 Lambda= 0.177869
 statistics sampled from 10076 (10085) to 10076 sequences
Algorithm: FASTA (3.7 Nov 2010) [optimized]
Parameters: BL50 matrix (15:-5), open/ext: -10/-2
 ktup: 2, E-join: 1 (0.712), E-opt: 0.2 (0.31), width:  16
 Scan time:  1.960

The best scores are:                                      opt bits E(32554)
CCDS3070.2 NUDT16 gene_id:131870|Hs108|chr3        ( 195) 1050 276.8 4.9e-75
CCDS54641.1 NUDT16 gene_id:131870|Hs108|chr3       ( 159)  828 219.8 5.8e-58
CCDS54640.1 NUDT16 gene_id:131870|Hs108|chr3       ( 227)  678 181.3 3.1e-46
CCDS10519.1 NUDT16L1 gene_id:84309|Hs108|chr16     ( 211)  573 154.4 3.7e-38
CCDS59257.1 NUDT16L1 gene_id:84309|Hs108|chr16     ( 192)  400 110.0 7.9e-25


>>CCDS3070.2 NUDT16 gene_id:131870|Hs108|chr3             (195 aa)
 initn: 1050 init1: 1050 opt: 1050  Z-score: 1460.6  bits: 276.8 E(32554): 4.9e-75
Smith-Waterman score: 1050; 100.0% identity (100.0% similar) in 162 aa overlap (1-162:34-195)

                                             10        20        30
pF1KE5                               MLFGRIPLRYAILMQMRFDGRLGFPGGFVD
                                     ::::::::::::::::::::::::::::::
CCDS30 ARRLELGEALALGSGWRHACHALLYAPDPGMLFGRIPLRYAILMQMRFDGRLGFPGGFVD
            10        20        30        40        50        60   

               40        50        60        70        80        90
pF1KE5 TQDRSLEDGLNRELREELGEAAAAFRVERTDYRSSHVGSGPRVVAHFYAKRLTLEELLAV
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS30 TQDRSLEDGLNRELREELGEAAAAFRVERTDYRSSHVGSGPRVVAHFYAKRLTLEELLAV
            70        80        90       100       110       120   

              100       110       120       130       140       150
pF1KE5 EAGATRAKDHGLEVLGLVRVPLYTLRDGVGGLPTFLENSFIGSAREQLLEALQDLGLLQS
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS30 EAGATRAKDHGLEVLGLVRVPLYTLRDGVGGLPTFLENSFIGSAREQLLEALQDLGLLQS
           130       140       150       160       170       180   

              160  
pF1KE5 GSISGLKIPAHH
       ::::::::::::
CCDS30 GSISGLKIPAHH
           190     

>>CCDS54641.1 NUDT16 gene_id:131870|Hs108|chr3            (159 aa)
 initn: 828 init1: 828 opt: 828  Z-score: 1154.1  bits: 219.8 E(32554): 5.8e-58
Smith-Waterman score: 828; 100.0% identity (100.0% similar) in 128 aa overlap (14-141:1-128)

               10        20        30        40        50        60
pF1KE5 MLFGRIPLRYAILMQMRFDGRLGFPGGFVDTQDRSLEDGLNRELREELGEAAAAFRVERT
                    :::::::::::::::::::::::::::::::::::::::::::::::
CCDS54              MQMRFDGRLGFPGGFVDTQDRSLEDGLNRELREELGEAAAAFRVERT
                            10        20        30        40       

               70        80        90       100       110       120
pF1KE5 DYRSSHVGSGPRVVAHFYAKRLTLEELLAVEAGATRAKDHGLEVLGLVRVPLYTLRDGVG
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS54 DYRSSHVGSGPRVVAHFYAKRLTLEELLAVEAGATRAKDHGLEVLGLVRVPLYTLRDGVG
        50        60        70        80        90       100       

              130       140       150       160            
pF1KE5 GLPTFLENSFIGSAREQLLEALQDLGLLQSGSISGLKIPAHH          
       :::::::::::::::::::::                               
CCDS54 GLPTFLENSFIGSAREQLLEAALHGPMKTEMRTLVLGREGRTWECFLIGSER
       110       120       130       140       150         

>>CCDS54640.1 NUDT16 gene_id:131870|Hs108|chr3            (227 aa)
 initn: 678 init1: 678 opt: 678  Z-score: 943.6  bits: 181.3 E(32554): 3.1e-46
Smith-Waterman score: 678; 100.0% identity (100.0% similar) in 104 aa overlap (1-104:34-137)

                                             10        20        30
pF1KE5                               MLFGRIPLRYAILMQMRFDGRLGFPGGFVD
                                     ::::::::::::::::::::::::::::::
CCDS54 ARRLELGEALALGSGWRHACHALLYAPDPGMLFGRIPLRYAILMQMRFDGRLGFPGGFVD
            10        20        30        40        50        60   

               40        50        60        70        80        90
pF1KE5 TQDRSLEDGLNRELREELGEAAAAFRVERTDYRSSHVGSGPRVVAHFYAKRLTLEELLAV
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS54 TQDRSLEDGLNRELREELGEAAAAFRVERTDYRSSHVGSGPRVVAHFYAKRLTLEELLAV
            70        80        90       100       110       120   

              100       110       120       130       140       150
pF1KE5 EAGATRAKDHGLEVLGLVRVPLYTLRDGVGGLPTFLENSFIGSAREQLLEALQDLGLLQS
       ::::::::::::::                                              
CCDS54 EAGATRAKDHGLEVGPAWDSVPFPISSSPKAFSPPRKHPWRKVFAPLTLPSPQLSWWSWD
           130       140       150       160       170       180   

>>CCDS10519.1 NUDT16L1 gene_id:84309|Hs108|chr16          (211 aa)
 initn: 586 init1: 314 opt: 573  Z-score: 798.5  bits: 154.4 E(32554): 3.7e-38
Smith-Waterman score: 573; 61.5% identity (82.4% similar) in 148 aa overlap (2-148:40-183)

                                            10        20        30 
pF1KE5                              MLFGRIPLRYAILMQMRFDGRLGFPGGFVDT
                                     ::::::.:...:::::::: ::::::::: 
CCDS10 KQISRVEAMRLGPGWSHSCHAMLYAANPGQLFGRIPMRFSVLMQMRFDGLLGFPGGFVDR
      10        20        30        40        50        60         

              40        50        60        70         80        90
pF1KE5 QDRSLEDGLNRELREELGEAAAAFRVERTDYRSSHVGSGP-RVVAHFYAKRLTLEELLAV
       .  :::::::: :   ::     .:. ..:: :::.  :: :::::.::..::::.: ::
CCDS10 RFWSLEDGLNRVLGLGLG----CLRLTEADYLSSHLTEGPHRVVAHLYARQLTLEQLHAV
      70        80            90       100       110       120     

              100       110       120       130       140       150
pF1KE5 EAGATRAKDHGLEVLGLVRVPLYTLRDGVGGLPTFLENSFIGSAREQLLEALQDLGLLQS
       : .:....:::::::::::::::: .: :::.:.:: :.:...:. ::: ::. :...  
CCDS10 EISAVHSRDHGLEVLGLVRVPLYTQKDRVGGFPNFLSNAFVSTAKCQLLFALKVLNMMPE
         130       140       150       160       170       180     

              160                
pF1KE5 GSISGLKIPAHH              
                                 
CCDS10 EKLVEALAAATEKQKKALEKLLPASS
         190       200       210 

>>CCDS59257.1 NUDT16L1 gene_id:84309|Hs108|chr16          (192 aa)
 initn: 360 init1: 233 opt: 400  Z-score: 559.2  bits: 110.0 E(32554): 7.9e-25
Smith-Waterman score: 400; 62.5% identity (80.8% similar) in 104 aa overlap (2-104:40-139)

                                            10        20        30 
pF1KE5                              MLFGRIPLRYAILMQMRFDGRLGFPGGFVDT
                                     ::::::.:...:::::::: ::::::::: 
CCDS59 KQISRVEAMRLGPGWSHSCHAMLYAANPGQLFGRIPMRFSVLMQMRFDGLLGFPGGFVDR
      10        20        30        40        50        60         

              40        50        60        70         80        90
pF1KE5 QDRSLEDGLNRELREELGEAAAAFRVERTDYRSSHVGSGP-RVVAHFYAKRLTLEELLAV
       .  :::::::: :   ::     .:. ..:: :::.  :: :::::.::..::::.: ::
CCDS59 RFWSLEDGLNRVLGLGLG----CLRLTEADYLSSHLTEGPHRVVAHLYARQLTLEQLHAV
      70        80            90       100       110       120     

              100       110       120       130       140       150
pF1KE5 EAGATRAKDHGLEVLGLVRVPLYTLRDGVGGLPTFLENSFIGSAREQLLEALQDLGLLQS
       : .:....::::::                                              
CCDS59 EISAVHSRDHGLEVGPPPGPRPPPRGLALAPWKAPMGAGPRAGPAVHPEGPSRRLPQLPE
         130       140       150       160       170       180     




162 residues in 1 query   sequences
18511270 residues in 32554 library sequences
 Tcomplib [36.3.4 Apr, 2011] (8 proc)
 start: Mon Nov  7 22:15:25 2016 done: Mon Nov  7 22:15:25 2016
 Total Scan time:  1.960 Total Display time: -0.020

Function used was FASTA [36.3.4 Apr, 2011]