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/W02B12.10
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  W02B12.10  CE03770    status:Partially_confirmed TR:Q23126 protein_id:CAA91400.1

Scores for sequence family classification (score includes all domains):
Model          Description                              Score    E-value  N
--------       -----------                              -----    ------- ---
Methyltransf_4 Putative methyltransferase               347.6   3.3e-112   1
heme_1         Heme/Steroid binding domain                6.3          3   1

Parsed for domains:
Model          Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------       ------- ----- -----    ----- -----      -----  -------
heme_1           1/1     127   147 ..    57    78 .]     6.3        3
Methyltransf_4   1/1      49   251 ..     1   254 []   347.6 3.3e-112

Alignments of top-scoring domains:
heme_1: domain 1 of 1, from 127 to 147: score 6.3, E = 3
                   *->aahreeaaeelleklykvGrvd<-*
                      a++r   a++++++++++G++
   W02B12.10   127    AVLR-SNAMKYMPNYFHKGQLS    147

Methyltransf_4: domain 1 of 1, from 49 to 251: score 347.6, E = 3.3e-112
                   *->wpengiairdypgfetlDwqalfgneqppylfsaekfksgpavlvlE
                       p+      ++      Dw +++g++++           g++v +++
   W02B12.10    49    TPN------HM------DWTKYYGDYTK-----------GRQVDFAD 72

                   IGcGmGrflvamAkknPdklfiGiEirvpgvdkalkkiialr.ltAlkgg
                   IGcG+G++l+ +++k+Pd l+iG+Eirv+++d++++ki+alr+++A++g+
   W02B12.10    73 IGCGYGGLLMRLSPKYPDNLMIGMEIRVKVSDYVNEKIQALRkHHAEAGH 122

                   lqNlrilcgdAekvlpnlfpdgsLqkifinFPDPWpKkRHhKRRllqpef
                   ++N+++l+++A+k++pn+f++g+L+k+f++FPDP++K+++hK+R+++p++
   W02B12.10   123 YRNVAVLRSNAMKYMPNYFHKGQLSKMFFLFPDPHFKNKKHKWRIITPTL 172

                   lkeyarvLkpGGvlhlaTDveeyfeemlevlsenfqfEvEikklaqhplf
                   l+eya+vL++GG ++++TDv++++e+m+++lse             hplf
   W02B12.10   173 LSEYAYVLREGGIIYTITDVKDLHEWMVKHLSE-------------HPLF 209

                   erihakvkklKkdiydlriestdlnesPisdelrpsTeyEqkvqrlggpi
                   er+                ++++++++Pi+++l  sTe++qkv r++g++
   W02B12.10   210 ERL----------------TEEEMKKDPIVEMLFESTEEGQKVTRNDGGK 243

                   feliFkkl<-*
                   +++iF++l
   W02B12.10   244 WPAIFRRL    251



[DB home][top]