bac95a147f49cd331052e597006e04b3deee40fc max Wed Apr 22 10:43:20 2026 -0700 lrSv/srSv: human-readable SV type filter labels, script cleanups Add human-readable labels to the supertrack-level svType filter on both the lrSv and srSv supertracks using the "CODE|CODE (Long name)" filterValues syntax: DEL -> "DEL (Deletion)", INS -> "INS (Insertion)", etc. Labels keep the short code up front so users can match what hgTracks shows next to each feature. Also sweep in the in-progress converter/as-file cleanups under scripts/lrSv/ and scripts/srSv/ (introduction of lrSvCommon.py helpers, consistent insLen / svLen / AC column naming, tightened field-description text) that had been piling up as an unstaged working tree. refs #36258 diff --git src/hg/makeDb/trackDb/human/hprc2Sv.html src/hg/makeDb/trackDb/human/hprc2Sv.html index ac1ba979642..11a002d914b 100644 --- src/hg/makeDb/trackDb/human/hprc2Sv.html +++ src/hg/makeDb/trackDb/human/hprc2Sv.html @@ -1,111 +1,118 @@

Description

This track shows structural variants (SVs) derived from the Human Pangenome Reference Consortium (HPRC) release-2 pangenome graph. The graph was built with minigraph-cactus from PacBio HiFi haplotype-resolved assemblies of 233 samples (including T2T-CHM13 and the diverse 1000 Genomes Project sample set). HPRC releases one VCF per reference path (GRCh38 and T2T-CHM13); we display both natively on the corresponding UCSC assembly (hg38 and hs1). Variants were extracted from the graph with vg deconstruct and decomposed into atomic alleles with vcfwave (WFA2-lib).

The hg38 track contains 1,483,114 SV-sized alleles (length ≥ 50 bp) split by type: 1,106,190 insertions, 192,597 deletions, 178,178 complex alleles and 6,149 inversions. The hs1 track is built from the parallel T2T-CHM13 wave VCF. Each row carries the allele count, allele frequency, number of samples with data and the snarl-nesting level of the variant in the pangenome decomposition tree.

Display Conventions and Configuration

Items are colored by SV type:

Insertions (INS) - blue
Deletions (DEL) - red
Complex alleles (COMPLEX) - purple
Inversions (INV) - orange

Insertions are placed at the insertion site with a width of 1 bp; deletions, complex alleles and inversions span the affected reference interval. Filters are available for SV type, SV length, allele frequency and snarl level (0 = top-level bubble; higher values are nested within parent bubbles).

Methods

-The HPRC v2.0 minigraph-cactus pangenome was downloaded from the HPRC S3 -release bucket: -

hg38: -hprc-v2.0-mc-grch38.sv.gfa.gz (graph) and -hprc-v2.0-mc-grch38.wave.vcf.gz (wave-decomposed VCF)
hs1: -hprc-v2.0-mc-chm13.wave.vcf.gz

-The VCF is the result of running vg deconstruct on the graph with -the corresponding reference path (GRCh38 or T2T-CHM13) and then -vcfwave / WFA2-lib to split complex multi-allelic records into -atomic alleles with per-allele TYPE and LEN fields. +HPRC release-2 is an open data release (not yet accompanied by a formal +peer-reviewed publication) built from PacBio HiFi haplotype-resolved +assemblies of 233 samples, including T2T-CHM13 and a diverse 1000 Genomes +Project panel. The pangenome graph was built with Minigraph-Cactus against +both GRCh38 and T2T-CHM13 reference paths; variants were extracted from +the graph with vg deconstruct and then decomposed into atomic +alleles with vcfwave / WFA2-lib, yielding per-allele TYPE and LEN +fields. For this track, each ALT in the wave VCF was emitted as its own +BED row, retaining alleles with |LEN| ≥ 50 bp or the INV flag; +allele counts, frequencies, sample counts and snarl levels are taken +directly from the per-allele INFO fields. On hg38 this yields 1,483,114 +SV-sized alleles (1,106,190 insertions, 192,597 deletions, 178,178 complex +alleles and 6,149 inversions); the hs1 track is built from the parallel +T2T-CHM13 wave VCF. Sample-list and assembly provenance for the graph are +maintained at HPRC in + +hprc_intermediate_assembly/alignments_v2.0.csv.

-For display here, the wave VCF was streamed and each ALT was emitted as -its own BED row. Alleles were retained if their absolute length was -≥ 50 bp or if the record carried the INV flag (inversions may -be shorter). Allele counts, frequencies, and sample counts are taken -directly from the per-allele INFO fields. +The HPRC v2.0 Minigraph-Cactus graph and wave-decomposed VCFs were +downloaded from the HPRC S3 release bucket: + +hprc-v2.0-mc-grch38.wave.vcf.gz (hg38) and + +hprc-v2.0-mc-chm13.wave.vcf.gz (hs1).

-A pointer to both the GRCh38 and CHM13 pangenome files (and the list of -assemblies that went into the graph) is maintained by HPRC at -human-pangenomics/hprc_intermediate_assembly -alignments_v2.0.csv, which links to both the hg38 and -CHM13/hs1 VCFs (and the underlying graph files) used for this track. +The step-by-step build commands (download, format conversion, bigBed build) +are recorded in the UCSC makeDoc for this track container: + +doc/hg38/lrSv.txt and + +doc/hs1/lrSv.txt. The conversion scripts and autoSql schemas live in + +makeDb/scripts/lrSv.

Data Access

The data can be explored interactively in table format with the Table Browser or the Data Integrator, and accessed programmatically through our API, track=hprc2Sv.

The bigBed is available from our download server for both assemblies:

GRCh38: hg38 hprc2.bb
T2T-CHM13: hs1 hprc2.bb

Example: bigBedToBed http://hgdownload.soe.ucsc.edu/gbdb/hg38/lrSv/hprc2.bb -chrom=chr21 -start=0 -end=100000000 stdout.

The original pangenome graph and the wave-decomposed VCF are available from the HPRC public S3 bucket, as linked from the HPRC release-2 announcement.

Credits

Thanks to the Human Pangenome Reference Consortium for building and publicly releasing the release-2 minigraph-cactus pangenome.

References

HPRC release-2 data is not yet described in a formal peer-reviewed publication. See the Human Pangenome Project release announcement for background and data-access details: HPRC data release 2.