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/K02G10.7
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  K02G10.7  CE30087   transmembrane channel protein status:Partially_confirmed TR:Q21159 protein_id:AAK39253.2

Scores for sequence family classification (score includes all domains):
Model           Description                             Score    E-value  N
--------        -----------                             -----    ------- ---
MIP             Major intrinsic protein                 118.2    9.6e-38   2
Transgly_assoc  Transglycosylase associated protein       5.9        8.6   1
Glyco_transf_28 Glycosyltransferase family 28 N-termi     4.4        3.9   1

Parsed for domains:
Model           Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------        ------- ----- -----    ----- -----      -----  -------
Glyco_transf_28   1/1      67    86 ..     1    20 [.     4.4      3.9
MIP               1/2      13   117 ..     1   117 [.    85.1  5.9e-27
MIP               2/2     170   220 ..   216   268 .]    38.4    1e-11
Transgly_assoc    1/1     202   225 ..    28    54 .]     5.9      8.6

Alignments of top-scoring domains:
Glyco_transf_28: domain 1 of 1, from 67 to 86: score 4.4, E = 3.9
                   *->vllaggGTrGHvePavALak<-*
                      v la+  ++GH+ Pa+++a+
    K02G10.7    67    VYLAASVSGGHLNPAISVAQ    86

MIP: domain 1 of 2, from 13 to 117: score 85.1, E = 5.9e-27
                   *->elrsvsfYlraviAEFlaTlLfvFfGvGavsalglkkayasaasvsc
                       + +  f +r+++AE+++T++++ +G++a  ++++++ + s
    K02G10.7    13    RIEDQQF-TRELLAECIGTFFLLLIGNAANIQAAVAVGGNS------ 52

                   agvgllavAlAfGlalfvlVyatghISGgHlNPAVTlallvgrrfsllra
                      +   +A+  G + + +Vy+++ +SGgHlNPA+ +a  + ++ ++ +
    K02G10.7    53 ---TSCHIAW--GIGFMFAVYLAASVSGGHLNPAISVAQSILGNLPPWKI 97

                   iyYiiAQlLGAivGallLyg<-*
                   i+Y iAQ++GA++Ga++ y+
    K02G10.7    98 IPYAIAQVIGAFLGAAVAYF    117

MIP: domain 2 of 2, from 170 to 220: score 38.4, E = 1e-11
                   *->GagMNPARsFGPavvtwkaGwgnlakt.GgrewddhWVywvGPliGA
                      G+++NPAR+FGP+v+   aG    +    + +  ++W+++vG liG+
    K02G10.7   170    GFAINPARDFGPRVFCLCAG---YGWEvFSAHGYYFWIPIVGALIGS 213

                   alaalvY<-*
                   +++a++Y
    K02G10.7   214 IIGAWIY    220

Transgly_assoc: domain 1 of 1, from 202 to 225: score 5.9, E = 8.6
                   *->fWVwiiSlivAviGAiilLaivrlvrg<-*
                      f  wi  ++ A+iG+ii  +i++ + g
    K02G10.7   202    F--WIP-IVGALIGSIIGAWIYKIFVG    225

[DB home][top]