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
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
|
********************************************************************************
MAST - Motif Alignment and Search Tool
********************************************************************************
MAST version 3.0 (Release date: 2004/08/18 09:07:01)
For further information on how to interpret these results or to get
a copy of the MAST software please access http://meme.sdsc.edu.
********************************************************************************
********************************************************************************
REFERENCE
********************************************************************************
If you use this program in your research, please cite:
Timothy L. Bailey and Michael Gribskov,
"Combining evidence using p-values: application to sequence homology
searches", Bioinformatics, 14(48-54), 1998.
********************************************************************************
********************************************************************************
DATABASE AND MOTIFS
********************************************************************************
DATABASE INO_up800.s (nucleotide)
Last updated on Mon Aug 16 21:19:59 2004
Database contains 7 sequences, 5600 residues
Scores for positive and reverse complement strands are combined.
MOTIFS meme.INO_up800.oops.txt (nucleotide)
MOTIF WIDTH BEST POSSIBLE MATCH
----- ----- -------------------
1 12 TTCACATGCCGC
2 10 TCTGGCACAG
PAIRWISE MOTIF CORRELATIONS:
MOTIF 1
----- -----
2 0.32
No overly similar pairs (correlation > 0.60) found.
Random model letter frequencies (from non-redundant database):
A 0.281 C 0.222 G 0.229 T 0.267
********************************************************************************
********************************************************************************
SECTION I: HIGH-SCORING SEQUENCES
********************************************************************************
- Each of the following 7 sequences has E-value less than 10.
- The E-value of a sequence is the expected number of sequences
in a random database of the same size that would match the motifs as
well as the sequence does and is equal to the combined p-value of the
sequence times the number of sequences in the database.
- The combined p-value of a sequence measures the strength of the
match of the sequence to all the motifs and is calculated by
o finding the score of the single best match of each motif
to the sequence (best matches may overlap),
o calculating the sequence p-value of each score,
o forming the product of the p-values,
o taking the p-value of the product.
- The sequence p-value of a score is defined as the
probability of a random sequence of the same length containing
some match with as good or better a score.
- The score for the match of a position in a sequence to a motif
is computed by by summing the appropriate entry from each column of
the position-dependent scoring matrix that represents the motif.
- Sequences shorter than one or more of the motifs are skipped.
- The table is sorted by increasing E-value.
********************************************************************************
SEQUENCE NAME DESCRIPTION E-VALUE LENGTH
------------- ----------- -------- ------
ACC1 sequence of the region up... 6.1e-05 800
CHO1 sequence of the region up... 0.00016 800
INO1 sequence of the region up... 0.00019 800
FAS1 sequence of the region up... 0.00022 800
OPI3 sequence of the region up... 0.00092 800
CHO2 sequence of the region up... 0.0029 800
FAS2 sequence of the region up... 0.0093 800
********************************************************************************
********************************************************************************
SECTION II: MOTIF DIAGRAMS
********************************************************************************
- The ordering and spacing of all non-overlapping motif occurrences
are shown for each high-scoring sequence listed in Section I.
- A motif occurrence is defined as a position in the sequence whose
match to the motif has POSITION p-value less than 0.0001.
- The POSITION p-value of a match is the probability of
a single random subsequence of the length of the motif
scoring at least as well as the observed match.
- For each sequence, all motif occurrences are shown unless there
are overlaps. In that case, a motif occurrence is shown only if its
p-value is less than the product of the p-values of the other
(lower-numbered) motif occurrences that it overlaps.
- The table also shows the E-value of each sequence.
- Spacers and motif occurences are indicated by
o -d- `d' residues separate the end of the preceding motif
occurrence and the start of the following motif occurrence
o [sn] occurrence of motif `n' with p-value less than 0.0001.
A minus sign indicates that the occurrence is on the
reverse complement strand.
********************************************************************************
SEQUENCE NAME E-VALUE MOTIF DIAGRAM
------------- -------- -------------
ACC1 6.1e-05 82_[+1]_137_[+2]_559
CHO1 0.00016 152_[+2]_396_[-2]_42_[+1]_17_
[+1]_149
INO1 0.00019 282_[-2]_327_[-1]_55_[+1]_102
FAS1 0.00022 43_[+2]_41_[+1]_694
OPI3 0.00092 185_[-2]_144_[+1]_449
CHO2 0.0029 353_[+1]_47_[-2]_378
FAS2 0.0093 184_[-2]_372_[+1]_222
********************************************************************************
********************************************************************************
SECTION III: ANNOTATED SEQUENCES
********************************************************************************
- The positions and p-values of the non-overlapping motif occurrences
are shown above the actual sequence for each of the high-scoring
sequences from Section I.
- A motif occurrence is defined as a position in the sequence whose
match to the motif has POSITION p-value less than 0.0001 as
defined in Section II.
- For each sequence, the first line specifies the name of the sequence.
- The second (and possibly more) lines give a description of the
sequence.
- Following the description line(s) is a line giving the length,
combined p-value, and E-value of the sequence as defined in Section I.
- The next line reproduces the motif diagram from Section II.
- The entire sequence is printed on the following lines.
- Motif occurrences are indicated directly above their positions in the
sequence on lines showing
o the motif number of the occurrence (a minus sign indicates that
the occurrence is on the reverse complement strand),
o the position p-value of the occurrence,
o the best possible match to the motif (or its reverse complement), and
o columns whose match to the motif has a positive score (indicated
by a plus sign).
********************************************************************************
ACC1
sequence of the region upstream from YNR016C
LENGTH = 800 COMBINED P-VALUE = 8.78e-06 E-VALUE = 6.1e-05
DIAGRAM: 82_[+1]_137_[+2]_559
[+1]
3.1e-07
TTCACATGCCGC
++++++++++ +
76 TAAAATCTTCACATGGCCCGGCCGCGCGCGCGTTGTGCCAACAAGTCGCAGTCGAAATTCAACCGCTCATTGCCA
[+2]
7.4e-07
TCTGGCACAG
++++++++++
226 TCGTATTCTGGCACAGTATAGCCTAGCACAATCACTGTCACAATTGTTATCGGTTCTACAATTGTTCTGCTCTCT
CHO1
sequence of the region upstream from YER026C
LENGTH = 800 COMBINED P-VALUE = 2.30e-05 E-VALUE = 0.00016
DIAGRAM: 152_[+2]_396_[-2]_42_[+1]_17_[+1]_149
[+2]
3.9e-05
TCTGGCACAG
++++++ +
151 CGTCTGGCGCCCTTCCCATTCCGAACCATGTTATATTGAACCATCTGGCGACAAGCAGTATTAAGCATAATACAT
[-2]
7.4e-07
CTGTGCCAGA
++++++++++
526 CAATCCCCACTCCTTCTCAATGTGTGCAGACTTCTGTGCCAGACACTGAATATATATCAGTAATTGGTCAAAATC
[+1] [+1]
8.7e-07 2.2e-05
TTCACATGCCGC TTCACATGCCGC
++++++ +++ + +++++++++ +
601 ACTTTGAACGTTCACACGGCACCCTCACGCCTTTGAGCTTTCACATGGACCCATCTAAAGATGAAGATCCGTATT
INO1
sequence of the region upstream from YJL153C
LENGTH = 800 COMBINED P-VALUE = 2.71e-05 E-VALUE = 0.00019
DIAGRAM: 282_[-2]_327_[-1]_55_[+1]_102
[-2]
1.8e-05
CTGTGCCAGA
+++ +++ ++
226 ACGTTGTATATGAAACGAGTAGTGAACGTTCGTACGATCTTTCACGCAGACATGCGACTGCGCCCGCCGTAGACC
[-1]
4.2e-08
GCGGCATGTGAA
++++++++++++
601 TGCGCTTCGGCGGCTAAATGCGGCATGTGAAAAGTATTGTCTATTTTATCTTCATCCTTCTTTCCCAGAATATTG
[+1]
1.3e-05
TTCACATGCCGC
+++++++++ ++
676 AACTTATTTAATTCACATGGAGCAGAGAAAGCGCACCTCTGCGTTGGCGGCAATGTTAATTTGAGACGTATATAA
FAS1
sequence of the region upstream from YKL182W
LENGTH = 800 COMBINED P-VALUE = 3.19e-05 E-VALUE = 0.00022
DIAGRAM: 43_[+2]_41_[+1]_694
[+2]
2.2e-05
TCTGGCACAG
++ +++++ +
1 CCGGGTTATAGCAGCGTCTGCTCCGCATCACGATACACGAGGTGCAGGCACGGTTCACTACTCCCCTGGCCTCCA
[+1]
4.2e-08
TTCACATGCCGC
++++++++++++
76 ACAAACGACGGCCAAAAACTTCACATGCCGCCCAGCCAAGCATAATTACGCAACAGCGATCTTTCCGTCGCACAA
OPI3
sequence of the region upstream from YJR073C
LENGTH = 800 COMBINED P-VALUE = 1.32e-04 E-VALUE = 0.00092
DIAGRAM: 185_[-2]_144_[+1]_449
[-2]
7.4e-07
CTGTGCCAGA
++++++++++
151 GTTAATCTGATCAACGCTACGCCGATGACAACGGTCTGTGCCAGATCTGGTTTTCCCCACTTATTTGCTACTTCC
[+1]
5.8e-06
TTCACATGCCGC
++++ ++ ++ +
301 AACTCCGTCAGGTCTTCCACGTGGAACTGCCAAGCCTCCTTCAGATCGCTCTTGTCGACCGTCTCCAAGAGATCC
CHO2
sequence of the region upstream from YGR157W
LENGTH = 800 COMBINED P-VALUE = 4.18e-04 E-VALUE = 0.0029
DIAGRAM: 353_[+1]_47_[-2]_378
[+1]
5.2e-07
TTCACATGCCGC
+++ ++++++++
301 ATATATATTTTTGCCTTGGTTTAAATTGGTCAAGACAGTCAATTGCCACACTTTTCTCATGCCGCATTCATTATT
[-2]
2.9e-05
CTGTGCCAGA
+ +++++++
376 CGCGAAGTTTTCCACACAAAACTGTGAAAATGAACGGCGATGCCAGAAACGGCAAAACCTCAAATGTTAGATAAC
FAS2
sequence of the region upstream from YPL231W
LENGTH = 800 COMBINED P-VALUE = 1.33e-03 E-VALUE = 0.0093
DIAGRAM: 184_[-2]_372_[+1]_222
[-2]
2.9e-05
CTGTGCCAGA
++++++++
151 AACAGGGTGTCGGTCATACCGATAAAGCCGTCAAGAGTGCCAGAAAAGCAAGAAAGAACAAGATTAGATGTTGGT
[+1]
1.9e-06
TTCACATGCCGC
+++++++++ +
526 GCTTAGCAAAATCCAACCATTTTTTTTTTATCTCCCGCGTTTTCACATGCTACCTCATTCGCCTCGTAACGTTAC
********************************************************************************
CPU: pmgm2
Time 0.030000 secs.
mast meme.INO_up800.oops.txt -text -stdout
|
