Given a population of permutations, this function mutates all individuals by randomly selecting two indices, and reversing the respective sub-permutation.
mutationPermutationReversal(population, parameters = list())
mutated population
List of permutations
list of parameters, currently only uses parameters$mutationRate, which should be between 0 and 1 (but can be larger than 1). The mutation rate determines the number of reversals performed, relative to the permutation length (N). 0 means none. 1 means N reversals. The default is 1/N.