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/T21E8.1
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: T21E8.1 CE23998 locus:pgp-6 p-glycoprotein status:Partially_confirmed TR:Q22656 protein_id:CAA94220.2
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
ABC_tran ABC transporter 402.8 3.4e-117 2
ABC_membrane ABC transporter transmembrane region 197.7 1.8e-55 2
SKI Shikimate kinase 23.8 1.2e-05 2
PRK Phosphoribulokinase / Uridine kinase fa 13.7 0.0074 2
AIG1 AIG1 family 11.9 0.04 2
APS_kinase Adenylylsulfate kinase 9.2 0.12 2
DUF258 Protein of unknown function, DUF258 8.0 0.26 1
CoaE Dephospho-CoA kinase 7.6 0.29 1
ATP-bind Conserved hypothetical ATP binding prot 7.1 0.44 1
AAA ATPase family associated with various c 7.1 0.43 2
Ribosomal_L20 Ribosomal protein L20 6.3 0.79 1
Clp_N Clp amino terminal domain 5.7 6.4 1
RNA_pol_Rpb4 RNA polymerase Rpb4 4.3 5.8 1
ras Ras family 3.9 4.4 1
SLAP Bacterial surface layer protein 3.9 1.4 1
7tm_6 7tm Odorant receptor 3.1 9.6 1
GTP_CDC Cell division protein 2.9 8 2
Rad17 Rad17 cell cycle checkpoint protein 2.6 5.4 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
SLAP 1/1 72 106 .. 105 144 .. 3.9 1.4
ABC_membrane 1/2 47 314 .. 1 266 [. 89.5 6.8e-23
GTP_CDC 1/2 362 386 .. 1 19 [. 2.5 10
AIG1 1/2 372 386 .. 1 15 [. 1.4 30
APS_kinase 1/2 370 387 .. 1 18 [. 3.4 6.6
CoaE 1/1 372 391 .. 1 20 [. 7.6 0.29
SKI 1/2 372 391 .. 1 20 [. 13.4 0.012
PRK 1/2 373 391 .. 1 19 [. 11.2 0.039
ATP-bind 1/1 373 396 .. 1 24 [. 7.1 0.44
AAA 1/2 373 409 .. 1 35 [. 2.3 11
ABC_tran 1/2 371 555 .. 1 199 [] 218.4 4.3e-62
RNA_pol_Rpb4 1/1 622 635 .. 101 115 .. 4.3 5.8
Ribosomal_L20 1/1 673 688 .. 96 111 .] 6.3 0.79
Clp_N 1/1 746 778 .. 21 53 .] 5.7 6.4
7tm_6 1/1 882 914 .. 1 33 [. 3.1 9.6
ABC_membrane 2/2 657 932 .. 1 285 [] 116.2 6.2e-31
AIG1 2/2 1007 1021 .. 1 15 [. 12.6 0.026
ras 1/1 1009 1022 .. 1 15 [. 3.9 4.4
SKI 2/2 1007 1022 .. 1 16 [. 11.4 0.043
GTP_CDC 2/2 997 1023 .. 1 21 [. 2.2 13
DUF258 1/1 1011 1026 .. 33 48 .. 8.0 0.26
AAA 2/2 1008 1026 .. 1 19 [. 3.8 4.3
Rad17 1/1 1011 1026 .. 53 68 .. 2.6 5.4
PRK 2/2 1010 1027 .. 1 20 [. 3.8 4.6
APS_kinase 2/2 1005 1028 .. 1 24 [. 8.7 0.18
ABC_tran 2/2 1006 1191 .. 1 199 [] 185.5 1.4e-52
Alignments of top-scoring domains:
SLAP: domain 1 of 1, from 72 to 106: score 3.9, E = 1.4
*->atitknekdkvvkpaeLqqvtssnAgaaYtVtvtdVslNF<-*
+tit++ ++k++ pa+L+ A Y +++V + F
T21E8.1 72 VTITNAINNKTIDPADLE-----KAYEEYERGMNQVVFHF 106
ABC_membrane: domain 1 of 2, from 47 to 314: score 89.5, E = 6.8e-23
*->lliailllilagatalvtfplllgrlldsgfplsdgnddh.......
+l ++++++ +++++ f ++ +l+ + + ++++ +d + ++
T21E8.1 47 VLFSLMTGFCQPFMSYT-FGEVSQVLVTITNAINNKTIDPadlekay 92
....arsslislailslfavfvliglllqgsfyllagerlgqrlrkrlfr
++ + ++ + + l++ + +i+++lq+ + ++g +++ r+rk+ ++
T21E8.1 93 eeyeRGMNQVVFHFF-LCGCAYFIFASLQHAIMKYVGDNTTYRVRKQYIS 141
allrqilglfdsffdtnsvGeltsRltnDvekirdglgeklgllfqslat
+ll+++ ++fds s+G l + l + +e+ r+ ++ek+ l++ ++
T21E8.1 142 RLLKKDAQYFDS----VSTGHLSTVLNDNLERFREVFNEKIALIIAFVTD 187
vvgglivmfyyswkLtLillailpllillsavlakklrklsrkeqkayak
+ g +++fy w+L+ +++++ ++++ ++ + +k + k++++ya+
T21E8.1 188 FTLGTALAFYTDWRLASYGIFFSLGIAFSGFINSAGVMKTTGKQNTHYAN 237
agsvaeEslsgirTVkafgrEeyelerfdkaledaekagikkaiiagllf
ags+a ++l++ +TV +++++++e+er++++l ek gi++a+ ++++
T21E8.1 238 AGSIAFQTLGAYKTVCSLNGQKTEIERYTEELKAGEKYGIRRALMYSISR 287
gitqlisylsyalalwfGgylvasvisggl<-*
g+t+++ + +++++l++G+++ + sg l
T21E8.1 288 GVTYFFCNSLNTVVLYVGATM---IYSGTL 314
GTP_CDC: domain 1 of 2, from 362 to 386: score 2.5, E = 10
*->kGfdF......TLMVVGeSGLGKtT<-*
kG++F+ ++++ + +VG+SG GK+T
T21E8.1 362 KGISFdvqngeCIALVGASGSGKST 386
AIG1: domain 1 of 2, from 372 to 386: score 1.4, E = 30
*->rnIVLVGRTGNGKSA<-*
+ I LVG +G GKS+
T21E8.1 372 ECIALVGASGSGKST 386
APS_kinase: domain 1 of 2, from 370 to 387: score 3.4, E = 6.6
*->rGctvWfTGLSGsGKSTI<-*
G ++ + G SGsGKST+
T21E8.1 370 NGECIALVGASGSGKSTV 387
CoaE: domain 1 of 1, from 372 to 391: score 7.6, E = 0.29
*->lvvGLTGGIGSGKSTvAnlF<-*
+++L G +GSGKSTv +l+
T21E8.1 372 ECIALVGASGSGKSTVVQLL 391
SKI: domain 1 of 2, from 372 to 391: score 13.4, E = 0.012
*->rnifLIGmmGaGKsTvGrlL<-*
++i L+G G+GKsTv +lL
T21E8.1 372 ECIALVGASGSGKSTVVQLL 391
PRK: domain 1 of 2, from 373 to 391: score 11.2, E = 0.039
*->vIGvaGsSGaGKTtvarri<-*
I++ G+SG+GK+tv + +
T21E8.1 373 CIALVGASGSGKSTVVQLL 391
ATP-bind: domain 1 of 1, from 373 to 396: score 7.1, E = 0.44
*->ailViGpAGSGKSTfcgalyqhle<-*
+i +G GSGKST ++ l +++
T21E8.1 373 CIALVGASGSGKSTVVQLLLHYYN 396
AAA: domain 1 of 2, from 373 to 409: score 2.3, E = 11
*->gvLLyGPPGtGKTlLAkavAkelg..vpfisisgsel<-*
++ L+G G+GK++++ + + ++ + i i+g +l
T21E8.1 373 CIALVGASGSGKSTVVQLLLHYYNidSGNIFIDGMDL 409
[output cut off at A = 10 top alignments]