Question: Using FeatureCounts for ChIP-seq normalised files?
gravatar for a.rex
15 months ago by
a.rex200 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 • 650 views
ADD COMMENTlink modified 15 months ago by Devon Ryan94k • written 15 months ago by a.rex200
gravatar for Devon Ryan
15 months ago by
Devon Ryan94k
Freiburg, Germany
Devon Ryan94k 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 15 months ago by Devon Ryan94k
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: 1176 users visited in the last hour