MEME version 3.5.0 (Release date: 2005/07/29 20:02:31)
For further information on how to interpret these results or to get a copy of the MEME software please access http://meme.nbcr.net.
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.nbcr.net.
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.
DATAFILE= C:\Documents_and_Settings\juhur\Desktop\ProtSeq.fas ALPHABET= ACDEFGHIKLMNPQRSTVWY Sequence name Weight Length Sequence name Weight Length ------------- ------ ------ ------------- ------ ------ At1g01140.1_4-2-4_SnRK3. 1.0000 447 At1g01140.2_SnRK3.12 1.0000 449 At1g01140.3_SnRK3.12 1.0000 451 At1g01450.1_2-1-1 1.0000 470 At1g01540.1_1-6-3 1.0000 386 At1g01540.2_Putative 1.0000 472 At1g01560.1_4-5-1_MPK11 1.0000 275 At1g01740.1_1-16-1 1.0000 483 At1g02970.1_4-3-1 1.0000 500 At1g03740.1_4-5-2 1.0000 542 At1g03920.1_4-2-6 1.0000 569 At1g03930.1_3-1-1-1_ADK1 1.0000 471 At1g04210.1_protein 1.0000 1112 At1g04440.1_3-1-1-1 1.0000 468 At1g04700.1_2-1-4-1_Raf1 1.0000 1042 At1g05100.1_4-4-1_MAPKKK 1.0000 339 At1g05700.1_1-8-1 1.0000 843 At1g06390.1_4-5-4_ASK-io 1.0000 407 At1g06390.2_ASK-iota 1.0000 407
This information can also be useful in the event you wish to report a problem with the MEME software. command: meme /home/meme/meme-rocks/LOGS/meme.5659.data -protein -mod zoops -nmotifs 3 -minw 6 -maxw 20 -evt 1e100 -time 7200 -maxsize 60000 -nostatus -maxiter 20 model: mod= zoops nmotifs= 3 evt= 1e+100 object function= E-value of product of p-values width: minw= 6 maxw= 20 minic= 0.00 width: wg= 11 ws= 1 endgaps= yes nsites: minsites= 2 maxsites= 19 wnsites= 0.8 theta: prob= 1 spmap= pam spfuzz= 120 em: prior= megap b= 50665 maxiter= 20 distance= 1e-05 data: n= 10133 N= 19 sample: seed= 0 seqfrac= 1 Dirichlet mixture priors file: prior30.plib Letter frequencies in dataset: A 0.053 C 0.017 D 0.055 E 0.071 F 0.041 G 0.066 H 0.024 I 0.051 K 0.065 L 0.099 M 0.023 N 0.048 P 0.046 Q 0.034 R 0.061 S 0.091 T 0.050 V 0.062 W 0.011 Y 0.032 Background letter frequencies (from dataset with add-one prior applied): A 0.053 C 0.017 D 0.055 E 0.071 F 0.041 G 0.065 H 0.024 I 0.051 K 0.065 L 0.099 M 0.023 N 0.048 P 0.046 Q 0.034 R 0.061 S 0.091 T 0.050 V 0.062 W 0.011 Y 0.032
Time 5.34 secs.
Time 10.58 secs.
Time 15.71 secs.
CPU: compute-0-1.local
MOTIFS
For each motif that it discovers in the training set, MEME prints the following information:
J. Kyte and R. Doolittle, 1982. "A Simple Method for Displaying the Hydropathic Character of a Protein", J. Mol Biol. 157, 105-132.
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 TTATGTGAACGACGTCACACT consensus AA T A G A GA AA sequence T C TT T
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.
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.