1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
|
hmmsearch - search a sequence database with a profile HMM
HMMER 2.2g (August 2001)
Copyright (C) 1992-2001 HHMI/Washington University School of Medicine
Freely distributed under the GNU General Public License (GPL)
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
HMM file: Peptidase_C1.hmm [Peptidase_C1]
Sequence database: cysprot1b.fa
per-sequence score cutoff: [none]
per-domain score cutoff: [none]
per-sequence Eval cutoff: <= 10
per-domain Eval cutoff: [none]
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query HMM: Peptidase_C1
Accession: PF00112
Description: Papain family cysteine protease
[HMM has been calibrated; E-values are empirical estimates]
Scores for complete sequences (score includes all domains):
Sequence Description Score E-value N
-------- ----------- ----- ------- ---
CATL_RAT 449.4 2e-135 1
CATL_HUMAN 444.5 6.1e-134 1
CATH_RAT 381.8 4.8e-115 1
PAPA_CARPA 337.7 9e-102 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CATL_RAT 1/1 114 332 .. 1 337 [] 449.4 2e-135
CATL_HUMAN 1/1 114 332 .. 1 337 [] 444.5 6.1e-134
CATH_RAT 1/1 114 330 .. 1 337 [] 381.8 4.8e-115
PAPA_CARPA 1/1 134 343 .. 1 337 [] 337.7 9e-102
Alignments of top-scoring domains:
CATL_RAT: domain 1 of 1, from 114 to 332: score 449.4, E = 2e-135
*->lPesfDWReWkggaVtpVKdQGiqCGSCWAFSavgalEgryciktgt
+P+++DWRe kg VtpVK+QG qCGSCWAFSa g lEg+ ++kt
CATL_RAT 114 IPKTVDWRE-KG-CVTPVKNQG-QCGSCWAFSASGCLEGQMFLKT-- 155
kawggklvsLSEQqLvDCdgedygnngesCGyGCnGGGlmdnAfeYikke
gkl+sLSEQ+LvDC++ d gn+ GCnG Glmd Af+Yik+
CATL_RAT 156 ----GKLISLSEQNLVDCSH-DQGNQ------GCNG-GLMDFAFQYIKE- 192
qIsnNgGlvtEsdYekgCkPYtdfPCgkdggndtyypCpgkaydpndTgt
NgGl++E++Y PY+ +kd g+
CATL_RAT 193 ----NGGLDSEESY-----PYE----AKD-------------------GS 210
CkynckknskypktyakikgygdvpynvsTydEealqkalaknGPvsVai
Cky+ + ++ a+++g++d+p++ E+al+ka+a++GP+sVa+
CATL_RAT 211 CKYR-AEYAV-----ANDTGFVDIPQQ-----EKALMKAVATVGPISVAM 249
dasedskgDFqlYksGendvgyGvYkhtsageCggtpfteLdHAVliVGY
das+ s q+Y+sG +Y+++ C+++ +LdH+Vl+VGY
CATL_RAT 250 DASHPS---LQFYSSG-------IYYEP---NCSSK---DLDHGVLVVGY 283
GteneggtfdetssskksesgiqvssgsngssgSSgssgapiedkgkdYW
G e+ ++++ +YW
CATL_RAT 284 GYEG-T------------------------------------DSNKDKYW 296
IVKNSWGtdWGEnGYfriaRgknksgkneCGIaseasypi<-*
+VKNSWG++WG++GY++ia+++n n+CG+a+ asypi
CATL_RAT 297 LVKNSWGKEWGMDGYIKIAKDRN----NHCGLATAASYPI 332
CATL_HUMAN: domain 1 of 1, from 114 to 332: score 444.5, E = 6.1e-134
*->lPesfDWReWkggaVtpVKdQGiqCGSCWAFSavgalEgryciktgt
+P s+DWRe kg +VtpVK+QG qCGSCWAFSa+galEg+ ++kt
CATL_HUMAN 114 APRSVDWRE-KG-YVTPVKNQG-QCGSCWAFSATGALEGQMFRKT-- 155
kawggklvsLSEQqLvDCdgedygnngesCGyGCnGGGlmdnAfeYikke
g l+sLSEQ+LvDC+g + gn+ GCnG Glmd+Af+Y+++
CATL_HUMAN 156 ----GRLISLSEQNLVDCSG-PQGNE------GCNG-GLMDYAFQYVQD- 192
qIsnNgGlvtEsdYekgCkPYtdfPCgkdggndtyypCpgkaydpndTgt
NgGl++E++Y PY+ +++ +
CATL_HUMAN 193 ----NGGLDSEESY-----PYE----ATE-------------------ES 210
CkynckknskypktyakikgygdvpynvsTydEealqkalaknGPvsVai
Ckyn +k s+ a+++g++d+p + E+al+ka+a++GP+sVai
CATL_HUMAN 211 CKYN-PKYSV-----ANDTGFVDIPKQ-----EKALMKAVATVGPISVAI 249
dasedskgDFqlYksGendvgyGvYkhtsageCggtpfteLdHAVliVGY
da++ s F +Yk G +Y ++ +C+++ + dH+Vl+VGY
CATL_HUMAN 250 DAGHES---FLFYKEG-------IYFEP---DCSSE---DMDHGVLVVGY 283
GteneggtfdetssskksesgiqvssgsngssgSSgssgapiedkgkdYW
G e+ e+++ +YW
CATL_HUMAN 284 GFES-T------------------------------------ESDNNKYW 296
IVKNSWGtdWGEnGYfriaRgknksgkneCGIaseasypi<-*
+VKNSWG++WG+ GY+++a+++ n+CGIas asyp+
CATL_HUMAN 297 LVKNSWGEEWGMGGYVKMAKDRR----NHCGIASAASYPT 332
CATH_RAT: domain 1 of 1, from 114 to 330: score 381.8, E = 4.8e-115
*->lPesfDWReWkggaVtpVKdQGiqCGSCWAFSavgalEgryciktgt
P s+DWR+ kg V+pVK+QG CGSCW FS++galE++ +i++
CATH_RAT 114 YPSSMDWRK-KGNVVSPVKNQG-ACGSCWTFSTTGALESAVAIAS-- 156
kawggklvsLSEQqLvDCdgedygnngesCGyGCnGGGlmdnAfeYikke
gk L EQqLvDC +++n+ GC+G Gl+++AfeYi++
CATH_RAT 157 ----GKMMTLAEQQLVDCAQ-NFNNH------GCQG-GLPSQAFEYILY- 193
qIsnNgGlvtEsdYekgCkPYtdfPCgkdggndtyypCpgkaydpndTgt
N+G++ E++Y PY gk+ g+
CATH_RAT 194 ----NKGIMGEDSY-----PYI----GKN-------------------GQ 211
CkynckknskypktyakikgygdvpynvsTydEealqkalaknGPvsVai
Ck+n +++++ a++k+ ++++ n dE+a+ +a+a + Pvs a+
CATH_RAT 212 CKFN-PEKAV-----AFVKNVVNITLN----DEAAMVEAVALYNPVSFAF 251
dasedskgDFqlYksGendvgyGvYkhtsageCggtpfteLdHAVliVGY
+++e DF++YksG vY++ +C +tp + ++HAVl+VGY
CATH_RAT 252 EVTE----DFMMYKSG-------VYSSN---SCHKTP-DKVNHAVLAVGY 286
GteneggtfdetssskksesgiqvssgsngssgSSgssgapiedkgkdYW
G +n g YW
CATH_RAT 287 GEQN-GLL----------------------------------------YW 295
IVKNSWGtdWGEnGYfriaRgknksgkneCGIaseasypi<-*
IVKNSWG++WG nGYf i+Rgkn +CG+a +asypi
CATH_RAT 296 IVKNSWGSNWGNNGYFLIERGKN-----MCGLAACASYPI 330
PAPA_CARPA: domain 1 of 1, from 134 to 343: score 337.7, E = 9e-102
*->lPesfDWReWkggaVtpVKdQGiqCGSCWAFSavgalEgryciktgt
+Pe +DWR+ kg aVtpVK+QG +CGSCWAFSav ++Eg+++i+t
PAPA_CARPA 134 IPEYVDWRQ-KG-AVTPVKNQG-SCGSCWAFSAVVTIEGIIKIRT-- 175
kawggklvsLSEQqLvDCdgedygnngesCGyGCnGGGlmdnAfeYikke
g+l +SEQ+L+DCd+ ++ GCnG G+++ A++ + +
PAPA_CARPA 176 ----GNLNEYSEQELLDCDR---RSY------GCNG-GYPWSALQLVAQ- 210
qIsnNgGlvtEsdYekgCkPYtdfPCgkdggndtyypCpgkaydpndTgt
G++ Y PY+ g++
PAPA_CARPA 211 -----YGIHYRNTY-----PYE----GVQ-------------------RY 227
CkynckknskypktyakikgygdvpynvsTydEealqkalaknGPvsVai
C+++ +k+ + +ak +g ++v+++ +E al + +a+ +PvsV
PAPA_CARPA 228 CRSR-EKGPY----AAKTDGVRQVQPY----NEGALLYSIAN-QPVSVVL 267
dasedskgDFqlYksGendvgyGvYkhtsageCggtpfteLdHAVliVGY
+a + DFqlY++G ++++ +Cg+ +dHAV++VGY
PAPA_CARPA 268 EAAGK---DFQLYRGG-------IFVG----PCGN----KVDHAVAAVGY 299
GteneggtfdetssskksesgiqvssgsngssgSSgssgapiedkgkdYW
G +Y
PAPA_CARPA 300 G---------------------------------------------PNYI 304
IVKNSWGtdWGEnGYfriaRgknksgkneCGIaseasypi<-*
++KNSWGt WGEnGY+ri+Rg+++s ++ CG+ ++ yp+
PAPA_CARPA 305 LIKNSWGTGWGENGYIRIKRGTGNS-YGVCGLYTSSFYPV 343
Histogram of all scores:
score obs exp (one = represents 1 sequences)
----- --- ---
> 337 4 -|====
% Statistical details of theoretical EVD fit:
mu = -195.8384
lambda = 0.1423
chi-sq statistic = 0.0000
P(chi-square) = 0
Total sequences searched: 4
Whole sequence top hits:
tophits_s report:
Total hits: 4
Satisfying E cutoff: 4
Total memory: 16K
Domain top hits:
tophits_s report:
Total hits: 4
Satisfying E cutoff: 4
Total memory: 20K
|
