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

Handle = in synonymous variant names #111

Merged
merged 1 commit into from
Feb 5, 2021
Merged

Conversation

susannasiebert
Copy link
Collaborator

CIViCpy check VCF fields for containing = and fails when such a variant is encountered because having a field value contain this character would break VCF specs since = has a special meaning for parsing. However, synonymous variants will legitimately contain this characters in their name (e.g. E55= (c.165G>A)). This PR implements a solution for this problem by replacing the = in the variant with the wildtype amino acid(s).

Closes #110.

Copy link
Member

@ahwagner ahwagner left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This seems like a good implementation.

Do we need to do any special handling for hgvs_p() or hgvs_g()? Synonymous variants are handled with = in HGVS as well.

@susannasiebert
Copy link
Collaborator Author

hgvs expressions are being encoded so they don't have this problem but we can switch to this implementation to make them look "nicer".

@ahwagner
Copy link
Member

ahwagner commented Feb 5, 2021

Okay. I'm fine with it either way–if you are happy with the current solution, I think it is good to merge!

@susannasiebert susannasiebert merged commit de055bd into master Feb 5, 2021
@susannasiebert susannasiebert deleted the synonymous_variants branch February 5, 2021 22:03
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

Successfully merging this pull request may close these issues.

synonymous CIViC variants with = in the name (e.g. E55= (c.165G>A)) cause VCF creation to fail
2 participants