Works! Plot is mutational spectra control/dosed. Pasted from Ipython... A lot of cleanup needed