This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to mask motif in HOMER?

Why when I do this, to mask motif known.1:

findMotifsGenome.pl $PEAK $GENOME $OUT -p 4 -mask $MOTIF

it says: motif not recognized.

I checked the path is correct, so what could else be happening?

mask motif chip-seq homer
• 2,955 views • link
updated by Kevin Blighe Kevin Blighe • written by salamandra salamandra

Please, show the contents of your bash variables, or show the actual command as it was executed.

Please, show some lines from the motifs file.

• 1 views • link
written by h.mon h.mon

Hello again

-mask does not take any arguments. It just instructs HOMER to mask repeats and masked (lower case) sequence.

Just run:

findMotifsGenome.pl $PEAK $GENOME $OUT -p 4 -mask

In your output directory will then be a list of motifs present in your genome.

-------------------------------------------------

Edit: I'm aware that the manual states:

Mask motifs ("-mask <motif file="">")

Mask the motif(s) in the supplied motif file before starting motif finding. Multiple motifs can be in the motif file.

The documentation for Homer is not great and often confusing.

Kevin

• 1 views • link
written by Kevin Blighe Kevin Blighe

1 answer

Maybe salamandra needs -maskMotif instead of -motif:

 -maskMotif <motif file1> [motif file 2]... (motifs to mask before motif finding)
• 1 views • link
written by h.mon h.mon

You are right. It's running! Thanks!

• 1 views • link
written by salamandra salamandra

Log in to answer this question.

More posts like this