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: hmm1/seq/C25A1.5
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: C25A1.5 CE08370 Yeast hypothetical 44.9 KD protein like status:Confirmed TR:Q9XVS0 protein_id:CAB02759.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
FA_hydroxylase Fatty acid hydroxylase 348.9 8.6e-102 1
heme_1 Heme/Steroid binding domain 58.2 4.5e-15 1
PqiA Paraquat-inducible protein A 4.1 4.7 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
heme_1 1/1 6 78 .. 12 78 .] 58.2 4.5e-15
PqiA 1/1 224 245 .. 111 132 .. 4.1 4.7
FA_hydroxylase 1/1 167 313 .. 1 155 [] 348.9 8.6e-102
Alignments of top-scoring domains:
heme_1: domain 1 of 1, from 6 to 78: score 58.2, E = 4.5e-15
*->kkdgrilvainGkVYDltrflafhpgGasallafaGrDaTaaf....
++++++lv ++G++YD+ +f+++hpgGa++l+++aG+++ +++++++
C25A1.5 6 DEKKPLLVQLDGQYYDIADFAPKHPGGAKVLNRLAGEEIGEFIrgek 52
....klaahreeaaeelleklykvGrvd<-*
+ + ++ h+ eaa+++le+ y v +
C25A1.5 53 rimgVRHEHS-EAAYNMLER-YNVNAIQ 78
PqiA: domain 1 of 1, from 224 to 245: score 4.1, E = 4.7
*->vPvlkillilyLllaarfkvkl<-*
vP+++i++i+yL++++ f++++
C25A1.5 224 VPATLIVGIFYLIYSNTFQWPV 245
FA_hydroxylase: domain 1 of 1, from 167 to 313: score 348.9, E = 8.6e-102
*->lfliGvflWTLiEYtLHRFLFHikpyspdSyiaitlHfllHGcHHKl
+f+iGv++WTL+EY+LHR++FH+kp spdS+++i+lHfl+HG+HHK+
C25A1.5 167 WFVIGVLTWTLTEYSLHRWVFHWKP-SPDSPNQILLHFLAHGLHHKT 212
PmDglRLVfPPtlasilvgpFyklvqlvlpDpavapAvfaGgLLGYVmYD
PmDg+RLVfPP++a+++vg+Fy++++++++ ++v++A++aG+L+GYV+YD
C25A1.5 213 PMDGDRLVFPPVPATLIVGIFYLIYSNTFQ-WPVFCAFGAGKLFGYVTYD 261
mtHYYLHHgsPpkpsylkkLKkyHlnHHFkdqdkGGshPFGiTSsLWDrV
m+HYYLHHgsP+++s+l+++K+yH+nHHFk++d+G FGi++sLWD+V
C25A1.5 262 MVHYYLHHGSPRPRSNLHFRKVYHHNHHFKNFDVG----FGISTSLWDYV 307
FGGttLpl<-*
F+ tL++
C25A1.5 308 FH--TLGM 313