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/C02A12.1
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  C02A12.1  CE07825   glutathione S-transferase status:Confirmed TR:O16198 protein_id:AAB65264.1

Scores for sequence family classification (score includes all domains):
Model    Description                                    Score    E-value  N
-------- -----------                                    -----    ------- ---
GST_N    Glutathione S-transferase, N-terminal domain    75.5    1.6e-20   1
GST_C    Glutathione S-transferase, C-terminal domain    42.1    1.3e-10   1

Parsed for domains:
Model    Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
-------- ------- ----- -----    ----- -----      -----  -------
GST_N      1/1       3    77 ..     1    88 []    75.5  1.6e-20
GST_C      1/1      83   199 ..     1   121 []    42.1  1.3e-10

Alignments of top-scoring domains:
GST_N: domain 1 of 1, from 3 to 77: score 75.5, E = 1.6e-20
                   *->sstlklfngrgspraervrillalekGveyeivrvdleagehksnkr
                      ++ l++fn rg  +ae+ r++    +Gve+e+vr+++
    C02A12.1     3    KYRLHYFNARG--YAEASRAMFH-MAGVEFEDVRYEI---------- 36

                   peflklnPlm...LpfgqvPaLvdgdGdgkkltESrAIarYlar<-*
                   +++ k    ++++ pfgq+P+L+++   g+k+ +S+AIar++a+
    C02A12.1    37 DDWIKEENTLkneMPFGQMPVLEVD---GEKIPQSVAIARFVAN    77

GST_C: domain 1 of 1, from 83 to 199: score 42.1, E = 1.3e-10
                   *->Gkdpderaaaerarwdmiedsrfadlrldlawlpygkpergaekdee
                      G++p+e +a+  a  d + d  f   + ++a+  +g p + + +de
    C02A12.1    83    GQTPVE-KAWADAFTDLYKD--FLIDMKPWAMIAFGYPGAAGDRDEL 126

                   eveealekllkvlealeeeLkghgkkylvGdkltlADialapaldwlerl
                    +    +++ k+++ l ++L+++++ +l  + +++ D+++++ +  l +l
    C02A12.1   127 KKTSLDPAKEKYFKLLSKRLEKSKSGFLLDSGISFPDLFFFETTTSLIEL 176

                   dpdplfreeyPnLkawrkrvaarp<-*
                   +  +l  +++P ++a++kr+a++p
    C02A12.1   177 EKGFL-GTDFPVVNAYFKRIAEHP    199



[DB home][top]