List of Nucleotides for each Node #1890
-
|
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, |
Beta Was this translation helpful? Give feedback.
Answered by
jeromekelleher
Nov 8, 2021
Replies: 1 comment
-
|
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! |
Beta Was this translation helpful? Give feedback.
0 replies
Answer selected by
AlirezaPa
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
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!