MEME

Search sequence databases with these motifs using MAST.
Submit these motifs to BLOCKS multiple alignment processor.
Build and use a motif-based hidden Markov model (HMM) using Meta-MEME.

MEME - Motif discovery tool

MEME version 3.0 (Release date: 2001/03/05 14:24:28)

For further information on how to interpret these results or to get a copy of the MEME software please access http://meme.sdsc.edu.

This file may be used as input to the MAST algorithm for searching sequence databases for matches to groups of motifs. MAST is available for interactive use and downloading at http://meme.sdsc.edu.

REFERENCE

If you use this program in your research, please cite:

Timothy L. Bailey and Charles Elkan, "Fitting a mixture model by expectation maximization to discover motifs in biopolymers", Proceedings of the Second International Conference on Intelligent Systems for Molecular Biology, pp. 28-36, AAAI Press, Menlo Park, California, 1994.

TRAINING SET

DATAFILE= ../positive_negative_promoters_3.0_200.fa
ALPHABET= ACGT
Sequence name           Weight Length  Sequence name           Weight Length  
-------------           ------ ------  -------------           ------ ------  
CG6640:                 1.0000    201  CG5150:                 1.0000    201  
CG4653:                 1.0000    201  CG3868:                 1.0000    201  
CG7874:                 1.0000    201  CG13095:                1.0000    201  
Cyp6a8:                 1.0000    201  Cyt-b5:                 1.0000    201  
CG5839:                 1.0000    201  CG5932:                 1.0000    201  
Chit:                   1.0000    201  CG6981:                 1.0000    201  
CG7298:                 1.0000    201  CG6296:                 1.0000    201  
CG18522:                1.0000    201  CG10433:                1.0000    201  
mfas:                   1.0000    201  glob1:                  1.0000    201  
Adh:                    1.0000    201  Pu:                     1.0000    201  
Hsp27:                  1.0000    201  CG4471:                 1.0000    201  
CG5171:                 1.0000    201  CG2471:                 1.0000    201  
CG8502:                 1.0000    201  CG11892:                1.0000    201  
ImpE2:                  1.0000    201  CG9021:                 1.0000    201  
CG11459:                1.0000    201

COMMAND LINE SUMMARY

This information can also be useful in the event you wish to report a
problem with the MEME software.

command: meme -dna -nmotifs 10 ../positive_negative_promoters_3.0_200.fa 

model:  mod=         zoops    nmotifs=        10    evt=           inf
object function=  E-value of product of p-values
width:  minw=            8    maxw=           50    minic=        0.00
width:  wg=             11    ws=              1    endgaps=       yes
nsites: minsites=        2    maxsites=       29    wnsites=       0.8
theta:  prob=            1    spmap=         uni    spfuzz=        0.5
em:     prior=   dirichlet    b=            0.01    maxiter=        50
        distance=    1e-05
data:   n=            5829    N=              29
strands: +
sample: seed=            0    seqfrac=         1
Letter frequencies in dataset:
A 0.276 C 0.232 G 0.215 T 0.277 
Background letter frequencies (from dataset with add-one prior applied):
A 0.276 C 0.232 G 0.215 T 0.277

P N
MOTIF 1 width = 16 sites = 10 llr = 131 E-value = 1.3e-001

NAME	START	P-VALUE				SITES
			Simplified	A	7	:	:	:	5	:	:	:	:	:	3	:	:	:	3	4
			pos.-specific	C	:	:	8	:	1	6	1	:	5	7	1	7	:	2	3	:
			probability	G	3	1	2	a	4	:	9	6	5	3	6	2	:	7	3	6
			matrix	T	:	9	:	:	:	4	:	4	:	:	:	1	a	1	1	:
.
			bits	2.2
				2.0
				1.8
				1.6
			Information	1.3
			content	1.1
			(18.9 bits)	0.9
				0.7
				0.4
				0.2
				0.0
.
			Multilevel		`A`	`T`	`C`	`G`	`A`	`C`	`G`	`G`	`C`	`C`	`G`	`C`	`T`	`G`	`A`	`G`
			consensus		`G`		`G`		`G`	`T`		`T`	`G`	`G`	`A`	`G`		`C`	`C`	`A`
			sequence																`G`

.

CG5171:	66	7.46e-09	`GAGTCGTTGA`		`A`	`T`	`C`	`G`	`G`	`C`	`G`	`G`	`G`	`C`	`G`	`G`	`T`	`G`	`G`	`G`	`TGGCGTACTG`
ImpE2:	143	1.85e-08	`TGCCAACAAA`		`G`	`T`	`C`	`G`	`A`	`C`	`G`	`G`	`C`	`G`	`G`	`C`	`T`	`G`	`C`	`G`	`CAGCCGGTGG`
mfas:	104	3.00e-08	`CGGCAGAGGA`		`A`	`T`	`C`	`G`	`A`	`C`	`G`	`G`	`C`	`G`	`A`	`C`	`T`	`G`	`C`	`G`	`GCAGAGCAGC`
Adh:	132	4.47e-07	`TGCTGCATCC`		`G`	`T`	`C`	`G`	`A`	`C`	`G`	`T`	`C`	`G`	`A`	`C`	`T`	`G`	`C`	`A`	`CTCGCCCCCA`
CG4471:	123	5.29e-07	`CTGACTTCTT`		`A`	`T`	`C`	`G`	`C`	`T`	`G`	`G`	`C`	`C`	`C`	`C`	`T`	`G`	`G`	`G`	`AAAGGTCGCA`
CG6640:	107	5.29e-07	`AGTTATGTTT`		`A`	`T`	`G`	`G`	`A`	`T`	`G`	`T`	`G`	`C`	`G`	`G`	`T`	`G`	`A`	`G`	`ATAAGCTGAT`
CG2471:	49	6.28e-07	`GTATCTGCGA`		`A`	`T`	`C`	`G`	`G`	`T`	`G`	`T`	`G`	`C`	`G`	`C`	`T`	`T`	`G`	`A`	`AAATGCCAAA`
CG5839:	174	7.45e-07	`GCCGTGTTCA`		`A`	`G`	`C`	`G`	`G`	`C`	`G`	`G`	`C`	`C`	`G`	`C`	`T`	`C`	`A`	`A`	`TAACAAAGTG`
CG4653:	44	2.22e-06	`TAGCGATCTA`		`A`	`T`	`C`	`G`	`A`	`T`	`G`	`G`	`G`	`C`	`A`	`T`	`T`	`G`	`T`	`A`	`TGGCGATAAG`
Pu:	29	3.61e-06	`ATGCCAAACC`		`G`	`T`	`G`	`G`	`G`	`C`	`C`	`T`	`G`	`C`	`G`	`C`	`T`	`C`	`A`	`G`	`CTCGGTTTAG`

Motif 1 block diagrams

Name

Lowest
p-value

Motifs

CG5171:

7.5e-09

ImpE2:

1.9e-08

mfas:

3e-08

Adh:

4.5e-07

CG4471:

5.3e-07

CG6640:

5.3e-07

CG2471:

6.3e-07

CG5839:

7.4e-07

CG4653:

2.2e-06

Pu:

3.6e-06

SCALE

\|	\|	\|	\|	\|	\|	\|	\|
1	25	50	75	100	125	150	175

Motif 1 in BLOCKS format

to BLOCKS multiple alignment processor.

Motif 1 position-specific scoring matrix

Motif 1 position-specific probability matrix






Time 33.81 secs.

P N
MOTIF 2 width = 13 sites = 20 llr = 185 E-value = 1.1e-002

NAME	START	P-VALUE				SITES
			Simplified	A	:	:	a	:	a	7	9	7	6	:	1	3	3
			pos.-specific	C	6	3	:	1	:	:	:	:	:	3	6	3	1
			probability	G	3	:	:	1	:	:	:	1	3	3	3	4	6
			matrix	T	1	8	:	9	:	3	1	3	1	3	:	:	:
.
			bits	2.2
				2.0
				1.8
				1.6
			Information	1.3
			content	1.1
			(13.4 bits)	0.9
				0.7
				0.4
				0.2
				0.0
.
			Multilevel		`C`	`T`	`A`	`T`	`A`	`A`	`A`	`A`	`A`	`G`	`C`	`G`	`G`
			consensus		`G`	`C`				`T`		`T`	`G`	`T`	`G`	`A`	`A`
			sequence											`C`		`C`

.

CG7298:	153	3.01e-07	`AATCGTAGAC`		`C`	`T`	`A`	`T`	`A`	`A`	`A`	`A`	`A`	`C`	`C`	`A`	`G`	`GCAATGGGAA`
CG2471:	155	1.04e-06	`CTTCTTCGTG`		`C`	`T`	`A`	`T`	`A`	`A`	`A`	`A`	`A`	`C`	`C`	`G`	`A`	`AAGTAAAAAG`
CG8502:	149	1.25e-06	`TGGCACATTT`		`C`	`T`	`A`	`T`	`A`	`T`	`A`	`A`	`A`	`T`	`C`	`C`	`G`	`GTGTACTGGA`
Hsp27:	158	1.25e-06	`CCCAGCGTCA`		`G`	`T`	`A`	`T`	`A`	`A`	`A`	`A`	`G`	`C`	`C`	`G`	`G`	`CGTCAACGTC`
CG6296:	156	5.55e-06	`CAGACTGAAG`		`G`	`T`	`A`	`T`	`A`	`A`	`A`	`T`	`A`	`G`	`C`	`A`	`G`	`TACCACGAAC`
CG5932:	154	1.14e-05	`TAAGAACGGG`		`C`	`T`	`A`	`T`	`A`	`T`	`A`	`A`	`A`	`G`	`A`	`A`	`G`	`GCCCCACGTC`
CG18522:	157	1.27e-05	`CGGAGAACCA`		`C`	`T`	`A`	`T`	`A`	`T`	`A`	`T`	`A`	`T`	`G`	`G`	`G`	`GCTACGGAAG`
CG9021:	159	1.40e-05	`GCTTCTCGTG`		`C`	`T`	`A`	`T`	`A`	`A`	`A`	`A`	`G`	`G`	`A`	`G`	`A`	`AAACAACGAC`
glob1:	153	1.57e-05	`GAGTTCGCTC`		`C`	`C`	`A`	`T`	`A`	`A`	`A`	`A`	`G`	`G`	`C`	`C`	`C`	`CACAATTTGC`
mfas:	130	1.90e-05	`GCAGAGCAGC`		`T`	`T`	`A`	`T`	`A`	`A`	`A`	`T`	`A`	`C`	`C`	`A`	`G`	`ATGCTCGACG`
CG4653:	147	2.10e-05	`ACGGGGTGTT`		`C`	`T`	`A`	`T`	`A`	`A`	`A`	`G`	`G`	`T`	`C`	`C`	`G`	`CCGGCCGGTA`
ImpE2:	172	2.74e-05	`CCGGTGGTTT`		`G`	`T`	`A`	`T`	`A`	`A`	`A`	`A`	`T`	`G`	`C`	`G`	`A`	`TGGGCCTAAA`
CG5839:	154	2.98e-05	`ATGGCCAAGT`		`C`	`C`	`A`	`T`	`A`	`A`	`A`	`A`	`A`	`T`	`G`	`C`	`C`	`GTGTTCAAGC`
CG11892:	148	3.51e-05	`TGAAATTTTG`		`C`	`C`	`A`	`T`	`A`	`A`	`A`	`A`	`G`	`C`	`G`	`C`	`A`	`GGTACTCTGA`
CG13095:	46	3.51e-05	`ATATGCCAGA`		`T`	`T`	`A`	`T`	`A`	`A`	`A`	`T`	`A`	`T`	`C`	`G`	`A`	`ATGCCGGACG`
CG6640:	154	5.95e-05	`TCTCGAACCC`		`T`	`T`	`A`	`T`	`A`	`A`	`A`	`T`	`A`	`G`	`A`	`A`	`G`	`TGGATAGTTT`
CG5171:	105	8.75e-05	`TTGTTCAAAA`		`G`	`C`	`A`	`C`	`A`	`A`	`A`	`A`	`A`	`C`	`C`	`C`	`G`	`GCGACAATTT`
CG7874:	154	8.75e-05	`GATCGGCCAG`		`C`	`T`	`A`	`T`	`A`	`T`	`A`	`A`	`T`	`T`	`G`	`G`	`C`	`CGTATTCTCG`
CG3868:	139	1.11e-04	`GAGTCAGTGC`		`C`	`C`	`A`	`G`	`A`	`T`	`A`	`A`	`G`	`G`	`C`	`A`	`G`	`CTATATAACC`
CG5150:	129	1.24e-04	`GCACCTGATG`		`G`	`T`	`A`	`T`	`A`	`T`	`T`	`A`	`A`	`T`	`G`	`G`	`G`	`TCGCATGGAA`

Motif 2 block diagrams

Name

Lowest
p-value

Motifs

CG7298:

3e-07

CG2471:

1e-06

CG8502:

1.3e-06

Hsp27:

1.3e-06

CG6296:

5.6e-06

CG5932:

1.1e-05

CG18522:

1.3e-05

CG9021:

1.4e-05

glob1:

1.6e-05

mfas:

1.9e-05

CG4653:

2.1e-05

ImpE2:

2.7e-05

CG5839:

3e-05

CG11892:

3.5e-05

CG13095:

3.5e-05

CG6640:

5.9e-05

CG5171:

8.8e-05

CG7874:

8.8e-05

CG3868:

0.00011

CG5150:

0.00012

SCALE

\|	\|	\|	\|	\|	\|	\|	\|
1	25	50	75	100	125	150	175

Motif 2 in BLOCKS format

to BLOCKS multiple alignment processor.

Motif 2 position-specific scoring matrix

Motif 2 position-specific probability matrix






Time 65.91 secs.

P N
MOTIF 3 width = 15 sites = 11 llr = 129 E-value = 3.1e+001

NAME	START	P-VALUE				SITES
			Simplified	A	:	:	:	:	:	:	3	:	4	:	:	:	5	5	1
			pos.-specific	C	:	1	:	1	:	4	3	a	:	:	4	a	3	4	3
			probability	G	1	9	4	:	:	4	:	:	:	5	1	:	:	1	:
			matrix	T	9	:	6	9	a	3	5	:	6	5	5	:	3	:	6
.
			bits	2.2
				2.0
				1.8
				1.6
			Information	1.3
			content	1.1
			(17.0 bits)	0.9
				0.7
				0.4
				0.2
				0.0
.
			Multilevel		`T`	`G`	`T`	`T`	`T`	`C`	`T`	`C`	`T`	`T`	`T`	`C`	`A`	`A`	`T`
			consensus				`G`			`G`	`A`		`A`	`G`	`C`		`C`	`C`	`C`
			sequence							`T`	`C`						`T`

.

CG4471:	100	1.13e-07	`CGTTGCAATC`		`T`	`G`	`T`	`T`	`T`	`C`	`T`	`C`	`T`	`T`	`T`	`C`	`C`	`C`	`T`	`GACTTCTTAT`
ImpE2:	125	1.93e-07	`CCTGTGTGCT`		`T`	`G`	`T`	`T`	`T`	`G`	`T`	`C`	`T`	`G`	`C`	`C`	`A`	`A`	`C`	`AAAGTCGACG`
CG7298:	68	3.01e-07	`AATGGAACAA`		`T`	`G`	`T`	`T`	`T`	`C`	`A`	`C`	`T`	`G`	`T`	`C`	`T`	`A`	`T`	`AATTCCTGAG`
CG2471:	91	3.69e-07	`TATTATTTTT`		`T`	`G`	`T`	`T`	`T`	`T`	`T`	`C`	`T`	`T`	`C`	`C`	`C`	`A`	`T`	`TTTATTTTTC`
CG6296:	126	1.00e-06	`TGCGCAGATT`		`T`	`G`	`G`	`T`	`T`	`C`	`T`	`C`	`A`	`T`	`T`	`C`	`A`	`A`	`C`	`ATTTGCAGAC`
CG6981:	94	1.69e-06	`ACAGGTTGGT`		`T`	`G`	`G`	`T`	`T`	`G`	`C`	`C`	`T`	`T`	`T`	`C`	`A`	`G`	`T`	`GATATTTTTC`
Hsp27:	121	3.59e-06	`GTTCCGTCCC`		`T`	`G`	`G`	`T`	`T`	`G`	`C`	`C`	`A`	`T`	`G`	`C`	`A`	`C`	`T`	`AGTGTGTGTG`
CG7874:	34	4.20e-06	`GGAACGTATC`		`T`	`G`	`G`	`C`	`T`	`G`	`C`	`C`	`T`	`G`	`T`	`C`	`T`	`A`	`T`	`CTGCCAGATA`
CG10433:	36	7.16e-06	`TCGCTCCGTC`		`T`	`C`	`T`	`T`	`T`	`C`	`T`	`C`	`T`	`T`	`T`	`C`	`T`	`C`	`C`	`GCGTCTAGTT`
Pu:	2	7.95e-06	`C`		`G`	`G`	`T`	`T`	`T`	`T`	`A`	`C`	`A`	`G`	`C`	`C`	`A`	`C`	`T`	`GGATGCCAAA`
CG5839:	14	7.95e-06	`GTAGACTAAT`		`T`	`G`	`T`	`T`	`T`	`T`	`A`	`C`	`A`	`G`	`C`	`C`	`C`	`A`	`A`	`TCATACTTCA`

Motif 3 block diagrams

Name

Lowest
p-value

Motifs

CG4471:

1.1e-07

ImpE2:

1.9e-07

CG7298:

3e-07

CG2471:

3.7e-07

CG6296:

1e-06

CG6981:

1.7e-06

Hsp27:

3.6e-06

CG7874:

4.2e-06

CG10433:

7.2e-06

Pu:

7.9e-06

CG5839:

7.9e-06

SCALE

\|	\|	\|	\|	\|	\|	\|	\|
1	25	50	75	100	125	150	175

Motif 3 in BLOCKS format

to BLOCKS multiple alignment processor.

Motif 3 position-specific scoring matrix

Motif 3 position-specific probability matrix






Time 95.21 secs.

P N
MOTIF 4 width = 34 sites = 3 llr = 100 E-value = 1.5e+003

NAME	START	P-VALUE				SITES
			Simplified	A	:	7	:	3	a	:	a	:	:	7	3	:	:	7	3	:	:	7	:	3	:	:	:	:	a	3	:	:	:	:	3	a	a	:
			pos.-specific	C	3	:	:	3	:	a	:	:	:	:	:	3	7	:	3	:	:	:	:	7	:	:	7	3	:	3	:	:	:	7	3	:	:	:
			probability	G	7	3	a	:	:	:	:	a	:	:	:	:	:	3	3	7	7	3	a	:	:	a	3	7	:	3	a	7	3	3	:	:	:	a
			matrix	T	:	:	:	3	:	:	:	:	a	3	7	7	3	:	:	3	3	:	:	:	a	:	:	:	:	:	:	3	7	:	3	:	:	:
.
			bits	2.2
				2.0
				1.8
				1.6
			Information	1.3
			content	1.1
			(47.9 bits)	0.9
				0.7
				0.4
				0.2
				0.0
.
			Multilevel		`G`	`A`	`G`	`A`	`A`	`C`	`A`	`G`	`T`	`A`	`T`	`T`	`C`	`A`	`A`	`G`	`G`	`A`	`G`	`C`	`T`	`G`	`C`	`G`	`A`	`A`	`G`	`G`	`T`	`C`	`A`	`A`	`A`	`G`
			consensus		`C`	`G`		`C`						`T`	`A`	`C`	`T`	`G`	`C`	`T`	`T`	`G`		`A`			`G`	`C`		`C`		`T`	`G`	`G`	`C`
			sequence					`T`											`G`											`G`					`T`

.

Adh:	161	6.33e-18	`GCCCCCACGA`		`G`	`A`	`G`	`A`	`A`	`C`	`A`	`G`	`T`	`A`	`T`	`T`	`T`	`A`	`A`	`G`	`G`	`A`	`G`	`C`	`T`	`G`	`C`	`G`	`A`	`A`	`G`	`G`	`T`	`C`	`C`	`A`	`A`	`G`	`TCATGCA`
CG5150:	157	1.81e-15	`TGGAACTACG`		`G`	`A`	`G`	`T`	`A`	`C`	`A`	`G`	`T`	`A`	`T`	`T`	`C`	`A`	`C`	`T`	`T`	`G`	`G`	`C`	`T`	`G`	`G`	`C`	`A`	`C`	`G`	`G`	`T`	`G`	`A`	`A`	`A`	`G`	`TATCTCGAAA`
CG7874:	73	4.77e-15	`TACGGACATA`		`C`	`G`	`G`	`C`	`A`	`C`	`A`	`G`	`T`	`T`	`A`	`C`	`C`	`G`	`G`	`G`	`G`	`A`	`G`	`A`	`T`	`G`	`C`	`G`	`A`	`G`	`G`	`T`	`G`	`C`	`T`	`A`	`A`	`G`	`CCCGATAAGG`

Motif 4 block diagrams

Name

Lowest
p-value

Motifs

Adh:

6.3e-18

CG5150:

1.8e-15

CG7874:

4.8e-15

SCALE

\|	\|	\|	\|	\|	\|	\|	\|
1	25	50	75	100	125	150	175

Motif 4 in BLOCKS format

to BLOCKS multiple alignment processor.

Motif 4 position-specific scoring matrix

Motif 4 position-specific probability matrix






Time 121.56 secs.

P N
MOTIF 5 width = 11 sites = 15 llr = 135 E-value = 3.7e+003

NAME	START	P-VALUE				SITES
			Simplified	A	:	9	:	1	:	:	9	3	7	1	1
			pos.-specific	C	a	1	4	6	:	7	:	7	1	3	5
			probability	G	:	:	5	:	:	:	:	:	:	6	:
			matrix	T	:	:	1	3	a	3	1	:	2	:	3
.
			bits	2.2
				2.0
				1.8
				1.6
			Information	1.3
			content	1.1
			(13.0 bits)	0.9
				0.7
				0.4
				0.2
				0.0
.
			Multilevel		`C`	`A`	`G`	`C`	`T`	`C`	`A`	`C`	`A`	`G`	`C`
			consensus				`C`	`T`		`T`		`A`	`T`	`C`	`T`
			sequence

.

CG13095:	15	1.87e-06	`ATCGAAATAT`		`C`	`A`	`C`	`C`	`T`	`C`	`A`	`C`	`A`	`C`	`C`	`TGAAATGATA`
glob1:	7	2.52e-06	`ATGACA`		`C`	`A`	`C`	`C`	`T`	`C`	`A`	`A`	`A`	`G`	`C`	`CCCCGTTGAA`
Chit:	79	3.14e-06	`CTTTGTTTAT`		`C`	`A`	`G`	`C`	`T`	`C`	`A`	`C`	`A`	`C`	`T`	`TTGCTCTCAA`
CG6640:	3	4.31e-06	`TT`		`C`	`A`	`C`	`C`	`T`	`C`	`A`	`C`	`T`	`G`	`C`	`GTGGAGAAGC`
CG11892:	124	5.79e-06	`GCTGTGGGAA`		`C`	`A`	`G`	`T`	`T`	`T`	`A`	`C`	`A`	`G`	`C`	`GTCTGAAATT`
CG7874:	182	1.34e-05	`TCTCGGCTAG`		`C`	`A`	`G`	`C`	`T`	`C`	`A`	`C`	`A`	`C`	`A`	`ACGATCGAC`
CG10433:	166	2.90e-05	`TAGCTGGGAT`		`C`	`A`	`G`	`T`	`T`	`C`	`A`	`A`	`C`	`G`	`C`	`AAACCGCGCC`
Cyp6a8:	57	3.27e-05	`CAATTCTTGC`		`C`	`A`	`T`	`C`	`T`	`C`	`A`	`C`	`A`	`G`	`T`	`GCTTGGCATT`
Pu:	141	3.92e-05	`GAGCTTTCCG`		`C`	`C`	`G`	`C`	`T`	`C`	`A`	`A`	`A`	`G`	`C`	`TCCGCTCGAG`
CG6296:	102	6.62e-05	`TTGATCTAAT`		`C`	`A`	`C`	`T`	`T`	`T`	`A`	`C`	`T`	`G`	`T`	`TGTTGCGCAG`
CG11459:	154	7.19e-05	`CTAGTTTTTC`		`C`	`A`	`G`	`C`	`T`	`C`	`A`	`A`	`A`	`A`	`T`	`CATTTACAAA`
CG4471:	181	7.19e-05	`GCATCTGCCG`		`C`	`A`	`G`	`A`	`T`	`T`	`A`	`C`	`A`	`C`	`C`	`ACTCTCGCGA`
Adh:	36	8.10e-05	`TTCTATTGGG`		`C`	`A`	`C`	`T`	`T`	`T`	`A`	`C`	`C`	`C`	`C`	`GCATTTGTTT`
ImpE2:	25	9.32e-05	`CAACTTGGTC`		`C`	`A`	`G`	`T`	`T`	`C`	`A`	`A`	`T`	`G`	`A`	`CACTCAGTGT`
CG4653:	174	1.11e-04	`CCGGTAGGCG`		`C`	`A`	`C`	`C`	`T`	`T`	`T`	`C`	`A`	`G`	`T`	`TGAGATCCCA`

Motif 5 block diagrams

Name

Lowest
p-value

Motifs

CG13095:

1.9e-06

glob1:

2.5e-06

Chit:

3.1e-06

CG6640:

4.3e-06

CG11892:

5.8e-06

CG7874:

1.3e-05

CG10433:

2.9e-05

Cyp6a8:

3.3e-05

Pu:

3.9e-05

CG6296:

6.6e-05

CG11459:

7.2e-05

CG4471:

7.2e-05

Adh:

8.1e-05

ImpE2:

9.3e-05

CG4653:

0.00011

SCALE

\|	\|	\|	\|	\|	\|	\|	\|
1	25	50	75	100	125	150	175

Motif 5 in BLOCKS format

to BLOCKS multiple alignment processor.

Motif 5 position-specific scoring matrix

Motif 5 position-specific probability matrix






Time 147.12 secs.

P N
MOTIF 6 width = 19 sites = 2 llr = 49 E-value = 3.8e+003

NAME	START	P-VALUE				SITES
			Simplified	A	:	:	5	:	5	:	a	:	:	:	:	:	a	a	:	:	:	:	:
			pos.-specific	C	:	a	:	:	5	:	:	5	:	:	a	:	:	:	5	a	:	:	:
			probability	G	a	:	:	a	:	a	:	5	:	a	:	a	:	:	:	:	:	a	a
			matrix	T	:	:	5	:	:	:	:	:	a	:	:	:	:	:	5	:	a	:	:
.
			bits	2.2
				2.0
				1.8
				1.6
			Information	1.3
			content	1.1
			(35.1 bits)	0.9
				0.7
				0.4
				0.2
				0.0
.
			Multilevel		`G`	`C`	`A`	`G`	`A`	`G`	`A`	`C`	`T`	`G`	`C`	`G`	`A`	`A`	`C`	`C`	`T`	`G`	`G`
			consensus				`T`		`C`			`G`							`T`
			sequence

.

CG8502:	123	5.68e-12	`TGAGAGCAGC`		`G`	`C`	`T`	`G`	`C`	`G`	`A`	`C`	`T`	`G`	`C`	`G`	`A`	`A`	`C`	`C`	`T`	`G`	`G`	`CACATTTCTA`
glob1:	123	2.31e-11	`ATCTTTGCCG`		`G`	`C`	`A`	`G`	`A`	`G`	`A`	`G`	`T`	`G`	`C`	`G`	`A`	`A`	`T`	`C`	`T`	`G`	`G`	`TGAGTTCGCT`

Motif 6 block diagrams

Name

Lowest
p-value

Motifs

CG8502:

5.7e-12

glob1:

2.3e-11

SCALE

\|	\|	\|	\|	\|	\|	\|	\|
1	25	50	75	100	125	150	175

Motif 6 in BLOCKS format

to BLOCKS multiple alignment processor.

Motif 6 position-specific scoring matrix

Motif 6 position-specific probability matrix






Time 169.76 secs.

P N
MOTIF 7 width = 11 sites = 8 llr = 85 E-value = 5.0e+003

NAME	START	P-VALUE				SITES
			Simplified	A	:	a	3	3	a	3	:	:	:	:	:
			pos.-specific	C	:	:	6	4	:	:	6	:	3	:	:
			probability	G	a	:	1	3	:	5	:	:	5	a	a
			matrix	T	:	:	:	1	:	3	4	a	3	:	:
.
			bits	2.2
				2.0
				1.8
				1.6
			Information	1.3
			content	1.1
			(15.3 bits)	0.9
				0.7
				0.4
				0.2
				0.0
.
			Multilevel		`G`	`A`	`C`	`C`	`A`	`G`	`C`	`T`	`G`	`G`	`G`
			consensus				`A`	`A`		`A`	`T`		`C`
			sequence					`G`		`T`			`T`

.

Chit:	163	3.76e-07	`CTCAACCTTC`		`G`	`A`	`C`	`A`	`A`	`G`	`C`	`T`	`G`	`G`	`G`	`ATGCGAAGTA`
CG6981:	69	2.26e-06	`TTGGCGACAC`		`G`	`A`	`C`	`G`	`A`	`A`	`C`	`T`	`G`	`G`	`G`	`CGCGACAGGT`
ImpE2:	12	6.11e-06	`TGCAATGCAA`		`G`	`A`	`C`	`C`	`A`	`A`	`C`	`T`	`T`	`G`	`G`	`TCCAGTTCAA`
CG13095:	124	9.10e-06	`CTTATCACCG`		`G`	`A`	`C`	`A`	`A`	`T`	`T`	`T`	`G`	`G`	`G`	`CTGTGATTGA`
CG8502:	17	9.97e-06	`AATCCAAAGG`		`G`	`A`	`G`	`C`	`A`	`G`	`C`	`T`	`C`	`G`	`G`	`GGTACAAGTC`
CG7298:	92	9.97e-06	`TAATTCCTGA`		`G`	`A`	`A`	`T`	`A`	`G`	`C`	`T`	`G`	`G`	`G`	`CTTATCTACA`
CG2471:	15	1.20e-05	`CAAGGCTAAT`		`G`	`A`	`C`	`C`	`A`	`T`	`T`	`T`	`C`	`G`	`G`	`AGCCAAGGAA`
Cyt-b5:	13	2.05e-05	`CGGCCACCTT`		`G`	`A`	`A`	`G`	`A`	`G`	`T`	`T`	`T`	`G`	`G`	`CTAAAATTAC`

Motif 7 block diagrams

Name

Lowest
p-value

Motifs

Chit:

3.8e-07

CG6981:

2.3e-06

ImpE2:

6.1e-06

CG13095:

9.1e-06

CG8502:

1e-05

CG7298:

1e-05

CG2471:

1.2e-05

Cyt-b5:

2.1e-05

SCALE

\|	\|	\|	\|	\|	\|	\|	\|
1	25	50	75	100	125	150	175

Motif 7 in BLOCKS format

to BLOCKS multiple alignment processor.

Motif 7 position-specific scoring matrix

Motif 7 position-specific probability matrix






Time 192.72 secs.

P N
MOTIF 8 width = 8 sites = 5 llr = 51 E-value = 5.7e+003

NAME	START	P-VALUE				SITES
			Simplified	A	:	:	:	a	:	a	2	:
			pos.-specific	C	:	4	:	:	:	:	:	:
			probability	G	:	6	a	:	a	:	8	a
			matrix	T	a	:	:	:	:	:	:	:
.
			bits	2.2
				2.0
				1.8
				1.6
			Information	1.3
			content	1.1
			(14.8 bits)	0.9
				0.7
				0.4
				0.2
				0.0
.
			Multilevel		`T`	`G`	`G`	`A`	`G`	`A`	`G`	`G`
			consensus			`C`					`A`
			sequence

.

CG18522:	119	9.68e-06	`AACTTGAACT`		`T`	`G`	`G`	`A`	`G`	`A`	`G`	`G`	`CGGGAGTCTC`
CG3868:	94	9.68e-06	`GAATCTAATC`		`T`	`G`	`G`	`A`	`G`	`A`	`G`	`G`	`TCTACCGCGA`
Pu:	157	2.01e-05	`AAAGCTCCGC`		`T`	`C`	`G`	`A`	`G`	`A`	`G`	`G`	`TTCTTATTGC`
CG10433:	188	2.01e-05	`AACCGCGCCA`		`T`	`C`	`G`	`A`	`G`	`A`	`G`	`G`	`AAGTGG`
CG6640:	15	3.25e-05	`CCTCACTGCG`		`T`	`G`	`G`	`A`	`G`	`A`	`A`	`G`	`CATATTGATT`

Motif 8 block diagrams

Name

Lowest
p-value

Motifs

CG18522:

9.7e-06

CG3868:

9.7e-06

Pu:

2e-05

CG10433:

2e-05

CG6640:

3.3e-05

SCALE

\|	\|	\|	\|	\|	\|	\|	\|
1	25	50	75	100	125	150	175

Motif 8 in BLOCKS format

to BLOCKS multiple alignment processor.

Motif 8 position-specific scoring matrix

Motif 8 position-specific probability matrix






Time 214.81 secs.

P N
MOTIF 9 width = 11 sites = 2 llr = 31 E-value = 6.1e+003

NAME	START	P-VALUE				SITES
			Simplified	A	a	:	:	:	a	5	:	:	:	:	:
			pos.-specific	C	:	a	:	:	:	:	:	:	a	a	a
			probability	G	:	:	a	a	:	5	a	a	:	:	:
			matrix	T	:	:	:	:	:	:	:	:	:	:	:
.
			bits	2.2
				2.0
				1.8
				1.6
			Information	1.3
			content	1.1
			(22.1 bits)	0.9
				0.7
				0.4
				0.2
				0.0
.
			Multilevel		`A`	`C`	`G`	`G`	`A`	`A`	`G`	`G`	`C`	`C`	`C`
			consensus							`G`
			sequence

.

CG11892:	75	1.01e-07	`CACAGTTGTG`		`A`	`C`	`G`	`G`	`A`	`G`	`G`	`G`	`C`	`C`	`C`	`TCGGATATAT`
CG18522:	173	2.31e-07	`ATATGGGGCT`		`A`	`C`	`G`	`G`	`A`	`A`	`G`	`G`	`C`	`C`	`C`	`ATACCCTCAG`

Motif 9 block diagrams

Name

Lowest
p-value

Motifs

CG11892:

1e-07

CG18522:

2.3e-07

SCALE

\|	\|	\|	\|	\|	\|	\|	\|
1	25	50	75	100	125	150	175

Motif 9 in BLOCKS format

to BLOCKS multiple alignment processor.

Motif 9 position-specific scoring matrix

Motif 9 position-specific probability matrix






Time 236.43 secs.

P N
MOTIF 10 width = 20 sites = 4 llr = 78 E-value = 6.9e+003

NAME	START	P-VALUE				SITES
			Simplified	A	:	a	3	a	:	:	5	:	:	:	:	:	3	8	3	3	:	8	:	:
			pos.-specific	C	8	:	5	:	:	a	:	3	8	8	a	a	:	:	5	:	:	:	:	8
			probability	G	3	:	3	:	8	:	3	:	3	3	:	:	8	:	:	8	a	:	:	3
			matrix	T	:	:	:	:	3	:	3	8	:	:	:	:	:	3	3	:	:	3	a	:
.
			bits	2.2
				2.0
				1.8
				1.6
			Information	1.3
			content	1.1
			(28.0 bits)	0.9
				0.7
				0.4
				0.2
				0.0
.
			Multilevel		`C`	`A`	`C`	`A`	`G`	`C`	`A`	`T`	`C`	`C`	`C`	`C`	`G`	`A`	`C`	`G`	`G`	`A`	`T`	`C`
			consensus		`G`		`A`		`T`		`G`	`C`	`G`	`G`			`A`	`T`	`A`	`A`		`T`		`G`
			sequence				`G`				`T`								`T`

.

CG7874:	128	5.09e-13	`CCGCCACCAT`		`C`	`A`	`C`	`A`	`G`	`C`	`A`	`T`	`C`	`C`	`C`	`C`	`G`	`A`	`C`	`G`	`G`	`A`	`T`	`C`	`GGCCAGCTAT`
CG5171:	135	2.46e-09	`TTTGTTTGAG`		`G`	`A`	`A`	`A`	`T`	`C`	`A`	`C`	`C`	`G`	`C`	`C`	`G`	`A`	`C`	`G`	`G`	`A`	`T`	`C`	`ATTATAAAAC`
mfas:	27	2.61e-09	`TTACCGATTT`		`C`	`A`	`C`	`A`	`G`	`C`	`G`	`T`	`C`	`C`	`C`	`C`	`A`	`A`	`A`	`A`	`G`	`T`	`T`	`C`	`TCTCTGTTCT`
Pu:	92	2.99e-09	`AAACTCGGAG`		`C`	`A`	`G`	`A`	`G`	`C`	`T`	`T`	`G`	`C`	`C`	`C`	`G`	`T`	`T`	`G`	`G`	`A`	`T`	`G`	`GTCTGTTGGC`

Motif 10 block diagrams

Name

Lowest
p-value

Motifs

CG7874:

5.1e-13

CG5171:

2.5e-09

mfas:

2.6e-09

Pu:

3e-09

SCALE

\|	\|	\|	\|	\|	\|	\|	\|
1	25	50	75	100	125	150	175

Motif 10 in BLOCKS format

to BLOCKS multiple alignment processor.

Motif 10 position-specific scoring matrix

Motif 10 position-specific probability matrix






Time 257.90 secs.

P N
SUMMARY OF MOTIFS

Combined block diagrams: non-overlapping sites with p-value < 0.0001

Name

Combined
p-value

Motifs

CG6640:

3.16e-05

CG5150:

1.78e-08

CG4653:

9.84e-04

CG3868:

1.06e-02

CG7874:

1.95e-21

CG13095:

8.63e-05

Cyp6a8:

3.52e-01

Cyt-b5:

2.28e-01

CG5839:

1.67e-05

CG5932:

1.59e-01

Chit:

2.44e-04

CG6981:

2.46e-03

CG7298:

5.58e-06

CG6296:

1.92e-03

CG18522:

8.22e-06

CG10433:

2.67e-03

mfas:

2.57e-10

glob1:

7.99e-09

Adh:

1.42e-14

Pu:

1.11e-10

Hsp27:

4.66e-04

CG4471:

2.69e-06

CG5171:

9.70e-09

CG2471:

4.85e-09

CG8502:

2.62e-10

CG11892:

7.10e-06

ImpE2:

2.51e-11

CG9021:

7.25e-02

CG11459:

4.44e-01

SCALE

\|	\|	\|	\|	\|	\|	\|	\|	\|
1	25	50	75	100	125	150	175	200

Motif summary in machine readable format.

Stopped because nmotifs = 10 reached.


CPU: monstre2

EXPLANATION OF MEME RESULTS

The MEME results consist of:

The version of MEME and the date it was released.
The reference to cite if you use MEME in your research.
A description of the sequences you submitted (the "training set") showing the name, "weight" and length of each sequence.
The command line summary detailing the parameters with which you ran MEME.
Information on each of the motifs MEME discovered, including:
1. A summary line showing the width, number of occurrences, log likelihood ratio and statistical significance of the motif.
2. A simplified position-specific probability matrix.
3. A diagram showing the degree of conservation at each motif position.
4. A multilevel consensus sequence showing the most conserved letter(s) at each motif position.
5. The occurrences of the motif sorted by p-value and aligned with each other.
6. Block diagrams of the occurrences of the motif within each sequence in the training set.
7. The motif in BLOCKS or FASTA format.
8. A position-specific scoring matrix (PSSM) for use by the MAST database search program.
9. The position specific probability matrix (PSPM) describing the motif.
A summary of motifs showing an optimized (non-overlapping) tiling of all of the motifs onto each of the sequences in the training set.
The reason why MEME stopped and the name of the CPU on which it ran.
This explanation of how to interpret MEME results.

MOTIFS

For each motif that it discovers in the training set, MEME prints the following information:

Summary Line
This line gives the width (`width'), number of occurrences in the training set (`sites'), log likelihood ratio (`llr') and E-value of the motif. Each motif describes a pattern of a fixed width--no gaps are allowed in MEME motifs. MEME numbers the motifs consecutively from one as it finds them. MEME usually finds the most statistically significant (low E-value) motifs first. The statistical significance of a motif is based on its log likelihood ratio, its width and number of occurrences, the background letter frequencies (given in the command line summary), and the size of the training set. The E-value is an estimate of the expected number of motifs with the given log likelihood ratio (or higher), and with the same width and number of occurrences, that one would find in a similarly sized set of random sequences. (In random sequences each position is independent with letters chosen according to the background letter frequencies.) The log likelihood ratio is the logarithm of the ratio of the probability of the occurrences of the motif given the motif model (likelihood given the motif) versus their probability given the background model (likelihood given the null model). (Normally the background model is a 0-order Markov model using the background letter frequencies, but higher order Markov models may be specified via the -bfile option to MEME.) Clicking on the buttons to the left of the motif summary line takes you to the previous motif (P) or next motif (N).
Simplified Position-Specific Probability Matrix
MEME motifs are represented by position-specific probability matrices that specify the probability of each possible letter appearing at each possible position in an occurrence of the motif. In order to make it easier to see which letters are most likely in each of the columns of the motif, the simplified motif shows the letter probabilities multiplied by 10 rounded to the nearest integer. Zeros are replaced by ":" (the colon) for readability.

Information Content Diagram

The information content diagram provides an idea of which positions in the motif are most highly conserved. Each column (position) in a motif can be characterized by the amount of information it contains (measured in bits). Highly conserved positions in the motif have high information; positions where all letters are equally likely have low information. (The information content is relative to the background letter frequencies which are given in the command line summary section.) The diagram is printed so that each column lines up with the same column in the simplified position-specific probability matrix above it. Columns in the information content diagram are colored according to the majority category of the letters occurring in that column of the alignment. If no letter category has frequency above 0.5, the column in the diagram is colored black. For DNA sequences, the letter categories contain one letter each. For proteins, the categories are based on the biochemical properties of the various amino acids. The categories and their colors are:

NUCLEIC ACIDS	COLOR
A	RED
C	BLUE
G	ORANGE
T	GREEN

AMINO ACIDS	COLOR
ACFILVM	BLUE
NQST	GREEN
DE	MAGENTA
KR	RED
H	PINK
G	ORANGE
P	YELLOW
Y	TURQUOISE

Summing the information content for each position in the motif gives the total information content of the motif (shown in parentheses to the left of the diagram). The total information content is approximately equal to the log likelihood ratio divided by the number of occurrences times ln(2). The total information content gives a measure of the usefulness of the motif for database searches. For a motif to be useful for database searches, it must as a rule contain at least log_2(N) bits of information where N is the number of sequences in the database being searched. For example, to effectively search a database containing 100,000 sequences for occurrences of a single motif, the motif should have an IC of at least 16.6 bits. Motifs with lower information content are still useful when a family of sequences shares more than one motif since they can be combined in multiple motif searches (using MAST).

Multilevel Consensus Sequence
The multilevel consensus sequence corresponding to the motif is an aid in remembering and understanding the motif. It is calculated from the motif position-specific probability matrix as follows. Separately for each column of the motif, the letters in the alphabet are sorted in decreasing order by the probability with which they are expected to occur in that position of motif occurrences. The sorted letters are then printed vertically with the most probable letter on top. Only letters with probabilities of 0.2 or higher at that position in the motif are printed. As an example, the multilevel consensus sequence of motif 1 in the sample output is:
```
Multilevel TTATGTGAACGACGTCACACT
consensus  AA  T A G A GA     AA
sequence           T C TT     T
```
This multilevel consensus sequence says several things about the motif. First, the most likely form of the motif can be read from the top line as TTATGTGAACGACGTCACACT. Second, that only letter A has probability more than 0.2 in position 3 of the motif, both T and A have probability greater than 0.2 in position 1, etc. Third, a rough approximation of the motif can be made by converting the multilevel consensus sequence into the Prosite signature
Occurrences of the Motif
MEME displays the occurrences (sites) of the motif in the training set. The sites are shown aligned with each other, and the ten sequence positions preceding and following each site are also shown. Each site is identified by the name of the sequence where it occurs, the strand (if both strands of DNA sequences are being used), and the position in the sequence where the site begins. When the DNA strand is specified, `+' means the sequence in the training set, and `-' means the reverse complement of the training set sequence. (For `-' strands, the `start' position is actually the position on the positive strand where the site ends.) The sites are listed in order of increasing statistical significance (p-value). The p-value of a site is computed from the the match score of the site with the position specific scoring matrix for the motif. The p-value gives the probability of a random string (generated from the background letter frequencies) having the same match score or higher. (This is referred to as the position p-value by the MAST algorithm.)
Block Diagrams of Motif Occurrences
The occurrences of the motif in the training set sequences are shown with MAST-style block diagrams. One diagram is printed for each sequence showing all the occurrences of the motif in that sequence. The sequences are sorted by the lowest p-value among all occurrences of the motif in a given sequence. (The p-value of an occurrence is the probability of a single random subsequence the length of the motif, generated according to the 0-order background model, having a score at least as high as the score of the occurrence.) When the DNA strand is specified, `+' means the motif appears from left to right on the sequence, and `-' means the motif appears from right to left on the complementary strand. A sequence position scale is shown at the end of each table of block diagrams. Very long sequences are shown with thick lines connecting the motifs and are not drawn to scale.
Motif in BLOCKS format or FASTA format
For use with BLOCKS tools, MEME prints the occurrences of the motif in BLOCKS format.
You can convert these blocks to PSSMs (position-specific scoring matrices), LOGOS (color representations of the motifs), phylogeny trees and search them against a database of other blocks by pasting everything from the "BL" line to the "//" line (inclusive) into the Multiple Alignment Processor. If you include the -print_fasta switch on the command line, MEME prints the motif sites in FASTA format instead of BLOCKS format.
Position-Specific Scoring Matrix
The position-specific scoring matrix corresponding to the motif is printed for use by database search programs such as MAST. This matrix is a log-odds matrix calculated by taking the log (base 2) of the ratio p/f at each position in the motif where p is the probability of a particular letter at that position in the motif, and f is the background frequency of the letter (given in the command line summary section.) This is the same matrix that is used above in computing the p-values of the occurrences of the motif in the Occurrences of the Motif and Block Diagrams of Motif Occurrences sections. The scoring matrix is printed "sideways"--columns correspond to the letters in the alphabet (in the same order as shown in the simplified motif) and rows corresponding to the positions of the motif, position one first. The scoring matrix is preceded by a line starting with "log-odds matrix:" and containing the length of the alphabet, width of the motif, number of characters in the training set and the scoring threshold used in the list of possible motif examples.
Position-Specific Probability Matrix
The motif itself is a position-specific probability matrix giving, for each position in the pattern, the probabilities of each possible letter occurring there. The probability matrix is printed "sideways"--columns correspond to the letters in the alphabet (in the same order as shown in the simplified motif) and rows corresponding to the positions of the motif, position one first. The motif is preceded by a line starting with "letter-probability matrix:" and containing the length of the alphabet, width of the motif and number of characters in the training set.

Go to top

P N MOTIF 1 width = 16 sites = 10 llr = 131 E-value = 1.3e-001

P N MOTIF 2 width = 13 sites = 20 llr = 185 E-value = 1.1e-002

P N MOTIF 3 width = 15 sites = 11 llr = 129 E-value = 3.1e+001

P N MOTIF 4 width = 34 sites = 3 llr = 100 E-value = 1.5e+003

P N MOTIF 5 width = 11 sites = 15 llr = 135 E-value = 3.7e+003

P N MOTIF 6 width = 19 sites = 2 llr = 49 E-value = 3.8e+003

P N MOTIF 7 width = 11 sites = 8 llr = 85 E-value = 5.0e+003

P N MOTIF 8 width = 8 sites = 5 llr = 51 E-value = 5.7e+003

P N MOTIF 9 width = 11 sites = 2 llr = 31 E-value = 6.1e+003

P N MOTIF 10 width = 20 sites = 4 llr = 78 E-value = 6.9e+003

P NSUMMARY OF MOTIFS

Combined block diagrams: non-overlapping sites with p-value < 0.0001

Stopped because nmotifs = 10 reached.

EXPLANATION OF MEME RESULTS

The MEME results consist of:

Summary Line

Simplified Position-Specific Probability Matrix

Information Content Diagram

Multilevel Consensus Sequence

Occurrences of the Motif

Block Diagrams of Motif Occurrences

Motif in BLOCKS format or FASTA format

Position-Specific Scoring Matrix

Position-Specific Probability Matrix

P N
MOTIF 1 width = 16 sites = 10 llr = 131 E-value = 1.3e-001

P N
MOTIF 2 width = 13 sites = 20 llr = 185 E-value = 1.1e-002

P N
MOTIF 3 width = 15 sites = 11 llr = 129 E-value = 3.1e+001

P N
MOTIF 4 width = 34 sites = 3 llr = 100 E-value = 1.5e+003

P N
MOTIF 5 width = 11 sites = 15 llr = 135 E-value = 3.7e+003

P N
MOTIF 6 width = 19 sites = 2 llr = 49 E-value = 3.8e+003

P N
MOTIF 7 width = 11 sites = 8 llr = 85 E-value = 5.0e+003

P N
MOTIF 8 width = 8 sites = 5 llr = 51 E-value = 5.7e+003

P N
MOTIF 9 width = 11 sites = 2 llr = 31 E-value = 6.1e+003

P N
MOTIF 10 width = 20 sites = 4 llr = 78 E-value = 6.9e+003

P N
SUMMARY OF MOTIFS