Freq := CreateArray(1..20,0):
Tot := sum(sum(CountMatrix)):
for i to 20 do
  Freq[i] := sum(CountMatrix[i])/Tot:
od:
print(Freq);
[0.08477395, 0.05103334, 0.03837665, 0.05740129, 0.01256165, 0.03471746, 
0.06883297, 0.07907659, 0.02077239, 0.06813598, 0.08906677, 0.06501537, 
0.02318017, 0.03823936, 0.04036729, 0.05668318, 0.05721648, 0.00871765, 
0.02785317, 0.07797831]
PercentageMutation := 100 * sum(Freq[i] * (1 - MutationMatrix[i,i]),i=1..20);
PercentageMutation := 44.5302