Convert missing genotype ./. to homozygous 0/0 in vcf

2

Dear all,

Is there a way to change missing genotype ./. to homozygous 0/0 in a vcf file? Fo example:

Input:

chr1    69849   .   G   A   100.59  PASS    AC=1;AF=0.500;AN=2;ClippingRankSum=0;DP=16;VQSLOD=1.66;culprit=NULL;set=variant-variant2    GT:AD:DP:GQ:PL  ./.:13,3:16:41:41,0,402

Output:

chr1    69849   .   G   A   100.59  PASS    AC=1;AF=0.500;AN=2;ClippingRankSum=0;DP=16;VQSLOD=1.66;culprit=NULL;set=variant-variant2    GT:AD:DP:GQ:PL  0/0:13,3:16:41:41,0,402


vcf

• 156 views



Source link