script to count reads from genomic features in a number of cell types