-
Notifications
You must be signed in to change notification settings - Fork 597
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
Fix RAW_MQ header inconsistencies after reblocking #6276
Conversation
I can confirm that this fixes the production WGS reblocking issue. |
Lots of files are changed because I changed the key name for RAW_MQ to include |
@tlangs once we cut a release with this PR we can update and merge your reblocking workflow in dsde-pipelines. |
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.
Looks good to me. Or at least as good as nasty vcf header shenanigans can look.
@@ -158,6 +158,7 @@ public void onTraversalStart() { | |||
final Set<VCFHeaderLine> headerLines = new HashSet<>(inputHeaders); | |||
// Remove GCVFBlocks, legacy headers, and annotations that aren't informative for single samples | |||
headerLines.removeIf(vcfHeaderLine -> vcfHeaderLine.getKey().startsWith(GVCFWriter.GVCF_BLOCK) || | |||
(vcfHeaderLine.getKey().equals("INFO")) && ((VCFInfoHeaderLine)vcfHeaderLine).getID().equals(GATKVCFConstants.RAW_RMS_MAPPING_QUALITY_DEPRECATED) || //remove old (maybe wrong type) and add new with deprecated note |
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.
I hate these classes so much.
I had a copy-paste guffaw when I deprecated the RAW_MQ header that changed the type to be incorrect. After reblocking, the new header is correct, but conflicting and GenomicsDB is understandably not pleased. Now there is a single, correct RAW_MQ header after reblocking.