Finally, I ran the script
evaluate.pywhich does what we said we were going to do last time. We count non-synonymous sites, and synonomous and non-synonomous differences between each pair of sequences.
There isn't a lot to say about how it works. Except of course that the whole point of the exercise is to notice the huge difference in non-synonymous changes allowed in IL6, say, compared to the others.
Here are the results: