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/E01G6.1
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: E01G6.1 CE03112 pancreatic trypsin inhibitor (kunitz) protein status:Partially_confirmed TR:Q19021 protein_id:CAA93531.1
Scores for sequence family classification (score includes all domains):
Model Description Score E-value N
-------- ----------- ----- ------- ---
Kunitz_BPTI Kunitz/Bovine pancreatic trypsin inhibi 203.2 4.3e-65 3
EB EB module 62.7 3.4e-15 3
Chitin_bind_1 Chitin recognition protein 5.6 3.8 1
Kelch Kelch motif 4.9 9.6 1
WSC WSC domain 4.7 9.5 1
Herpes_PAP Herpesvirus polymerase accessory protei 1.4 6.3 1
Parsed for domains:
Model Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
EB 1/3 188 238 .. 26 63 .] -8.7 3.8e+04
Kunitz_BPTI 1/3 255 309 .. 1 62 [] 75.1 2.8e-23
Kunitz_BPTI 2/3 362 416 .. 1 62 [] 79.7 8.5e-25
WSC 1/1 400 416 .. 24 40 .. 4.7 9.5
Herpes_PAP 1/1 432 441 .. 312 321 .] 1.4 6.3
Kunitz_BPTI 3/3 466 519 .. 1 62 [] 64.8 6.3e-20
Kelch 1/1 881 908 .. 19 45 .] 4.9 9.6
Chitin_bind_1 1/1 1215 1231 .. 1 18 [. 5.6 3.8
EB 2/3 1189 1245 .. 1 63 [] 78.3 1.6e-19
EB 3/3 1248 1260 .. 48 63 .] 6.7 4.4
Alignments of top-scoring domains:
EB: domain 1 of 3, from 188 to 238: score -8.7, E = 3.8e+04
*->gCeaseQCqggwpgSvCingmC..................eCpegvp
+C+a+ C++g +C + +C + + + + + ++++ C e+
E01G6.1 188 QCSATKPCPNG---FSCNDKKCcpyasacpagaplnngpmKCSEN-- 229
kttrvngvC<-*
+ ++++ C
E01G6.1 230 APCPEGHSC 238
Kunitz_BPTI: domain 1 of 3, from 255 to 309: score 75.1, E = 2.8e-23
*->CslPpdtGptkcpCkggnfarywRwYYnpatgktTkCepFtYgGCgG
C++P+ +G + C ++ +++ R+Y++ t + C++F+++ CgG
E01G6.1 255 CNMPRNSGVQ---CASS-MKSSNRYYFDIGTAT---CRSFKFTQCGG 294
NaNNFetleeCertC<-*
NaNNF +leeCe +C
E01G6.1 295 NANNFGSLEECEGFC 309
Kunitz_BPTI: domain 2 of 3, from 362 to 416: score 79.7, E = 8.5e-25
*->CslPpdtGptkcpCkggnfarywRwYYnpatgktTkCepFtYgGCgG
C++ + +G pC g++ ++ R+Y+np+t+k C+pF+Y GC G
E01G6.1 362 CNEMVSAGT---PCFGRG-LTIQRFYFNPSTQK---CQPFHYYGCSG 401
NaNNFetleeCertC<-*
N+NNFet ++C+++C
E01G6.1 402 NGNNFETVDQCQNFC 416
WSC: domain 1 of 1, from 400 to 416: score 4.7, E = 9.5
*->vadsysfqsveaCqsaC<-*
++++++f++v +Cq++C
E01G6.1 400 SGNGNNFETVDQCQNFC 416
Herpes_PAP: domain 1 of 1, from 432 to 441: score 1.4, E = 6.3
*->nPNtKkQkCg<-*
nPN+ Q+C
E01G6.1 432 NPNSELQRCS 441
Kunitz_BPTI: domain 3 of 3, from 466 to 519: score 64.8, E = 6.3e-20
*->CslPpdtGptkcpCkggnfarywRwYYnpatgktTkCepFtYgGCgG
C+ + +G+ Ckg ++ + +w+Y+++++k C +F+Y+GCgG
E01G6.1 466 CNANMSRGN---GCKG--STQRSMWFYDKSKKK---CSQFVYNGCGG 504
NaNNFetleeCertC<-*
N+F t+ +C + C
E01G6.1 505 TPNRFTTKVACTESC 519
Kelch: domain 1 of 1, from 881 to 908: score 4.9, E = 9.6
*->Gydg.qslnsvevYDpetntWeklpsmp<-*
G++ +q+l+ v + +p tn+ +++++m+
E01G6.1 881 GVSQpQCLRGVAYVNPATNQRQFCSPMK 908
Chitin_bind_1: domain 1 of 1, from 1215 to 1231: score 5.6, E = 3.8
*->aqqCGrQaggalCpnnlC<-*
+qC + + + +C+n++C
E01G6.1 1215 -EQCSSREANSTCDNGYC 1231
EB: domain 2 of 3, from 1189 to 1245: score 78.3, E = 1.6e-19
*->ClsgqVlingqCvakelnkvnpGesgCeaseQCqggwpgSvCingmC
C++g+ +++g+C l+ + G++gC+++eQC++ +++S+C+ng+C
E01G6.1 1189 CPPGFFNMDGKC----LKILFAGQKGCISDEQCSSREANSTCDNGYC 1231
eCpegvpkttrvngvC<-*
Cp++ k ++g+C
E01G6.1 1232 ICPTD--KPLVHGGKC 1245
EB: domain 3 of 3, from 1248 to 1260: score 6.7, E = 4.4
*->eCpegvpkttrvngvC<-*
+Cpeg + + g+C
E01G6.1 1248 TCPEG---FSNIAGRC 1260
[output cut off at A = 10 top alignments]