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.