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:            hmm2/seq/F26A3.8
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  F26A3.8  CE27141  locus:rrf-1  status:Confirmed TR:Q9NDH1 protein_id:CAC42290.1

Scores for sequence family classification (score includes all domains):
Model        Description                                Score    E-value  N
--------     -----------                                -----    ------- ---
RdRP         RNA dependent RNA polymerase               803.4   2.4e-246   1
Saccharop_dh Saccharopine dehydrogenase                   2.3          5   1
Caulimo_VI   Caulimovirus viroplasmin                     1.5        7.1   1

Parsed for domains:
Model        Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------     ------- ----- -----    ----- -----      -----  -------
Saccharop_dh   1/1     687   703 ..   475   492 .]     2.3        5
RdRP           1/1     628  1042 ..     1   450 []   803.4 2.4e-246
Caulimo_VI     1/1    1298  1330 ..   471   508 .]     1.5      7.1

Alignments of top-scoring domains:
Saccharop_dh: domain 1 of 1, from 687 to 703: score 2.3, E = 5
                   *->pYdPfLerLnkyGLpwke<-*
                      + +PfL++++k+ L++ +
     F26A3.8   687    I-EPFLDEIRKWALVNGV    703

RdRP: domain 1 of 1, from 628 to 1042: score 803.4, E = 2.4e-246
                   *->eipDiegkpcIhsdGtgYcFsDGvGmiSlnklArkvarklnldrgGs
                       ++D++g++  +++G++Y+FsDGvGm+S+ ++A++v++ ++++
     F26A3.8   628    STYDLTGGR--ATNGSEYTFSDGVGMMSY-QFAQEVSQAMQFG---- 667

                   rvkevPsAfQiRygydGyAyKGVlavdpnlddirnW.........kgrkv
                     k vPs+fQiR++  G   KGV+a++p ld+ir+W   ++ ++ k+ k+
     F26A3.8   668 --KAVPSCFQIRFR--G--NKGVIAIEPFLDEIRKWalvngvtsmKMAKC 711

                   qlRpSMiKFeadskntlleqlEvvkfsrpspvtLNRqlItlLsqlGVeds
                    +RpS+iKF+  +k + ++q+E+vkfs+ + v+LN+++I++L+q+++++s
     F26A3.8   712 LFRPSQIKFQ--AKAISGDQIEMVKFSSAVLVALNKPFINILDQVSEMQS 759

                   vfeemRiTnqeellldlnalleskqkaqevlaanklgenyrltidilalm
                   +++++RiT+++e+l+d++ l+++kq+ +e++++nkl+e++r+ idi++l+
     F26A3.8   760 LDCHKRITSRIEELMDRQILSFAKQMNEETFCRNKLKEFPRR-IDIDNLR 808

                   liaGfelsdEPflrsmLsllrksklreLkeKarIfvPideGrsLlGvvDE
                   +++Gf+ls+EPf+rs++++++k++++++++K++I++P+++Grs+lGvvDE
     F26A3.8   809 TMWGFTLSSEPFFRSLIKASIKFSITKQLCKEQIQIPSELGRSMLGVVDE 858

                   tGtLehygqVfvqvskldiekllpktsnrkksgksklevitGtVlVtKNP
                   tG+L+ ygq+fvq++k++++kl+p++sn+k++g s  e++tGtVl+tKNP
     F26A3.8   859 TGRLQ-YGQIFVQYTKNYKKKLPPRDSNNKVHG-S--EIVTGTVLLTKNP 904

                   cLHPGDVRVlkAVdvPaLehyvGNlkdvVvFPQkGpRphpdecSGsDLDG
                   c++PGDVR+++AVd+P+L+h    ++dvVvFPQ+GpRphpde++GsDLDG
     F26A3.8   905 CIVPGDVRIFEAVDIPELHH----MCDVVVFPQHGPRPHPDEMAGSDLDG 950

                   DeYFVIWDqkLlpernpeePmlsesekskiyiGrepseldreetdeelie
                   DeY+VIWDq+Ll+ern eeP+++++ek        +  +dre++d +++e
     F26A3.8   951 DEYSVIWDQELLLERN-EEPFDFAVEKI-------KVPYDREKLDVLMRE 992

                   fFvstkYlksdslGlIsnAHlaiaDrfltLGDerepekarskkCleLAkl
                   f+v+  Ylk+ds+G+Isn+Hl+++D        ++++++r  +C++LAk+
     F26A3.8   993 FYVT--YLKLDSVGQISNSHLHNSD--------QYGLNSR--VCMDLAKK 1030

                   hsdAVDfpKtGv<-*
                   +++AVDf+K+G+
     F26A3.8  1031 NCQAVDFTKSGQ    1042

Caulimo_VI: domain 1 of 1, from 1298 to 1330: score 1.5, E = 7.1
                   *->eshlCqyCGdhPrdeknvpdkskavitsyNneekSkkd<-*
                      e+h +qy  dh  d kn++++    it     ekSk++
     F26A3.8  1298    EEHRKQYLIDHDADFKNFCEL-DHLITG----EKSKEA    1330



[DB home][top]