README.md

VCF positions for precise SV calls

Status: In progress

Overview

This is guide to clarify the usage of positions for POS, END, and MATE_POS fields for structural variants inferred at basepair resolution and rendered to VCF 4.1 (or higher). The VCF spec leaves many clues on this topic but there is considerable room for off-by-one errors at VCF generation and parse stages.

Conventions used in this guide

'Left' and 'Right'

All discussions below are based on the forward DNA strand, with "left" used to denote positions with a lower reference coordinate relative to "right" positions.

Breakend Homology

Most SVs will have some amount of breakpoint homology, its representation is reasonably well defined in the VCF spec already. The guide below follows the NGS community convention of standarizing variants by describing the local breakend standardized to the left-most position within its homology range.

Examples

All examples are drawn from public Platinum Genomes data for NA12878

Exact VCF position examples for common SV types:

Deletions (using VCF symbolic <DEL> allele)

Left end

POS indicates right-most position before crossing the left-shifted deletion left breakend:

Right end

END indicates right-most position before crossing the left-shifted deletion right breakend