Using FeatureCounts for ChIP-seq normalised files?
1
1
Entering edit mode
2.4 years ago
a.rex ▴ 250

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.

deeptools ChIP-Seq • 1.2k views
ADD COMMENT
3
Entering edit mode
2.4 years ago

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 COMMENT

Login before adding your answer.

Traffic: 2525 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6