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/M03A1.3
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: M03A1.3 CE07390 locus:dao-2 status:Partially_confirmed TR:Q21476 protein_id:AAA93434.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
DB DB module 162.1 9.2e-45 1
Aa_trans Transmembrane amino acid transporter protein 5.0 1.1 2
DX DX module 3.7 5.4 1
GASA Gibberellin regulated protein 2.0 1.6 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
Aa_trans 1/2 245 274 .. 437 466 .] 2.8 5.2
Aa_trans 2/2 548 577 .. 437 466 .] 3.0 4.4
GASA 1/1 727 740 .. 59 72 .. 2.0 1.6
DX 1/1 740 763 .. 8 31 .. 3.7 5.4
DB 1/1 670 768 .. 1 124 [] 162.1 9.2e-45
Alignments of top-scoring domains:
Aa_trans: domain 1 of 2, from 245 to 274: score 2.8, E = 5.2
*->fkkildviciviglllmiygviglilslvd<-*
++ i+++++ ++g ++ ++++ ++i+ +v+
M03A1.3 245 QRIIFNWAHRCVGIVAYTLALTAIIIAAVQ 274
Aa_trans: domain 2 of 2, from 548 to 577: score 3.0, E = 4.4
*->fkkildviciviglllmiygviglilslvd<-*
++ i+++++ ++gll +++v+++++ + +
M03A1.3 548 QRSIFNWAHRIVGLLSYTFAVAAIYVAAAN 577
GASA: domain 1 of 1, from 727 to 740: score 2.0, E = 1.6
*->CskssrkklClrfC<-*
C+k++ +Cl++C
M03A1.3 727 CKKAGVSGECLKYC 740
DX: domain 1 of 1, from 740 to 763: score 3.7, E = 5.4
*->CdakekiptkyatfaFCdpkTGrv<-*
C+a + +p++y +aFC+ + +
M03A1.3 740 CTAHKGAPSNYLDYAFCTENFNEI 763
DB: domain 1 of 1, from 670 to 768: score 162.1, E = 9.2e-45
*->CCesdrkvpeeClskkCsfntysksalnmvvqkFkqlmffkqddsCp
CC++++++++eC++++C+fn++s++++ l+f+++ C+
M03A1.3 670 CCATLKDADKECVDRFCDFNALSQANI---------LNFMST---CG 704
lkgkalpemlfCAAqgrDHtdCClrnGVttPqTlaGkkCLtfCdqspgPR
++g++++em++CA+ +++H+dCC+++GV++ +CL++C++++g
M03A1.3 705 ERGPTVGEMWDCASLRHNHEDCCKKAGVSG-------ECLKYCTAHKG-- 745
RFWfvtqldksYvpCldrfesIkqCFy<-*
++++++ +Y++C+++f++I++CFy
M03A1.3 746 ---APSNYL-DYAFCTENFNEIRDCFY 768