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/W09H1.6b
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: W09H1.6b CE16577 locus:lec-1 Vertebrate galactoside-binding lectins (2 domains) status:Partially_confirmed TR:O45904 protein_id:CAB0496
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
Gal-bind_lectin Galactoside-binding lectin 347.4 1.5e-100 2
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
Gal-bind_lectin 1/2 18 149 .. 1 139 [] 122.0 3.6e-34
Gal-bind_lectin 2/2 157 284 .. 1 139 [] 232.9 2.2e-66
Alignments of top-scoring domains:
Gal-bind_lectin: domain 1 of 2, from 18 to 149: score 122.0, E = 3.6e-34
*->pglvalnlglkpGktltVkGtVapknakrFavNlgkgskE..Endlv
p+ ++l++ + pG+tl+VkG+ + ++rF +Nl++++ + + nd+
W09H1.6b 18 PYRSVLQEKFEPGQTLIVKGSTID-ESQRFTINLHSKTADfsGNDVP 63
LHfNPRFneaHGDqntvVcNSkenGDNeWGtEqReaafPFqaGqpFeisI
LH+ +RF+ ++ +V+NS+ nG eWG+E+R +++P ++G+ F+i+I
W09H1.6b 64 LHVSVRFD-EG----KIVLNSFSNG--EWGKEER-KSNPIKKGDSFDIRI 105
sveedkfkVkvndghef.eFphRlkleavqyLgikGDikltsikf<-*
+ + d+f+++v ef+++ hRl+l++++ L+i+GD l + +
W09H1.6b 106 RAHDDRFQIIVD-HKEFkDYEHRLPLSSISHLSIDGDLYLNHVHW 149
Gal-bind_lectin: domain 2 of 2, from 157 to 284: score 232.9, E = 2.2e-66
*->pglvalnlglkpGktltVkGtVapknakrFavNlgkgskEEndlvLH
p++++l++gl++Gk+l V+GtV++ +akrF+vNl++++ +d+++H
W09H1.6b 157 PYESGLANGLPVGKSLLVFGTVEK-KAKRFHVNLLRKN---GDISFH 199
fNPRFneaHGDqntvVcNSkenGDNeWGtEqReaafPFqaGqpFeisIsv
fNPRF+e+H v++NS++++ eWG+E+Re+++PF++G++F+++I++
W09H1.6b 200 FNPRFDEKH-----VIRNSLAAN--EWGNEEREGKNPFEKGVGFDLVIQN 242
eedkfkVkvndghefeFphRlkleavqyLgikGDikltsikf<-*
ee++f+V+vn++++++F+hR++++++++L+i+GDi+l++i++
W09H1.6b 243 EEYAFQVFVNGERYISFAHRADPHDIAGLQISGDIELSGIQI 284