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:            hmm2/seq/D1053.1
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  D1053.1  CE03099   glutathione S-transferase status:Predicted SW:Q18938 protein_id:CAA91449.1

Scores for sequence family classification (score includes all domains):
Model    Description                                    Score    E-value  N
-------- -----------                                    -----    ------- ---
GST_N    Glutathione S-transferase, N-terminal domain    81.7    2.7e-22   1
GST_C    Glutathione S-transferase, C-terminal domain    23.6    2.2e-05   1

Parsed for domains:
Model    Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
-------- ------- ----- -----    ----- -----      -----  -------
GST_N      1/1       3    78 ..     1    88 []    81.7  2.7e-22
GST_C      1/1      87   197 ..     1   121 []    23.6  2.2e-05

Alignments of top-scoring domains:
GST_N: domain 1 of 1, from 3 to 78: score 81.7, E = 2.7e-22
                   *->sstlklfngrgspraervrillalekGveyeivrvdleagehksnkr
                       ++  l++++ s++++rvri+la  k+v+ye++ vdl   e+ +
     D1053.1     3    NQKPVLYSYWRSSCSWRVRIALA-LKNVDYEYKTVDL-LSEEAK--- 44

                   peflklnPlmLpfgqvPaLvdgdGdgkkltESrAIarYlar<-*
                   ++++++nP      +vP++v++   g ++tES+AI++Yl++
     D1053.1    45 SKLKEINPA----AKVPTFVVD---GQVITESLAIIEYLEE    78

GST_C: domain 1 of 1, from 87 to 197: score 23.6, E = 2.2e-05
                   *->Gkdpderaaaerarwdmiedsrfadlrldlawlpygkperg.aekde
                      +kdp +r  a    ++ +++s++++l + +++ ++ k+e+g +++ +
     D1053.1    87    PKDPIKR--AHARAISLLVASGIQPLHNLKVLQLLNKKEAGfGGQFA 131

                   eeveealekllkvlealeeeLkghgkkylvGdkltlADialapaldwler
                         +  + + l ale  Lk+h +ky+vGd +t+AD+   p+++ + r
     D1053.1   132 ------KQFVVEGLTALEILLKQHSGKYAVGDDVTIADLSIPPLIYSANR 175

                   ldpdplfreeyPnLkawrkrvaarp<-*
                   +++d+   ++yP +++  +++a  p
     D1053.1   176 FNLDL---SPYPTVNRINETLADIP    197



[DB home][top]