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: hmm3/seq/F52H3.7a
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: F52H3.7a CE03403 locus:lec-2 galactoside-binding lectin status:Partially_confirmed TR:Q20684 protein_id:CAA91327.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
Gal-bind_lectin Galactoside-binding lectin 281.8 8.7e-81 2
pro_isomerase Cyclophilin type peptidyl-prolyl cis- 4.7 3.1 1
isoamylase_N Isoamylase N-terminal domain 4.2 5.4 1
cytochr_b558a Cytochrome Cytochrome b558 alpha-subu 3.1 2 1
Cl_perox Vanadium chloroperoxidase 1.6 6.5 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
Cl_perox 1/1 304 314 .. 603 613 .] 1.6 6.5
cytochr_b558a 1/1 374 396 .. 169 191 .] 3.1 2
pro_isomerase 1/1 864 890 .. 149 180 .] 4.7 3.1
Gal-bind_lectin 1/2 994 1125 .. 1 139 [] 95.1 2.3e-26
isoamylase_N 1/1 1221 1230 .. 112 121 .] 4.2 5.4
Gal-bind_lectin 2/2 1133 1261 .. 1 139 [] 191.4 2.5e-54
Alignments of top-scoring domains:
Cl_perox: domain 1 of 1, from 304 to 314: score 1.6, E = 6.5
*->EEqAPvvdvAP<-*
EEq Pvvd P
F52H3.7a 304 EEQTPVVDASP 314
cytochr_b558a: domain 1 of 1, from 374 to 396: score 3.1, E = 2
*->eeaaavpaGgPqvnPiPvtDEvv<-*
ee+aa+++ q++P P DE
F52H3.7a 374 EEVAAAATSEDQKEPAPASDEHT 396
pro_isomerase: domain 1 of 1, from 864 to 890: score 4.7, E = 3.1
*->VkkIenvgtdeyGksgdkpskeVkIadcgqle<-*
V++Ie+v+td ++ + ++V+ +g +e
F52H3.7a 864 VESIESVPTD----AQ-SVEEPVVQQHVGTVE 890
Gal-bind_lectin: domain 1 of 2, from 994 to 1125: score 95.1, E = 2.3e-26
*->pglvalnlglkpGktltVkGtVapknakrFavNlgkgskE..Endlv
p+ + l + + pG+tl+VkG+ ++ rF vNl++ + + + nd+
F52H3.7a 994 PYRSKLVEPFEPGQTLIVKGKTGE-DSIRFTVNLHTNTADfsGNDVP 1039
LHfNPRFneaHGDqntvVcNSkenGDNeWGtEqReaafPFqaGqpFeisI
LH+ +RF+ ++ V N++ +G eWG+E+R +++P ++G++++i+I
F52H3.7a 1040 LHVSVRFD-EG----KLVFNTFSKG--EWGKEER-KSNPYKKGDDIDIRI 1081
sveedkfkVkvndghef.eFphRlkleavqyLgikGDikltsikf<-*
+ + +kf+++v + e++e+ hR +l+++ + i+GD+ + +i +
F52H3.7a 1082 RAHDSKFQIFVD-QKEVkEYEHRVPLSSITHFTIDGDVLVNYIHW 1125
isoamylase_N: domain 1 of 1, from 1221 to 1230: score 4.2, E = 5.4
*->PYAkavvgrp<-*
PYA+av++++
F52H3.7a 1221 PYAFAVTVNG 1230
Gal-bind_lectin: domain 2 of 2, from 1133 to 1261: score 191.4, E = 2.5e-54
*->pglvaln.lglkpGktltVkGtVapknakrFavNlgkgskEEndlvL
p++++l+++gl pGktltV+G++++ +akrF++Nl+k++ +d++L
F52H3.7a 1133 PYESGLAgEGLAPGKTLTVFGIPEK-KAKRFHINLLKKN---GDIAL 1175
HfNPRFneaHGDqntvVcNSkenGDNeWGtEqReaafPFqaGqpFeisIs
H+N RF+e+H vV+NS++n+ +WG+E+Re+++PF++ ++F++ I+
F52H3.7a 1176 HLNARFDEKH-----VVRNSLINS--AWGNEEREGKMPFEKAVGFDLEIH 1218
veedkfkVkvndghefeFphRlkleavqyLgikGDikltsikf<-*
+e+++f V+vn+++ ++++hRl++ +v++L+i GD+++t+i++
F52H3.7a 1219 NEPYAFAVTVNGERFASYAHRLSPDEVNGLQIGGDVEITGIQL 1261