This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Fimo output for classes of TFs from specific database

Currently if getting a fimo database for specific TFs can do the following:

fimo -max-stored-scores 10000000 --thresh 1e-6 -oc /path/to/output/directory -motif /path/to/motif/file \
  /path/to/HOCOMOCOv11_HUMAN_mono_meme_format.meme /path/to/whole_genome.fa

Is there a good way to do this for classes of TFs all at once, for instance if I wanted to look at all Jaspar non redundant vertebrates Forkhead motifs and where they are in mm10 and get an output for the whole subclass instead of one at a time as above?

Thanks!

chip-seq fimo meme
• 1,877 views • link
updated by Alex Reynolds Alex Reynolds • written by rbronste rbronste

I may point your attention to the FAQs of FIMO. You might read the section on why it is not preferrable to run it on a whole genome fasta, due to the multiple testing problem. Has nothng to do with your question directly, but should be read IMHO.

• 1 views • link
written by ATpoint ATpoint

1 answer

Perhaps take out the --motif argument and pass in a concatenation of the MEME-formatted models. It should treat each model independently. Via: http://meme-suite.org/doc/fimo.html

• 1 views • link
written by Alex Reynolds Alex Reynolds

Log in to answer this question.

More posts like this