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/Y65B4A.9
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  Y65B4A.9  CE22732    status:Partially_confirmed TR:Q9BL55 protein_id:AAK29951.1

Scores for sequence family classification (score includes all domains):
Model    Description                                    Score    E-value  N
-------- -----------                                    -----    ------- ---
pkinase  Protein kinase domain                           60.8    9.2e-16   2
MT-A70   MT-A70                                           2.9        6.9   1

Parsed for domains:
Model    Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
-------- ------- ----- -----    ----- -----      -----  -------
pkinase    1/2      21    51 ..     1    32 [.    23.1  2.5e-05
MT-A70     1/1     143   155 ..   181   195 .]     2.9      6.9
pkinase    2/2      89   172 ..    71   149 ..    36.3  5.8e-09

Alignments of top-scoring domains:
pkinase: domain 1 of 2, from 21 to 51: score 23.1, E = 2.5e-05
                   *->yelleklGeGsfGkVykakhkdktgkiVAvKi<-*
                      ye+ ++l  G+fG++yk+k+  ++gk +AvK+
    Y65B4A.9    21    YEVTKPLATGTFGSIYKVKRE-SDGKFFAVKC    51

MT-A70: domain 1 of 1, from 143 to 155: score 2.9, E = 6.9
                   *->RdhnlrPgWltwGne<-*
                      Rd  + P+ +++G+e
    Y65B4A.9   143    RD--IKPSHFCIGRE    155

pkinase: domain 2 of 2, from 89 to 172: score 36.3, E = 5.8e-09
                   *->ddhlylvmEymegGdLfdylrrng..gplsekeakkialQilrGleY
                      ++ l +vm++ +   L++++++ ++++++s+   +++a Q+l ++
    Y65B4A.9    89    NRFLFIVMPLYGE-NLYELMMNTNkdRKFSMATGLHLAEQTLAAIRD 134

                   LHsngivHRDLKpeNILlden....dgtvKiaDFGLAr<-*
                   LH+ng++HRD+Kp+ + + ++ ++++++v + DFGL++
    Y65B4A.9   135 LHRNGFIHRDIKPSHFCIGREidgqHHQVYLLDFGLCK    172



[DB home][top]