MEME output: HTML plain text XML Motifs discovered in 600 (randomly chosen) trimmed (central 100bp) input sequences. TOMTOM output: HTML plain text XML Motifs from JASPAR CORE 2009 that match motifs MEME discovers. MAST output: HTML plain text XML Predicted locations of all MEME motifs (p < 0.0001) in the input sequences. AMA output: plain text XML Estimated binding affinity of each MEME motif to each input sequence. DREME output: plain text Motifs discovered in the trimmed (central 100bp) input sequences. TOMTOM output: HTML plain text XML Motifs from JASPAR CORE 2009 that match motifs DREME discovers. MAST output: HTML plain text XML Predicted locations of all DREME motifs (p < 0.0005) in the input sequences. AME output: HTML plain text JASPAR CORE 2009 motifs enriched in the trimmed (central 100bp) input sequences.
input: sequences your original untrimmed sequences fasta-center output: seqs-centered your sequences centered and trimed to width 100 fasta-dinucleotide-shuffle output: seqs-shuffled the centered sequences randomly shuffled maintaining their dinucleotide frequency, used as a background for AME AME input: seqs-centered_w_bg the centered sequences followed by the same sequences after dinucleotide shuffling fasta-subsample (used) output: seqs-sampled a random sample of 600 of the centered sequences, used by MEME fasta-subsample (discarded) output: seqs-discarded the centered sequences omitted from the sample used by MEME
fasta-center -len 100 < sequences > seqs-centered fasta-dinucleotide-shuffle -f seqs-centered -t -dinuc > seqs-shuffled cat seqs-centered seqs-shuffled > seqs-centered_w_bg fasta-subsample seqs-centered 600 -rest seqs-discarded > seqs-sampled meme -oc meme_out -nostatus seqs-sampled -sf RM2_SCL-width500-nogene.fasta_1278655789.masked -dna -mod zoops -nmotifs 3 -minw 6 -maxw 30 -time 7200 -revcomp tomtom -verbosity 1 -oc meme_tomtom_out -min-overlap 5 -dist pearson -evalue -thresh 0.1 -no-ssc meme_out/meme.html databases/motif_databases/JASPAR_CORE_2009.meme mast -nostatus -oc meme_mast_out meme_out/meme.html sequences -ev 2400 ama --oc ama_out --verbosity 1 --sdbg 0 meme_out/meme.html sequences dreme -v 1 -p seqs-centered > dreme.txt tomtom -verbosity 1 -oc dreme_tomtom_out -min-overlap 5 -dist pearson -evalue -thresh 0.1 -no-ssc ./dreme.txt databases/motif_databases/JASPAR_CORE_2009.meme mast -nostatus -oc dreme_mast_out ./dreme.txt sequences -ev 2400 -mt 0.0005 ame --oc ame_out --verbose 1 --fix-partition 2386 --bgformat 0 seqs-centered_w_bg databases/motif_databases/JASPAR_CORE_2009.meme