TAD and DIST definition should use absolute values when calculating differences between chromosome coordinates

avoiding PLOT_CHR for now, commented output and call to plotting script; increased RE input reads from 250000 to 750000

if two families are not observed in an interaction, that should give NA and not (log) 0 in the heatmap

rewrote to ggplot, unified treatment of NAs (when either of obs/exp is 0); only 12 iterations of sample are now part of the output; modified code for when insufficient amount of data is present

Matrix ordering fixes for heatmaps, now they are always square, no need to test for it when plotting a matrix created with reorder_mtx