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/ZK131.6
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  ZK131.6  CE04501  locus:his-12 histone H2A status:Partially_confirmed TN:CAB05836

Scores for sequence family classification (score includes all domains):
Model         Description                               Score    E-value  N
--------      -----------                               -----    ------- ---
histone       Core histone H2A/H2B/H3/H4                188.1    1.4e-52   1
CBFD_NFYB_HMF Histone-like transcription factor (CBF/     7.6        3.1   1

Parsed for domains:
Model         Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------      ------- ----- -----    ----- -----      -----  -------
CBFD_NFYB_HMF   1/1      27    92 ..     1    67 []     7.6      3.1
histone         1/1       4   121 ..     1   117 []   188.1  1.4e-52

Alignments of top-scoring domains:
CBFD_NFYB_HMF: domain 1 of 1, from 27 to 92: score 7.6, E = 3.1
                   *->LPiAnveRImKkaLPgaarvSkDAkellaecveEfisfiaseAveic
                       P+ +  RI++k    a+rv + A+++la ++e  + ++   A   +
     ZK131.6    27    FPVGRLHRILRKG-NYAQRVGAGAPVYLAAVLEYLAAEVLELAGNAA 72

                   khekRKTikaEDIllAlesL<-*
                    ++k+  i++ +++lA+++
     ZK131.6    73 RDNKKTRIAPRHLQLAVRND    92

histone: domain 1 of 1, from 4 to 121: score 188.1, E = 1.4e-52
                   *->agkkgagkk..kkrsrsrkaglqfyvgRvlkqlhpdtgissragaia
                       gk+g++k+++k++srs++aglqf+vgR ++ l+ ++++++r+ga+a
     ZK131.6     4    RGKGGKAKTggKAKSRSSRAGLQFPVGRLHRILR-KGNYAQRVGAGA 49

                   nsyvadvlErlaaEaleLagnakkdtKkkRitpreiqlAvrldlelelak
                   ++y+a+vlE+laaE+leLagna++d+Kk+Ri pr++qlAvr+d+el++++
     ZK131.6    50 PVYLAAVLEYLAAEVLELAGNAARDNKKTRIAPRHLQLAVRNDEELNKLL 99

                   hAvtegtkaVtpniskrllgkk<-*
                   ++vt+++++V+pni+++ll+kk
     ZK131.6   100 AGVTIAQGGVLPNIQAVLLPKK    121



[DB home][top]