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: hmm4/seq/M116.5
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: M116.5 CE25957 status:Partially_confirmed TR:Q965L0 protein_id:AAK71387.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
CH Calponin homology (CH) domain 273.2 3.4e-78 4
GAS2 Growth-Arrest-Specific Protein 2 Domain 146.1 3.8e-42 2
CENP-B_N CENP-B N-terminal DNA-binding domain 5.3 6.4 1
Guanylin Guanylin precursor 2.2 2.8 1
V1R Vomeronasal organ pheromone receptor family, 2.2 8 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CH 1/4 40 156 .. 1 106 [] 83.2 1.3e-21
CH 2/4 193 299 .. 1 106 [] 112.5 7.8e-30
CH 3/4 399 466 .. 48 106 .] 31.2 3.3e-07
CH 4/4 587 652 .. 48 106 .] 49.6 2.6e-12
Guanylin 1/1 678 693 .. 92 107 .. 2.2 2.8
GAS2 1/2 791 821 .. 1 33 [. 44.8 5.6e-12
V1R 1/1 891 899 .. 257 265 .] 2.2 8
GAS2 2/2 885 926 .. 34 75 .] 99.1 3.7e-28
CENP-B_N 1/1 1051 1062 .. 1 12 [. 5.3 6.4
Alignments of top-scoring domains:
CH: domain 1 of 4, from 40 to 156: score 83.2, E = 1.3e-21
*->sqekellrWinsvlagypgkkstlPippvpvtnfftdwlkDGlaLca
q+k++++W+n++l + ++ ++ ++f d++ DG aL+a
M116.5 40 MQKKTFTKWVNKHLTK-TDH---------TIVDLFVDLR-DGFALIA 75
Llnklr............................eNlnlaldaaekeG.l
Ll++l +++ +++++ ++ ++ ++ + + + N++ aek+ ++
M116.5 76 LLEALTgeriqrengytrfhriqnvqycldflkmKNVSFLERQAEKKTnI 125
givlvlldpeDiveaegnpdPDeklvltylwqLirrfq<-*
++ lv++ peDive gn kl+l+++w++i+ fq
M116.5 126 IK-LVNIRPEDIVE--GNG----KLTLGLIWTIILNFQ 156
CH: domain 2 of 4, from 193 to 299: score 112.5, E = 7.8e-30
*->sqekellrWinsvlagypgkkstlPippvpvtnfftdwlkDGlaLca
s++++ll+W+++++a+yp+ +v nf+ +w+ DGla++a
M116.5 193 SARDALLQWAKRTTAEYPR---------GNVINFSNSWR-DGLAFNA 229
Llnklr.................eNlnlaldaaekeGlgivlvlldpeDi
+++++r++ + ++ ++++ +++e+ln+a+ aa++e +g+ +lld+eD+
M116.5 230 IIHRYRssaidwnkinsdsvsntERLNIAFAAADRE-FGV-ERLLDAEDV 277
veaegnpdPDeklvltylwqLirrfq<-*
n+ PDek+++ty+++L+++ +
M116.5 278 D---TNN-PDEKSIITYVSSLYNALP 299
CH: domain 3 of 4, from 399 to 466: score 31.2, E = 3.3e-07
*->Llnklr..............eNlnlaldaaekeGlgivlvlldpeDi
L+ ++ + ++ ++++ +++e lnla+ aa+ e +g +lld+eD+
M116.5 399 LMKHVEdaysdinsdsvpntESLNLAFAAADPE-FGA-GRLLDAEDV 443
veaegnpdPDeklvltylwqLirrfq<-*
+n+ PDek v+ty+++L+++ q
M116.5 444 DT--NNN-PDEKPVITYVSSLKNAPQ 466
CH: domain 4 of 4, from 587 to 652: score 49.6, E = 2.6e-12
*->Llnklr............eNlnlaldaaekeGlgivlvlldpeDive
L++++ ++ ++++ +++e+ln+al aa+ e +g +lld+eD+
M116.5 587 LIDRCNasninsdsvsntERLNIALAAADAE-FGA-GRLLDAEDVDT 631
aegnpdPDeklvltylwqLirrfq<-*
+n+ PDek+++ty+++L ++f+
M116.5 632 --NNN-PDEKSIITYVSSLNNAFP 652
Guanylin: domain 1 of 1, from 678 to 693: score 2.2, E = 2.8
*->eAasIlqaLrtIaqDd<-*
+A I++aL++I++D+
M116.5 678 DADHIFKALEEISKDE 693
GAS2: domain 1 of 2, from 791 to 821: score 44.8, E = 5.6e-12
*->dklddeVkrivedPPCkCpnkFqvekvSEGkYR<-*
+k++ e+++++e+ C+C++++q+ek+SE+ YR
M116.5 791 EKITNEIELLKET--CSCCTPYQIEKISENYYR 821
V1R: domain 1 of 1, from 891 to 899: score 2.2, E = 8
*->KrIikflrS<-*
Kr++++lrS
M116.5 891 KRMVRILRS 899
GAS2: domain 2 of 2, from 885 to 926: score 99.1, E = 3.7e-28
*->fGDsqiLrlVRiLRStVMVRVGGGWetLdefLlKhDPCRakg<-*
fGD++i+r+VRiLRStVMVRVGGGWe L+e L+KhDPCR kg
M116.5 885 FGDTHIKRMVRILRSTVMVRVGGGWEPLEESLQKHDPCRTKG 926
CENP-B_N: domain 1 of 1, from 1051 to 1062: score 5.3, E = 6.4
*->rkkRrsLTLeEK<-*
r kR sL+L+EK
M116.5 1051 RQKRISLSLKEK 1062