-
Notifications
You must be signed in to change notification settings - Fork 45
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
Split variant flags #1623
Split variant flags #1623
Conversation
@rileyhgrant I've decided to put the VA/VRS-related fixes in the other VA topic branch I've been working on, which I'll rebase on top of this, hence this PR doesn't include VA fixes. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
LGTM, very nicely done, I especially like the care taken to the testing, and to a neatly organized series of commits.
Just one minor comment about the naming of a common region flag (based on the docs, it might be named "non_par", rather than "par", but I wasn't quite sure and wanted to double check).
Nice one!
const baseExomeFlags = variant.exome?.flags || [] | ||
const baseGenomeFlags = variant.genome?.flags || [] | ||
|
||
const regionalFlags = ['par', 'segdup', 'lcr'].filter( |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Based on the hail tables help text, it appears as though the par flag might be named "non_par"?
I did a quick grep and didn't see any place where we rename this flag in our pipeline, so just wanted to check, this flag seems to potentially differ in meaning from in the SV datasets where there is a "par" flag.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
The v4 data pipeline does address this,
gnomad-browser/data-pipeline/src/data_pipeline/datasets/gnomad_v4/gnomad_v4_variants.py
Line 286 in b804dbe
hl.or_missing( |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Ah whoops, I naively only grepped for "non_par", my bad for not finding that, and nice catch!
6ccd21e
to
9447737
Compare
This updates the API and frontend to support flags both on a variant generally, or on the exome and/or genome specifically, and fixes a number of related bugs related to display of flags.
fixes #1614