Why ranges of narrowPeak file and that of bed file are not identical

Question

[Wang]

Hi everyone!

I did peak call using MACS2 to analysis a Chip-seq data and I found the ranges in the narrowPeak file is not same as that in the bed file.what dose this mean?and why?Thanks!

more ff_peaks.narrowPeak

chr1    3044565 3044674 ff_peak_1       97      .       3.97001 12.89295        9.76606 21
chr1    3050086 3050252 ff_peak_2       332     .       13.90061        38.02962        33.29285        77
chr1    3050434 3050592 ff_peak_3       91      .       6.33852 12.18861        9.14707 69
..

 

more ff_summits.bed 

chr1    3044586 3044587 ff_peak_1       9.76606
chr1    3050163 3050164 ff_peak_2       33.29285
chr1    3050503 3050504 ff_peak_3       9.14707
..

[Ram]

Please use the formatting bar (especially the code option) to present your post better. You can use backticks for inline code (`text` becomes text), or select a chunk of text and use the highlighted button to format it as a code block. I've done it for you this time.

[Wang]

ok,I got it!.thank you!

[Wang]

I run macs2 this way

macs2 callpeak -t ff.bam -c input.bam -n ff -f BAM -g mm

ff.bam and input.bam was sorted by samtools

samtools sort -@ 32 1.bam  >ff.bam
samtools sort -@ 32 2.bam  >input.bam

[GenoMax]

Please use ADD COMMENT/ADD REPLY when responding to existing posts to keep the threads logically organized. SUBMIT ANSWER is for new answers to original questions.

[Wang]

Sorry! I'm newbie here. Maybe I should edit the original post and append these informations.

Answer 1

[ATpoint]

The narrowPeak file contains the entire interval of the peak, the ff_summits.bed the peak summits so the base with the highest read pipeup.

[Wang]

Thank you,I really appreciate it! Is there any detailed materials or article I can read? I failed to find these information in the official README.md

[ATpoint]

Check the data format definitions at UCSC and the macs documentation, the latter explicitely explains the outputs.

[Wang]

OK,thank you.May be I missed some useful information.I'll read these again.

[Jingyue]

I have been using narrowPeak files for downstream analysis (annotation, motif finding, differential peak calling), should I use the summits.bed or narrowPeak files?

I'm asking this because narrowPeak files have many peaks start-end ranges if use --call-summits in MACS2. So I only keep 1 unique peak as the entire interval of the peak and plot Venn diagram among samples. The number will be different when I use summits.bed because it has more lines (1 peak with multiple sub-summit intervals).

[buffealo]

Hello, have you got any answer for your question from any different source? When I try to visualize overlapping peaks with the summits.bed results do not seem good