simulations from T. cristinae GBS data