Hello
I have the peaks file from macs2 output, I want to quantify how many reads of each peak from the bam file. Is there a package or bash/R/Python code that can do this quantification?
Thanks
•
written
by
Dan •
2 answers
Lots of great information on this site and even a few published best-practices depending on your ultimate goal.
Best practices for differential accessibility with ATACseq here, which includes some R scripts for using csaw and lots of great background.
You could also quantify peaks directly from bam files directly.
In any case, if you are assessing differential accessibility between conditions, ensure you at least consider normalization to account for differences in library output. Deciding how to define consensus peaks across conditions also dramatically impacts DA calls (see Reske et al 2020).
Log in to answer this question.
-
Overestimate of ATAC-seq differential peaks
written by Darryl •
-
Are there pipelines to analyze the Quality Controls for ATAC-seq and ChIP-seq data?
written by Dan •
-
Chipseq data peak calling issue
written by Adyasha •
-
How to quantify the reads of each peak for ChIP-seq data?
written by Dan •
-
ATAC DiffBind: use peak or summit file?
written by a_bis •
-
MACS: peak calling from snATAC data
written by abiuma •
-
Normalizing ChIP-seq and ATAC-seq data?
written by star
-
How to visualize selected few ChIPseq peaks
written by Mike
-
How can I increase the number of peaks using macs2 peak caller?
written by GK1610
-
macs and bedtools
written by marina-orlova
