Compute mutation densities for subgenic regions

[FerriolCalvet]

This can take advantage of the files that are computed for Omega particularly in the expandregions step.

See the expandregions output for some examples on these files.

• get the mutdensity modules run with the output of the expanded regions used by omega