MEME

Search sequence databases for the best combined matches with these motifs using MAST.
Search sequence databases for all matches with these motifs using FIMO.
Submit these motifs to BLOCKS multiple alignment processor.

MEME - Motif discovery tool

MEME version 4.3.0 (Release date: Sat Sep 26 01:51:56 PDT 2009)

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= memeout/meme.fasta
ALPHABET= ACDEFGHIKLMNPQRSTVWY
Sequence name	Weight	Length	Sequence name	Weight	Length
Bacillus_subtilis_strain_168	1.0000	236	Natronobacterium	1.0000	225
Natrialba	1.0000	219	Sulfolobus	1.0000	220
Pyrobaculum	1.0000	212	Microlunatus	1.0000	235
Mycobacterium	1.0000	227	Bifidobacterium	1.0000	267
Scardovia	1.0000	226	Zunongwangia	1.0000	215
Bacteroides	1.0000	215	Stanieria	1.0000	222
Chloroflexus	1.0000	218	Leptospira	1.0000	231
Isosphaera	1.0000	233	Candidatus	1.0000	222
Helicobacter	1.0000	218	Octadecabacter	1.0000	228
Acinetobacter	1.0000	224	Comamonas	1.0000	223
Achromobacter	1.0000	220	Acanthamoeba	1.0000	601
Magnaporthe	1.0000	509	Torulaspora	1.0000	552
Vanderwaltozyma	1.0000	549	Gibberella	1.0000	510
Penicillium	1.0000	510	Tuber	1.0000	510
Ustilago	1.0000	642	Verticillium	1.0000	368
Schizosaccharomyces	1.0000	551	Coprinopsis	1.0000	556
Puccinia	1.0000	515	Oryza	1.0000	609
Solanum	1.0000	603	Physcomitrella	1.0000	546
Thalassiosira	1.0000	218	Phaeodactylum	1.0000	226
Perkinsus	1.0000	504	Naegleria	1.0000	538

COMMAND LINE SUMMARY

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

command: meme memeout/meme.fasta -mod zoops -nmotifs 3 -prior dirichlet -nostatus -protein -oc memeout/

model:	mod=	zoops	nmotifs=	3	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=	40	wnsites=	0.8
theta:	prob=	1	spmap=	pam	spfuzz=	120
em:	prior=	dirichlet	b=	0.01	maxiter=	50
	distance=	1e-05
data:	n=	14353	N=	40
sample:	seed=	0	seqfrac=	1
Letter frequencies in dataset:
A 0.106 C 0.017 D 0.051 E 0.069 F 0.041 G 0.061 H 0.028 I 0.052 K 0.048 L 0.099 M 0.024 N 0.031 P 0.042 Q 0.040 R 0.042 S 0.065 T 0.057 V 0.067 W 0.017 Y 0.043
Background letter frequencies (from dataset with add-one prior applied): A 0.106 C 0.017 D 0.051 E 0.069 F 0.041 G 0.061 H 0.029 I 0.052 K 0.048 L 0.099 M 0.024 N 0.031 P 0.042 Q 0.040 R 0.042 S 0.065 T 0.057 V 0.067 W 0.017 Y 0.043

MOTIF 1 width = 46 sites = 40 llr = 2760 E-value = 8.9e-702

SEQUENCE LOGO

PNG LOGOS require CONVERT from ImageMagick; see MEME installation guide

Information Content
108.3 (bits)

Relative Entropy
99.5 (bits)

Download LOGO
Without SSC:[EPS][PNG]
With SSC:[EPS][PNG]

NAME	START	P-VALUE	SITES
Verticillium	143	5.01e-40	`RPDVRDVWKT`	`FVYHPFVMGLGNGTLPLASFKGYLVQDYLYLVHFARANALASYKAT`	`SMEDIVAGAT`
Gibberella	308	6.65e-39	`RPDVRDVWKE`	`FVYHPFVMAMGDGTLPLESFKGYIIQDYLYLIHFSRANALAAYKAQ`	`NVEDISRATQ`
Magnaporthe	309	1.50e-37	`RPDVAPIWER`	`YIKHPFVLGLGDGTLPLSSFKSYLVQDYLYLIQFARANALASYKSK`	`NLADISAAAA`
Octadecabacter	23	5.01e-36	`RAGCPEPWRA`	`YTHHPFVQGLGDGTLPQASFLHYLVQDYVFLVHFARAWSLGVVKAE`	`TLDEMKICAS`
Pyrobaculum	16	5.87e-36	`RAWADPIWQE`	`ILRHPFVVEVYRGSLPMDKFRYYLLQDYNYLVNFAKALSLAAARAP`	`GVEQMKTALE`
Tuber	311	8.03e-36	`HPSVKPIWHA`	`YTHHPFTTQLGAGTLRRKSFQNYIVQDYLFLIQFARANSLAGYKST`	`STAAISKTAS`
Comamonas	17	1.28e-35	`WNANQALFQS`	`TLELPFNQELASGTLSRERFRHYMIQDAHYLVAYGRALAVTAAKSD`	`NAEGVVQFAN`
Helicobacter	17	1.28e-35	`YQNAQSIWGD`	`CISHPFVQGIGRGTLERDKFRFYIIQDYLFLLEYAKVFALGVVKAY`	`DEAVMREFSN`
Achromobacter	17	2.36e-35	`WARNAALYET`	`TRDMPFNQELASGQLDENAFKHYMIQDAHYLVAFGRALAIASAKAD`	`DADGVVQFAD`
Acanthamoeba	388	3.71e-35	`WTSIEGVFEK`	`IIAHPFIQGLTSGELPEDVFKFYIVQDMHYLREYSRILALLAAKAD`	`CNAGMLLFLE`
Microlunatus	32	3.71e-35	`WESGAATYAE`	`ILRHPFLTGLTDGTLSKDAFRYFVIQDSHYLRAYARALSLVSARAT`	`DPEAVALFAG`
Chloroflexus	18	6.72e-35	`WSSITDIYQA`	`IIYHPFNEELARGILPREKFAFYMQQDALYLADFARALAIMAGRAP`	`DEAAIIQFVR`
Torulaspora	348	2.13e-34	`HPKVKPHWES`	`YVHHDFVKQVAQGTLKREKFRFFIEQDYAYLIDYARVHCIAGSKAP`	`TLEDIEKELN`
Sulfolobus	18	2.13e-34	`WNSIQDIYSS`	`ILKHPFILELVEGTLSRNKFEYYIIQDYLYLREFSKALALLSAKAE`	`DEENALLFAI`
Penicillium	307	4.31e-34	`RQDVQQVWHQ`	`FTHHEFVEGMGQGTLPIERFKAYLVQDYLYLVQFARSNALASYKAK`	`SMSSIAASAE`
Natrialba	17	5.68e-34	`LDAGADIWDA`	`QKRHPFVTELAAGTLEEASFRHWVKQDYRYLLDYARLFSIAGTKAR`	`DEATMTHLLG`
Mycobacterium	24	1.47e-33	`WESIEPTFTA`	`ILEHPFLTGLTDGSLDDAAFAHYVAQDVHYLRDYARALAIVGAKAP`	`TLADTAMFSR`
Stanieria	17	3.70e-33	`WNKIDNIYNS`	`ILIMPFVQELKSGLLKREIFQHYMIQDAIYLGEFARVLAIISAKAP`	`EPDLQLQFAN`
Acinetobacter	17	4.21e-33	`WQRNLNLYQK`	`ILDLPFNQELANGTLDKEAFCHYVIQDAQYLVAYGRVLAVAAAKAF`	`DADDIMQFSD`
Natronobacterium	17	6.40e-32	`LEDGESVWEA`	`QKRHPFVRELAAGSLDEEAFLTWVRQDYRYLLDYARVFAIGGAKAR`	`DEETMTHLLG`
Zunongwangia	17	9.10e-32	`WQTIRPIYNR`	`IIEMPFITELMDGDLALEKFQFYMSQDSAYLENFGRALALIAARAY`	`DTENTLAFLQ`
Puccinia	307	1.02e-31	`AATGHSWDVF`	`TSSHPFLQGVRDGTLPIRCFQHFLRQDYIFLTHYARIHSLMALKSD`	`EMNEIEATAT`
Bacillus_subtilis_strain_168	17	2.56e-31	`RSAAAEWWEG`	`SFVHPFVQGIGDGTLPIDRFKYYVLQDSYYLTHFAKVQSFGAAYAK`	`DLYTTGRMAS`
Vanderwaltozyma	346	8.59e-31	`HPAVAPHWNS`	`YVNHDFVKQVAQGTIAPKKFQFFIEQDYSYLVDYGKVYCIAASKAP`	`NLDEVEKALA`
Perkinsus	28	1.07e-30	`WKENQDLALM`	`SLHNPFVQGLGDGTLDPEAFKTYLAQDTLYLNGYIRALSYCISKSD`	`VTATGGELLT`
Schizosaccharomyces	339	1.32e-30	`HPQVVPAWKE`	`YINHKFTNMLAKGTLPLPAFQDYLKQDYLYLVNFARAYSLKGYKEN`	`TFPNILEAAQ`
Bifidobacterium	52	2.47e-30	`RQAADPVWEE`	`GYRQPFIQELGAGTLDRSKFAFYLMQDELYLGAYAKVHALAVTKTD`	`DREVMAWMAS`
Ustilago	429	2.74e-30	`ISNSIDVWKA`	`FVQHAFVLGLADGSLPRESFEWFMKQDYLFLRHYARIWAQAAASPD`	`NTFQEVSTLA`
Isosphaera	30	3.04e-30	`WRAIEPINDA`	`LKSHPFLVGLTTGDLERDRFAFYIVQDALYLVEFAQALSALAAFAP`	`DPADSAMFNR`
Bacteroides	17	6.83e-30	`WNSIQPIFKK`	`TISHPFIQELMNGNLPLEKFNFYIQQDSIYLSSYGKLLTRLASKFK`	`DPEHSRAFIH`
Candidatus	17	8.33e-30	`RKLAQPVWDA`	`QLTHPFVVALGKGTLAERKFRYYILQDARFLADLARVFAAGSLRAP`	`DSESALRFAK`
Coprinopsis	347	2.20e-29	`IRNTKEVWKG`	`YVEHEFVKQLGKGVLAKEKFVHFIKQDYHYLKYYARAYGLLASKSS`	`TFPEIHAAAT`
Leptospira	24	5.09e-29	`KSSAKQSFGA`	`SFSHPFVLALADGSLDPKVFRFYQIQDAKYLESFSDACAILSTRVT`	`DPEDKLWLVD`
Naegleria	32	2.30e-28	`SMMEDQIMHG`	`ILGLPFIREMMEGTLNQNIFRFYIIQDTIYLQTFSKIFVLAATKID`	`LLKEVNGEAQ`
Physcomitrella	25	5.79e-28	`WASCKREAYA`	`SLYHPFVVALAAGTLPKQTFQRYMAQDAYFLEAFKNAYQLAMETTT`	`DEEAKAIIES`
Phaeodactylum	16	8.03e-28	`WDVAAPLIAV`	`TEKHPFLVAMVDGSLDMDSFQYYVVQDALYLHDFAYALRRLGDNDG`	`ISRQDSERLH`
Solanum	71	1.65e-27	`WINSKKESTF`	`ALYTPFVVCLASGTLNIDTFRHYIAQDVHFLKAFVQGCEAAEVCTD`	`DDDAKVGISE`
Scardovia	23	1.92e-27	`HQAAAPLWQK`	`SIHHPFITELASGKLPLATFRYYLTQDTKYLTAFEKLHEKTAALLP`	`QSQGALLLKL`
Oryza	72	8.63e-27	`IAAASREAAF`	`AAYTPFLVSLAAGALRLDSFRQYIAQDAYFLHAFARAYEMAEECAD`	`DDDDKATIVV`
Thalassiosira	19	1.39e-24	`RTSVAARWER`	`VISHKFTVELAKGEIDRDVLKKYLVQDFRFLDAFVILLASVIANAR`	`SLEDRIPGCQ`

Motif 1 block diagrams

Name

Lowest
p-value

Motifs

Verticillium

5.01e-40

Gibberella

6.65e-39

Magnaporthe

1.50e-37

Octadecabacter

5.01e-36

Pyrobaculum

5.87e-36

Tuber

8.03e-36

Comamonas

1.28e-35

Helicobacter

1.28e-35

Achromobacter

2.36e-35

Acanthamoeba

3.71e-35

Microlunatus

3.71e-35

Chloroflexus

6.72e-35

Torulaspora

2.13e-34

Sulfolobus

2.13e-34

Penicillium

4.31e-34

Natrialba

5.68e-34

Mycobacterium

1.47e-33

Stanieria

3.70e-33

Acinetobacter

4.21e-33

Natronobacterium

6.40e-32

Zunongwangia

9.10e-32

Puccinia

1.02e-31

Bacillus_subtilis_strain_168

2.56e-31

Vanderwaltozyma

8.59e-31

Perkinsus

1.07e-30

Schizosaccharomyces

1.32e-30

Bifidobacterium

2.47e-30

Ustilago

2.74e-30

Isosphaera

3.04e-30

Bacteroides

6.83e-30

Candidatus

8.33e-30

Coprinopsis

2.20e-29

Leptospira

5.09e-29

Naegleria

2.30e-28

Physcomitrella

5.79e-28

Phaeodactylum

8.03e-28

Solanum

1.65e-27

Scardovia

1.92e-27

Oryza

8.63e-27

Thalassiosira

1.39e-24

SCALE

\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|
1	25	50	75	100	125	150	175	200	225	250	275	300	325	350	375	400	425	450	475	500	525	550	575	600	625

Motif 1 in BLOCKS format

to BLOCKS multiple alignment processor.

Motif 1 position-specific scoring matrix

Scan sequence databases for the best match in each sequence using MAST.

Motif 1 position-specific probability matrix

Scan sequence databases for all matches with this motif using FIMO.

Motif 1 regular-expression

I[LI]xHPFV[QV][EG]L[AG]DGTL[PD][RL][ED][KS]F[RKQ][FH]Y[IL][IV]QD[YA]L[YF]LV[AD][FY]A[RK][AV]L[AS][LI]A[AS]AKA[DP]

Time 17.81 secs.

MOTIF 2 width = 37 sites = 31 llr = 1692 E-value = 6.0e-354

SEQUENCE LOGO

Information Content
84.8 (bits)

Relative Entropy
78.7 (bits)

Download LOGO
Without SSC:[EPS][PNG]
With SSC:[EPS][PNG]

NAME	START	P-VALUE	SITES
Chloroflexus	109	2.10e-31	`TRQPAPACFT`	`YTNFLLATAACRSYEEGMAALLPCFWIYREVGHDIYR`	`RAAPNNPYQK`
Sulfolobus	109	1.32e-30	`DYEMSPTNLA`	`YTSYLLAVAYSRPFHEVISAVLPCYWIYMEVGKELLK`	`KGSRDKYYQK`
Comamonas	111	1.95e-30	`ATPLTPACHH`	`YTSYLLATAWSTTYPVAVAALLPCFWIYAEVGRDIHA`	`RSAKDNPYQA`
Achromobacter	111	8.87e-30	`ATPLTPACHH`	`YTSYLIATAWSAPYPVALAALLPCFWIYAEIGREIHA`	`RAAQPNPYAA`
Acinetobacter	111	3.71e-29	`NTPLTLAGHH`	`YTSFLTATAWSESYPVVLAALLPCFWIYAEVGKDIVS`	`KSVPNNPYQA`
Stanieria	111	3.71e-29	`ATEPSPTCLN`	`YTNFLIATAYRDSFAVTVAAVLPCFWIYSEVGKHIYQ`	`TAKIDNNPYQ`
Acanthamoeba	492	2.36e-28	`AHPASPTNML`	`YTSYMLRIGFDRPYYEGVAAVLPCAWIYNEVGHYLKS`	`KGSPHPQYAR`
Mycobacterium	118	3.26e-28	`STPVAPTTRA`	`YTSYLLAAAYSGTFAEGLAAVLPCYWIYARVGTALLE`	`RGSPDRRYQR`
Isosphaera	124	3.25e-27	`ASELAPTCRA`	`YVSHLKSTVLGRPFHEGLAAVLPCYWIYQDVGQRLAA`	`QGSPDPLYAR`
Bacillus_subtilis_strain_168	111	5.79e-27	`AFKPSPTAYS`	`YTSHMYRSVLSGNFAEILAALLPCYWLYYEVGEKLLH`	`CDPGHPIYQK`
Verticillium	236	1.77e-26	`ATEEHQACTA`	`YTRYVLDVGMSGDWIGLQVALAPCLLGYGAVAKQLHG`	`DAASVRGEGN`
Magnaporthe	402	2.66e-26	`QSEEHVACTA`	`YSRYVLDVGQSSDYVALQMALAPCLLGYGAIAKHLHA`	`DPKSKRDGSP`
Phaeodactylum	108	6.70e-26	`GVEQMPHSLL`	`YTSYMKQVVATRPHAEGLAVLLPCFWVYMHVGQCMLK`	`LRQELGDSVQ`
Gibberella	401	7.62e-26	`ATPEKQACTA`	`YTRYVLDVGQNGDWLGLQMALAPCLLGYGAAAKMLHD`	`HENTVREDNT`
Octadecabacter	117	8.66e-26	`AASEEFENLA`	`YTRYVMDAGLQGDFLDLMAALAPCCFGYGEIGARLAQ`	`SAAPDTPYRD`
Microlunatus	126	1.84e-25	`AAGSGPTTTA`	`YASYLLGVCATGTYAEALAAVLPCYWVYRDVGRELLT`	`RSSPDPLYAQ`
Tuber	404	2.35e-25	`ETEESLACTA`	`YTRFVLDIGMSQDLFALQVAMAPCLLGYQEVALRLAK`	`EGKKDDNIYW`
Penicillium	400	3.38e-25	`GTPETIACTA`	`YSRYILDIGQSEDWLALQMALAPCLIGYGAIAKRLYT`	`GETTLRDGNT`
Bifidobacterium	146	2.13e-24	`HGQQSAFARA`	`YTTNMLTTAYSGDLLDILVSVLPCAWVYADYGERLAR`	`DFAQGLEHNP`
Candidatus	111	2.37e-24	`NEPMAPTNYA`	`YTRHMLAVAQGGTATEIAVVALPCAWIYCVVGKHFLK`	`NGPPPAKHPY`
Ustilago	525	5.02e-24	`ETMESAATLA`	`YTRFVLDTARSSDSLDLLVAVSPCMVGYAQIGLWLAQ`	`NRKEGIDKDY`
Zunongwangia	107	6.19e-24	`KAIIQPACHH`	`YIHFLKSMAALENLEVAMAAVLPCFWIYKEVGDYIYR`	`NQKSGNNPYQ`
Natronobacterium	111	9.33e-24	`TVQKSPTCVA`	`YTNFLVRTAHEGSLAEIAAAIYPCGQGYLDIAEHIAD`	`NLEEEDNPYL`
Pyrobaculum	110	1.03e-23	`GAEPNFVNVS`	`YMSYLLSVCSLEGFYQCMAALLPCFWSYAEIAQRHRD`	`ALERNPVPLY`
Natrialba	111	1.39e-23	`SVEKAPTCVA`	`YTNFLVRTAYEGSLAEIAGALYPCMQGYLDVAEHMAE`	`LADEEHRYTP`
Naegleria	130	2.76e-23	`DAKASKSCLL`	`YTSFLLASANSDTFLQTLCSCLPCYYIYYRVAQYITN`	`KMMNESGEVL`
Puccinia	400	4.02e-23	`QTPESNQTIA`	`YTRYVMDIGHSKSLLHLRVATAPCLLGYGEMAMKLKS`	`DPLTVKDESN`
Bacteroides	105	1.30e-22	`DIEPSPACEL`	`YTSYLLKHLAVSPVEVALAAVMPCFVIYQKVGSYIYK`	`HQVKGTNPYQ`
Leptospira	118	2.17e-22	`ETEPTPNNLA`	`YQNHMIASVVKGSVVEGFAAIAPCPWLYIDLGQHLLK`	`VKGSIPEDHP`
Coprinopsis	440	5.76e-22	`NTPESTATTA`	`YGAFIINSGVQGDACALMMAVLSCLLGYGEVGLWLKK`	`ESELPNSWVV`
Helicobacter	111	1.44e-21	`NACPTLANKS`	`YTSYMLAEGIKGSIKEVTVAVLACGWSYLVIAQNLSQ`	`IPNALEHAFY`

Motif 2 block diagrams

Name

Lowest
p-value

Motifs

Chloroflexus

2.10e-31

Sulfolobus

1.32e-30

Comamonas

1.95e-30

Achromobacter

8.87e-30

Acinetobacter

3.71e-29

Stanieria

3.71e-29

Acanthamoeba

2.36e-28

Mycobacterium

3.26e-28

Isosphaera

3.25e-27

Bacillus_subtilis_strain_168

5.79e-27

Verticillium

1.77e-26

Magnaporthe

2.66e-26

Phaeodactylum

6.70e-26

Gibberella

7.62e-26

Octadecabacter

8.66e-26

Microlunatus

1.84e-25

Tuber

2.35e-25

Penicillium

3.38e-25

Bifidobacterium

2.13e-24

Candidatus

2.37e-24

Ustilago

5.02e-24

Zunongwangia

6.19e-24

Natronobacterium

9.33e-24

Pyrobaculum

1.03e-23

Natrialba

1.39e-23

Naegleria

2.76e-23

Puccinia

4.02e-23

Bacteroides

1.30e-22

Leptospira

2.17e-22

Coprinopsis

5.76e-22

Helicobacter

1.44e-21

SCALE

\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|
1	25	50	75	100	125	150	175	200	225	250	275	300	325	350	375	400	425	450	475	500	525	550	575	600	625

Motif 2 in BLOCKS format

to BLOCKS multiple alignment processor.

Motif 2 position-specific scoring matrix

Scan sequence databases for the best match in each sequence using MAST.

Motif 2 position-specific probability matrix

Scan sequence databases for all matches with this motif using FIMO.

Motif 2 regular-expression

YT[SR][YF][LMV]L[AD][TV][AG]xSG[DS][FY][AL]ELL[AV]A[LV][LA]PC[FL]W[IG]Y[AG]E[VI][GA]Kx[LI][AL]K

Time 32.69 secs.

MOTIF 3 width = 41 sites = 40 llr = 1845 E-value = 2.8e-347

SEQUENCE LOGO

Information Content
71.0 (bits)

Relative Entropy
66.5 (bits)

Download LOGO
Without SSC:[EPS][PNG]
With SSC:[EPS][PNG]

NAME	START	P-VALUE	SITES
Acinetobacter	151	3.33e-31	`VGKDIVSKSV`	`PNNPYQAWIDTYAGEEFHTAVRNVIATVDKVAARCDADTLE`	`KMHAAYTMGA`
Chloroflexus	149	4.67e-28	`VGHDIYRRAA`	`PNNPYQKWIDTYAGQEFDEVVTQAIALTDTIAEQASATQQQ`	`RMRDAFIYSA`
Verticillium	280	2.13e-27	`LHGDAASVRG`	`EGNVYWKWIENYIAEDYQQALRTGSELLERHAARQSPEGIE`	`RLAKIFIHGT`
Comamonas	151	1.68e-25	`VGRDIHARSA`	`KDNPYQAWVDTYASEEFHAAVRGVCATADRLAEEATEATRA`	`AMHAAYKDAA`
Isosphaera	164	1.68e-25	`VGQRLAAQGS`	`PDPLYARWIATYGGAEFEAVTRQVLDLTDRLAPTLTPDQRA`	`AMTQQFVLSS`
Zunongwangia	148	1.98e-25	`GDYIYRNQKS`	`GNNPYQAWIDTYAGEEFGEAVQKAISICDDHAINTTPQIRT`	`KMTEAFVSAS`
Stanieria	152	2.33e-25	`GKHIYQTAKI`	`DNNPYQKWIETYADPDFEASVNYIIQVVDQQAKIASVQELK`	`LMEQAFYRAS`
Achromobacter	151	3.23e-25	`IGREIHARAA`	`QPNPYAAWIEAYACEDFHALVRRVIATVDRVAEQASAETLA`	`QMHAAYTRAA`
Gibberella	444	1.34e-24	`MLHDHENTVR`	`EDNTYWAWIKNYNEEDYTDAVKLGSALLEKHIQLQSPSRIE`	`ELVQIFIHAL`
Oryza	223	2.12e-24	`KELAVFLKQD`	`ENHPYKKWIETYASSDFENNALQIEELLDKLSVSLTGEELE`	`IIGKLYQQAM`
Magnaporthe	445	1.67e-23	`HLHADPKSKR`	`DGSPYWTWIQNYVADDYVQAVKTGSELLERHAMLQSPSRIN`	`ELVEVFKHAT`
Penicillium	443	3.38e-23	`RLYTGETTLR`	`DGNTYWKWIENYVADDYTKAVQIGSALLEDHMQHVSPSRME`	`ELIKIFVRAT`
Natronobacterium	152	6.75e-23	`AEHIADNLEE`	`EDNPYLPFVEKYTSDEFHEAVDWMRSFVDRCGERYPGEHEA`	`MREAFLTSAR`
Naegleria	179	8.85e-23	`MNESGEVLLT`	`DSHPYADWIKSYSSDEFKKGVEQFSAFINTYAPTASESEKE`	`KCFENVKMAS`
Physcomitrella	166	1.73e-22	`QELKKHVDHV`	`ADHPYQEWIDTYSAAEFEAAASKIEQLLDKLTATLTGKHEI`	`AFLESLYLQA`
Solanum	221	1.73e-22	`KELQVCLEGE`	`KNHPYEKWIDSYAAEGFQASALQTEDLLDKLSVTLTGEELD`	`IIEKLYHQAM`
Vanderwaltozyma	483	1.73e-22	`INHKDYVVKN`	`EGNVYYEWCQVYSSTWYHEGMEKGKQLLNQIASTYPKENIE`	`VLVNIYAQVC`
Schizosaccharomyces	474	2.25e-22	`ARLMKEPFRN`	`PQGPYQKWVDNYFCEDYLSAVRRGCRQIEEIVLKLSPERIQ`	`ELIEIFIRAT`
Mycobacterium	158	2.56e-22	`VGTALLERGS`	`PDRRYQRWIDSYGGDEFAATVAEVLGVTDRVGPTLGEEDDA`	`AARAHFEVTA`
Bacillus_subtilis_strain_168	151	2.56e-22	`VGEKLLHCDP`	`GHPIYQKWIGTYGGDWFRQQVEEQINRFDELAENSTEEVRA`	`KMKENFVISS`
Natrialba	151	2.91e-22	`VAEHMAELAD`	`EEHRYTPFIEMYTGEEFREATAWCREFVDRCGERYPGEREA`	`MCEAFLTSAK`
Microlunatus	166	4.87e-22	`VGRELLTRSS`	`PDPLYAQWIASYGSPEFDAVVESVLAVTDGLGSDIGPSERV`	`RCHEHFAVTT`
Puccinia	444	6.28e-22	`LKSDPLTVKD`	`ESNPYWEWILHYADPGFQSAVALGREVLEQMVLRDPPSVTS`	`FQTLKQIFKK`
Coprinopsis	488	7.13e-22	`SELPNSWVVM`	`EGNPYKKWIEDYSGEGYQSAVRVGLATIEKMLVADHPSPAR`	`LKEWLGIFQR`
Acanthamoeba	532	1.17e-21	`VGHYLKSKGS`	`PHPQYARWIETYGSEEFSATTKSLIEITTGVAAGLGDEQRR`	`RMRECFVQTS`
Bacteroides	146	1.33e-21	`GSYIYKHQVK`	`GTNPYQAWIDTYAGSEHVESVITATRICDEIAKSCTIEQQE`	`AMLNAYKTSS`
Bifidobacterium	189	1.69e-21	`RLARDFAQGL`	`EHNPCRQWIDMYKTDEFWQSSVWLLAQIERRVAGVGEERRR`	`GLVDSFVRGV`
Candidatus	153	3.09e-21	`KHFLKNGPPP`	`AKHPYRDWLMLYASPEFEAVQEWMRARVDTWARTAGHEEKK`	`RMEQAFLISS`
Ustilago	567	1.39e-20	`LWLAQNRKEG`	`IDKDYAAWIEAYSGQDFQAVVRRAMRLVEAKAAKDPPSPQR`	`LRKLQTTWNA`
Tuber	445	3.35e-20	`ALRLAKEGKK`	`DDNIYWEWVENYSGKGYGEAVIKSRELLEASAQDIGVKKIG`	`ELVAIFAKAT`
Torulaspora	485	7.06e-20	`INQKKFITVE`	`EDDVYYEWCQTYLSVWYQDAMARGKELLNQIALTGSKAELD`	`TLVTIFADVC`
Sulfolobus	149	7.83e-20	`VGKELLKKGS`	`RDKYYQKWIETYGGEDYERGVRAVLGIVNGLKVSEEEFNKM`	`KIHFRTASIY`
Helicobacter	153	3.54e-19	`QNLSQIPNAL`	`EHAFYGHWIKGYSSKEFQACVSWNINLLDSLTHASSKQEIE`	`KLKDIFITTS`
Leptospira	161	5.21e-19	`HLLKVKGSIP`	`EDHPYASWLLMYSDPGFNEYMTNLLSRLQKYADLADENSKK`	`RAVAAFQQSC`
Phaeodactylum	155	6.30e-19	`LRQELGDSVQ`	`RPPQFDAWIDMYGGEDFEHEVKDYIAMNNKACQSADAETLK`	`NMEKHFLTCC`
Octadecabacter	157	1.01e-18	`IGARLAQSAA`	`PDTPYRDWIATYADDDYQNVMVSVGHMIDRAVTRRLGGDPA`	`GSARWLHLQA`
Scardovia	153	3.22e-18	`IADYWKDQTS`	`PVPVYQAFFDSYQTVGFDTSTQQLIGLVNELAANSNETVKQ`	`QMMTAFLRSS`
Thalassiosira	154	1.33e-17	`GEKVLPFTNR`	`DDFTTYEWVDLHSGTYFSEVVTYLRGLLDKEGTLIDEEGKK`	`KCKDAFLRAI`
Perkinsus	158	1.18e-16	`RELTADREVP`	`ETHPFRRWLLSHADEPINTSTRVLESLLDKHIRPGEFKEVA`	`QAYRRAMELE`
Pyrobaculum	152	2.35e-16	`QRHRDALERN`	`PVPLYRRWASVYLTPEYNALVEKLRALLESSGAPPEALWPY`	`FRQASIYELK`

Motif 3 block diagrams

Name

Lowest
p-value

Motifs

Acinetobacter

3.33e-31

Chloroflexus

4.67e-28

Verticillium

2.13e-27

Comamonas

1.68e-25

Isosphaera

1.68e-25

Zunongwangia

1.98e-25

Stanieria

2.33e-25

Achromobacter

3.23e-25

Gibberella

1.34e-24

Oryza

2.12e-24

Magnaporthe

1.67e-23

Penicillium

3.38e-23

Natronobacterium

6.75e-23

Naegleria

8.85e-23

Physcomitrella

1.73e-22

Solanum

1.73e-22

Vanderwaltozyma

1.73e-22

Schizosaccharomyces

2.25e-22

Mycobacterium

2.56e-22

Bacillus_subtilis_strain_168

2.56e-22

Natrialba

2.91e-22

Microlunatus

4.87e-22

Puccinia

6.28e-22

Coprinopsis

7.13e-22

Acanthamoeba

1.17e-21

Bacteroides

1.33e-21

Bifidobacterium

1.69e-21

Candidatus

3.09e-21

Ustilago

1.39e-20

Tuber

3.35e-20

Torulaspora

7.06e-20

Sulfolobus

7.83e-20

Helicobacter

3.54e-19

Leptospira

5.21e-19

Phaeodactylum

6.30e-19

Octadecabacter

1.01e-18

Scardovia

3.22e-18

Thalassiosira

1.33e-17

Perkinsus

1.18e-16

Pyrobaculum

2.35e-16

SCALE

\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|
1	25	50	75	100	125	150	175	200	225	250	275	300	325	350	375	400	425	450	475	500	525	550	575	600	625

Motif 3 in BLOCKS format

to BLOCKS multiple alignment processor.

Motif 3 position-specific scoring matrix

Scan sequence databases for the best match in each sequence using MAST.

Motif 3 position-specific probability matrix

Scan sequence databases for all matches with this motif using FIMO.

Motif 3 regular-expression

[EP]D[NH]PYQ[KA]WI[DE]TY[AS][GS][ED][ED][FY]Q[AE]AVRx[GV]I[AE]L[LV][DE][KR]LAATA[ST][PE][ES]Ex[AE]

Time 45.53 secs.

SUMMARY OF MOTIFS

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

Name

Combined
p-value

Motifs

Bacillus_subtilis_strain_168

3.90e-68

Natronobacterium

3.29e-66

Natrialba

1.77e-67

Sulfolobus

1.98e-72

Pyrobaculum

8.72e-64

Microlunatus

3.57e-70

Mycobacterium

1.20e-71

Bifidobacterium

1.25e-63

Scardovia

3.62e-39

Zunongwangia

8.32e-69

Bacteroides

7.15e-62

Stanieria

3.19e-75

Chloroflexus

7.22e-82

Leptospira

4.03e-58

Isosphaera

1.74e-70

Candidatus

4.34e-63

Helicobacter

4.48e-64

Octadecabacter

3.95e-68

Acinetobacter

6.15e-81

Comamonas

4.69e-79

Achromobacter

6.96e-78

Acanthamoeba

2.71e-71

Magnaporthe

1.10e-73

Torulaspora

1.46e-50

Vanderwaltozyma

3.27e-49

Gibberella

1.18e-75

Penicillium

7.16e-69

Tuber

8.91e-68

Ustilago

4.64e-61

Verticillium

1.26e-80

Schizosaccharomyces

1.12e-52

Coprinopsis

1.31e-59

Puccinia

3.29e-63

Oryza

7.03e-42

Solanum

1.88e-40

Physcomitrella

1.49e-43

Thalassiosira

2.45e-34

Phaeodactylum

2.36e-60

Perkinsus

1.81e-38

Naegleria

7.65e-61

SCALE

\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|	\|
1	25	50	75	100	125	150	175	200	225	250	275	300	325	350	375	400	425	450	475	500	525	550	575	600	625

Motif summary in machine readable format.

Stopped because Stopped because nmotifs = 3 reached..

CPU: kodomo.fbb.msu.ru

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 sequence LOGO.
3. The information content of the motif.
4. The relative entropy of the motif.
5. Downloadable LOGO files suitable for publication.
6. The occurrences of the motif sorted by p-value and aligned with each other.
7. Block diagrams of the occurrences of the motif within each sequence in the training set.
8. The motif in BLOCKS or FASTA format.
9. A position-specific scoring matrix (PSSM) for use by the MAST database search program.
10. The position specific probability matrix (PSPM) describing the motif.
11. A regular expression 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).

Sequence LOGO

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. These are displayed as "sequence LOGOS", containing stacks of letters at each position in the motif. The total height of the stack is the "information content" of that position in the motif in bits. The height of the individual letters in a stack is the probability of the letter at that position multiplied by the total information content of the stack.

Note: The MEME LOGO differs from those produced by the Weblogo program because a small-sample correction is NOT applied. However, MEME LOGOs in PNG and encapsulated postscript (EPS) formats with small-sample correction (SSC) are available by clicking on one of the links named "With SSC" (EPS or PNG) under Download LOGO. The MEME LOGOs without small sample correction are similarly available. Error bars are included in the LOGOs with small-sample correction.

The information content of each motif position is computed as described in the paper by Schneider and Stephens, "Sequence Logos: A New Way to Display Consensus Sequences" but the small-sample correction, e(n), is set to zero for the LOGO displayed in the MEME output. The corrected information content of position i is given by

           R(i) for amino acids   = log2(20) - (H(i) + e(n))   (1a) 
           R(i) for nucleic acids =    2    - (H(i) + e(n))    (1b)

where H(i) is the entropy of position i,

           H(l) = - (Sum f(a,i) * log2[ f(a,i) ]).             (2)

Here, f(a,i) is the frequency of base or amino acid a at position i, and e(n) is the small-sample correction for an alignment of n letters. The height of letter a in column i is given by

           height = f(a,i) * R(i)                              (3)

The approximation for the small-sample correction, e(n), is given by:

           e(n) = (s-1) / (2 * ln(2) * n),                     (4)

where s is 4 for nucleotides, 20 for amino acids, and n is the number of sequences in the alignment.

The letters in the logos are colored as follows. 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	PROPERTIES
A, C, F, I, L, V, W and M	BLUE	Most hydrophobic[Kyte and Doolittle, 1982]
NQST	GREEN	Polar, non-charged, non-aliphatic residues
DE	MAGENTA	Acidic
KR	RED	Positively charged
H	PINK
G	ORANGE
P	YELLOW
Y	TURQUOISE

J. Kyte and R. Doolittle, 1982. "A Simple Method for Displaying the Hydropathic Character of a Protein", J. Mol Biol. 157, 105-132.

Note: the "text" output format of MEME preserves the historical MEME format where LOGOS are replaced by a simplified probability matrix, a relative entropy plot, and a multi-level consensus sequence.

Information Content
This is the information content of the motif in bits. It is equal to the sum of the uncorrected information content, R(), in the columns of the LOGO. This is equal relative entropy of the motif relative to a uniform background frequency model.
Relative Entropy
This is the relative entropy of the motif, computed in bits and relative to the background letter frequencies given in the command line summary. It is equal to the log-likelihood ratio (llr) divided by the number of occurrences (sites) of the motif times 1/ln(2),
```
               re = llr / (sites * ln(2)).
	    
```
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 100 times 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, the scoring threshold (obsolete) and the motif E-value.
Note: The probability p used to compute the PSSM is not exactly the same as the corresponding value in the Position Specific Probability Matrix (PSPM). The values of p used to compute the PSSM take into account the motif prior, whereas the values in the PSPM are just the observed frequencies of letters in the motif sites.
Position-Specific Probability Matrix
The motif itself is a position-specific probability matrix giving, for each position in the pattern, the observed frequency ("probability") of each possible letter. 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, number of occurrences of the motif, and the E-value of the motif.
Note: Earlier versions of MEME gave the posterior probabilities--the probability after applying a prior on letter frequencies--rather than the observed frequencies. These versions of MEME also gave the number of possible positions for the motif rather than the actual number of occurrences. The output from these earlier versions of MEME can be distinguished by "n=" rather than "nsites=" in the line preceding the matrix.
Regular Expression
This is a regular expression (RE) describing the motif. In each column, all letters with observed frequencies greater than 0.2 are shown; less-frequent letters are not included in the RE. MEME regular expressions are interpreted as follows: single letters match that letter; groups of letters in square brackets match any of the letters in the group. Regular expressions can be used for searching for the motif in sequences (using, for example, PatMatch ) but the search accuracy will usually be better with the PSSM (using, for example MAST.)
Motif Summary Tiling
The motif summary tiling is done using the same algorithm as used by MAST. The motif occurrences shown in the motif summary may not be exactly the same as those reported in each motif section because only motifs with a position p-value of 0.0001 that don't overlap other, more significant motif occurrences are shown. The format of the machine readable motif-summary is:
```
        [sequence_name combined_p-value number_of_motif_occurrences [motif_number start_of_motif position_p-value]+]+
        
```
See the documentation for MAST output for the definition of position and combined p-values.