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