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

Statistics:  Expectation_n fit: rho(ln(x))= 4.7656+/-0.000834; mu= 18.7176+/- 0.050
 mean_var=64.8781+/-13.431, 0's: 0 Z-trim(106.2): 21  B-trim: 13 in 2/46
 Lambda= 0.159230
 statistics sampled from 8820 (8824) to 8820 sequences
Algorithm: FASTA (3.7 Nov 2010) [optimized]
Parameters: BL50 matrix (15:-5), open/ext: -10/-2
 ktup: 2, E-join: 1 (0.642), E-opt: 0.2 (0.271), width:  16
 Scan time:  2.000

The best scores are:                                      opt bits E(32554)
CCDS31424.1 TMEM41B gene_id:440026|Hs108|chr11     ( 291) 1895 443.9 6.5e-125
CCDS53600.1 TMEM41B gene_id:440026|Hs108|chr11     ( 127)  785 188.6   2e-48
CCDS3271.1 TMEM41A gene_id:90407|Hs108|chr3        ( 264)  558 136.7 1.7e-32


>>CCDS31424.1 TMEM41B gene_id:440026|Hs108|chr11          (291 aa)
 initn: 1895 init1: 1895 opt: 1895  Z-score: 2356.0  bits: 443.9 E(32554): 6.5e-125
Smith-Waterman score: 1895; 100.0% identity (100.0% similar) in 291 aa overlap (1-291:1-291)

               10        20        30        40        50        60
pF1KSD MAKGRVAERSQLGAHHTTPVGDGAAGTRGLAAPGSRDHQKEKSWVEAGSARMSLLILVSI
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS31 MAKGRVAERSQLGAHHTTPVGDGAAGTRGLAAPGSRDHQKEKSWVEAGSARMSLLILVSI
               10        20        30        40        50        60

               70        80        90       100       110       120
pF1KSD FLSAAFVMFLVYKNFPQLSEEERVNMKVPRDMDDAKALGKVLSKYKDTFYVQVLVAYFAT
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS31 FLSAAFVMFLVYKNFPQLSEEERVNMKVPRDMDDAKALGKVLSKYKDTFYVQVLVAYFAT
               70        80        90       100       110       120

              130       140       150       160       170       180
pF1KSD YIFLQTFAIPGSIFLSILSGFLYPFPLALFLVCLCSGLGASFCYMLSYLVGRPVVYKYLT
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS31 YIFLQTFAIPGSIFLSILSGFLYPFPLALFLVCLCSGLGASFCYMLSYLVGRPVVYKYLT
              130       140       150       160       170       180

              190       200       210       220       230       240
pF1KSD EKAVKWSQQVERHREHLINYIIFLRITPFLPNWFINITSPVINVPLKVFFIGTFLGVAPP
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS31 EKAVKWSQQVERHREHLINYIIFLRITPFLPNWFINITSPVINVPLKVFFIGTFLGVAPP
              190       200       210       220       230       240

              250       260       270       280       290 
pF1KSD SFVAIKAGTTLYQLTTAGEAVSWNSIFILMILAVLSILPAIFQKKLKQKFE
       :::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS31 SFVAIKAGTTLYQLTTAGEAVSWNSIFILMILAVLSILPAIFQKKLKQKFE
              250       260       270       280       290 

>>CCDS53600.1 TMEM41B gene_id:440026|Hs108|chr11          (127 aa)
 initn: 785 init1: 785 opt: 785  Z-score: 982.9  bits: 188.6 E(32554): 2e-48
Smith-Waterman score: 785; 99.2% identity (100.0% similar) in 123 aa overlap (1-123:1-123)

               10        20        30        40        50        60
pF1KSD MAKGRVAERSQLGAHHTTPVGDGAAGTRGLAAPGSRDHQKEKSWVEAGSARMSLLILVSI
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS53 MAKGRVAERSQLGAHHTTPVGDGAAGTRGLAAPGSRDHQKEKSWVEAGSARMSLLILVSI
               10        20        30        40        50        60

               70        80        90       100       110       120
pF1KSD FLSAAFVMFLVYKNFPQLSEEERVNMKVPRDMDDAKALGKVLSKYKDTFYVQVLVAYFAT
       ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS53 FLSAAFVMFLVYKNFPQLSEEERVNMKVPRDMDDAKALGKVLSKYKDTFYVQVLVAYFAT
               70        80        90       100       110       120

              130       140       150       160       170       180
pF1KSD YIFLQTFAIPGSIFLSILSGFLYPFPLALFLVCLCSGLGASFCYMLSYLVGRPVVYKYLT
       ::.                                                         
CCDS53 YILYPFN                                                     
                                                                   

>>CCDS3271.1 TMEM41A gene_id:90407|Hs108|chr3             (264 aa)
 initn: 492 init1: 492 opt: 558  Z-score: 696.7  bits: 136.7 E(32554): 1.7e-32
Smith-Waterman score: 558; 32.2% identity (72.3% similar) in 242 aa overlap (54-289:4-245)

            30        40        50        60        70           80
pF1KSD AAGTRGLAAPGSRDHQKEKSWVEAGSARMSLLILVSIFLSAAFVMFLVYKNFPQ---LSE
                                     :: :. .: . .:...:.   .:.   :. 
CCDS32                            MRPLLGLLLVFAGCTFALYLLSTRLPRGRRLGS
                                          10        20        30   

                  90       100       110       120       130       
pF1KSD EERV---NMKVPRDMDDAKALGKVLSKYKDTFYVQVLVAYFATYIFLQTFAIPGSIFLSI
        :..   ..  : :. . . :..:: .:.    . :.. . ..:.. : :::::: ::..
CCDS32 TEEAGGRSLWFPSDLAELRELSEVLREYRKEHQAYVFLLFCGAYLYKQGFAIPGSSFLNV
            40        50        60        70        80        90   

       140       150       160       170       180       190       
pF1KSD LSGFLYPFPLALFLVCLCSGLGASFCYMLSYLVGRPVVYKYLTEKAVKWSQQVERHREHL
       :.: :.   :.:.: :. ...::. ::.:: . :. .: .:. .:..  ...::..:. :
CCDS32 LAGALFGPWLGLLLCCVLTSVGATCCYLLSSIFGKQLVVSYFPDKVALLQRKVEENRNSL
           100       110       120       130       140       150   

       200       210       220       230       240       250       
pF1KSD INYIIFLRITPFLPNWFINITSPVINVPLKVFFIGTFLGVAPPSFVAIKAGTTLYQLTTA
       . ...:::. :. ::::.:...:..:.:.  ::.....:. : .:. ...:. :  ::. 
CCDS32 FFFLLFLRLFPMTPNWFLNLSAPILNIPIVQFFFSVLIGLIPYNFICVQTGSILSTLTSL
           160       170       180       190       200       210   

       260       270       280       290                  
pF1KSD GEAVSWNSIFILMILAVLSILPAIFQKKLKQKFE                 
           ::...: :. .:.....:. . ::..::                   
CCDS32 DALFSWDTVFKLLAIAMVALIPGTLIKKFSQKHLQLNETSTANHIHSRKDT
           220       230       240       250       260    




291 residues in 1 query   sequences
18511270 residues in 32554 library sequences
 Tcomplib [36.3.4 Apr, 2011] (8 proc)
 start: Wed Nov  2 23:30:48 2016 done: Wed Nov  2 23:30:49 2016
 Total Scan time:  2.000 Total Display time: -0.010

Function used was FASTA [36.3.4 Apr, 2011]