code for getting per-bp coverage from read alignment files