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:            hmm5/seq/Y53F4B.31
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  Y53F4B.31  CE22418   glutathione-S-transferase status:Predicted TR:Q9NAB2 protein_id:CAB61105.1

Scores for sequence family classification (score includes all domains):
Model    Description                                    Score    E-value  N
-------- -----------                                    -----    ------- ---
GST_N    Glutathione S-transferase, N-terminal domain    61.4    2.1e-16   1
GST_C    Glutathione S-transferase, C-terminal domain    49.1    1.4e-12   1

Parsed for domains:
Model    Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
-------- ------- ----- -----    ----- -----      -----  -------
GST_N      1/1       3    75 ..     1    88 []    61.4  2.1e-16
GST_C      1/1      81   195 ..     1   121 []    49.1  1.4e-12

Alignments of top-scoring domains:
GST_N: domain 1 of 1, from 3 to 75: score 61.4, E = 2.1e-16
                   *->sstlklfngrgspraervrillalekGveyeivrvdleagehksnkr
                      +++l++fn rg  +ae++ril    ++v++e+ r   ++g
   Y53F4B.31     3    TYKLTYFNARG--YAEAARILFH-LADVPFEDFRMTHGDGTW----- 41

                   peflklnPlmLpfgqvPaLvdgdGdgkkltESrAIarYlar<-*
                   ++++   P+    gqvP+L ++   g+++ +S AI+rYla+
   Y53F4B.31    42 EKLKDGAPF----GQVPYLTVD---GFDIPQSAAIIRYLAN    75

GST_C: domain 1 of 1, from 81 to 195: score 49.1, E = 1.4e-12
                   *->Gkdpderaaaerarwdmiedsrfadlrldlawlpygkperg...aek
                      Gk+p     +e    d i+d +f d ++ +  + ++ +++g+++ e
   Y53F4B.31    81    GKTP-----EEQVWADAIVD-QFKDFMGSFR-ERIMAHFAGksqEEI 120

                   deeeveealekllkvlealeeeLkghgkkylvGdkltlADialapaldwl
                    +   e + ++ + ++++l++ L+++++ +lvGd++t+ADi+  + +  l
   Y53F4B.31   121 AKITSEVVIPARDSYFKILNGLLEKSKSGFLVGDSITFADIVISENITTL 170

                   erldpdplfreeyPnLkawrkrvaarp<-*
                   e  ++ +   +e PnL+a+r++v a p
   Y53F4B.31   171 EKNQLFVA--SEQPNLAALREKVYAIP    195



[DB home][top]