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

Statistics:  Expectation_n fit: rho(ln(x))= 5.0534+/-0.000513; mu= 13.3440+/- 0.031
 mean_var=64.0734+/-12.628, 0's: 0 Z-trim(114.8): 19  B-trim: 0 in 0/50
 Lambda= 0.160227
 statistics sampled from 15322 (15340) to 15322 sequences
Algorithm: FASTA (3.7 Nov 2010) [optimized]
Parameters: BL50 matrix (15:-5), open/ext: -10/-2
 ktup: 2, E-join: 1 (0.834), E-opt: 0.2 (0.471), width:  16
 Scan time:  1.220

The best scores are:                                      opt bits E(32554)
CCDS58345.1 TRPM1 gene_id:4308|Hs108|chr15         ( 136)  950 226.9 3.1e-60
CCDS10024.2 TRPM1 gene_id:4308|Hs108|chr15         (1603)  518 127.7 2.6e-29
CCDS58346.1 TRPM1 gene_id:4308|Hs108|chr15         (1625)  518 127.7 2.6e-29
CCDS58347.1 TRPM1 gene_id:4308|Hs108|chr15         (1642)  518 127.7 2.6e-29
CCDS43835.1 TRPM3 gene_id:80036|Hs108|chr9         (1707)  301 77.5 3.4e-14


>>CCDS58345.1 TRPM1 gene_id:4308|Hs108|chr15              (136 aa)
 initn: 950 init1: 950 opt: 950  Z-score: 1195.0  bits: 226.9 E(32554): 3.1e-60
Smith-Waterman score: 950; 99.3% identity (100.0% similar) in 136 aa overlap (1-136:1-136)

               10        20        30        40        50        60
pF1KE5 MKDSNRCCCGQFTNQHIPPLPSATPSKNEEENKQVETQPEKWSVAKHTQSYPTDSYGVLE
       :::::::::::::::::::::::::::::::.::::::::::::::::::::::::::::
CCDS58 MKDSNRCCCGQFTNQHIPPLPSATPSKNEEESKQVETQPEKWSVAKHTQSYPTDSYGVLE
               10        20        30        40        50        60

               70        80        90       100       110       120
pF1KE5 FQGGGYSNKAMVRKAFRHGATRITAFIGGQSPSPKLQIPGLLHGCGSIFLDISLKNQEIY
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS58 FQGGGYSNKAMVRKAFRHGATRITAFIGGQSPSPKLQIPGLLHGCGSIFLDISLKNQEIY
               70        80        90       100       110       120

              130      
pF1KE5 LCTWLLAMRLGNWTPL
       ::::::::::::::::
CCDS58 LCTWLLAMRLGNWTPL
              130      

>>CCDS10024.2 TRPM1 gene_id:4308|Hs108|chr15              (1603 aa)
 initn: 522 init1: 499 opt: 518  Z-score: 639.7  bits: 127.7 E(32554): 2.6e-29
Smith-Waterman score: 518; 60.9% identity (76.8% similar) in 138 aa overlap (1-133:1-134)

               10        20        30        40        50        60
pF1KE5 MKDSNRCCCGQFTNQHIPPLPSATPSKNEEENKQVETQPEKWSVAKHTQSYPTDSYGVLE
       :::::::::::::::::::::::::::::::.::::::::::::::::::::::::::::
CCDS10 MKDSNRCCCGQFTNQHIPPLPSATPSKNEEESKQVETQPEKWSVAKHTQSYPTDSYGVLE
               10        20        30        40        50        60

               70         80        90       100          110      
pF1KE5 FQGGGYSNKAM-VRKAFRHGATRITAFIGGQSPSPKLQIPGLL---HGCGSIFLDISLKN
       ::::::::::: .: ..    :.  ...  .  . .:..: ::   :: :   .... : 
CCDS10 FQGGGYSNKAMYIRVSYD---TKPDSLLHLMVKDWQLELPKLLISVHG-GLQNFEMQPKL
               70           80        90       100        110      

        120        130                                             
pF1KE5 QEIYLCTWL-LAMRLGNWTPL                                       
       ....    .  ::  : :                                          
CCDS10 KQVFGKGLIKAAMTTGAWIFTGGVSTGVISHVGDALKDHSSKSRGRVCAIGIAPWGIVEN
        120       130       140       150       160       170      

>>CCDS58346.1 TRPM1 gene_id:4308|Hs108|chr15              (1625 aa)
 initn: 522 init1: 499 opt: 518  Z-score: 639.7  bits: 127.7 E(32554): 2.6e-29
Smith-Waterman score: 518; 60.9% identity (76.8% similar) in 138 aa overlap (1-133:23-156)

                                     10        20        30        
pF1KE5                       MKDSNRCCCGQFTNQHIPPLPSATPSKNEEENKQVETQ
                             :::::::::::::::::::::::::::::::.::::::
CCDS58 MGQKSWIEKTFCKRECIFVIPSMKDSNRCCCGQFTNQHIPPLPSATPSKNEEESKQVETQ
               10        20        30        40        50        60

       40        50        60        70         80        90       
pF1KE5 PEKWSVAKHTQSYPTDSYGVLEFQGGGYSNKAM-VRKAFRHGATRITAFIGGQSPSPKLQ
       ::::::::::::::::::::::::::::::::: .: ..    :.  ...  .  . .:.
CCDS58 PEKWSVAKHTQSYPTDSYGVLEFQGGGYSNKAMYIRVSYD---TKPDSLLHLMVKDWQLE
               70        80        90       100          110       

       100          110       120        130                       
pF1KE5 IPGLL---HGCGSIFLDISLKNQEIYLCTWL-LAMRLGNWTPL                 
       .: ::   :: :   .... : ....    .  ::  : :                    
CCDS58 LPKLLISVHG-GLQNFEMQPKLKQVFGKGLIKAAMTTGAWIFTGGVSTGVISHVGDALKD
       120        130       140       150       160       170      

CCDS58 HSSKSRGRVCAIGIAPWGIVENKEDLVGKDVTRVYQTMSNPLSKLSVLNNSHTHFILADN
        180       190       200       210       220       230      

>>CCDS58347.1 TRPM1 gene_id:4308|Hs108|chr15              (1642 aa)
 initn: 522 init1: 499 opt: 518  Z-score: 639.6  bits: 127.7 E(32554): 2.6e-29
Smith-Waterman score: 518; 60.9% identity (76.8% similar) in 138 aa overlap (1-133:40-173)

                                             10        20        30
pF1KE5                               MKDSNRCCCGQFTNQHIPPLPSATPSKNEE
                                     ::::::::::::::::::::::::::::::
CCDS58 KSSTSGSQKGQKSWIEKTFCKRECIFVIPSMKDSNRCCCGQFTNQHIPPLPSATPSKNEE
      10        20        30        40        50        60         

               40        50        60        70         80         
pF1KE5 ENKQVETQPEKWSVAKHTQSYPTDSYGVLEFQGGGYSNKAM-VRKAFRHGATRITAFIGG
       :.::::::::::::::::::::::::::::::::::::::: .: ..    :.  ...  
CCDS58 ESKQVETQPEKWSVAKHTQSYPTDSYGVLEFQGGGYSNKAMYIRVSYD---TKPDSLLHL
      70        80        90       100       110          120      

      90       100          110       120        130               
pF1KE5 QSPSPKLQIPGLL---HGCGSIFLDISLKNQEIYLCTWL-LAMRLGNWTPL         
       .  . .:..: ::   :: :   .... : ....    .  ::  : :            
CCDS58 MVKDWQLELPKLLISVHG-GLQNFEMQPKLKQVFGKGLIKAAMTTGAWIFTGGVSTGVIS
        130       140        150       160       170       180     

CCDS58 HVGDALKDHSSKSRGRVCAIGIAPWGIVENKEDLVGKDVTRVYQTMSNPLSKLSVLNNSH
         190       200       210       220       230       240     

>>CCDS43835.1 TRPM3 gene_id:80036|Hs108|chr9              (1707 aa)
 initn: 217 init1: 194 opt: 301  Z-score: 368.2  bits: 77.5 E(32554): 3.4e-14
Smith-Waterman score: 301; 37.9% identity (65.7% similar) in 140 aa overlap (2-133:82-217)

                                            10        20        30 
pF1KE5                              MKDSNRCCCGQFTNQHIPPLPSATPSKNEEE
                                     :: .:::::.. .::.   :: .  .::..
CCDS43 LPSVRLLKAQKSWIERAFYKRECVHIIPSTKDPHRCCCGRLIGQHVGLTPSISVLQNEKN
              60        70        80        90       100       110 

                 40        50        60        70         80       
pF1KE5 NKQV---ETQPEKWSVAKHTQSYPTDSYGVLEFQGGGYSNKAM-VRKAFRHGATRITAFI
       ....   . : ::::..::::  :::..:..::::::.::::: :: .:    :.   ..
CCDS43 ESRLSRNDIQSEKWSISKHTQLSPTDAFGTIEFQGGGHSNKAMYVRVSF---DTKPDLLL
             120       130       140       150       160           

        90       100          110       120        130             
pF1KE5 GGQSPSPKLQIPGLL---HGCGSIFLDISLKNQEIYLCTWL-LAMRLGNWTPL       
         ..   .:..: ::   :: :   .... : ....    .  ::  : :          
CCDS43 HLMTKEWQLELPKLLISVHG-GLQNFELQPKLKQVFGKGLIKAAMTTGAWIFTGGVNTGV
      170       180        190       200       210       220       

CCDS43 IRHVGDALKDHASKSRGKICTIGIAPWGIVENQEDLIGRDVVRPYQTMSNPMSKLTVLNS
       230       240       250       260       270       280       




136 residues in 1 query   sequences
18511270 residues in 32554 library sequences
 Tcomplib [36.3.4 Apr, 2011] (8 proc)
 start: Mon Nov  7 22:24:30 2016 done: Mon Nov  7 22:24:30 2016
 Total Scan time:  1.220 Total Display time: -0.010

Function used was FASTA [36.3.4 Apr, 2011]