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/F59G1.3
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: F59G1.3 CE26372 S. cerevisiae vacuolar sorting protein status:Partially_confirmed TR:Q21053 protein_id:AAK31535.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
Vps35 Vacuolar protein sorting-associated prot 1553.7 0 1
Anthrax_toxA Anthrax toxin LF subunit 3.6 7.5 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
Anthrax_toxA 1/1 378 396 .. 1 20 [. 3.6 7.5
Vps35 1/1 24 774 .. 1 1065 [] 1553.7 0
Alignments of top-scoring domains:
Anthrax_toxA: domain 1 of 1, from 378 to 396: score 3.6, E = 7.5
*->DslkavARelntklqiRPVe<-*
++ +v+Rel + l i P++
F59G1.3 378 EAFGKVGRELTKLLNI-PID 396
Vps35: domain 1 of 1, from 24 to 774: score 1553.7, E = 0
*->LRTskLSPkkYYeLYMkvfdaleLrvLeiyLkeehrkkghkLadLYE
+RT++LSPk+YY+LYM++++ eL++Le++L++e++++++kL++LYE
F59G1.3 24 MRTAELSPKFYYRLYMDSMH--ELQCLEVNLVQEYAQEPAKLGNLYE 68
LVQYaGnIvPRLYLLITVGsVYIkskeapakeILKDLVEMCRGVQHPLRG
+VQYa++I+PRLYLL+T+G+V+Ik++ +++keILKDLVEMCRGVQHPLRG
F59G1.3 69 CVQYASAIIPRLYLLVTIGGVFIKCGLGSRKEILKDLVEMCRGVQHPLRG 118
LFLRnYLlQcTRdkLPdigsRRvsdvgngpageLgkeepedeedgGtvsD
LFLRnYL+QcTR++LPd+++ +++++++++++L+k++p++++++Gtv+D
F59G1.3 119 LFLRNYLMQCTRSVLPDFPE--TEEMLVAHNDNLSKGTPKLKPRDGTVDD 166
avdFVLqNFtEMNKLWVRmQaGsYvnHQGPstnsvtaCsPttaaLPatsL
++dFVL+NF+EMNKLWVRmQ HQGP
F59G1.3 167 TIDFVLINFAEMNKLWVRMQ------HQGP-------------------- 190
PPaastaassRssDaatRPGsPsRssrekekRekERnELriLVGtNLVRL
s+ekekRek+R+ELriLVGtNLVRL
F59G1.3 191 -------------------------SKEKEKREKDRMELRILVGTNLVRL 215
SQLegvdlemYkdtVLPrILEQvVnCRDalAQeYLmEcIIQVFPDEFHLq
+QLe++++emY+++VLP+ILEQ+V+CRD+++QeYLmEc+IQVF+D+FHL+
F59G1.3 216 AQLEALTEEMYVKDVLPSILEQIVSCRDPISQEYLMECVIQVFADDFHLA 265
TLDtlLgAcaqLqpsVnvknlliaLmdRLadYaareqEsdsevestseev
TL+++L+Ac+qLq++Vn+k+lliaL+dRLa+Y+++ ++e+
F59G1.3 266 TLTEFLNACGQLQQDVNIKILLIALVDRLALYTTS-----------YNEG 304
lpklvdikineevqqkDeqeCPGDkviplqvelFevlsYsasnqvatVIq
+p +p++++lFe++ s+q++t+I+
F59G1.3 305 QP------------------------APTKMQLFEIF----SEQATTLIK 326
srsDmPletPmmLGdiVtLlvsLlnlalkvYPdrlDYvdkvLsnivdvva
+r+DmPl+ diV+L+vsL++la+k+YPdr+DY+++++
F59G1.3 327 NRPDMPLD------DIVALHVSLVSLAVKCYPDRQDYANMTF-------- 362
skfpsapadsseskessgvveTavasaetlthVsigtcvqiikqlslkni
++++q+i+++++++i
F59G1.3 363 -----------------------------------QGLRQVIEEKGVTDI 377
leesehLfLvLLsfqavkkeLtasttaPPnsPvtskklLsiPLisqlekY
++++ +v++eLt+ lL+iP +++Y
F59G1.3 378 EAFG-----------KVGRELTK--------------LLNIP----IDEY 398
nniltllslkLqnyppVmeylDyelrksmAayivqsilknassGaGGGGW
+n+l+l ++L++y++Vm+y+Dy++++++A+y++q++l++
F59G1.3 399 KNVLRL--SQLPEYIKVMNYFDYRGQCNIASYMIQNMLEE---------- 436
kPtvvvtvDseskLhPPpsntlivtaDkLqVdsLfeLispLIkdqddqps
+t+++++D +Vds+f+Lis+L+kdq++q+s
F59G1.3 437 -------------------ETVFRNQD--DVDSAFSLISSLLKDQEKQSS 465
derliykypdeedFaDEvffnrhDeeesWiLDPiQnlVARliHlLrnDDP
d+ +++e+FaDE QnlVARl+Hl+r+DD+
F59G1.3 466 DS------HETEEFADE-----------------QNLVARLLHLIRADDV 492
rkqtmknkiqfsLeeqlkiLnsvrkhlikGGprRlkyTlPpLVfLhRRaa
++q+++Lns+rk+l++GG++Rl+yTlPp++f ++
F59G1.3 493 -------------DSQFLLLNSARKTLGEGGRHRLRYTLPPIIF----EL 525
lrLarqyrrisesllKkrPdngkddgdeasaaAekarkifkfvnqCfshv
+rL++q++++ kd++d+++a k+rk+f+++++
F59G1.3 526 YRLVLQFSDM------------KDEDDKWDA---KIRKMFVCAMG----- 555
hsgtInaLvkvfnveapeLaLkLYLqcAeaAdqcgFenletvAYEFftqA
tI+aLv+ ++e++eL++kLYL++A++Ad+++Fe+++tv+YEF+++A
F59G1.3 556 ---TIGALVS--TAELAELPMKLYLNGAITADRVPFEDNHTVVYEFVSKA 600
FsLYEeeisDSkaQvqaLtlIigtLqkmrvfseENydtLitkctlyASKL
+s++E++++DS+++v++L+l++gtL k+++++eEN+++L+++++l+A+K+
F59G1.3 601 LSILEDDVVDSRDRVRCLHLTVGTLLKTTHLPEENWQPLANQTVLAAAKM 650
LKKPDQCraVylCahLFWsgqlsdedgEeGitDnslkDgkRVleCLkraL
+KKPDQ+r+++++a+L+W+gq+++++gE+ +k+gk+V+++L++a+
F59G1.3 651 FKKPDQVRSLVTVAALYWHGQTLETNGEK------MKNGKKVVDILRKAA 694
KIAdqcldmAnaaRGslvqvqLfveiLNrYlYffihtekgnslitvkhin
KIA++cl+ +lvq+qLf+++L++Y+Y++ e+++s+++v+hi+
F59G1.3 695 KIARECLE-------PLVQQQLFIQLLSAYTYYY---EDNCSEVNVDHIE 734
eLIelIkeeselkslevseetdliaaDdddsdnaanmfasiekhqLkYiP
eLI++++++ +++l+vs+e+d s+ek
F59G1.3 735 ELIARTQDN--AVQLDVSAEAD-----------------SLEK------- 758
ihhfrntleyirsrrese<-*
+++++++++++++ ++
F59G1.3 759 --QLGEAIRRLQLAKLDV 774