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/R05D3.11
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  R05D3.11  CE25069    status:Partially_confirmed SW:P34544 protein_id:AAK21437.1

Scores for sequence family classification (score includes all domains):
Model     Description                                   Score    E-value  N
--------  -----------                                   -----    ------- ---
Pre-SET   Pre-SET motif                                 161.1    8.4e-51   1
SET       SET domain                                    130.2    3.8e-35   1
MBD       Methyl-CpG binding domain                      74.6    5.2e-19   1
pkinase_C Protein kinase C terminal domain                4.3        4.6   1

Parsed for domains:
Model     Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------  ------- ----- -----    ----- -----      -----  -------
pkinase_C   1/1     165   178 ..    55    70 .]     4.3      4.6
MBD         1/1     830   885 ..     1    55 []    74.6  5.2e-19
Pre-SET     1/1     913  1051 ..     1   124 []   161.1  8.4e-51
SET         1/1    1077  1307 ..     9   265 .]   130.2  3.8e-35

Alignments of top-scoring domains:
pkinase_C: domain 1 of 1, from 165 to 178: score 4.3, E = 4.6
                   *->eeFrGFSYvnsPersf<-*
                      e+F++F++ n P+ ++
    R05D3.11   165    EQFKNFEFTN-PN-YL    178

MBD: domain 1 of 1, from 830 to 885: score 74.6, E = 5.2e-19
                   *->drlpLppGWrRetvqRksG..rgagkvDvyYysPvcGKklRSksevi
                       ++pL++GWrR ++++++G++rga k++++Y+sP cG +l+++s+v
    R05D3.11   830    IHTPLLCGWRRIMYTMSTGkkRGAVKKNIIYFSP-CGAALHQISDVS 875

                   rYLgangdls<-*
                   +Y++ +++l+
    R05D3.11   876 EYIHVTRSLL    885

Pre-SET: domain 1 of 1, from 913 to 1051: score 161.1, E = 8.4e-51
                   *->DisngkEsvpIplVENeVDldkpiPedFtYIneyiptegviipqseW
                      D+s g E++pIplV N+VD+d p P + +Y++ +++++++++++s+
    R05D3.11   913    DFSLGTEGIPIPLV-NSVDNDEP-P-SLEYSKRRFQYNDQVDISSV- 955

                   WSleflsGCdCtddkgsCspAdsskCaCarlsgnggaypph.k.ipe...
                    S++f+sGC+C +   +Cs  d+skC+C++ls+++++ +ph+++++++++
    R05D3.11   956 -SRDFCSGCSCDG---DCS--DASKCECQQLSIEAMKRLPHnLqFDGhde 999

                   ...................egylnnrilvepgplIYECNsrCkCgp.sCk
                     ++++++++  +    + ++y+n+ +  ++ +++YECN++C+C+++sC+
    R05D3.11  1000 lyessekqnkflklfffrvPHYQNRLLSSKVISGLYECNDQCSCHRkSCY 1049

                   NR<-*
                   NR
    R05D3.11  1050 NR    1051

SET: domain 1 of 1, from 1077 to 1307: score 130.2, E = 3.8e-35
                   *->legfkt...GwGlrAtedIpkGefiiEYpGelitseeaekreddllg
                      l++fkt ++GwG+rA+ dIp  +fi+ Y+G  i++++++ +e
    R05D3.11  1077    LQIFKTaqsGWGVRALTDIPQSTFICTYVG-AILTDDLA-DE----- 1116

                   pkyydddgdksldeYlfdldalfLlrEkleSWkpYiiLPktydsPlfWse
                    ++ +d+       Y+ dld
    R05D3.11  1117 -LRNADQ-------YFADLDL----------------------------- 1129

                   eellLkGTllevdeeelkeykdllaeeievflqlklFpdveteefsyedy
                                       kd
    R05D3.11  1130 --------------------KDT--------------------------- 1132

                   lWaasivSRafsl.....................................
                             ++l++++++++++ + ++++++ +++++++++++++  ++
    R05D3.11  1133 ----------VELekgredhetdfgyggdesdyddeegsdgdsgddvmnk 1172

                   ..................................................
                     ++++++++++++++ ++++++++++++++++ ++++++++++ ++++
    R05D3.11  1173 mvkrqdssesgeetkrltrqkrkqskksgkggsvekddttprdsmekdni 1222

                   ..................vegvnreeenseeelglviDatrkGNlaRFiN
                   ++++++  + ++  ++ ++               +viDa+ +GNl+RF+N
    R05D3.11  1223 eskdepvfnwdkyfepfpL---------------YVIDAKQRGNLGRFLN 1257

                   HSCsPNcevqvvfvdgedlgtsseeetDAeekaglepgdprivifAlrdI
                   HSC+PN+ vq+v++d++dl                  + p++++f+ + +
    R05D3.11  1258 HSCDPNVHVQHVMYDTHDL------------------RLPWVAFFTRKYV 1289

                   kaGeELtydYgydlsnge<-*
                   kaG+ELt+dY+y  +  +
    R05D3.11  1290 KAGDELTWDYQYTQDQTA    1307



[DB home][top]