I know that
bedtools genomecov -ibam xx.bam -bg
would report genome coverage in BEDGRAPH format
What if I need the coverage not for the whole genome, but just for specific regions (specific exomes) ? how can I achieve this ? my input files are a bam file and a .bed file containing the start and end positions of exomes.