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: hmm1/seq/C05A9.1
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: C05A9.1 CE05221 locus:pgp-5 p-glycoprotein status:Partially_confirmed TR:Q17645 protein_id:CAA94202.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
ABC_tran ABC transporter 412.2 4.8e-120 2
ABC_membrane ABC transporter transmembrane region 188.2 1.3e-52 2
SKI Shikimate kinase 22.3 3.2e-05 2
PRK Phosphoribulokinase / Uridine kinase fa 11.3 0.035 2
AIG1 AIG1 family 7.7 0.56 2
Auxin_eff Auxin Efflux Carrier 6.8 0.26 1
CoaE Dephospho-CoA kinase 6.6 0.58 2
ATP-bind Conserved hypothetical ATP binding prot 6.4 0.74 2
GTP_CDC Cell division protein 6.0 1 2
APS_kinase Adenylylsulfate kinase 5.9 1.2 2
AAA ATPase family associated with various c 5.3 1.4 2
ATP-synt_B ATP synthase B/B' CF(0) 4.7 4.4 1
DNA_primase_S DNA primase small subunit 4.7 1.9 1
Art_glycop Arterivirus glycoprotein 2.6 5.4 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
ATP-synt_B 1/1 316 329 .. 1 14 [. 4.7 4.4
ABC_membrane 1/2 68 358 .. 1 285 [] 89.5 7e-23
AIG1 1/2 434 448 .. 1 15 [. 1.4 30
CoaE 1/2 434 449 .. 1 16 [. 6.4 0.66
APS_kinase 1/2 432 449 .. 1 18 [. 1.5 23
GTP_CDC 1/2 424 450 .. 1 21 [. 5.0 2
SKI 1/2 434 453 .. 1 20 [. 12.0 0.029
PRK 1/2 435 453 .. 1 19 [. 10.9 0.047
ATP-bind 1/2 435 458 .. 1 24 [. 4.8 2.2
AAA 1/2 435 468 .. 1 32 [. 3.2 6.3
ABC_tran 1/2 433 617 .. 1 199 [] 222.2 3.5e-63
Art_glycop 1/1 834 849 .. 168 183 .] 2.6 5.4
Auxin_eff 1/1 831 855 .. 497 521 .] 6.8 0.26
ABC_membrane 2/2 715 993 .. 1 285 [] 106.6 4.7e-28
AIG1 2/2 1065 1079 .. 1 15 [. 8.7 0.31
SKI 2/2 1065 1080 .. 1 16 [. 11.1 0.053
CoaE 2/2 1065 1080 .. 1 16 [. 1.8 17
GTP_CDC 2/2 1066 1081 .. 6 21 .. 3.5 5.3
AAA 2/2 1066 1084 .. 1 19 [. 1.5 21
PRK 2/2 1068 1085 .. 1 20 [. 2.1 14
APS_kinase 2/2 1063 1086 .. 1 24 [. 10.2 0.066
ATP-bind 2/2 1066 1091 .. 1 26 [. 2.0 15
DNA_primase_S 1/1 1082 1117 .. 315 350 .] 4.7 1.9
ABC_tran 2/2 1064 1249 .. 1 199 [] 191.9 2e-54
Alignments of top-scoring domains:
ATP-synt_B: domain 1 of 1, from 316 to 329: score 4.7, E = 4.4
*->tLetnlInfliLvg<-*
++tn++nf+iL++
C05A9.1 316 YFFTNALNFVILYF 329
ABC_membrane: domain 1 of 2, from 68 to 358: score 89.5, E = 7e-23
*->lliailllilagatalvtfplllgrlldsgfplsdgnddh.......
++++ +++il g ++ ++ lg ++ ++++++ + ++ +
C05A9.1 68 FFLGVVFSILTGMCQPF-ESYTLGETSQVLVKVTNAINNKtidpvdl 113
.........arsslislailslfavfvliglllqgsfyllagerlgqrlr
+ + ++ ++ + l++l++fa f++ +l ++ ++ g +++ ++r
C05A9.1 114 ahayklfesDMNRVVLLFFLVGFAYFTFGFLQFSIMKFV--GDNTAYNVR 161
krlfrallrqilglfdsffdtnsvGeltsRltnDvekirdglgeklgllf
++ +++llr+++++fd +s+G l l + +e+ r+ ++ek+ l++
C05A9.1 162 RQYISRLLRKDISYFDG----MSTGHLSIVLNDNMERFREVFNEKIALII 207
qslatvvgglivmfyyswkLtLillailpllillsavlakklrklsrkeq
l+ +v+g i++fy w+L++ ++ i+ls l ++ + +++++q
C05A9.1 208 ALLTDFVIGTILAFYTDWRLACYGT-VFSFGIVLSGLLDSWGKMKNNEKQ 256
kayak.agsvaeEslsgirTVkafgrEeyelerfdkaledaekagikkai
+ + +ags+a ++l+ +TV +++++++e+er++++l + ek ++ +a+
C05A9.1 257 NEHISnAGSIAFQALGCYKTVSSLNGQQQEVERYTEELKNGEKYALNRAF 306
iagllfgitqlisylsyalalwfGgylvasvisgglsvgtlfaflslgnq
+ +l+ ++ +++++++++++l+fG+ + + g + +g ++ l++ ++
C05A9.1 307 VFSLSRSADYFFTNALNFVILYFGANM---IYEGTIEPGVVVRILYYILF 353
ligpl<-*
++ l
C05A9.1 354 GSYCL 358
AIG1: domain 1 of 2, from 434 to 448: score 1.4, E = 30
*->rnIVLVGRTGNGKSA<-*
+ I LVG +G GKS+
C05A9.1 434 ECIALVGASGSGKST 448
CoaE: domain 1 of 2, from 434 to 449: score 6.4, E = 0.66
*->lvvGLTGGIGSGKSTv<-*
+++L G +GSGKSTv
C05A9.1 434 ECIALVGASGSGKSTV 449
APS_kinase: domain 1 of 2, from 432 to 449: score 1.5, E = 23
*->rGctvWfTGLSGsGKSTI<-*
G ++ + G SGsGKST+
C05A9.1 432 GGECIALVGASGSGKSTV 449
GTP_CDC: domain 1 of 2, from 424 to 450: score 5.0, E = 2
*->kGfdFT......LMVVGeSGLGKtTlI<-*
k+++F+ ++++ + +VG+SG GK+T I
C05A9.1 424 KEISFNvqggecIALVGASGSGKSTVI 450
SKI: domain 1 of 2, from 434 to 453: score 12.0, E = 0.029
*->rnifLIGmmGaGKsTvGrlL<-*
++i L+G G+GKsTv +lL
C05A9.1 434 ECIALVGASGSGKSTVIQLL 453
PRK: domain 1 of 2, from 435 to 453: score 10.9, E = 0.047
*->vIGvaGsSGaGKTtvarri<-*
I++ G+SG+GK+tv + +
C05A9.1 435 CIALVGASGSGKSTVIQLL 453
ATP-bind: domain 1 of 2, from 435 to 458: score 4.8, E = 2.2
*->ailViGpAGSGKSTfcgalyqhle<-*
+i +G GSGKST + l +++
C05A9.1 435 CIALVGASGSGKSTVIQLLLHYYN 458
AAA: domain 1 of 2, from 435 to 468: score 3.2, E = 6.3
*->gvLLyGPPGtGKTlLAkavAkelg..vpfisisg<-*
++ L+G G+GK+++ + + ++ + +isi+g
C05A9.1 435 CIALVGASGSGKSTVIQLLLHYYNidSGRISIDG 468
[output cut off at A = 10 top alignments]