You would think NIST/GiaB would explicitly state the reference build for their putative gold standard SV calls. But I can't see it anywhere
I'm assuming it's in GRCh37???
Anyone else knows? I can't find any information in READMEs or in the published svclassify paper.
2 answers
If you are referring to this (svclassify: a method to establish benchmark structural variant calls), then —yes— they aligned to NCBI's GRCh37 reference genome:
...raw reads were mapped to the National Center for Biotechnology Information (NCBI) build 37 using the Burrows-Wheeler Aligner (BWA) “bwa mem” v.0.7.5a with default parameters
Variants were mapped to human reference coordinates (NCBI build 37) by walking the read overlap graph in both directions until an “anchor” read, where a continuous 65 bps matches the reference, denoted the beginning and end of each variant.
If you are referring to the original published works (Extensive sequencing of seven human genomes to characterize benchmark reference materials), then, the same:
The sequencing data were aligned by bwa mem6 against b37 human decoy reference genome.
It looks like it's GRCh37: On p. 11 of Parikh et al. (2014) in the first paragraph of "Methods", they say that they mapped the Platinum Genomes 2x100bp HiSeq data to NCBI "build 37" using bwa mem v.0.7.5a with default parameters, and that aligned (meaning, presumably, to the same reference) BAM files were publicly available for the Illumina 250bp, PacBio and Moleculo data. Also near the top of the next page they write that the Spiral Genetics variants in category C were mapped to NCBI build 37, though it's not yet clear to me how the subheading this falls under corresponds to the 11 rows of Table 2.
Log in to answer this question.