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/Y47D9A.1a
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: Y47D9A.1a CE22073 mannose-1-phosphate gaunyl transferase status:Partially_confirmed TR:Q9N4V3 protein_id:AAF60647.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
NTP_transferase Nucleotidyl transferase 52.9 5.5e-14 1
hexapep Bacterial transferase hexapeptide (th 26.2 2.9e-05 3
GatB_N PET112 family, N terminal region 5.2 0.98 1
CTP_transf_3 Cytidylyltransferase 3.1 8.9 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CTP_transf_3 1/1 27 54 .. 16 44 .. 3.1 8.9
NTP_transferase 1/1 4 191 .. 1 235 [. 52.9 5.5e-14
GatB_N 1/1 186 200 .. 275 289 .] 5.2 0.98
hexapep 1/3 266 283 .. 1 18 [] 6.7 7.6
hexapep 2/3 284 301 .. 1 18 [] 13.8 0.08
hexapep 3/3 318 335 .. 1 18 [] 5.7 14
Alignments of top-scoring domains:
CTP_transf_3: domain 1 of 1, from 27 to 54: score 3.1, E = 8.9
*->KnLlklgGkPLiahvieaAlqSklfdkdv<-*
K+L +++G+PLi+h i+++ q +++ + +
Y47D9A.1a 27 KPLFPIAGVPLIEHHIDQLCQLSGLSE-I 54
NTP_transferase: domain 1 of 1, from 4 to 191: score 52.9, E = 5.5e-14
*->kaiILAGG..sGTRLwPLTrtlaKplvpvlgkypmiqytlsrlanlg
ka++L+GG+++GTR +PL++ l+Kpl+p++g p+i + + l +
Y47D9A.1a 4 KAVVLVGGpqKGTRFRPLSLQLPKPLFPIAGV-PLIEHHIDQLCQ-L 48
agirepivictqehrflvaeqlgdgskfglgvdGfVevtyalqpeprpag
+g e i+ ++ + ++ + + + + + V++ y+ +p p+
Y47D9A.1a 49 SGLSE-ILLLGFFPSDVFTDFISRCQQTY-----RVSIKYLEEPNPL--- 89
kkWFqGTApAvalaadflgdddsknedelvLVLggDhvigmeedfeqalq
GTA+ + + + d ++ v+V ++D+ + +e++ +
Y47D9A.1a 90 -----GTAGGLVSFKKQILAGD----PDAVFVINADVCGDL--PIEDMGA 128
karakaadgktvtfgivpvedptg.....YGviefdelesveegqdtevy
k ++ + +++ + ++t++++ ++G +++d +eg
Y47D9A.1a 129 KLDSLSG-----SSMLMLTTEATRqqsinFGSVVTD-----SEG------ 162
rVirFvEKPdlpkAkmqVDTtLGlpeakssnyasmGlYlfnsgvf<-*
rVi+ v KP s+ +s G+Yl + +v+
Y47D9A.1a 163 RVIHYVDKPTTFV----------------STNISCGVYLIKAEVI 191
GatB_N: domain 1 of 1, from 186 to 200: score 5.2, E = 0.98
*->IeyEiqRQvnllerG<-*
I++E++RQ++l +G
Y47D9A.1a 186 IKAEVIRQLDLPLNG 200
hexapep: domain 1 of 3, from 266 to 283: score 6.7, E = 7.6
*->gviIGdnvvIgpgvvIgd<-*
+v+I + ++++p ++Ig+
Y47D9A.1a 266 DVFIDPSAKVHPTAKIGP 283
hexapep: domain 2 of 3, from 284 to 301: score 13.8, E = 0.08
*->gviIGdnvvIgpgvvIgd<-*
+v IG+ +vIg gv+I +
Y47D9A.1a 284 NVSIGPKSVIGKGVRIKE 301
hexapep: domain 3 of 3, from 318 to 335: score 5.7, E = 14
*->gviIGdnvvIgpgvvIgd<-*
++IG ++v+g ++I +
Y47D9A.1a 318 QSVIGWRSVVGMWARIEG 335