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:            hmm3/seq/F56B3.10
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  F56B3.10  CE27168    status:Partially_confirmed TR:O45115 protein_id:AAK29835.1

Scores for sequence family classification (score includes all domains):
Model    Description                                    Score    E-value  N
-------- -----------                                    -----    ------- ---
GST_N    Glutathione S-transferase, N-terminal domain    57.3    3.3e-15   1
GST_C    Glutathione S-transferase, C-terminal domain    36.6    4.7e-09   1

Parsed for domains:
Model    Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
-------- ------- ----- -----    ----- -----      -----  -------
GST_N      1/1       8    81 ..     1    88 []    57.3  3.3e-15
GST_C      1/1      91   199 ..     1   121 []    36.6  4.7e-09

Alignments of top-scoring domains:
GST_N: domain 1 of 1, from 8 to 81: score 57.3, E = 3.3e-15
                   *->sstlklfngrgspraervrillalekGveyeivrvdleagehksnkr
                      +++l++ n+rg  rae++r++    +Gv++e+ r d +++      +
    F56B3.10     8    RYKLYYVNMRG--RAEPIRLVFH-FLGVDFEDYRMDKGDF------T 45

                   peflklnPlmLpfgqvPaLvdgdGdgkkltESrAIarYlar<-*
                     ++ + P+    +q+P  +++ G +++l++++ I+rYla+
    F56B3.10    46 GAMKDKAPM----KQLPFIEIDGG-KTTLCQTVSICRYLAK    81

GST_C: domain 1 of 1, from 91 to 199: score 36.6, E = 4.7e-09
                   *->Gkdpderaaaerarwdmiedsrfadlrldlawlpygkpergaekdee
                      G++       + a +dm +d +fadl+  +a+  y ++     kd
    F56B3.10    91    GATK-----TDSAKVDMMAD-GFADLFQLAAMAKYAPDP---IKDS- 127

                   eveealekllkvlealeeeLkghgkkylvGdkltlADialapaldwlerl
                   +    +e + + le   + Lk+++++y+vG ++ ++D++ + +l+ l +
    F56B3.10   128 MMGTYKESIGPKLENMQDLLKKSKGEYFVGKSIHWCDVYILGVLQALDES 177

                   dpdplfreeyPnLkawrkrvaarp<-*
                   d  ++  +  P+L  ++ r+ ++p
    F56B3.10   178 DDGVF--DNLPELREYYLRMRNLP    199



[DB home][top]