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/ZK40.1
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: ZK40.1 CE28186 status:Partially_confirmed TR:Q23087 protein_id:AAB04844.2
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
Acyltransferase Acyltransferase 159.7 4.8e-44 1
HGWP HGWP repeat 5.2 4.4 1
CO_dh CO dehydrogenase beta subunit/acetyl- 3.6 5.3 1
GASA Gibberellin regulated protein 0.5 7.6 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
HGWP 1/1 61 82 .. 10 33 .] 5.2 4.4
GASA 1/1 183 189 .. 103 109 .] 0.5 7.6
Acyltransferase 1/1 80 225 .. 1 167 [] 159.7 4.8e-44
CO_dh 1/1 287 303 .. 164 180 .] 3.6 5.3
Alignments of top-scoring domains:
HGWP: domain 1 of 1, from 61 to 82: score 5.2, E = 4.4
*->rLhGcWpimPpllGvyrlvftagl<-*
rL G W imP l y +f a +
ZK40.1 61 RLVGIWVIMPGALCNY--IFGANI 82
GASA: domain 1 of 1, from 183 to 189: score 0.5, E = 7.6
*->gGgpKCP<-*
+G+ KCP
ZK40.1 183 EGTDKCP 189
Acyltransferase: domain 1 of 1, from 80 to 225: score 159.7, E = 4.8e-44
*->vevrglenlpe.kgpaiivsNHqSylDplvlslllykrgpwprplvf
++r +++++ +++pa+i++NH+++lD+l+++ +lyk++pw++ + +
ZK40.1 80 ANIRIKGDFINhEEPALIIMNHRTRLDWLFFWNALYKMDPWLCTTEK 126
iakkellggisvPllGwlmrlagaifidRnnhkdk.pllnel.yl.g.er
i+ k +l ++vP++Gw+m+ a +if+dR++++dk++l+n+l+y+ ++e+
ZK40.1 127 ISLKGML--KYVPGAGWAMQAASYIFLDRSFDTDKtKLDNILnYYaEtEY 174
.yqriwkrkaygaslleyavelLrkGelvviFPEGTRsrgkdetrygell
+yq +++FPEGT ++ ++t++++++
ZK40.1 175 kYQ-------------------------LLLFPEGTDKCP-KATERSRIH 198
pfKkGafrlAaellph.tgagvpIvPv<-*
+ KkG++++ ++l+p+ tg+ + ++ +
ZK40.1 199 SEKKGLVHYQYVLHPRvTGFVHIVQAM 225
CO_dh: domain 1 of 1, from 287 to 303: score 3.6, E = 5.3
*->NlskeeeeylemLeEml<-*
Nl+ +ee+l+ ++Em
ZK40.1 287 NLWRNKEEKLKRFYEMP 303