How to mask motif in HOMER?
1
0
Entering edit mode
5.9 years ago
salamandra ▴ 550

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.2k views
ADD COMMENT
1
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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.

ADD REPLY
3
Entering edit mode
5.9 years ago
h.mon 35k

Maybe salamandra needs -maskMotif instead of -motif:

 -maskMotif <motif file1> [motif file 2]... (motifs to mask before motif finding)
ADD COMMENT
0
Entering edit mode

You are right. It's running! Thanks!

ADD REPLY

Login before adding your answer.

Traffic: 2803 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6