vcftools --minGQ not filtering

[AP38]

Hello,

I am trying to filter based on GQ < 15. I do the following:

vcftools --vcf infile.vcf --minGQ 15 --recode --out filtered

However, this filtering does not work, nothing is being removed: "After filtering, kept 1287174 out of a possible 1287174 Site”. I confirm that the GQ tag is present in my VCF file. Other filters such as min/maxDP or minQ work just fine. I am using VCFtools - v0.1.13

Any thoughts on this would be greatly appreciated.

Thanks!

[AP38]

I cross-posted on biostars since I received no answers here:

https://www.biostars.org/p/194918/

[yzzhang]

Previously I tried to use vcftools v0.1.13 to filter my vcf file with --minGQ option, and it worked fine. So it supposed to be some problem with the vcf file than the software. Have you tried to use a different minGO score? which pipeline was used for the vcf generation. I used bowtie2-samtools v0.1.19, and the generated vcf file was good to feed vcftools. hope it helps.

[AP38]

Yes, I tried to use a different minGQ, and I had the same issue. The vcf file was created following this workflow:

- samtools mpileup -C 50 -E -t SP -t DP -u -I -f genome -b bam_list.txt > out.bcf
- bcftools call -v -c -f gq out.bcf > out.vcf

Maybe something went wrong with bcftools call?

[AP38]

Just posted an answer here:

https://www.biostars.org/p/194918/