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:            hmm5c/seq/ZK863.3
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  ZK863.3  CE15439   Yeast LPG22P protein like status:Partially_confirmed TR:Q23651 protein_id:CAB01454.1

Scores for sequence family classification (score includes all domains):
Model        Description                                Score    E-value  N
--------     -----------                                -----    ------- ---
Radical_SAM  Radical SAM superfamily                     48.7    1.6e-12   1
Acetyltransf Acetyltransferase (GNAT) family             43.5    4.8e-14   1

Parsed for domains:
Model        Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------     ------- ----- -----    ----- -----      -----  -------
Radical_SAM    1/1      96   308 ..     1   163 []    48.7  1.6e-12
Acetyltransf   1/1     442   539 ..     1    80 []    43.5  4.8e-14

Alignments of top-scoring domains:
Radical_SAM: domain 1 of 1, from 96 to 308: score 48.7, E = 1.6e-12
                   *->iiitrgCnlr......C.yCafp........................
                      ++  ++C++ + +++ C yC+ +++++ + ++++ ++ ++++ +  +
     ZK863.3    96    MSKPHRCPHInftgniCvYCPGGpdsdfeystqsytgyeptsmrair 142

                   .gk.relspeeileeakelkarlGve.viltGGePlllpdlvellerllk
                    ++++ l  +  l+++ +l+++  + ++i+ GG+++ lp+ +++ ++  +
     ZK863.3   143 aRYnPYLQTRGRLNQLMQLGHSVDKVeFIVMGGTFMSLPEDYRDFFIRNL 192

                   l.....................girisletngtlld..pivrldeellee
                   ++  +++ + + ++    ++++++++  +t +t++d        +++l +
     ZK863.3   193 HdalsghtsasveeavaysersKMKCIGITIETRPDycL-----PRHLND 237

                   LkeagldrvsislqsgddekvlkiinrkhtfeevleaiellreaGipvvi
                   +  +g++r+ i++qs  ++ v + +nr+ht++ v e++ ++++ G++v
     ZK863.3   238 MLLYGCTRLEIGVQSTYED-VARDTNRGHTVKSVCETFHMAKDTGYKV-- 284

                   vvdnivglpGend.edleellell<-*
                   v ++++ lp++  + d e++lel+
     ZK863.3   285 VIHMMPDLPNVGLeRDKEQFLELF    308

Acetyltransf: domain 1 of 1, from 442 to 539: score 43.5, E = 4.8e-14
                   *->lvaeedgelvGfaglspide...............vaeieglaVd..
                      + + + + l+G  +l++i ++ ++++ +++ +  ++++++g +V+
     ZK863.3   442    YEDPKQDILIGLLRLRKISDkahrpelkgnvsvvrELHVYGSVVSva 488

                   ....peyRgkGiGtaLlealeeyareelGlkrieleveedNeaAialYek
                   +++++ ++ +G+G+ L+e++e++aree+G ++i  +  ++ +   ++Y+k
     ZK863.3   489 drdpKKFQHQGYGSLLMEEAERIAREEHGSDKIAVISGVGTR---EYYRK 535

                   lGFk<-*
                   lG++
     ZK863.3   536 LGYE    539



[DB home][top]