VCF (Variant Call Format)
The Variant Call Format (VCF) was designed as part of the 1000 Genomes Project as a standardized means to report genetic variation coming from SNP, INDEL, genomic rearrangements and structural variant detection programs and ownership has been subsequently transferred to Global Alliance for Genomics and Health Data Working group file format team.
The format can be used to show information about all kinds of genomic variation. The Variant Call Format (VCF) specifies the format of a text file used in bioinformatics for storing gene sequence variations. The format has been developed with the advent of large-scale genotyping and DNA sequencing projects.
VCF is a text file format most likely stored in a compressed manner. It contains meta-information lines, a header line, and then data lines each containing information about a position in the genome. The format also has the ability to contain genotype information on samples for each position.
The text files in VCF are saved with .vcf or .vcf.gz extensions.
Every VCF file has three parts in the following order:
- Meta-information lines
- One header line
- Data lines contain marker and genotype data (one variant per line). A data line is called a VCF record.
Each VCF record has the same number of tab-separated fields as the header line. The symbol “.” is used to denote missing data.
Each meta-information line must have the form ##KEY=VALUE and cannot contain white-space. The first meta-information line must specify the VCF version number. Additional meta-information lines are optional, but are often included to describe terms used in the FILTER, INFO, and FORMAT fields.
The header line syntax
The header line names the 8 fixed, mandatory columns. These columns are as follows:
There are 8 fixed fields per record. All data lines are tab-delimited. In all cases, missing values are specified with a dot (“.”). Fixed fields are:
- CHROM chromosome: an identifier from the reference genome. All entries for a specific CHROM should form a contiguous block within the VCF file.(Alphanumeric String, Required)
- POS position: The reference position, with the 1st base having position 1. Positions are sorted numerically, in increasing order, within each reference sequence CHROM. (Integer, Required)
- ID semi-colon separated list of unique identifiers where available. If this is a dbSNP variant it is encouraged to use the rs number(s). (Alphanumeric String)
- REF reference base(s): Each base must be one of A,C,G,T,N. Bases should be in uppercase. Multiple bases are permitted. The value in the POS field shows the position of the first base in the String. For InDels, the reference String must include the base before the event. (String, Required).
- ALT comma separated list of alternate non-reference alleles called on at least one of the samples. Options are base Strings made up of the bases A,C,G,T,N, or an angle-bracketed ID String (”<ID>“). If there are no alternative alleles, then the missing value should be used. Bases should be in uppercase. (Alphanumeric String; no whitespace, commas, or angle-brackets are permitted in the ID String itself)
- QUAL phred-scaled quality score for the assertion made in ALT. If ALT is ”.” (no variant) then this is -10log_10 p(variant), and if ALT is not ”.” this is -10log_10 p(no variant). (Numeric)
- FILTER filter: PASS if this position has passed all filters, i.e. a call is made at this position. Otherwise, if the site has not passed all filters, a semicolon-separated list of codes for filters that fail. (Alphanumeric String)
- INFO additional information: INFO fields are encoded as a semicolon-separated series of short keys with optional values in the format: <key>=<data>[,data]. (Alphanumeric String)
Common INFO fields:
- AA ancestral allele
- AC allele count
- AF allele frequency
- AN total number of alleles in called genotypes
- BQ RMS base quality
- CIGAR cigar string describing how to align an alternate allele to the reference allele
- DB dbSNP membership
- DP combined depth across samples
- END end position of the variant
- H2 membership in hapmap2
- MQ RMS mapping quality
- MQ0 Number of MAPQ == 0
- NS Number of samples with data
- SB strand bias at this position
- SOMATIC indicates somatic mutation
- VALIDATED describes validated data
Common FORMAT fields:
- AD Read depth for each allele
- ADF Read depth for each allele on the forward strand
- ADR Read depth for each allele on the reverse strand
- DP Read depth
- EC Expected alternate allele counts
- FT Filter
- GL Genotype likelihoods
- GP Conditional genotype quality
- GQ Conditional genotype quality
- GT Genotype
- HQ Haplotype quality
- MQ RMS mapping quality
- PL Phred-scaled genotype likelihoods
- PQ Phasing quality
- PS Phase set