Entering edit mode
5 weeks ago
Марта
•
0
Hello, all. I'm trying to add FORMAT/TAG annotation to all samples in my vcf file and I want it to be a missing value (".") or 0 for all of them. Is there a way to do this with bcftools?
Thank you for the answer! This works, although I noticed that if add a second TAG to the edited vcf with the same method it replaces the first one I added. Could you help me figure out how to add more than one TAG?
run the script twice by changing the name "TAG" to "TAG2" ...
I did, I ran
and then
In the resulting out2.vcf, both ##FORMAT=<ID=TAG1..> and ##FORMAT=<ID=TAG2..> are present, but all of the samples only have TAG2 in their FORMAT fields. I'm not familiar with Java so I'm not sure what can be done to correct this. Could you help please?
works on my machine: