added kmer substitution analysis - heatmap commented out due to failing on sparse matrices as found in test outputs