Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Classifying reversion mutations #94

Open
hyanwong opened this issue Feb 8, 2023 · 1 comment
Open

Classifying reversion mutations #94

hyanwong opened this issue Feb 8, 2023 · 1 comment

Comments

@hyanwong
Copy link
Contributor

hyanwong commented Feb 8, 2023

We thought it would be a good idea to get a handle on revertants (i.e. mutations whose derived state is equivalent to the state of their parent mutation). We think these are mostly erroneous, so it would be good to identify what's going wrong. We could classify them by the nodes they are above (e.g are they above a recombinant, or a reversion push)

@jeromekelleher
Copy link
Owner

jeromekelleher commented Feb 8, 2023

We have code for classifying mutations as reversions here so I think we want two things:

  1. A plot to summarise the path length and time back to the thing its reverting (i.e., either the parent mutation, or the root. All mutations have a time in sc2ts). Both are useful, as we'd like to know if the majority of reversions revert a mutation one node above the immediate parent, e.g. So, I guess a method giving a histogram of path length and time difference would be a useful thing to add to the TreeInfo
  2. A plot to classify the types of nodes that these reversions happen on. A simple first thing to do would be to do a bar chart with the categorial axis is the node flags. This would also be a useful thing to add to the TreeInfo

I think it's handy to add these QC-type functions to the notebooks/qc-template as they are done, so that we have something that you can run on a new set of trees to get an overall view of the quality.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants