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: hmm5/seq/Y57A10A.18
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: Y57A10A.18 CE28715 locus:pqn-87 status:Partially_confirmed TR:Q9U207 protein_id:CAB60772.2
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
MATH MATH domain 67.1 3.4e-17 1
FlaG FlaG protein 7.5 1.5 1
TCL1_MTCP1 TCL1/MTCP1 family 6.0 4.6 1
Colicin Colicin pore forming domain 4.4 4 1
DUF691 Protein of unknown function (DUF691) 1.3 4.3 1
complex1_49Kd Respiratory-chain NADH dehydrogenase, 4 -0.2 9.5 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
complex1_49Kd 1/1 20 25 .. 1 6 [. -0.2 9.5
FlaG 1/1 110 148 .. 1 39 [. 7.5 1.5
MATH 1/1 181 306 .. 1 159 [] 67.1 3.4e-17
Colicin 1/1 316 326 .. 193 203 .] 4.4 4
DUF691 1/1 378 405 .. 53 80 .. 1.3 4.3
TCL1_MTCP1 1/1 1068 1086 .. 32 52 .] 6.0 4.6
Alignments of top-scoring domains:
complex1_49Kd: domain 1 of 1, from 20 to 25: score -0.2, E = 9.5
*->DvGAlT<-*
DvG++T
Y57A10A.18 20 DVGCQT 25
FlaG: domain 1 of 1, from 110 to 148: score 7.5, E = 1.5
*->atsksleqtiqqqkssaaqsvkeaaeenaaeeqkisked<-*
+ ++s+++t++qq+++ +q+ ++a ++ + +e k+++ +
Y57A10A.18 110 SLNSSTNTTQPQQAGQTTQKARAASTTGVVGEGKKQNSN 148
MATH: domain 1 of 1, from 181 to 306: score 67.1, E = 3.4e-17
*->tiknfSkikeeakeGregeeyytspyeerfnipWrlnvlriyrnggg
+i fS+++ e + sp ++ ++ W+ +++ ++ +
Y57A10A.18 181 VINRFSQLE----------ETMCSPPKSIEGVSWK---IMVMPKQHM 214
egrs.gkflglyLhClkeekdsptienlkWsieteftlklvsdnGksirr
++ ++k +g++L+C +e+ s +Ws+++ ++++s + +
Y57A10A.18 215 VQKKqQKCMGFFLQCAPERAYS-----DQWSVHAIADMRMISYKPN---V 256
mssttltkktkdaknnshvFekp.tgegwgksGfkkFisWddLeddyngy
++ ++++t+ h+++ ++++g+++ F+ W d+ d+ gy
Y57A10A.18 257 PH---FARRTT------HTYTSKeNDWGYSC-----FMTWADIIDEAQGY 292
lvdDsiiiEaeVkI<-*
+ dD+++ E+ Vk
Y57A10A.18 293 IRDDTVVLEIAVKA 306
Colicin: domain 1 of 1, from 316 to 326: score 4.4, E = 4
*->dkfveklNkli<-*
++f++k+Nk+i
Y57A10A.18 316 EDFLDKINKWI 326
DUF691: domain 1 of 1, from 378 to 405: score 1.3, E = 4.3
*->eaeeiekPPVPkavskepkatsLlakmr<-*
++e+iek PV ++ ++ k+tsL+++++
Y57A10A.18 378 SIERIEKGPVVCTDTTHGKPTSLRQALT 405
TCL1_MTCP1: domain 1 of 1, from 1068 to 1086: score 6.0, E = 4.6
*->vvihiksrlvedvLLrlldvd<-*
+v+ +k+r++ ++ +l+dv+
Y57A10A.18 1068 IVEREKERDR--ARKELEDVT 1086