adapt script for array of paired reads
use regular hg19 ref for script