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/C56G2.6
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: C56G2.6 CE30639 locus:dhs-10 status:Confirmed SW:Q09517 protein_id:AAA64333.2
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
adh_short short chain dehydrogenase 113.0 6.4e-31 1
3HCDH_N 3-hydroxyacyl-CoA dehydrogenase, NAD bi 4.6 3.1 1
Lipoprotein_9 NLPA lipoprotein 3.0 8.3 1
Saccharop_dh Saccharopine dehydrogenase 2.8 3.4 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
3HCDH_N 1/1 68 91 .. 24 47 .. 4.6 3.1
Saccharop_dh 1/1 5 94 .. 1 47 [. 2.8 3.4
Lipoprotein_9 1/1 85 103 .. 232 253 .] 3.0 8.3
adh_short 1/1 46 239 .. 1 197 [. 113.0 6.4e-31
Alignments of top-scoring domains:
3HCDH_N: domain 1 of 1, from 68 to 91: score 4.6, E = 3.1
*->AskGfpVvlkDinpkaLdkalkeI<-*
A+ Gf+V+l++++++ Ld +keI
C56G2.6 68 ARRGFNVLLVSRTQSKLDETKKEI 91
Saccharop_dh: domain 1 of 1, from 5 to 94: score 2.8, E = 3.4
*->VLvlGAGgVsqpvaeklarnsdvfs......................
+ ++GAG+V+ +++ + r ++fs+ ++ ++ + +++ + +
C56G2.6 5 CFLVGAGYVALAAVAY--RLLTIFSnilgpyvllspidlkkragasw 49
.......................eitvASrtkskadalaasiagk<-*
++ +++ ++ + +++ ++++ Srt+sk d+++++i +k
C56G2.6 50 avvtgatdgigkayafelarrgfNVLLVSRTQSKLDETKKEILEK 94
Lipoprotein_9: domain 1 of 1, from 85 to 103: score 3.0, E = 8.3
*->eeVkkfiekkyPYykGaviPaw<-*
+e+kk+i +k Y ++ +
C56G2.6 85 DETKKEILEK--Y-SSIEVRTA 103
adh_short: domain 1 of 1, from 46 to 239: score 113.0, E = 6.4e-31
*->tgKvaLvTGassGIGlaiAkrLakeGakVvvvdrreekaeqvaaelk
+a+vTGa+ GIG+a+A +La++G++V++v r+++k++++ +e+
C56G2.6 46 GASWAVVTGATDGIGKAYAFELARRGFNVLLVSRTQSKLDETKKEIL 92
aelG.dralfiqlDvtdee.qvkaavaqaverlGd.rlDvLVNNAGilgd
+ + + +++++++D t ++++ a++ ++ l + vL+NN+G++ +
C56G2.6 93 EKYSsIEVRTAAFDFTNAApS---AYKDLLATLNQvEIGVLINNVGMSYE 139
gppfe...elseedwervidvNltGvflltravlpamdhmlkrkgGrIvN
p+ ++ + e + +N + ll++ lp +m +rk+G+IvN
C56G2.6 140 YPDVLhkvDGGIERLANITTINTLPPTLLSAGILP---QMVARKAGVIVN 186
isSvaGlrvgvpglsaYsASKaavigltrsLAlElaphgtgIrVnavaPG
++S aG + + +++YsA+K v lt L++E+ g I+V +aP
C56G2.6 187 VGSSAGAN-QMALWAVYSATKKYVSWLTAILRKEYEHQG--ITVQTIAPM 233
gvdTdm<-*
+v+T m
C56G2.6 234 MVATKM 239