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

HGVS - reference validation now fails #491

Closed
davmlaw opened this issue Sep 21, 2021 · 9 comments
Closed

HGVS - reference validation now fails #491

davmlaw opened this issue Sep 21, 2021 · 9 comments
Assignees

Comments

@davmlaw
Copy link
Contributor

davmlaw commented Sep 21, 2021

Searching or submitting classification for these now fails validating that = not in GATC - issue with validating our internal representation

NM_005629.3(SLC6A8):c.1403A=
NM_018090.5:c.462=

@davmlaw davmlaw added the HGVS label Sep 21, 2021
davmlaw added a commit that referenced this issue Sep 21, 2021
@EmmaTudini
Copy link
Contributor

@davmlaw Not quite sure what the intended output is here. NM_005629.3(SLC6A8):c.1403A= can be imported, matched and searched for but NM_018090.5:c.462= fails import and searching (error about not having a base)

@davmlaw
Copy link
Contributor Author

davmlaw commented Nov 7, 2021

Fix was to allow creation of a reference HGVS, didn't work due to a bug.

The error about not having a base is due to our enforcing that we specify the reference in our HGVSs - so that's the right behavior the example I listed was wrong and forgot about that rule

@TheMadBug
Copy link
Member

NM_018090.5:c.462= in prod complains about
HGVS: Invalid HGVS name "NM_018090.5:c.462=": alt=NN contains non-standard (A,C,G,T) bases: NN (GRCh38)

Is this expected?

@davmlaw
Copy link
Contributor Author

davmlaw commented Nov 16, 2021

Yeah I think we should just take the reference base here, if we want to stop people from doing that can always give a warning / flag etc but should resolve it

Needs latest pyhgvs

@EmmaTudini
Copy link
Contributor

@davmlaw What has changed since your last comment?

@davmlaw davmlaw added the rework label Dec 2, 2021
@davmlaw
Copy link
Contributor Author

davmlaw commented Dec 2, 2021

Fails due to https://github.com/SACGF/variantgrid_private/issues/2142

"NM_018090.5:c.462=" now fails with "Coordinate span (1) not equal to ref length 2"

@davmlaw
Copy link
Contributor Author

davmlaw commented Dec 9, 2021

This was due to the HGVS library not handling reference variants without reference bases supplied - my validation change just exposed it

Raised as bug counsyl/hgvs#61

I have fixed it in our version of HGVS library, so upgrade via:

sudo python3 -m pip install --upgrade git+https://github.com/SACGF/hgvs#egg=pyhgvs

@davmlaw davmlaw assigned EmmaTudini and unassigned davmlaw Dec 9, 2021
@davmlaw
Copy link
Contributor Author

davmlaw commented Dec 9, 2021

James upgraded Shariant test - tested submitting "NM_018090.5:c.462=" works fine just raises a matching warning whne it inserts the gene symbol and reference base

Search also works (not sure if it was broken before)

@EmmaTudini
Copy link
Contributor

Testing:

  • imported and searched for NM_001354689.1(RAF1):c.1=
    Expected output: Should insert ref base and search correctly
    Actual output: passed

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

No branches or pull requests

3 participants