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: hmm4/seq/W02B12.9
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: W02B12.9 CE03769 mitochondrial RNA splicing MSR4 like protein status:Confirmed TR:Q23125 protein_id:CAA91399.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
mito_carr Mitochondrial carrier protein 240.3 2.7e-68 3
FliH Flagellar assembly protein FliH 5.1 1.7 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
FliH 1/1 39 49 .. 248 258 .] 5.1 1.7
mito_carr 1/3 16 110 .. 1 125 [] 91.8 1.1e-23
mito_carr 2/3 113 203 .. 1 125 [] 65.4 1.7e-16
mito_carr 3/3 204 310 .. 1 125 [] 83.3 2.3e-21
Alignments of top-scoring domains:
FliH: domain 1 of 1, from 39 to 49: score 5.1, E = 1.7
*->ervksvlqelC<-*
++vk+++q+lC
W02B12.9 39 DSVKTRMQSLC 49
mito_carr: domain 1 of 3, from 16 to 110: score 91.8, E = 1.1e-23
*->sflksllaGgiAGaiaatvtyPlDvvKtRLQtqgssslskssvnsdv
s++++l aG++AGa+ ++v++P+D vKtR+Q+ ++ +
W02B12.9 16 SVPVHLTAGALAGAVEHCVMFPFDSVKTRMQSLCPCPET-------- 54
gksksyrsakksnasasrkYkgildcfkkivkeEGliktrgLYkGlgpnl
k ++++ +++ivk EG++ r +G+ +
W02B12.9 55 ------------------KCPTPVHSLMSIVKREGWL--R-PLRGVNAVA 83
lrvapygaaiyFgtYEqlKklllkkkge<-*
++ +p+ a+yF++YE++K +l++++ +
W02B12.9 84 AGSMPA-HALYFTVYEKMKGYLTGNSAG 110
mito_carr: domain 2 of 3, from 113 to 203: score 65.4, E = 1.7e-16
*->sflksllaGgiAGaiaatvtyPlDvvKtRLQtqgssslskssvnsdv
++l++ ++G +A +i ++++P +vvK R+Q+ s
W02B12.9 113 NTLAYGASGVVATLIHDAIMNPAEVVKQRMQMAFS------------ 147
gksksyrsakksnasasrkYkgildcfkkivkeEGliktrgLYkGlgpnl
Y + l+c++ ++ EG + +++Y+ ++ l
W02B12.9 148 ------------------PYGSSLECARCVYNREG-V--AAFYRSYTTQL 176
lrvapygaaiyFgtYEqlKklllkkkge<-*
+ ++p+ +ai F++YE+ +++l +
W02B12.9 177 AMNVPF-QAIHFMSYEFWQHVLNPEHKY 203
mito_carr: domain 3 of 3, from 204 to 310: score 83.3, E = 2.3e-21
*->sflksllaGgiAGaiaatvtyPlDvvKtRLQtqgssslsk.ssvnsd
++ ++l+aGg+AG++aa++t+P+D vKt L tq++ +++ ++ +
W02B12.9 204 DPKSHLIAGGLAGGLAAALTTPMDCVKTVLNTQQAAEADPaNR---R 247
vgksksyrsakksnasasrkYkgildcfkkivkeEGliktrgLYkGlgpn
++ + y +Y+gi d++++i+ + G + +g+ Gl ++
W02B12.9 248 IFLQARY------------RYRGISDAVRTIYSQRG-L--SGFSCGLQAR 282
llrvapygaaiyFgtYEqlKklllkkkge<-*
++ ++p+ +a+ +++YE +K +l+ +g
W02B12.9 283 VIFQVPA-TALSWSVYELFKFMLSFEGGH 310