Hi, Jiseon! Well, here's an example of a chromosome with only selected mutations happening in a 10kb chunk in the middle:
initialize() {
initializeMutationRate(c(0, 1e-7, 0), ends=c(50000, 60000, 100000));
initializeMutationType("m1", 0.5, "f", 0.1);
initializeGenomicElementType("g1", m1, 1.0);
initializeGenomicElement(g1, 50000, 60000);
initializeRecombinationRate(1e-8);
}
But, maybe I didn't understand the question? If not, a short, simplified script might help?
-peter