ChIPQC read depth doesn't match BAM files

Hello,

I am using ChIPQC, and the report is listing read depths that are much lower than I can plainly see in the MACS2 output. For instance, one of my BAM files reports ~13 million fragments, but then only ~550,000 reads in the ChIPQC report. Any suggestions as to why this might be happening?

Thanks,

library(ChIPQC)
library(SerialParam())
register(SerialParam())
chipObj <- ChIPQC(B)
#B is my sample sheet
ChIPQCreport(chipObj, reportName="ChIP QC report: Rpb1 ChIPseq broadpeak 3 rep")

chipseq chipqc

• 1,018 views • link

updated 4.0 years ago by Rory Stark

• written 4.1 years ago by bertb •

1 answer

Which part of the report do you see the lower number of reads? By default, ChIPQC only checks a subset of the reads (those mapping to the first chromosome on which there is a peak) for its QC metrics rather than processing all of the reads. Generally, this runs much more quickly and gives indicative results.

• 0 views • link

written 4.0 years ago by Rory Stark