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/C46E1.2
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: C46E1.2 CE19736 locus:gcy-36 GUANYLATE CYCLASE SOLUBLE, BETA-2 CHAIN (EC 4.6.1.2) (76 KD CHAIN) status:Partially_confirmed TR:Q9XTE0 p
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
guanylate_cyc Adenylate and Guanylate cyclase catalyt 229.2 5.9e-65 1
VlpA_repeat Mycoplasma hyorhinis VlpA repeat 9.8 0.36 1
Hpr_kinase HPr Serine kinase 2.8 6.8 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
VlpA_repeat 1/1 202 216 .. 1 13 [] 9.8 0.36
Hpr_kinase 1/1 610 627 .. 43 60 .. 2.8 6.8
guanylate_cyc 1/1 451 635 .. 1 226 [] 229.2 5.9e-65
Alignments of top-scoring domains:
VlpA_repeat: domain 1 of 1, from 202 to 216: score 9.8, E = 0.36
*->KtEntqqs..EAlGt<-*
Kt nt qs++EAlGt
C46E1.2 202 KTLNTDQSneEALGT 216
Hpr_kinase: domain 1 of 1, from 610 to 627: score 2.8, E = 6.8
*->GlFtyvpigrVqlLGKrE<-*
G F + p+grVq+ GK+E
C46E1.2 610 GRFEFEPRGRVQIKGKGE 627
guanylate_cyc: domain 1 of 1, from 451 to 635: score 229.2, E = 5.9e-65
*->vyaerydeVtilFaDIvgFTalssrhssNngPeevvrlLNelftrFD
v+a +y+e t++F+D++ F+++ ++ P +v lLNelft+FD
C46E1.2 451 VEAREYEEATVMFTDVPTFQQIVPLCT----PKDIVHLLNELFTKFD 493
elvdehkfgvyKiKTiGDaYMaasGlppaakakqtssKlaeeldedsaaH
+l++ yK++T+GD+YM + G+p +H
C46E1.2 494 RLIGIQ--KAYKVETVGDSYMSVGGIPDL-----------------VDDH 524
aakladfALamvealqevnvghteggpPlvLyrglqlRiGIHtGpVvaGv
+ +++ AL+mv +++v+ + t+ +p l +R GIH+GpVvaGv
C46E1.2 525 CEVICHLALGMVMEARTVCDPITN-TP-------LHIRAGIHSGPVVAGV 566
iGgaskrprYdvwGdTVNvASRmEStGvpgekIhVsestyrlLkelesfq
+G +k+prY+++GdTVN+ SRmES++ g +Ih se ++ + + ++f+
C46E1.2 567 VG--AKMPRYCLFGDTVNTSSRMESHSPIG-RIHCSENAKKCAESTGRFE 613
fefaleprgevsvkGkgkpGaietYfLkg<-*
fe+ rg v kGkg++ tYfL +
C46E1.2 614 FEP----RGRVQIKGKGEM---NTYFLLR 635