The Variant Call Format (VCF) is a flexible file format specification that allows us to represent many different variant types ranging from SNPs, indels to copy number variations. However, variant representation in VCF is non-unique for variants that have explicitly expressed reference and alternate sequences.
On this wiki page, we describe a a variant classification system for VCF
= Definitions =
The definition of a variant is based on the definition of each allele with respect to the reference sequence. We consider 5 major types loosely
defined as follows.
: The reference and alternate sequences are of length 1 and the base nucleotide is different from one another.
a.The reference and alternate sequences are of the same length and have to be greater than 1 and all nucleotides in the sequences differ from one another.
b. all reference and alternate sequences have the same length.
a. The reference and alternate sequence are not the same length . : AND : b. The removal of a subsequence of the longer sequence would reduce the longer sequence to the smaller sequence.
a. A clumping of nearby SNPs, MNPs or Indels.
: The alternate sequence is represented by
a angled bracket tag .
= Examples =
We present the following examples to explain the
concepts explained earlier.
== Simple Biallelic Examples ==
#SNP, 1 ts
#2 ts ALT GC
REF AT #1 del
== Complex Biallelic Examples ==
REF AT #1 ts
REF ATT #1 ts, 1 tv, 1 del
#2 ts ALT GTTTC #since all the alleles are of the sample length, classified as MNP too.
REF ATTTTTTTT #2 ts, 1 del
== Simple Multiallelic Examples ==
#1 ts ALT C #1 tv
#1 ts, 1 tv ALT CT #2 tv
#1 del ALT ATTTT #1 ins
== Complex Multiallelic Examples ==
#SNP, 1 ts ALT AC #SNP, 1 ts #since all the alleles are of the sample length, classified as MNP too.
#CLUMPED, 1 ts, 1 tv ALT ATTTC #SNP, 1 tv #since all the alleles are of the sample length, classified as MNP too.
#SNP, 1 tv ALT AG #MNP, 2 tv ALT GTT #INDEL, 1 ins
#MNP |INDEL, 2 tv, 1 del ALT AG #MNP |INDEL, 1 ts, 1 tv ALT GTGTG #SNP |INDEL |CLUMPED, 1 tv
Weird Examples ==
Structured Variants Examples ==
no. of structural variants : 41217 2 alleles : 38079 deletion : 13135 <DEL> insertion : 16451 <INS> mobile element : 16253 <INS:ME> ALU : 12513 <INS:ME:ALU> LINE1 : 2911 <INS:ME:LINE1> SVA : 829 <INS:ME:SVA> numt : 198 <INS:MT> duplication : 664 <DUP> inversion : 100 <INV> copy number variation : 7729 <CN4> >=3 alleles : 3138 copy number variation : 3138 <CN4>,<CN8>
= Output =
3 alleles : 273 (0.89) [537/601]
4 alleles : 3 (1.00) [9/9] <br>
no. of Indel : 6600770
2 alleles : 6285861 (0.88) [2937096/3348765] #ins/del ratio and the respective counts
3 alleles : 280892 (8.72) [503977/57807]
4 alleles : 34 (1.16) [109/94]
>=5 alleles : 4 (0.76) [13/17] <br>
no. of Complex Substitutions : 159298 #equivalent to categories not including
simple SNPs, Block Substitutions and Simple Indels
2 alleles : 81508 (0.61) [60312/98113] (0.66) [32479/49029]
3 alleles : 71003 (0.69) [35811/51840] (0.34) [34268/100942]