Calculate coverage in plot_read_distribution