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.6a
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: W09H1.6a CE16576 locus:lec-1 Vertebrate galactoside-binding lectins (2 domains) status:Confirmed SW:P36573 protein_id:CAB04959.1
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 12 143 .. 1 139 [] 122.0 3.6e-34
Gal-bind_lectin 2/2 151 278 .. 1 139 [] 232.9 2.2e-66
Alignments of top-scoring domains:
Gal-bind_lectin: domain 1 of 2, from 12 to 143: score 122.0, E = 3.6e-34
*->pglvalnlglkpGktltVkGtVapknakrFavNlgkgskE..Endlv
p+ ++l++ + pG+tl+VkG+ + ++rF +Nl++++ + + nd+
W09H1.6a 12 PYRSVLQEKFEPGQTLIVKGSTID-ESQRFTINLHSKTADfsGNDVP 57
LHfNPRFneaHGDqntvVcNSkenGDNeWGtEqReaafPFqaGqpFeisI
LH+ +RF+ ++ +V+NS+ nG eWG+E+R +++P ++G+ F+i+I
W09H1.6a 58 LHVSVRFD-EG----KIVLNSFSNG--EWGKEER-KSNPIKKGDSFDIRI 99
sveedkfkVkvndghef.eFphRlkleavqyLgikGDikltsikf<-*
+ + d+f+++v ef+++ hRl+l++++ L+i+GD l + +
W09H1.6a 100 RAHDDRFQIIVD-HKEFkDYEHRLPLSSISHLSIDGDLYLNHVHW 143
Gal-bind_lectin: domain 2 of 2, from 151 to 278: score 232.9, E = 2.2e-66
*->pglvalnlglkpGktltVkGtVapknakrFavNlgkgskEEndlvLH
p++++l++gl++Gk+l V+GtV++ +akrF+vNl++++ +d+++H
W09H1.6a 151 PYESGLANGLPVGKSLLVFGTVEK-KAKRFHVNLLRKN---GDISFH 193
fNPRFneaHGDqntvVcNSkenGDNeWGtEqReaafPFqaGqpFeisIsv
fNPRF+e+H v++NS++++ eWG+E+Re+++PF++G++F+++I++
W09H1.6a 194 FNPRFDEKH-----VIRNSLAAN--EWGNEEREGKNPFEKGVGFDLVIQN 236
eedkfkVkvndghefeFphRlkleavqyLgikGDikltsikf<-*
ee++f+V+vn++++++F+hR++++++++L+i+GDi+l++i++
W09H1.6a 237 EEYAFQVFVNGERYISFAHRADPHDIAGLQISGDIELSGIQI 278