powered by
countBamInGRanges.exomeDepth(bam.file, index = bam.file, granges, min.mapq = 1, read.width = 1)