How to ignore "GT" from vcf file
I used medaka_haploid_variant caller to identify variants in my viral reads. After running medaka tools annotate, I got an annotated vcf file. Then filtering by bcftools,
bcftools complaining about GT not defined in the annotated vcf
bcftools filter -Ob -e 'DP<1000' -o filtered.vcf annotated.vcf
I got the following error
FORMAT 'GT' at NC_xxxxx:33 is not defined in the header, assuming Type=String
From the github post I made https://github.com/nanoporetech/medaka/issues/257, I came to that GT does not appear in the VCF header but does occur in the records which is a bug in medaka.
Now, how to run bcftools by ignoring the GT? Is it possible?
• 2,012 views
•
link
updated
by
Ram
•
written
by
SUDOsundu •
•
written
by
SUDOsundu •
1 answer
fix your VCF header.
awk '/^#CHROM/ {printf("##FORMAT=<ID=GT,Number=1,Type=String,Description=\"Genotype\">\n");} {print}' input.vcf
• 0 views
•
link
written
by
Pierre Lindenbaum
Log in to answer this question.
More posts like this
-
update FMT/GT in VCF file using bcftools annotate
written by karen2 •
-
Adding custom annotation to VEP output VCF
written by avelarbio46 •
-
add INFO tag to vcf with bcftools
written by from the mountains
-
bcftools view error: the tag "mis" is not defined in the VCF header
written by skylerz •
-
Medaka tools annotation missing QUAL values
written by SUDOsundu •
-
Confused about the AF tag in VCF from samtools/bcftools mpileup
written by Makplus T
-
Fix or remove format tag in .vcf file
written by j.lunger18 •
-
bcftools view error: "F_PASS" is not defined in the VCF header
written by Luther_Blisset •
-
Format in vcf files in Geuvadis data not recognized
written by jamespower •
-
does bcftools has anything like vcf-check?
written by roalva1
easiest would be to add one line to the header. You can do that with any text editor. The other option is to remove the GT fields from the format column, which can be a bit more complicated, depending on your prior experience.