minor update:
- changed ordering method to "position weighted summation" (PWS) and moved to Rscript (potentially fixed banding in the heatmap)
PWS weights the bins based on their location in the matrix, where the farthest bin from the matrix midpoint are weighted the highest and the weight decreases as it approaches the middle of the matrix in either direction.