R script for generation coverage plots for BAM files