initialize()
{
    setSeed(23);
	defineConstant("L", 1e2);
	initializeSLiMOptions(keepPedigrees=T, nucleotideBased=T);
    initializeTreeSeq();
	initializeAncestralNucleotides(randomNucleotides(L));
	initializeMutationTypeNuc("m1", 0.5, "f", 0.0);
	initializeGenomicElementType("g1", m1, 1.0, mmJukesCantor(4e-2));
	initializeGenomicElement(g1, 0, L-1);
	initializeRecombinationRate(1e-2);
}

1 late() { 
    // output immediately! Nothing should have changed.
    sim.readFromPopulationFile("restart_nucleotides.init.trees");
    sim.treeSeqOutput("restart_nucleotides.trees");
}

20 {
    catn("Done.");
    sim.simulationFinished();
}
