Skip to content

Add chr prefix to gnomAD.r2.1.1-sv.vcf.gz for svdb/query test - #2268

Merged
ramprasadn merged 1 commit into
nf-core:modulesfrom
ramprasadn:fix-svdb-gnomad-chr-prefix
Sep 11, 2026
Merged

ramprasadn merged 1 commit into
nf-core:modulesfrom
ramprasadn:fix-svdb-gnomad-chr-prefix

Conversation

@ramprasadn

@ramprasadn ramprasadn commented Sep 11, 2026 •

Copy link
Copy Markdown

Summary

gnomAD.r2.1.1-sv.vcf.gz used bare chromosome names (22) while sv_query.vcf.gz (the query VCF used against it in nf-core/modules' svdb/query test) uses chr-prefixed names (chr22). This mismatch means the two files describe the same genomic positions but never lexically match.

svdb <= 2.8.4 silently stripped chr prefixes on read, so the mismatch didn't matter in practice. svdb 2.9+ intentionally removed that normalization, so with svdb 2.12.0 the two files no longer overlap at all, and the svdb/query test's only annotated hit silently disappeared.

This PR fixes that by adding chr prefix to gnomAD.r2.1.1-sv.vcf.gz's CHROM column, ##contig header IDs, and CHR2 INFO values so it's consistent with sv_query.vcf.gz, restoring the intended annotation match under current and future svdb versions.

This file is only referenced by nf-core/modules' svdb/query test, so this change is isolated.

sv_query.vcf.gz (the query input in svdb/query's nf-test) uses chr-prefixed
contig names (chr22), while this gnomAD SV database used bare names (22).
svdb <= 2.8.4 silently normalized away chr prefixes on read, so the
mismatch never mattered, but svdb 2.9+ dropped that normalization
(J35P312/SVDB@a520684), so the two files no longer overlap and the
query test's only annotated hit was silently lost. Align the naming so
the fixture actually exercises a match.

@beatrizsavinhas beatrizsavinhas left a comment

Copy link
Copy Markdown

Choose a reason for hiding this comment

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

Very nice description of the investigation! And great that you already checked that the file is not used for anything else too! 🌟

@ramprasadn
ramprasadn merged commit 9e3d19a into nf-core:modules Sep 11, 2026
2 checks passed
@ramprasadn
ramprasadn deleted the fix-svdb-gnomad-chr-prefix branch September 11, 2026 13:46
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.

2 participants