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:            hmm1/seq/B0218.1
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  B0218.1  CE06684   amidase status:Partially_confirmed TR:Q17449 protein_id:AAB00667.1

Scores for sequence family classification (score includes all domains):
Model          Description                              Score    E-value  N
--------       -----------                              -----    ------- ---
Amidase        Amidase                                  595.6     3e-175   1
HAMP           HAMP domain                                5.2        7.4   1
Biotin_carb_C  Biotin carboxylase C-terminal domain       4.9        5.1   1
Polysacc_synt  Polysaccharide biosynthesis protein        4.0        6.4   1
EMP24_GP25L    emp24/gp25L/p24 family                     3.8        6.2   1
PS_pyruv_trans Polysaccharide pyruvyl transferase         3.3        4.8   1

Parsed for domains:
Model          Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------       ------- ----- -----    ----- -----      -----  -------
EMP24_GP25L      1/1       1    30 [.     1    37 [.     3.8      6.2
HAMP             1/1       1    50 [.     1    54 [.     5.2      7.4
PS_pyruv_trans   1/1      36    56 ..   392   412 .]     3.3      4.8
Biotin_carb_C    1/1     283   302 ..    86   105 ..     4.9      5.1
Polysacc_synt    1/1     407   431 ..   271   296 .]     4.0      6.4
Amidase          1/1      84   559 ..     1   521 []   595.6   3e-175

Alignments of top-scoring domains:
EMP24_GP25L: domain 1 of 1, from 1 to 30: score 3.8, E = 6.2
                   *->mlslllllalvtsvdsaltFylpagekkEqvCfyedl<-*
                         ++ l++lv+ + +a+++++ +++kk    ++e+l
     B0218.1     1    --MIFYLVLLVLGA-IAFYVHFSNNRKK----LIERL    30

HAMP: domain 1 of 1, from 1 to 50: score 5.2, E = 7.4
                   *->lllvlllalllalllalllarrltrPlrrLaeaaqriasGrgGdltv
                      +++ l+l++l+a+++++ ++++  + ++rL+ +aqr +    +dl++
     B0218.1     1    MIFYLVLLVLGAIAFYVHFSNNRKKLIERLEIVAQRRR----DDLSK 43

                   rlpvsgr<-*
                   +++  ++
     B0218.1    44 NVEQARK    50

PS_pyruv_trans: domain 1 of 1, from 36 to 56: score 3.3, E = 4.8
                   *->ekettlskrveevrksalknf<-*
                      +++++lsk ve++rk+a+k +
     B0218.1    36    RRRDDLSKNVEQARKAADKLD    56

Biotin_carb_C: domain 1 of 1, from 283 to 302: score 4.9, E = 5.1
                   *->evkTNIpFLialLedpdFls<-*
                      +vkTN++FL+++  d dF+s
     B0218.1   283    DVKTNVEFLRNVWGDIDFQS    302

Polysacc_synt: domain 1 of 1, from 407 to 431: score 4.0, E = 6.4
                   *->kllvillsliasivtavlfPrfsklq<-*
                      ++ v ++++++++v+  +fPr+ +l+
     B0218.1   407    MVPVWIQRILSYPVS-LVFPRMGMLM    431

Amidase: domain 1 of 1, from 84 to 559: score 595.6, E = 3e-175
                   *->ElveaflarieaanpklNvtafvtvffeeAlaaAkaadkrrarkrgg
                      E+++a+ +++  a++k+N   +vt+f+ +A+++A+++d++++  + +
     B0218.1    84    EAIRAYFHKAILAHEKTN---AVTCFILDAERQAEELDEQAK-LPYY 126

                   ekgPLhGvPialKDnidvkGvpTTaGskalegypppyDAtvVerLraaGA
                    k+PL+GvP++lK+++ vkG++TT+G+ + +++p+++D+  Ve ++++G
     B0218.1   127 VKPPLFGVPLSLKECLKVKGYDTTRGFVQDAYHPATEDSIQVEHYKKLGL 176

                   vilGKTNmdEfamgsttensafGpTrNPwdlsrtpGGSSGGSAAaVAagl
                   +++++TN++++++++ ++n+ fG+T NP+d +rt+GGSSGG++A+++ag+
     B0218.1   177 IPFCQTNVPQSLLSYNCSNPLFGTTTNPYDSTRTCGGSSGGEGALIGAGG 226

                   vplAiGtDTGGSIRiPAafCGlvGlKPTyGRvSryGvvgSveplssslDq
                   ++++iGtD+GGS+RiP++f+G++G+KP+  R+ ++G  +Sv++ + + +
     B0218.1   227 SLIGIGTDVGGSVRIPCHFTGTAGIKPSKMRFAHRGGGASVPGKPLIDAN 276

                   vGPlARsVeDaAllLdviaGydpgDptsadsppPtvpwdpvpdyaeglkk
                    GP+A++V+ ++ +L+ ++G+  + ++ ++++pP v+w+++      +++
     B0218.1   277 DGPMAKDVKTNVEFLRNVWGDIDFQSDRDPYCPP-VHWNES------VYS 319

                   dlkgLRiGvpredfgddLdhVaActIRLdgldfvEaDipialnyasldpe
                   ++k+LR+G++ +d+                             +++++p+
     B0218.1   320 SEKKLRVGYYIDDG-----------------------------WFTPTPA 340

                   vdeavekaaaaLeelGhhevvevslpslkealplyyliapaeassnlary
                   +++av++ ++ Le++G h+v+++ +p+l++ ++ly+  ++ +++++++++
     B0218.1   341 LQRAVLESKKHLEAAG-HTVIPFYPPRLPSVMQLYFRAVCLDGGQYVLNK 389

                   dGvryGergleddeasdldelyprtraeglgdevkrrielgayalsagys
                   +  +++e+ ++++   +  ++ p ++++ l+  v + +++++ +++  +
     B0218.1   390 LLKDIIEPTIRFQ---VTLWMVPVWIQRILSYPVSLVFPRMGMLMQSLTR 436

                   dayeylkaqkvrrllrrefdglFeehgvDvLisPTtptpAprigepdkli
                   d++e+++a+++++ +r+ef+gl+ ++++Dv+++P+  +pAp+++ p+k++
     B0218.1   437 DTFELREAYADIEAYREEFVGLMMKDNLDVILCPASIMPAPQHDIPSKVV 486

                   seaddytvlyllddftantvPaNlaGlPaiSvPvgfspeds.wdalvk.e
                    ++ +yt+ly+l+df a++vP++ + + ++   ++++pe+++w++++k++
     B0218.1   487 -SGVSYTCLYNLLDFGAGVVPVTAVSKSDEEKLINEYPETDkWYQITKkA 535

                   ylpegyvGLPvGlQiiGkpgdeetLL<-*
                   +l  g+vG+P+G+Q+ ++p++ee +L
     B0218.1   536 TL--GAVGMPIGVQVAAPPYREEAVL    559



[DB home][top]