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/F18A12.6
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: F18A12.6 CE09459 peptidase status:Partially_confirmed TR:O16790 protein_id:AAB66078.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
Peptidase_M13 Peptidase family M13 215.6 1.7e-66 1
PhoU PhoU family 6.5 1.8 1
Vault Major Vault Protein repeat 6.3 2.2 1
C1q C1q domain 4.3 9.2 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
Vault 1/1 361 385 .. 1 28 [. 6.3 2.2
C1q 1/1 478 490 .. 125 137 .] 4.3 9.2
PhoU 1/1 561 590 .. 17 47 .. 6.5 1.8
Peptidase_M13 1/1 451 648 .. 1 222 [] 215.6 1.7e-66
Alignments of top-scoring domains:
Vault: domain 1 of 1, from 361 to 385: score 6.3, E = 2.2
*->ewlvvgpgaetYlPqvpeeVvekPvpai<-*
+++++ + ++P v+ +V+ k v+ai
F18A12.6 361 LEMIRESD--SFTPSVKKRVLKK-VEAI 385
C1q: domain 1 of 1, from 478 to 490: score 4.3, E = 9.2
*->gegvhSvFSGFLl<-*
+e+v+ +F GFL+
F18A12.6 478 PEYVNLIFTGFLI 490
PhoU: domain 1 of 1, from 561 to 590: score 6.5, E = 1.8
*->eAlaerDenelAeevieelDdeiDdlyeqie<-*
+A+ ++D+ ++ +++i ++D+ +D+ly++i+
F18A12.6 561 RAFKKMDL-SKMQKLIGFEDKNLDQLYFRIQ 590
Peptidase_M13: domain 1 of 1, from 451 to 648: score 215.6, E = 1.7e-66
*->NAyynpefNsitfpagilvpPffdltyPravnskyGgiGfVlgHEig
NA y+ +N +t +a +l++Pffd tyP++vn + +Gf +gHE g
F18A12.6 451 NAHYDSKDNALTVLAPFLDDPFFDSTYPEYVN--LIFTGFLIGHEFG 495
HgFDdqGiqldeeGsLadcwnlrdWkwtdednkcfqdrvkclidqyseav
H D ++ + +W ++ed++++++r++clidqy++++
F18A12.6 496 HSID---PKILRRD---------GW-YKTEDMTEYGKRAQCLIDQYDNYD 532
kPddkgwfnncankGattlgEniADlgGlraAlkaykkkssdeereeqrl
+Pd g+++n G+++ gE+++D G + ++a+kk ++++ q+l
F18A12.6 533 DPDH-GKQMN----GTYCIGEIVGDWVGRDVTWRAFKKMDLSKM---QKL 574
pgLenftpdQLfFetyAyswCrksrpenslrpllvdpHsPnelRvNgavr
+g+e+ + dQL+F+ +++C +r +sl +ll+dpH++ +RvNg+++
F18A12.6 575 IGFEDKNLDQLYFRIQSLFFCG-PRSMKSLEQLLSDPHPTEVFRVNGIYS 623
nmPaFasaFnCkpGdkMfpepekrc<-*
nmP+Fa+aFnC+ G++M+pe+++++
F18A12.6 624 NMPQFAKAFNCPIGSPMNPEKKCKM 648