Motif confirmation by HMMer


hmmpfam - search a single seq against HMM database
HMMER 2.1.1 (Dec 1998)
Copyright (C) 1992-1998 Washington University School of Medicine
HMMER is freely distributed under the GNU General Public License (GPL).
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
HMM file:                 /home/niguts/usr02/tshini/ykclst/db/Pfam_fs
Sequence file:            hmm2/seq/D1046.3
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  D1046.3  CE03096   mitochondrial carrier protein status:Confirmed TR:Q18934 protein_id:CAA92291.1

Scores for sequence family classification (score includes all domains):
Model     Description                                   Score    E-value  N
--------  -----------                                   -----    ------- ---
mito_carr Mitochondrial carrier protein                 161.1    1.9e-44   3

Parsed for domains:
Model     Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------  ------- ----- -----    ----- -----      -----  -------
mito_carr   1/3       7    86 ..     1   125 []    41.9  4.6e-10
mito_carr   2/3      88   176 ..     1   125 []    37.9  5.8e-09
mito_carr   3/3     180   264 ..     1   118 [.    86.2  3.7e-22

Alignments of top-scoring domains:
mito_carr: domain 1 of 3, from 7 to 86: score 41.9, E = 4.6e-10
                   *->sflksllaGgiAGaiaatvtyPlDvvKtRLQtqgssslskssvnsdv
                      s+++ l++G++AG++  +  yPlD++K R+Q+ +
     D1046.3     7    SVVRWLVCGATAGLAVDIGLYPLDTIKSRMQSKQ------------- 40

                   gksksyrsakksnasasrkYkgildcfkkivkeEGliktrgLYkGlgpnl
                                                 ++  G +  +  Y+G +  l
     D1046.3    41 ----------------------------GFIAAGG-F--KDIYRGMISVL 59

                   lrvapygaaiyFgtYEqlKklllkkkge<-*
                   ++ ap+ aai+F+tY ++   ++   +e
     D1046.3    60 VGSAPG-AAIFFLTYKYINGQMKQVIEE    86

mito_carr: domain 2 of 3, from 88 to 176: score 37.9, E = 5.8e-09
                   *->sflksllaGgiAGaiaatvtyPlDvvKtRLQtqgssslskssvnsdv
                      + l++++++++A ++a +v  P ++ K R Q+ +   l+
     D1046.3    88    NALVDAVSASLAEIAACAVRVPTELCKQRGQVNKNERLTL------- 127

                   gksksyrsakksnasasrkYkgildcfkkivkeEGliktrgLYkGlgpnl
                                              k+i+ + G i  rg+Y+G+g ++
     D1046.3   128 -------------------------ICKEIMETKG-I--RGFYRGYGSTV 149

                   lrvapygaaiyFgtYEqlKklllkkkge<-*
                   +r +p+ + i+F + E lK+   +kk++
     D1046.3   150 AREIPF-SIIQFPIWEALKRAVANKKES    176

mito_carr: domain 3 of 3, from 180 to 264: score 86.2, E = 3.7e-22
                   *->sflksllaGgiAGaiaatvtyPlDvvKtRLQtqgssslskssvnsdv
                      s+l+ +++G++AG+iaa +t+PlDv+KtR+++ + ++++
     D1046.3   180    SPLEGAACGSVAGFIAAGLTTPLDVAKTRIMLTKNGPAP-------- 218

                   gksksyrsakksnasasrkYkgildcfkkivkeEGliktrgLYkGlgpnl
                                        gil+++k+++ ++G +  rgLY G++p++
     D1046.3   219 ---------------------GILSTLKEVYTSNG-V--RGLYSGVVPRV 244

                   lrvapygaaiyFgtYEqlKkl<-*
                   ++++ + + ++Fg+YE+ +++
     D1046.3   245 MWISGG-GFVFFGAYETAMHF    264



[DB home][top]