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: hmm5/seq/Y48C3A.7
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: Y48C3A.7 CE28132 locus:mac-1 status:Partially_confirmed TR:Q9NAG4 protein_id:CAB55106.2
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
AAA ATPase family associated with various 488.7 2.1e-146 2
RNA_helicase RNA helicase 19.3 0.00025 2
Parvo_NS1 Parvovirus non-structural protein NS1 12.4 0.0073 2
SKI Shikimate kinase 10.8 0.067 1
DUF668 Protein of unknown function (DUF668) 10.6 0.019 1
IstB IstB-like ATP binding protein 8.4 0.22 1
Viral_helicase1 Viral (Superfamily 1) RNA helicase 7.0 0.34 1
NB-ARC NB-ARC domain 6.3 0.64 2
SRP-alpha_N Signal recognition particle, alpha su 5.8 0.84 1
ABC_tran ABC transporter 5.4 2.1 2
URO-D Uroporphyrinogen decarboxylase (URO-D 3.8 3.4 1
ADK Adenylate kinase 3.5 9.6 2
G-alpha G-protein alpha subunit 2.7 9.6 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
SRP-alpha_N 1/1 111 133 .. 288 310 .] 5.8 0.84
G-alpha 1/1 121 147 .. 1 27 [. 2.7 9.6
ADK 1/2 244 252 .. 1 9 [. 4.9 3.8
Viral_helicase1 1/1 242 253 .. 1 12 [. 7.0 0.34
NB-ARC 1/2 244 259 .. 46 61 .. 2.5 8.1
Parvo_NS1 1/2 245 260 .. 175 194 .. 2.7 5.6
ABC_tran 1/2 244 263 .. 6 25 .. 0.2 67
RNA_helicase 1/2 243 267 .. 127 151 .. 11.3 0.041
AAA 1/2 241 453 .. 1 216 [] 216.0 4e-64
DUF668 1/1 469 485 .. 68 84 .. 10.6 0.019
URO-D 1/1 526 545 .. 339 358 .] 3.8 3.4
ADK 2/2 573 581 .. 1 9 [. 3.9 7.1
IstB 1/1 569 587 .. 50 68 .. 8.4 0.22
RNA_helicase 2/2 571 590 .. 126 145 .. 8.4 0.27
ABC_tran 2/2 571 592 .. 4 25 .. 5.1 2.6
Parvo_NS1 2/2 571 593 .. 172 198 .. 10.5 0.027
NB-ARC 2/2 575 596 .. 48 72 .. 3.9 3.2
SKI 1/1 569 597 .. 1 29 [. 10.8 0.067
AAA 2/2 570 759 .. 1 216 [] 280.4 1.4e-83
Alignments of top-scoring domains:
SRP-alpha_N: domain 1 of 1, from 111 to 133: score 5.8, E = 0.84
*->GgktkdgDvvsedeSDdEddden<-*
G++ ++vv++d+S+dE ++ +
Y48C3A.7 111 GQAAAAKEVVLSDDSEDERAARQ 133
G-alpha: domain 1 of 1, from 121 to 147: score 2.7, E = 9.6
*->seeekeqakrnkeIekqLkqekkkakr<-*
+++ e +++ +++ekq++ k+ +++
Y48C3A.7 121 LSDDSEDERAARQLEKQIESLKTNRAN 147
ADK: domain 1 of 2, from 244 to 252: score 4.9, E = 3.8
RF xxxxxxxxx
*->llGpPGaGK<-*
+ GpPG+GK
Y48C3A.7 244 VHGPPGCGK 252
Viral_helicase1: domain 1 of 1, from 242 to 253: score 7.0, E = 0.34
*->ivVhGVpGcGKS<-*
+VhG pGcGK+
Y48C3A.7 242 FIVHGPPGCGKT 253
NB-ARC: domain 1 of 2, from 244 to 259: score 2.5, E = 8.1
*->ivGMGGvGKTTLAkqi<-*
++G++G+GKT A+++
Y48C3A.7 244 VHGPPGCGKTMFAQAV 259
Parvo_NS1: domain 1 of 2, from 245 to 260: score 2.7, E = 5.6
*->yGpastGKtnlysGAAqaIa<-*
+Gp+++GKt+ Aqa a
Y48C3A.7 245 HGPPGCGKTMF----AQAVA 260
ABC_tran: domain 1 of 2, from 244 to 263: score 0.2, E = 67
*->lvGpNGaGKSTLLklisGll<-*
+ Gp G+GK+ + + +G+l
Y48C3A.7 244 VHGPPGCGKTMFAQAVAGEL 263
RNA_helicase: domain 1 of 2, from 243 to 267: score 11.3, E = 0.041
*->ylhGaPGqGKSllanvlaralskhl<-*
+ hG+PG+GK+++a+++a l+ +
Y48C3A.7 243 IVHGPPGCGKTMFAQAVAGELAIPM 267
AAA: domain 1 of 2, from 241 to 453: score 216.0, E = 4e-64
*->gvLLyGPPGtGKTlLAkavAkelgvpfisisgselvskyvGesekrv
g + +GPPG+GKT+ A avA+el +p++ + ++elvs+ Ge+e+++
Y48C3A.7 241 GFIVHGPPGCGKTMFAQAVAGELAIPMLQLAATELVSGVSGETEEKI 287
ralfelArkslkkaaPspiiFIDEiDalapkRgd.egdvssrvvnqLlte
r lf+ A+++ +P+ i+++D+iDa+ap+R+ ++++++rvv+qL
Y48C3A.7 288 RRLFDTAKQN----SPC-ILILDDIDAIAPRRETaQREMERRVVSQLCSS 332
mDLe.......................rigfekhylrgsdvvDlsnvivi
+D e ++++++ +++ + +++++ ig ++ +v+vi
Y48C3A.7 333 LD-ElvlpprekplkdqltfgddgsvaIIGDSPTAA-------GAGVLVI 374
aaTNrpdlldpaLlrpGRfdreievplPdeeerleIlkihlkkmplalcq
++T rpd d L+r GRf+ ei ++ Pde +r++Il++++k ++la
Y48C3A.7 375 GTTSRPDAVDGGLRRAGRFENEISLGIPDETAREKILEKICK-VNLA--- 420
erselakdvdldelakelArrtpgfsgadlaalcreAalralr<-*
+dv l++ +A++tpg++gadl+al+reAa a++
Y48C3A.7 421 ------GDVTLKQ----IAKLTPGYVGADLQALIREAAKVAID 453
DUF668: domain 1 of 1, from 469 to 485: score 10.6, E = 0.019
*->llAtewKaamekiLeWL<-*
l+++++K+++++ L WL
Y48C3A.7 469 LTVEQIKEELDRVLAWL 485
[output cut off at A = 10 top alignments]