custom header test
sim9_single(speciesData = matrix(rbinom(100, 1, 0.5), nrow = 10))
Function for a single iteration of the sequential swap