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/Y59C2A.1
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  Y59C2A.1  CE31658   carboxypeptidase status:Confirmed TR:Q9TZC6 protein_id:AAC68743.2

Scores for sequence family classification (score includes all domains):
Model        Description                                Score    E-value  N
--------     -----------                                -----    ------- ---
Zn_carbOpept Zinc carboxypeptidase                      241.4    1.3e-68   1
TSC22        TSC-22/dip/bun family                        3.6          6   1

Parsed for domains:
Model        Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------     ------- ----- -----    ----- -----      -----  -------
TSC22          1/1      85   102 ..    43    60 .]     3.6        6
Zn_carbOpept   1/1     145   350 ..     1   215 []   241.4  1.3e-68

Alignments of top-scoring domains:
TSC22: domain 1 of 1, from 85 to 102: score 3.6, E = 6
                   *->nvSpEqLaqlqaqlqlae<-*
                      +v pE L+++ + l  ae
    Y59C2A.1    85    LVAPEFLQKFKSLLTSAE    102

Zn_carbOpept: domain 1 of 1, from 145 to 350: score 241.4, E = 1.3e-68
                   *->ytYhnlqeilawLdsvaseypslvrkesiGkSyeGRelyvLkisadp
                      + Yh ++e+++++++++ ++ ++v+ +++ +S eGR++y +ki  +p
    Y59C2A.1   145    NAYHSYNEMVEFMKLLSEQKSDMVEMVKVATSSEGRSIYGVKIH-PP 190

                   gehgpgkPeifldariHarEwvsselmlglleellseygsdariTkLrds
                   g   p+kP+i +da++HarEw+++++ l+ ++++++eyg++++ T  +++
    Y59C2A.1   191 GPSPPEKPSIIVDAGVHAREWIAPAVGLFMIRKIVEEYGRNPQVTANLQK 240

                   lrfyIvPvlNPDGgEYevtwgsdRmWRKtRsgragSs..cvGtDLNRNwp
                   +++yI+P +NPDG  Ye+++++dR+WRKtRs++  +++ cvG+D NRNw+
    Y59C2A.1   241 FDWYIMPQVNPDG--YEYSRTTDRLWRKTRSKNVTVNrwCVGADANRNWG 288

                   dqWgceggasssPCSEtYrgpaaesePEtkavkdfirsrvlgekgqikvv
                   + Wg e ga + PCS++Y g+ + sePE++  k f   +     + +++
    Y59C2A.1   289 YRWG-EAGANRTPCSNIYMGSHPYSEPEIRGLKEFF--T----WQITNPM 331

                   lyidlHgySqlanlfPYGyt<-*
                   +yi lH+y ql+ l P+Gyt
    Y59C2A.1   332 VYISLHSYGQLL-LSPWGYT    350



[DB home][top]