lh3 / dipcall

Reference-based variant calling pipeline for a pair of phased haplotype assemblies
MIT License
96 stars 10 forks source link

Genotype in the vcf output #15

Open daggoo opened 5 months ago

daggoo commented 5 months ago

Hi, I have two calls at pos1 and pos2. The last columns are the following

pos1: . . GT:AD 1|0:1,1 pos2: . . GT:AD 1|0:1,1

The VAF of the call at pos1 from the read-based variant calling is 0.9; at pos2 -- 0.4.

Could you explain how to read the last column? I saw the header but was confused about why they had the same genotype but different zygosity.