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/D2013.8
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: D2013.8 CE17621 locus:scp-1 status:Partially_confirmed TR:Q18968 protein_id:CAA87777.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
WD40 WD domain, G-beta repeat 63.8 2.4e-16 3
IspA Intacellular septation protein A 6.6 1.2 1
Cation_ATPase_N Cation transporter/ATPase, N-terminus 4.3 9.5 1
Patched Patched family 0.7 7.9 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
IspA 1/1 266 276 .. 1 11 [. 6.6 1.2
Patched 1/1 392 428 .. 392 428 .. 0.7 7.9
Cation_ATPase_N 1/1 645 659 .. 70 87 .] 4.3 9.5
WD40 1/3 694 733 .. 1 39 [] 27.2 6.2e-06
WD40 2/3 813 849 .. 1 39 [] 26.0 1.3e-05
WD40 3/3 857 895 .. 1 39 [] 12.0 0.13
Alignments of top-scoring domains:
IspA: domain 1 of 1, from 266 to 276: score 6.6, E = 1.2
*->MKqllDliPLI<-*
+K+++D+iPLI
D2013.8 266 LKTFSDYIPLI 276
Patched: domain 1 of 1, from 392 to 428: score 0.7, E = 7.9
*->IGaiTptPeIqlFCiytavAiffdfiYQlTFFaAila<-*
IG++T + Iq FC ++ +++ +df Ql F+a++l+
D2013.8 392 IGFLTRISDIQEFCQFSVICVTVDFYMQLFFYAPCLT 428
Cation_ATPase_N: domain 1 of 1, from 645 to 659: score 4.3, E = 9.5
*->llLlaAAiiSALdfvlyg<-*
+lLl+ +iS f++y
D2013.8 645 WLLLITTVIS---FFAYV 659
WD40: domain 1 of 3, from 694 to 733: score 27.2, E = 6.2e-06
*->geclatl.gHtspvtcvafspdsgpylasgSrDgtiklWd<-*
++ +++gH p+++va++p++ ++++s++ g +++W+
D2013.8 694 DTLPIVYqGHRFPIESVAIDPEDTSTFVSCCQEGVVYVWN 733
WD40: domain 2 of 3, from 813 to 849: score 26.0, E = 1.3e-05
*->geclatl.gHtspvtcvafspdsgpylasgSrDgtiklWd<-*
e+l++ ++H++pv+ +a+ + +l+++S D++ik+W+
D2013.8 813 IELLKSVrAHQKPVCRIAIWKS---QLITSSFDRSIKMWN 849
WD40: domain 3 of 3, from 857 to 895: score 12.0, E = 0.13
*->geclatl.gHtspvtcvafspdsgpylasgSrDgtiklWd<-*
+ ++ +H+spv +a++ +++ s+++ g i+ W+
D2013.8 857 IDISNVFlAHNSPVVNLAVDEA-TSIMYSSCEEGVICWWN 895