Calculate the coverage from a bam file