Hello all, I'd like to modify my vcfs to contain several info only. There are "GT:AD:DP:FT:GQ:PL:PP" in format column, but I want to include "GT,DP,GQ" only.
I ran "--extract-FORMAT-info" of vcftools, but the result file is not vcf format. What I want to have finally, is vcf file with "GT,DP,GQ" fields only.
Does anyone know how to handle this? Thanks.
3 answers
I would go for using SNPSift program This is an instance of how I use it to extract: position, ref allele, alt allele, Allele Depth (AD from genotype field), and functional annotation by gene ID.
!/bin/bash
for i in *.vcf;
do java -jar /data/software/snpEff/snpEff/SnpSift.jar\
extractFields "$i"\
POS REF ALT GEN[*].AD ANN[*].GENEID > "filt_${i}";
done
Have you come across SnpSift in your search for tools to do this?
http://snpeff.sourceforge.net/SnpSift.html
Section 10 in the link above (SnpSift documentation) could potentially be useful: 10. SnpSift Extract Fields
Log in to answer this question.
Can you please provide a line having GT:AD:DP:FT:GQ:PL:PP from the vcf file?
Thanks, but I am looking for tools to handle this.
you can use bcftools annotate to keep/remove FORMAT fields
Have you tried AWK aleady?
Actually I extracted these fields with python. But when I tried to analyze this file using rare variant association tool like pseq, vtools, and rvtest, something is not working properly.
I wonder it there are validated tools to extract, not with linux or python.