List of Nucleotides for each Node

[AlirezaPa]

Hi msprime community,

I have a simple question. Is there any way to see the nucleotides (ACCGGTC...) for the following command and tree. I expect to see 10 nucleotides for each node. No mutation rate and no recombination rate.

ts = msprime.sim_ancestry(2, sequence_length=10)

Thank you,
Alireza

[jeromekelleher]

Hi @AlirezaPa, welcome 👋

The sim_ancestry function only simulates the ancestral trees for a sample (like shown above). To get sequence data we have to simulate mutations on top of this ancestry, using the sim_mutations function. Hope this helps!