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: hmm5c/seq/ZK896.8
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: ZK896.8 CE15463 locus:gcy-18 guanylate cyclase status:Partially_confirmed TR:O02344 protein_id:CAB05325.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
guanylate_cyc Adenylate and Guanylate cyclase catal 253.5 2.9e-72 1
ANF_receptor Receptor family ligand binding region 116.8 2.9e-31 1
pkinase Protein kinase domain 43.9 4.3e-11 1
Herpes_glycop_D Herpesvirus glycoprotein D 5.8 1.4 1
2-Hacid_DH_C D-isomer specific 2-hydroxyacid dehyd 5.8 0.052 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
2-Hacid_DH_C 1/1 169 182 .. 1 14 [. 5.8 0.052
ANF_receptor 1/1 47 442 .. 1 456 [] 116.8 2.9e-31
Herpes_glycop_D 1/1 499 533 .. 374 409 .] 5.8 1.4
pkinase 1/1 683 846 .. 131 294 .] 43.9 4.3e-11
guanylate_cyc 1/1 908 1095 .. 1 226 [] 253.5 2.9e-72
Alignments of top-scoring domains:
2-Hacid_DH_C: domain 1 of 1, from 169 to 182: score 5.8, E = 0.052
*->npvsVAEhvigliL<-*
n++s+AE +++l+L
ZK896.8 169 NTNSIAEATVALLL 182
ANF_receptor: domain 1 of 1, from 47 to 442: score 116.8, E = 2.9e-31
*->avrgfpllvavgfafelinvgpallpgialgaaildvnslpnyalnl
++ ++ +++ ++a vg + p+ a ++ i+ n ++
ZK896.8 47 KEHKRVIKIGHIGA-----VG--VMPNDARILNISKENLIEEGLVGD 86
glkflrgidydsacdcsfavavavdlllekgvdewvvAiIGPscssaaia
+ f+ + + ++ ++vava++l++ ++v A+IGP+c a
ZK896.8 87 DIEFEIVSRQAC-SESFEGVAVAAELYHVHQVR----AFIGPYCAAELEA 131
varlanewniPlisygatapelsdketrfptflRTlPsdtkq.gealvdl
v+ +a+ wniP+isy+++ +++sd+ + + t++R+ t +ea v+l
ZK896.8 132 VTKMATFWNIPIISYSSVPNAVSDR-SVYKTLARVSSKNTNSiAEATVAL 180
lkhfgWtwvalvysdddygeegvaaliealeerggicialkvsvpeidea
l h++W +va+ ++ e v+ ++e ++ r+g+ i kv++ e + a
ZK896.8 181 LLHYKWLKVAIATNTGSTAFERVSIFEEIMH-REGVTIVKKVMFDENTDA 229
deedyrallkrikrrkstaRvivvcgsSVLiGgedlrqlleqavrlGmtg
+e + + l ++ + aR+i++ +s + e ++++++ + m+
ZK896.8 230 NEMMNSGQLGDLAAN---ARIIICLFS---STKELSKEFMQATYTMRMNN 273
ggyvfisedwlsssldldpdlellapgaitlgtdgsdnvlgfayyavlvi
+y +i+ + s + dl+p +ga+ g + + v ++ + a v
ZK896.8 274 AEYAYIIPWLQSGTKDLTPW-----IGAD--G-EMLQRVKDHYANAIIVD 315
tldspdnpefnefweelfecaltpfnfnksesdrkCtgvecnsqalledt
+ + d++ + f+e+ +++ ++++ + ++
ZK896.8 316 DVNGFDDSVVSSFVEKIEKHGMQKSDI---------DTNING-------- 348
svnpyagavyDaVylyAhAlhnellrngggctlcaasdwrdGsqllrylr
++++D+ lyA+A+++ l + + + + +G +++++r
ZK896.8 349 -----YLHLFDSLKLYALAIRKVLNETDN------EAYVTNGQFIWNRMR 387
nvtFeGltGp.............eVqfDenRGDrnpdYdllnlqdsetgk
+ FeG+ +++++++++ + ++V +D+ Dr+p ++ +++ ++
ZK896.8 388 RMKFEGVVSRssseenkdagaigTVLMDDV-ADRAPIFSAFYI-----SP 431
levvgSGmPitfyydgvgklvvv<-*
+d+v k+v++
ZK896.8 432 ------------NRDKVMKMVNM 442
Herpes_glycop_D: domain 1 of 1, from 499 to 533: score 5.8, E = 1.4
*->ivvlalligaiivgvvyyilrRrRkeakkyvrlpps<-*
i+v++l+i++i++ vv++ l+R++ k + + p
ZK896.8 499 ISVGSLIILLILISVVFFFLFRYCE-NKQLEKMPWR 533
pkinase: domain 1 of 1, from 683 to 846: score 43.9, E = 4.3e-11
*->ILldendgtvKiaDFGLArlle................sssklttfv
+L+d+n vK++D G+A le+ ++++ + +++++ + +
ZK896.8 683 CLIDRN-WMVKLSDYGIANPLErwekqgaieiaaakdsDDKSQASQA 728
GTpwYmmAPEvi...................legrgysskvDvWSlGviL
++ Ym APE++++++++++++ +++ ++++++r + D++S+G+++
ZK896.8 729 TSIIYM-APELLknretnkrrgmdqswvkqsMLRRQAG---DIYSFGMVM 774
yElltggplfpgadlpaftggdevdqliifvlklPfsdelpktridplee
yE+l++ lPf+d + +++e
ZK896.8 775 YEILFR--------------------------SLPFRD------NTNISE 792
lfriikrp....glrlplpsncSeelkdLlkkcLnkDPskRpGsatakei
l++ + + +++ +++ ++ + ++l Ll++c+ +P+ Rp +++ +
ZK896.8 793 LVDYLADGsktvSPEIQNQMGLHPDLNALLRDCWSENPEIRP---SIRRV 839
.lnhpwf<-*
+ln +
ZK896.8 840 rLNTEMV 846
guanylate_cyc: domain 1 of 1, from 908 to 1095: score 253.5, E = 2.9e-72
*->vyaerydeVtilFaDIvgFTalssrhssNngPeevvrlLNelftrFD
v ++ y + tilF+DIvgFT+++s ++ P evv++LN l+t FD
ZK896.8 908 VAPKLYSSATILFSDIVGFTTICSGST----PLEVVNMLNGLYTGFD 950
elvdehkfgvyKiKTiGDaYMaasGlppaakakqtssKlaeeldedsaaH
e + ++ yK++TiGDaYM++sG+p+++ ++ H
ZK896.8 951 ECITRN--KSYKVETIGDAYMVVSGIPEEN----------------EYNH 982
aakladfALamvealqevnvghteggpPlvLyrglqlRiGIHtGpVvaGv
+ +a AL+m+++l+ +++h++ + ++ R G HtG+V aGv
ZK896.8 983 SRNIANTALDMRQYLTGYQIPHRPTHR-------VRCRWGFHTGSVAAGV 1025
iGgaskrprYdvwGdTVNvASRmEStGvpgekIhVsestyrlLkelesfq
+G + prY+++GdTVNv SRmEStG pg I++se+++++ +++ ++
ZK896.8 1026 VG--LTCPRYCLFGDTVNVSSRMESTGTPG-MIQMSEEAHMHIRAHH-PV 1071
fefaleprgevsvkGkgkpGaietYfLkg<-*
f++ ++rgev vkGk G+ t++L+
ZK896.8 1072 FTT--TERGEVQVKGK---GTCRTFWLED 1095