the construction of the muse-gaut codon rate matrix now seems to work