******************************************************************************** MEME - Motif discovery tool ******************************************************************************** MEME version 4.10.0 (Release date: Wed May 21 10:35:36 2014 +1000) 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. ******************************************************************************** ******************************************************************************** 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= MEME14.txt ALPHABET= ACGT Sequence name Weight Length Sequence name Weight Length ------------- ------ ------ ------------- ------ ------ codB 1.0000 200 purE 1.0000 200 pyrC 1.0000 200 purR 1.0000 200 cvpA 1.0000 200 purM 1.0000 200 guaB 1.0000 200 glnB 1.0000 200 purL 1.0000 200 purA 1.0000 200 folD 1.0000 200 rpiA 1.0000 200 carA 1.0000 200 pdhR 1.0000 200 ******************************************************************************** ******************************************************************************** COMMAND LINE SUMMARY ******************************************************************************** This information can also be useful in the event you wish to report a problem with the MEME software. command: meme MEME14.txt -dna -oc . -nostatus -time 18000 -maxsize 60000 -mod oops -nmotifs 1 -minw 16 -maxw 16 -revcomp -psp priors.psp model: mod= oops nmotifs= 1 evt= inf object function= E-value of product of p-values width: minw= 16 maxw= 16 minic= 0.00 width: wg= 11 ws= 1 endgaps= yes nsites: minsites= 14 maxsites= 14 wnsites= 0.8 theta: prob= 1 spmap= uni spfuzz= 0.5 global: substring= yes branching= no wbranch= no em: prior= dirichlet b= 0.01 maxiter= 50 distance= 1e-05 data: n= 2800 N= 14 strands: + - sample: seed= 0 seqfrac= 1 Letter frequencies in dataset: A 0.273 C 0.227 G 0.227 T 0.273 Background letter frequencies (from dataset with add-one prior applied): A 0.273 C 0.227 G 0.227 T 0.273 ******************************************************************************** ******************************************************************************** MOTIF 1 MEME width = 16 sites = 14 llr = 175 E-value = 7.4e-013 ******************************************************************************** -------------------------------------------------------------------------------- Motif 1 Description -------------------------------------------------------------------------------- Simplified A ::4a951131:::2:: pos.-specific C 7:4::29::11:9412 probability G 2a1:11192:1811:2 matrix T 1:1::1::5992:496 bits 2.1 * 1.9 * * 1.7 * * * 1.5 * ** * * * Relative 1.3 * ** ** ** * Entropy 1.1 ** ** ** **** * (18.1 bits) 0.9 ** ** ** **** * 0.6 ** ** ** **** ** 0.4 ** ** ******* ** 0.2 **************** 0.0 ---------------- Multilevel CGAAAACGTTTGCCTT consensus G C C A T T C sequence G A G -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 sites sorted by position p-value -------------------------------------------------------------------------------- Sequence name Strand Start P-value Site ------------- ------ ----- --------- ---------------- purM + 122 1.06e-09 AAAGCAGTCT CGCAAACGTTTGCTTT CCCTGTTAGA codB + 119 4.28e-09 TATTTCCCCA CGAAAACGATTGCTTT TTATCTTCAG cvpA + 130 2.12e-08 GAAATCCCTA CGCAAACGTTTTCTTT TTCTGTTAGA purE + 115 4.54e-08 TTCACAGCCA CGCAACCGTTTTCCTT GCTCTCTTTC purR - 139 5.16e-08 AAAATCGCAA GGTAAACGTTTGCCTT TACACACCTT purL - 108 1.44e-07 GATGCGCTGA CGAAACCGTTTGCGTG GAAATAAAAT pyrC - 133 1.44e-07 AAAGGATAAG CGGAAACGTTTTCCTT TGCACGAAAA guaB - 131 3.70e-07 TTATACAGAG CGTAACCGATTGCATC TACCCCTTTT purA + 80 6.45e-07 TACATGTTGA GGAAAACGATTGGCTG AACAAAAAAC glnB - 117 1.37e-06 ATTCATTCCT TGAAATCGTTTGCATC CAGCTCGTGT carA - 91 6.87e-06 ACGAATTTCT GGCAAACGGCGGCATT CTGGAGATAT rpiA + 74 8.21e-06 ATTTGCGGGG CGAAAGGGGATGCCTG CCATTGCGCG folD + 87 9.21e-06 GCCTCACCTT CGCAAGAGGTCGCTTC ACGCGATAAA pdhR - 160 6.08e-05 AGCCACTTGC CGAAGTCAATTGGTCT TACCAATTTC -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 block diagrams -------------------------------------------------------------------------------- SEQUENCE NAME POSITION P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- purM 1.1e-09 121_[+1]_63 codB 4.3e-09 118_[+1]_66 cvpA 2.1e-08 129_[+1]_55 purE 4.5e-08 114_[+1]_70 purR 5.2e-08 138_[-1]_46 purL 1.4e-07 107_[-1]_77 pyrC 1.4e-07 132_[-1]_52 guaB 3.7e-07 130_[-1]_54 purA 6.4e-07 79_[+1]_105 glnB 1.4e-06 116_[-1]_68 carA 6.9e-06 90_[-1]_94 rpiA 8.2e-06 73_[+1]_111 folD 9.2e-06 86_[+1]_98 pdhR 6.1e-05 159_[-1]_25 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 in BLOCKS format -------------------------------------------------------------------------------- BL MOTIF 1 width=16 seqs=14 purM ( 122) CGCAAACGTTTGCTTT 1 codB ( 119) CGAAAACGATTGCTTT 1 cvpA ( 130) CGCAAACGTTTTCTTT 1 purE ( 115) CGCAACCGTTTTCCTT 1 purR ( 139) GGTAAACGTTTGCCTT 1 purL ( 108) CGAAACCGTTTGCGTG 1 pyrC ( 133) CGGAAACGTTTTCCTT 1 guaB ( 131) CGTAACCGATTGCATC 1 purA ( 80) GGAAAACGATTGGCTG 1 glnB ( 117) TGAAATCGTTTGCATC 1 carA ( 91) GGCAAACGGCGGCATT 1 rpiA ( 74) CGAAAGGGGATGCCTG 1 folD ( 87) CGCAAGAGGTCGCTTC 1 pdhR ( 160) CGAAGTCAATTGGTCT 1 // -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific scoring matrix -------------------------------------------------------------------------------- log-odds matrix: alength= 4 w= 16 n= 2590 bayes= 7.52356 E= 7.4e-013 -1045 166 -8 -193 -1045 -1045 214 -1045 65 66 -166 -94 187 -1045 -1045 -1045 176 -1045 -166 -1045 87 -8 -67 -94 -193 192 -166 -1045 -193 -1045 203 -1045 6 -1045 -8 87 -193 -166 -1045 165 -1045 -166 -166 165 -1045 -1045 179 -35 -1045 192 -67 -1045 -35 66 -166 39 -1045 -166 -1045 176 -1045 -8 -8 106 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 position-specific probability matrix -------------------------------------------------------------------------------- letter-probability matrix: alength= 4 w= 16 nsites= 14 E= 7.4e-013 0.000000 0.714286 0.214286 0.071429 0.000000 0.000000 1.000000 0.000000 0.428571 0.357143 0.071429 0.142857 1.000000 0.000000 0.000000 0.000000 0.928571 0.000000 0.071429 0.000000 0.500000 0.214286 0.142857 0.142857 0.071429 0.857143 0.071429 0.000000 0.071429 0.000000 0.928571 0.000000 0.285714 0.000000 0.214286 0.500000 0.071429 0.071429 0.000000 0.857143 0.000000 0.071429 0.071429 0.857143 0.000000 0.000000 0.785714 0.214286 0.000000 0.857143 0.142857 0.000000 0.214286 0.357143 0.071429 0.357143 0.000000 0.071429 0.000000 0.928571 0.000000 0.214286 0.214286 0.571429 -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- Motif 1 regular expression -------------------------------------------------------------------------------- [CG]G[AC]AA[AC]CG[TAG]TT[GT]C[CTA]T[TCG] -------------------------------------------------------------------------------- Time 0.14 secs. ******************************************************************************** ******************************************************************************** SUMMARY OF MOTIFS ******************************************************************************** -------------------------------------------------------------------------------- Combined block diagrams: non-overlapping sites with p-value < 0.0001 -------------------------------------------------------------------------------- SEQUENCE NAME COMBINED P-VALUE MOTIF DIAGRAM ------------- ---------------- ------------- codB 1.58e-06 118_[+1(4.28e-09)]_66 purE 1.68e-05 114_[+1(4.54e-08)]_70 pyrC 5.35e-05 132_[-1(1.44e-07)]_52 purR 1.91e-05 138_[-1(5.16e-08)]_46 cvpA 7.83e-06 129_[+1(2.12e-08)]_55 purM 3.90e-07 121_[+1(1.06e-09)]_63 guaB 1.37e-04 130_[-1(3.70e-07)]_54 glnB 5.08e-04 116_[-1(1.37e-06)]_68 purL 5.35e-05 107_[-1(1.44e-07)]_21_\ [-1(7.80e-05)]_40 purA 2.39e-04 79_[+1(6.45e-07)]_105 folD 3.40e-03 86_[+1(9.21e-06)]_98 rpiA 3.03e-03 73_[+1(8.21e-06)]_111 carA 2.54e-03 90_[-1(6.87e-06)]_94 pdhR 2.22e-02 159_[-1(6.08e-05)]_25 -------------------------------------------------------------------------------- ******************************************************************************** ******************************************************************************** Stopped because nmotifs = 1 reached. ******************************************************************************** CPU: compute-0-2.local ********************************************************************************