homer findMotifsGenome.pl
0
0
Entering edit mode
11 months ago

Hi all,

I have an ATAC-seq data and from differential expressed peaks I would like to find any motif that can be useful for me. This is my bed file (first 5 row);

chr1 816940 818102 Interval_1 0 +
chr1 826720 828040 Interval_2 0 +
chr1 858479 858908 Interval_3 0 +
chr1 869659 870152 Interval_4 0 +
chr1 875146 875666 Interval_5 0 +

this is my command for running findMotifsGenome.pl;

findMotifsGenome.pl ./DE_bed.bed hg38 ./170124/ -size given -preparsedDir ./170124_preparsed

I downloaded my genome data using following command;

perl /path-to-homer/configureHomer.pl -install hg38

When I run the findMotifsGenome.pl, I get this error;

 !!! Something is wrong... are you sure you chose the right length for motif finding?
    !!! i.e. also check your sequence file!!!

    -blen automatically set to 2
    Scanning input files...

I couldn't find the solution of this. If you can help me I will be happy. Thanks

chip-seq atac-seq homer • 546 views
ADD COMMENT
0
Entering edit mode

Not a solution, but I don't see the flag for -preparsedDir in the manual. What happens if you remove that?

ADD REPLY

Login before adding your answer.

Traffic: 2206 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