Question: Using FeatureCounts for ChIP-seq normalised files?
gravatar for a.rex
23 months ago by
a.rex240 wrote:

I have a set of bed intervals (that correspond to genomic regions of open ATAC NFR regions). I also have ChIP-seq bigwig files for a histone mark that have been normalised to input with Deeptools2' bamCompare. I would like to use FeaturesCount to count the number of normalised reads in each interval. However, I need a bam file as an input. How can I convert from Bigwig or use another software.

chip-seq deeptools • 1.0k views
ADD COMMENTlink modified 23 months ago by Devon Ryan97k • written 23 months ago by a.rex240
gravatar for Devon Ryan
23 months ago by
Devon Ryan97k
Freiburg, Germany
Devon Ryan97k wrote:

I suggest that instead of trying to get counts from the bigWig files you instead use alignmentSieve to generate NFR BAM files. DESeq2 is then OK if you have particular regions you're interested in (e.g., peaks of NFR signal) but otherwise give CSAW a try. It's quite good at finding differentially accessible regions.

ADD COMMENTlink written 23 months ago by Devon Ryan97k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1315 users visited in the last hour