calculate coverage in samtools depth file