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/R153.1a
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: R153.1a CE02038 cyclic nucleotide phosphodiesterase status:Partially_confirmed SW:Q22000 protein_id:AAA68292.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
PDEase 3'5'-cyclic nucleotide phosphodiesterase 455.2 2.7e-158 1
NHL NHL repeat 7.6 2.5 1
S_100 S-100/ICaBP type calcium binding domain 5.9 5.8 1
ATP-synt ATP synthase 2.2 9.6 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
ATP-synt 1/1 40 56 .. 1 17 [. 2.2 9.6
S_100 1/1 117 131 .. 30 44 .] 5.9 5.8
NHL 1/1 416 432 .. 13 29 .] 7.6 2.5
PDEase 1/1 281 525 .. 1 279 [] 455.2 2.7e-158
Alignments of top-scoring domains:
ATP-synt: domain 1 of 1, from 40 to 56: score 2.2, E = 9.6
*->AsLkeIrdRIkSVKnTq<-*
AsL+ +r+ + S++n+q
R153.1a 40 ASLRNVRSNLISITNIQ 56
S_100: domain 1 of 1, from 117 to 131: score 5.9, E = 5.8
*->LKeLlekELpnfLkn<-*
++++l+kEL++f +
R153.1a 117 FRKMLNKELSHFAES 131
NHL: domain 1 of 1, from 416 to 432: score 7.6, E = 2.5
*->dGrvyVaDsenhrvqvf<-*
+++v+V D +n ++qv+
R153.1a 416 GNNVIVLDKYNDKIQVL 432
PDEase: domain 1 of 1, from 281 to 525: score 455.2, E = 2.7e-158
*->YHNwiHAfdVtQtthlLlltlaleryLtdlEvLalvfAAaiHDvDHr
YHN+iHA+dV Q++h+Ll + l+ ++tdlEvLa++fA a HDvDH+
R153.1a 281 YHNFIHAADVAQSMHVLLMSPVLTEVFTDLEVLAAIFAGAVHDVDHP 327
GTnNsFqinsLqkSeLAlLYndegSVLEnHHlaqafkLLqdEecnIfqNL
G +N+++ins ++eLA++Ynde SVLE HHla+afkLLqd +c+ + NL
R153.1a 328 GFTNQYLINS--NNELAIMYNDE-SVLEQHHLAVAFKLLQDSNCDFLANL 374
skkdfrtlrdlvieaILATDmslHlqklkdlktmveqkkvyetgvEWtqY
s+k++ + r+ vi+++LATDms+H+++l+dlktmve kkv++ v
R153.1a 375 SRKQRLQFRKIVIDMVLATDMSKHMSLLADLKTMVEAKKVAGNNV----I 420
lldnythkilllsllmtAADLSnpTKpwslskRwAelimeEFFeQGDlEr
ld y +ki++l+++++ ADLSnpTKp +l+ +w +imeE ++QGD+E+
R153.1a 421 VLDKYNDKIQVLQSMIHLADLSNPTKPIELYQQWNQRIMEEYWRQGDKEK 470
elGldrpspmcDRtsAayvpksQvgFidfIvePvfklLadvvekGRttse
elGl +spmcDR + +++ksQvgFid+Iv+P+++++ad+v +
R153.1a 471 ELGLE-ISPMCDRGN-VTIEKSQVGFIDYIVHPLYETWADLVYP------ 512
aiDanhLCWvaLDeevRnddiaplldriednR<-*
d++ +ld +e+nR
R153.1a 513 -------------------DAQNILDQLEENR 525