MACS2 Total Reads Within a Peak?
Hi,
I wanted to know if there is a way to determine the number of TOTAL reads within a peak being called by the MACS2 command callpeaks?
I am aware that it gives the pileup of reads for the summit peak but I want to know the total reads for each peak being outputted.
P.S I have already tried coveragebed from bedtools.
Thank you!
• 1,788 views
•
link
written
by
takamichael •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Get the strongest TF binding regions for Chip-Seq
written by Alexander Sun •
-
Recovering summit position from ENCODE narrowPeak (BED6+4) format
written by Ready2Rapture •
-
MACS2 Summit TSSPlot dividing BDG into bins
written by Seigfried
-
Why pileup values of MACS2 output have decimal values?
written by Aspire
-
Macs2 output shows only single basepair peak instead of a section of DNA
written by fdewards •
-
Handling replicates in Chip-seq
written by chaudharyc61 •
-
Data value in WIG file
written by Rishikesh Lotke •
-
macs2: how to retrieve normalized counts
written by olechnwin •
-
Discrepancy with MACS2.0 and MACS2.1 broad peaks
written by ivivek_ngs
-
Why does MACS2 report multiple records for the same peak region?
written by James Ashmore
and ?
I am not seeing the total reads in the .cnt.bed file I just get the following information.
chr4 93874012 93874062 HWI-ST926:425:C5WF1ACXX:4:1101:1118:1982 60 - 93874012 93874062 0,0,0 1 50, 0, 0