R scripts for sampling loci at random from across reference genome