Modify VCF to support GRC assembly model (GRCh37, GRCh38, GRCm38) about hts-specs HOT 2 OPEN

On the representation of alt contigs, I think we should develop a best practice before modifying the spec. What is the intended output from variant callers? Is it practical for callers to generate such output? How downstream tools are supposed to use the vcf?

Specifically, you proposed to add HT, but in my experience, alt contigs frequently recombine with each other, which makes the tag not applicable most of times. In addition, how are we supposed to use ALTLOCS? If we know a locus overlapping an alt contig, what can we do with it?

We will be clearer about the answers and then develop the right spec when more researchers get experiences on h38. Tools determine the adoption of alt contigs. It is not urgent to change the spec.

from hts-specs.

I think this is a bit of a chicken and egg problem. If we want variant callers to be able to use the Alt loci, we need to be able to express the variants in VCF. This doesn't work well with the current spec (see how dbSNP distributes data).

I think the issue is, there are multiple ways to use VCF- it is just a reporting tool. dbSNP uses it to dump data- so you want to report all genomic contexts for a given SNV. An argument could be made that in the context of an individual genome, you may only want to report one context for a SNP- but how do you handle that when you have multiple samples in the VCF? I fear that decision making will be hard.
I agree the trying to define some best practices is useful.
To attempt to address some specific issues:

• knowing an alt-locus is allelic with a region on the chromosome let's you put the alt-locus in chromsome context which is useful for reporting. Granted, this could be done as some sort of post-processing step, but then you'd have to convert the data to some other format (which may be OK, but likely inconvenient for everyone who wants VCF).
• I agree there are few loci with named haplotypes, but for the ones that exist, this is useful information.

This is really meant to start the discussion about how we want to represent variation on GRCh38. It will be good to have some concrete examples.

from hts-specs.