fixed erroneous command to calculate coverage on regions; updated R script to plot heatmap/region coverage