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/F08F8.2
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: F08F8.2 CE30649 Hydroxymethylglutaryl-coA reductase status:Confirmed TR:Q19207 protein_id:AAA68386.2
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
HMG-CoA_red Hydroxymethylglutaryl-coenzyme A redu 897.0 1.7e-272 1
CUT CUT domain 4.9 2.8 1
Glyco_tran_28_C Glycosyltransferase family 28 C-termi 3.6 5.9 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CUT 1/1 236 255 .. 69 88 .] 4.9 2.8
HMG-CoA_red 1/1 112 500 .. 1 426 [] 897.0 1.7e-272
Glyco_tran_28_C 1/1 473 523 .. 125 179 .] 3.6 5.9
Alignments of top-scoring domains:
CUT: domain 1 of 1, from 236 to 255: score 4.9, E = 2.8
*->MqnWLsdpnavrdlilqqek<-*
M++WL+ p+++ +++ + ++
F08F8.2 236 MKRWLEHPENQDRARQEFQS 255
HMG-CoA_red: domain 1 of 1, from 112 to 500: score 897.0, E = 1.7e-272
*->kiplyeLEsrLgadperAveiRReaieklatgkklearLkklPyegy
+ +++eLEsr++ +e+A+ iRR++i+k k+e+ lPy gy
F08F8.2 112 QAKSRELESRFP--AEQAIPIRRTFINK-----KFEN----LPYLGY 147
DyekvlgaCCENVIGYvpiPvGvAGPLllDGkyakGeyyvPMATTEGaLV
Dy++++++CCENVIGY p+PvGvAGPL+l+G++ e+yvPMATTEGaL+
F08F8.2 148 DYTLATECCCENVIGYTPVPVGVAGPLTLNGTS---EIYVPMATTEGALI 194
AStNRGckaitasGGatsvvlkdgMTRAPvVrfpsvvrAaelklwlespE
AStNRG+++i+a+GG++++++++gMTRAPvV+fp++++A+++k wle+pE
F08F8.2 195 ASTNRGMNVIRAAGGVETSIFNSGMTRAPVVKFPTARDAVSMKRWLEHPE 244
NfealkeaFdsLTSRFgrLqsIkiiviDtiaGrnlYlRFratTGDAMGMN
N+++++++F+s +SRF++L+sI+i ti+G+ +YlRF+a+TGDAMGMN
F08F8.2 245 NQDRARQEFQS-CSRFAKLKSIDI----TIDGNLAYLRFDAHTGDAMGMN 289
MVSKgvEkaLefleeitgegfpdMeiiSlSGNyCtDKKpAAiNwieGRGK
M+SK ++++++fl+e+ fp+M++++lSGN+C+DKK+AA Nw+eGRG+
F08F8.2 290 MISKSCDSTMRFLMEN----FPEMTVLALSGNLCVDKKAAAKNWTEGRGR 335
sVVaEatIpgdvVrkvLKTtvealvelNiaKNLiGSAmAGSiGgnGFNAH
sVVaE++Ip++vV+k+L+Tt+e+l++l+++K +iGS++AG++G G NAH
F08F8.2 336 SVVAECLIPREVVTKTLRTTPEQLAYLTTTKLHIGSSRAGAVG--GSNAH 383
aANIVtAiFiATGQDpAQnVESSnCiTlmEvveePMaVGlvGGAdggDLy
aANIV+AiFiATGQD+AQ+V SS+C+T mEv++ +++Ly
F08F8.2 384 AANIVAAIFIATGQDAAQVVSSSMCSTRMEVTA------------DKNLY 421
iSVTMPSlEVGTVGGGTrLapQsACLdlLGVkGanatsPGqNAhqLAriv
+S+T+P++EVGTVGGGT+LapQ+ACL+ LG+ G+n+++PGqNA++LA+++
F08F8.2 422 VSCTLPCVEVGTVGGGTILAPQRACLESLGCAGPNKEQPGQNAERLAEVI 471
aaaVLAGELSLmsALaaghLVksHmklNR<-*
aa+VLAGELSLm+AL++++LV+sHmklNR
F08F8.2 472 AATVLAGELSLMAALTTNELVSSHMKLNR 500
Glyco_tran_28_C: domain 1 of 1, from 473 to 523: score 3.6, E = 5.9
*->vkagaaevllqkeltllpekLaeallklllkplklyemakaagaakr
++++a+e l + lt +++L++ +kl++++++ly+ ++ +a
F08F8.2 473 ATVLAGELSLMAALT--TNELVSSHMKLNRSKQQLYADDSG--KATH 515
ikdaierl<-*
+ + +e+
F08F8.2 516 FEKEVEKA 523