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/R107.7
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: R107.7 CE00302 locus:gst-1 Glutathione S-transferase P subunit status:Confirmed SW:P10299 protein_id:CAA78471.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
GST_N Glutathione S-transferase, N-terminal 88.0 3.9e-24 1
GST_C Glutathione S-transferase, C-terminal 85.1 8.5e-23 1
Retrotrans_gag Retrotransposon gag protein 5.2 4.7 1
EFP Elongation factor P (EF-P) 3.1 4.6 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
GST_N 1/1 2 72 .. 1 88 [] 88.0 3.9e-24
Retrotrans_gag 1/1 77 96 .. 81 100 .] 5.2 4.7
EFP 1/1 133 149 .. 74 90 .. 3.1 4.6
GST_C 1/1 78 186 .. 1 121 [] 85.1 8.5e-23
Alignments of top-scoring domains:
GST_N: domain 1 of 1, from 2 to 72: score 88.0, E = 3.9e-24
*->sstlklfngrgspraervrillalekGveyeivrvdleagehksnkr
+ +l++f++ g +ae++r+lla +k v ye++rv
R107.7 2 TLKLTYFDIHG--LAEPIRLLLA-DKQVAYEDHRVTY---------- 35
peflklnPlmLpfgqvPaLvdgdGdgkkltESrAIarYlar<-*
++++ ++P+m fgqvP+L g ++ +++S+AI+r+lar
R107.7 36 EQWADIKPKM-IFGQVPCLLSG---DEEIVQSGAIIRHLAR 72
Retrotrans_gag: domain 1 of 1, from 77 to 96: score 5.2, E = 4.7
*->hgdvrdeealisaFlrGLre<-*
+g ++e i++F++GLr+
R107.7 77 NGSNETETTFIDMFYEGLRD 96
EFP: domain 1 of 1, from 133 to 149: score 3.1, E = 4.6
*->YLYkDGdeyvFMDqetY<-*
++Yk+G+ yv D+e+Y
R107.7 133 HTYKNGEHYVIGDKESY 149
GST_C: domain 1 of 1, from 78 to 186: score 85.1, E = 8.5e-23
*->Gkdpderaaaerarwdmiedsrfadlrldlawlpygkpergaekdee
G+++ e +dm+++ +++dl ++ ++y+++e +++
R107.7 78 GSNE-----TETTFIDMFYE-GLRDLHTKYTTMIYRNYE---DGKAP 115
eveealekllkvlealeeeLkghgkkylvGdkltlADialapaldwlerl
+ + +l+ l le+l ++k+ g++y+ Gdk ++AD++l+++ld + +l
R107.7 116 YIKDVLPGELARLEKLFHTYKN-GEHYVIGDKESYADYVLFEELDIHLIL 164
dpdplfreeyPnLkawrkrvaarp<-*
p++l ++ P+Lk++ +r+a+rp
R107.7 165 TPNAL--DGVPALKKFHERFAERP 186