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:            hmm3/seq/K02F3.2
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  K02F3.2  CE01348    status:Partially_confirmed SW:Q21153 protein_id:AAK21421.1

Scores for sequence family classification (score includes all domains):
Model          Description                              Score    E-value  N
--------       -----------                              -----    ------- ---
mito_carr      Mitochondrial carrier protein            301.0    1.5e-86   3
efhand         EF hand                                   53.2    8.4e-13   3
formyl_trans_C Formyl transferase, C-terminal domain      9.5       0.17   1
UPF0021        PP-loop family UPF0021                     5.4        2.1   1
DUF479         Protein of unknown function, DUF479        2.4        8.7   1
PFK            Phosphofructokinase                        1.0        5.6   1

Parsed for domains:
Model          Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------       ------- ----- -----    ----- -----      -----  -------
formyl_trans_C   1/1      71    85 ..    91   105 .]     9.5     0.17
efhand           1/3      97   125 ..     1    29 []    12.1     0.18
efhand           2/3     131   159 ..     1    29 []    17.6   0.0056
efhand           3/3     205   230 ..     4    29 .]    23.5  0.00014
DUF479           1/1     232   243 ..     1    12 [.     2.4      8.7
UPF0021          1/1     241   250 ..     1    10 [.     5.4      2.1
mito_carr        1/3     363   461 ..     1   125 []   112.4  8.9e-30
PFK              1/1     502   508 ..     1     7 [.     1.0      5.6
mito_carr        2/3     462   553 ..     1   125 []    82.3  4.3e-21
mito_carr        3/3     555   649 ..     1   125 []   113.1  5.4e-30

Alignments of top-scoring domains:
formyl_trans_C: domain 1 of 1, from 71 to 85: score 9.5, E = 0.17
                   *->GKKamkaedflnGrr<-*
                      GKK+m  edf++G+
     K02F3.2    71    GKKLMTPEDFIRGYL    85

efhand: domain 1 of 3, from 97 to 125: score 12.1, E = 0.18
                   *->elkeaFkefDkDgDGkIsfeEfkaalkkl<-*
                       ++ + + +D+ +DG IsfeEf+a+ + l
     K02F3.2    97    TVRLLASAADTTKDGDISFEEFCAFEALL    125

efhand: domain 2 of 3, from 131 to 159: score 17.6, E = 0.0056
                   *->elkeaFkefDkDgDGkIsfeEfkaalkkl<-*
                          aF++fD++ + +Is +Ef a+++++
     K02F3.2   131    LYLTAFELFDRNASDTISCDEFEAVIRHT    159

efhand: domain 3 of 3, from 205 to 230: score 23.5, E = 0.00014
                   *->eaFkefDkDgDGkIsfeEfkaalkkl<-*
                      +aFk +Dk+g+G+Is  +f+++++++
     K02F3.2   205    QAFKRYDKNGNGTISSLDFQQIMTTV    230

DUF479: domain 1 of 1, from 232 to 243: score 2.4, E = 8.7
                   *->GnLlgDFVkGdp<-*
                      G Ll DFV+ ++
     K02F3.2   232    GHLLTDFVRHNL    243

UPF0021: domain 1 of 1, from 241 to 250: score 5.4, E = 2.1
                   *->kiLvAvSGGl<-*
                      ++L+AvSGG+
     K02F3.2   241    HNLIAVSGGA    250

mito_carr: domain 1 of 3, from 363 to 461: score 112.4, E = 8.9e-30
                   *->sflksllaGgiAGaiaatvtyPlDvvKtRLQtqgssslskssvnsdv
                        ++++l G++AGa++at++yP+D+vKtR+Q+q++s+
     K02F3.2   363    ESAYRFLLGSVAGACGATAVYPIDLVKTRMQNQRTSG---------- 399

                   gksksyrsakksnasasrkYkgildcfkkivkeEGliktrgLYkGlgpnl
                    +  ++             Yk+ ldcfkk+vk EGl+   gLY+Gl+p +
     K02F3.2   400 -SFVGEV-----------MYKNSLDCFKKVVKFEGLL---GLYRGLLPQI 434

                   lrvapygaaiyFgtYEqlKklllkkkge<-*
                   ++vap   ai+++  +++++ ++k+++
     K02F3.2   435 VGVAPE-KAIKLTMNDYMRDKFTKDGKI    461

PFK: domain 1 of 1, from 502 to 508: score 1.0, E = 5.6
                   *->KkIgVLT<-*
                      KkIgV+T
     K02F3.2   502    KKIGVFT    508

mito_carr: domain 2 of 3, from 462 to 553: score 82.3, E = 4.3e-21
                   *->sflksllaGgiAGaiaatvtyPlDvvKtRLQtqgssslskssvnsdv
                      +++  ++aGg++G++ +++t+Pl+ vK+RLQt g+  +++
     K02F3.2   462    PLYGEIIAGGTGGMCQVVFTNPLEIVKIRLQTAGEVQQAG------- 501

                   gksksyrsakksnasasrkYkgildcfkkivkeEGliktrgLYkGlgpnl
                                     k+ g++    +++ke G++   gLYkG  + +
     K02F3.2   502 ------------------KKIGVF----TVLKELGFL---GLYKGSRACF 526

                   lrvapygaaiyFgtYEqlKklllkkkge<-*
                   lr++p+ +aiyF +Y   K    + +g+
     K02F3.2   527 LRDIPF-SAIYFPAYAHAKLASADEDGM    553

mito_carr: domain 3 of 3, from 555 to 649: score 113.1, E = 5.4e-30
                   *->sflksllaGgiAGaiaatvtyPlDvvKtRLQtqgssslskssvnsdv
                      s+++ ++++ iAG+ aa +++P Dv+KtRLQ+  + +++
     K02F3.2   555    SPGTLFASAFIAGVPAAGLVTPADVIKTRLQVAARAGQT-------- 593

                   gksksyrsakksnasasrkYkgildcfkkivkeEGliktrgLYkGlgpnl
                                      Y+g++dc++k++keEG++   +L+kG ++++
     K02F3.2   594 ------------------TYNGVIDCARKLIKEEGPM---SLWKGTAARV 622

                   lrvapygaaiyFgtYEqlKklllkkkge<-*
                    r +p+ +a++++tYE+l++l++ ++ +
     K02F3.2   623 CRSSPQ-FAVTLLTYEVLQRLFYVDFAG    649

	[output cut off at A = 10 top alignments]


[DB home][top]