I'd like to find the highest PL value (posterior likelihood for each possible genotype) for each sample in my VCf file.
From the documentation, I know that you can do
MAX(DP) which will find you the highest genotype depth over all samples
Since the PL field has a value for each possible genotype (e.g. 0/0, 0/1 and 1/1 for a biallelic variant, so 3 PL values), this mechanism is not applicable here.
If for instance
Then I want to find 0.767
Is this possible with bcftools v1.2?
Thanks for the input.