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]