If you can have a bedfile with regions where you would like to count, bedtools multicov doest the job fast and well. Required an indexed bam....