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/Y54G9A.3
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  Y54G9A.3  CE19231  locus:kqt-3 Ion transport proteins status:Partially_confirmed TR:Q9XWG9 protein_id:CAA21699.1

Scores for sequence family classification (score includes all domains):
Model          Description                              Score    E-value  N
--------       -----------                              -----    ------- ---
KCNQ_channel   KCNQ voltage-gated potassium channel     153.2      6e-43   1
ion_trans      Ion transport protein                     85.3    5.6e-22   1
IQ             IQ calmodulin-binding motif                9.7       0.87   2
PS_pyruv_trans Polysaccharide pyruvyl transferase         7.8       0.26   1

Parsed for domains:
Model          Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------       ------- ----- -----    ----- -----      -----  -------
ion_trans        1/1      52   325 ..     1   228 []    85.3  5.6e-22
IQ               1/2     347   367 ..     1    21 []    13.6    0.065
IQ               2/2     489   509 ..     1    21 []    -3.0  3.4e+03
PS_pyruv_trans   1/1     501   546 ..   367   412 .]     7.8     0.26
KCNQ_channel     1/1     479   583 ..    86   192 ..   153.2    6e-43

Alignments of top-scoring domains:
ion_trans: domain 1 of 1, from 52 to 325: score 85.3, E = 5.6e-22
                   *->ileildyvftviFtlEmllkiialgfklhsriprkYlrspwnilDfl
                      +l+il++v +v+F +E++++++++g++       kY++        +
    Y54G9A.3    52    LLYILEIVLVVFFSVEFIVRLWSAGCR------SKYIGVYGR--LKF 90

                   vvlpslislilflsgggsv.kl.ir.l.......................
                   v  p   +li ++   g+++++ ir++ ++++++  ++++ +++ + +++
    Y54G9A.3    91 VRKPI--TLIGRF---GGRhLQlIRdMlrhggpsfcgkrdagdtvsadse 135

                   ...res.l.flalse.vsygilsvlrllRllRvLRllrllrrapglf.e.
                   + +r+s   ++ ++  +  + + ++++l +     +++ +  +  ++ ++
    Y54G9A.3   136 dapRRSaGwHMVEIRiSKIPKIEIFKSLKIG----IFKNFANWNFQQfRk 181

                   ..................................fe.lrtlG.Wklqslg
                    + ++  +++  ++ +  ++ ++++++++ + + f++l +l+++  + lg
    Y54G9A.3   182 lefskilkknlagrvhywgsatpssgpdpldvgcFPkLHQLAsF-RRLLG 230

                   .rslkslfnlllllllllfifaiigvqlfggelnkccnpingntnFdnfg
                   ++ ++   +l+  l + +++ ++ +++++ +e+++   +++g+  F+++
    Y54G9A.3   231 sVVFIHRQELITTLYIGFLGLIFSSYFVYLAEKDHI--GVDGRQAFTSYA 278

                   eaflwlFqtlTteGwgdimpdtlagtvlgkiyfvifiilggflllNlfia
                   +a++w ++t+Tt+G+gd++p+    t+lg i++ +++ ++++ +++l+ +
    Y54G9A.3   279 DALWWGVITMTTIGYGDVVPQ----TWLGRIVA-SCFSIFAISFFALPAG 323

                   vi<-*
                   ++
    Y54G9A.3   324 IL    325

IQ: domain 1 of 2, from 347 to 367: score 13.6, E = 0.065
                   *->rraaikIQaawRGylaRkryk<-*
                       +aa++IQ+ wR++ a k+
    Y54G9A.3   347    PAAATLIQCLWRCHAAEKKVS    367

IQ: domain 2 of 2, from 489 to 509: score -3.0, E = 3.4e+03
                   *->rraaikIQaawRGylaRkryk<-*
                      +r ai+  ++++ ++aR+r++
    Y54G9A.3   489    HRNAIRAIRRVKYFVARRRFQ    509

PS_pyruv_trans: domain 1 of 1, from 501 to 546: score 7.8, E = 0.26
                   *->yaidaeelaagpllevvldlleqllekettlskrveevrksalknf<
                      y+++++++++ ++++ v+d  eq+++   ++  r++e+++++++
    Y54G9A.3   501    YFVARRRFQQARKPYDVRDVIEQYSQGHLNMMVRIKELQRRLDQTL  546

                   -*

    Y54G9A.3     -    -

KCNQ_channel: domain 1 of 1, from 479 to 583: score 153.2, E = 6e-43
                   *->itHvsqLResHRvtvKviRRmqYFvAKKKFqqARKPYDvRDviEqYs
                      i Hv  L + HR  ++ iRR++YFvA+++FqqARKPYDvRDviEqYs
    Y54G9A.3   479    ISHVCELTDAHRNAIRAIRRVKYFVARRRFQQARKPYDVRDVIEQYS 525

                   qGHLnLmvRiKELqRRLDqslGKPslFlsvsEKskDRGk.ntiGARLnRv
                   qGHLn+mvRiKELqRRLDq lGKP  +     K   +G++ tiG RL R
    Y54G9A.3   526 QGHLNMMVRIKELQRRLDQTLGKPGQY---DGKGSRKGHpVTIGSRLSRL 572

                   EeKvvqlDqkL<-*
                   E ++  lD k+
    Y54G9A.3   573 ELQMSSLDRKV    583



[DB home][top]