Can you support Illumina DRAGEN msVCF?

I see what you would like to achieve now. You would have to rerun the iterative gVCF genotyper with the parameter --gg-remove-nonref false, then you get three values for LPL (Local normalized, Phred-scaled likelihoods for genotypes as in original gVCF).

chr1    10109   .       AACCCT  A,<NON_REF>     .       PASS    AC=1,0;AN=6;NS=3;NS_GT=3;NS_NOGT=0;NS_NODATA=0;GAC=1,0;GAN=6;GNS=3;GNS_GT=3;GNS_NOGT=0;GNS_NODATA=0     GT:GQ:AD:FT:LPL:LAA     0/1:8:2,1,0:PASS:25,0,5,29,9,37:1,2     0/0:5:101,0,13:PASS:0,5,189:2   0/0:5:84,0,10:PASS:0,5,101:2
1 Like