How to get Motif positions by peak (not by Chromosome) using Homer ??
Entering edit mode
19 months ago
carina2817 ▴ 20


I would like to get positions of motifs (start and end positions) for the motifs found in every peak using of Homer. I tried the following: ad_close.bed Aque1.31 -m knownResults/*.motif -gff3 Aqu2.1_Genes.gff3  -mbed motifs.bed > motifs.tsv

The output file motifs.bed has positions for the peaks present in each contig:

Contig13470 239433 239445 MYB3R4(MYB)/Homer 8.288071 +

Contig13464 37847 37859 MYB3R4(MYB)/Homer 7.195835 +

Contig13462 264324 264336 MYB3R4(MYB)/Homer 9.530264 -

Contig13036 3430 3442 MYB3R4(MYB)/Homer 7.475983 -

but I want the output to have positions for the motifs found in each peak (by peak id) , I tried the option -mpeak instead of -mbed and only changes the format of the output, I also tried with a peak file as input instead of bed file. Any ideas?


Homer annotatePeaks Motif positions peakid • 556 views

Login before adding your answer.

Traffic: 1854 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6