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: hmm5a/seq/ZK524.2a
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: ZK524.2a CE15371 locus:unc-13 UNC-13 protein status:Partially_confirmed SW:P27715 protein_id:CAB07173.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
C2 C2 domain 188.0 8.6e-53 3
DAG_PE-bind Phorbol esters/diacylglycerol binding 86.8 4.5e-22 1
VPS9 Vacuolar sorting protein 9 (VPS9) dom 8.2 0.95 1
GerE Bacterial regulatory proteins, luxR f 6.3 4 1
CARD Caspase recruitment domain 4.5 8.3 1
Prothymosin Prothymosin/parathymosin family 4.1 4.1 1
LexA_DNA_bind LexA DNA binding domain 3.8 5.5 1
Glyco_hydro_3_C Glycosyl hydrolase family 3 C termina 3.3 6.4 1
DUF615 Protein of unknown function (DUF615) 2.9 6.2 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
C2 1/3 21 96 .. 1 88 [] 23.8 2e-05
Glyco_hydro_3_C 1/1 289 316 .. 256 285 .] 3.3 6.4
DAG_PE-bind 1/1 693 742 .. 1 51 [] 86.8 4.5e-22
DUF615 1/1 755 791 .. 142 182 .] 2.9 6.2
Prothymosin 1/1 883 901 .. 94 111 .] 4.1 4.1
C2 2/3 817 908 .. 1 88 [] 110.4 2.1e-30
CARD 1/1 989 1005 .. 74 90 .] 4.5 8.3
VPS9 1/1 1310 1325 .. 103 118 .] 8.2 0.95
GerE 1/1 1551 1562 .. 1 12 [. 6.3 4
LexA_DNA_bind 1/1 1551 1563 .. 1 13 [. 3.8 5.5
C2 3/3 1648 1738 .. 1 88 [] 59.3 1.2e-15
Alignments of top-scoring domains:
C2: domain 1 of 3, from 21 to 96: score 23.8, E = 2e-05
*->LrVtvieAknLppkDkkngksDPYVkvslggdpkdqkkkTkvikktl
L +t+++A+ ++D + + YV+v+l +k T ++++
ZK524.2a 21 LCITIKKARLQGAVD-EF---NSYVTVKLQT----VKSTTVAVRGN- 58
NPvWNEtFtFevpppdlqetLeitVyDydrfgrddflGevs<-*
P W+++F Fe + +++ +++++ ++ ++ d+++G +
ZK524.2a 59 LPCWEQEFIFE-TNRPDDG-MVLELWAKGVLW-DKLIGVHY 96
Glyco_hydro_3_C: domain 1 of 1, from 289 to 316: score 3.3, E = 6.4
*->dedgfpqddfYeegifvdYryFdkgnitpr<-*
d d+++q+ +Ye+ + +++r+F++ p+
ZK524.2a 289 DQDQQHQN-IYED-TVTPVREFGESAVPPA 316
DAG_PE-bind: domain 1 of 1, from 693 to 742: score 86.8, E = 4.5e-22
*->HrFkrttfyksptfCdhCgellwglakQGlkCsnCglnvHkrChekV
H+F ttf ++ptfC C+ llwgla+QGl+C++C+++vH++C e++
ZK524.2a 693 HNFATTTF-QTPTFCYECEGLLWGLARQGLRCTQCQVKVHDKCRELL 738
ptnC<-*
++C
ZK524.2a 739 SADC 742
DUF615: domain 1 of 1, from 755 to 791: score 2.9, E = 6.2
*->eADR.QkLRtLiRnAkKEkaqNKPPkafReiFQyLkelmetq<-*
eADR+Q L + iR K ++qNKP e+FQ+++++ +
ZK524.2a 755 EADRtQSLVNVIRDRMKIQEQNKP-----EVFQMIRTVFDVD 791
Prothymosin: domain 1 of 1, from 883 to 901: score 4.1, E = 4.1
*->DEDdDvetK.KqKktDEDD<-*
DED+D K +qK t E D
ZK524.2a 883 DEDNDLKSKlRQKLTRESD 901
C2: domain 2 of 3, from 817 to 908: score 110.4, E = 2.1e-30
*->LrVtvieAknLppkDkkngksDPYVkvslggdpkdqkkkTkvikktl
+++tv++A++L +kD k+gksDPYV+ ++g+ +k++T++i+ l
ZK524.2a 817 ITLTVLCAQGLIAKD-KTGKSDPYVTAQVGK----TKRRTRTIHQEL 858
NPvWNEtFtFevpppdlqetLeitVyDydrfgr...........ddflGe
NPvWNE+F+Fe ++++ ++ +++ V+D+d +++ +++ ++++ddflG+
ZK524.2a 859 NPVWNEKFHFE-CHNSTDR-IKVRVWDEDNDLKsklrqkltresDDFLGQ 906
vs<-*
++
ZK524.2a 907 TV 908
CARD: domain 1 of 1, from 989 to 1005: score 4.5, E = 8.3
*->aLretdqelladlllde<-*
+++et+qe ++++++++
ZK524.2a 989 CFQETGQEIAEEFAMRY 1005
VPS9: domain 1 of 1, from 1310 to 1325: score 8.2, E = 0.95
*->LTsLeaAleFiktlde<-*
+T+L aAl+++k +d+
ZK524.2a 1310 FTQLNAALKLLKQMDC 1325
GerE: domain 1 of 1, from 1551 to 1562: score 6.3, E = 4
*->ldsLspREreVL<-*
+sL+pR+++VL
ZK524.2a 1551 EKSLTPRQCTVL 1562
LexA_DNA_bind: domain 1 of 1, from 1551 to 1563: score 3.8, E = 5.5
*->MkpLTeRQrevLd<-*
k LT+RQ +vLd
ZK524.2a 1551 EKSLTPRQCTVLD 1563
[output cut off at A = 10 top alignments]