Substantial update to matrix construction, specifically scaling procedure.
All nucleotide, amino-acid, empirical codon, and mechanistic codon model matrices are scaled as usual ("persite") scaling.
Heterogeneous codon model matrices are scaled such that the mean substitution rate for the _mean_ dN/dS across all categories is equal to 1.
Mutation-selection model matrices are scaled such that the mean _neutral_ substitution rate is equal to 1.