-
Notifications
You must be signed in to change notification settings - Fork 260
Description
If the tag with the same name exists both in the INFO and the FORMAT --rename-annots as well as -c INFO/OLD:=INFO/NEW renames both tags
Example:
chr1 10250 rs199706086 A C 231.64 . AC=1;AF=0.5;AN=2;ASSEMBLED_HAPS=13;BaseQRankSum=0.328;DB;DP=95;ExcessHet=0;FILTERED_HAPS=9;FS=0;HAPCOMP=0;HAPDOM=0.647;HEC=45;MLEAC=1;MLEAF=0.5;MQ=13.23;MQRankSum=-1.858;QD=8.58;ReadPosRankSum=0.642;SOR=0.582;VARIANT_TYPE=snp;X_CSS=non-skip;X_GCC=0.4;X_HIL=0;X_HIN=.;X_IC=NA;X_IL=.;X_LM=CCTAA;X_RM=CCCTA;LCR=TRUE;MAP_UNIQUE=TRUE;UG_HCR=TRUE;GNOMAD_AF=0.011 GT:AD:DP:GQ:MQ0C:PL:SB:SCL:SCR 0/1:11,16:27:99:20,34:239,0,250:3,8,5,11:51,47:24,39
after running
bcftools annotate -c INFO/DP_RAW:=INFO/DP -Oz -o new.vcf.gz old.vcf.gz
becomes
chr1 10250 rs199706086 A C 231.64 . AC=1;AF=0.5;AN=2;ASSEMBLED_HAPS=13;BaseQRankSum=0.328;DB;DP_RAW=95;ExcessHet=0;FILTERED_HAPS=9;FS=0;HAPCOMP=0;HAPDOM=0.647;HEC=45;MLEAC=1;MLEAF=0.5;MQ=13.23;MQRankSum=-1.858;QD=8.58;ReadPosRankSum=0.642;SOR=0.582;VARIANT_TYPE=snp;X_CSS=non-skip;X_GCC=0.4;X_HIL=0;X_HIN=.;X_IC=NA;X_IL=.;X_LM=CCTAA;X_RM=CCCTA;LCR=TRUE;MAP_UNIQUE=TRUE;UG_HCR=TRUE;GNOMAD_AF=0.011 GT:AD:DP_RAW:GQ:MQ0C:PL:SB:SCL:SCR 0/1:11,16:27:99:20,34:239,0,250:3,8,5,11:51,47:24,39