Multi-sample variant calling files made with GATK4 have both genotype, genotype phase and genotype phase block information.
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
##FORMAT=<ID=PGT,Number=1,Type=String,Description="Physical phasing haplotype information, describing how the alternate alleles are phased
##FORMAT=<ID=PID,Number=1,Type=String,Description="Physical phasing ID information, where each unique ID within a given sample (but not across samples) connects records within a phasing group">
If I run the following bcftools query I can get the 3 above fields for each sample.
bcftools query -f '%CHROM\t%POS\t%REF\t%ALT[\t%GT\t%PGT\t%PID]\n' my_samples.vcf.gz
Combinations of GT, PGT and PID make sense to make. Heterozygous genotypes (0/1) can be both on the paternal (e.g. 0|1) or maternal haplotype (e.g. 1|0).
0/1 0|1 748_A_C
0/1 1|0 627_C_T
But how can a homozygous genotype (0/0) be phased on one of the parental haplotypes (e.g. 0|1)?
There is no alternative allele (=1) in the homozygous genotype (0/0)?
0/0 0|1 627_C_T
The above 0/0 genotype combined with 0|1 genotype phase does not make sense to me. Does anyone know what this means? How to interpret this. To me it just seems nonsens, the 0|1 627_C_T information can be discarded?
0 answers
No answers yet.
Log in to answer this question.
Not sure how this is possible. Perhaps you should bring it up on the GATK forums or search to see if it has already been brought up? https://gatkforums.broadinstitute.org/gatk/categories/ask-the-team