Produce a featureCounts matrix for consensus peak sets

[Tomrrr1]

I spoke with Maria, and she said it would be helpful to output a featureCounts matrix for each consensus peak set. A featureCounts matrix will contain, for each consensus peak, the number of reads in the constituent BAM files that overlap with that peak. The matrix should be like this:

	bam1	bam2	bam3
peak1	20	30	40
peak2	10	40	40
peak3	80	90	50