This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Vcf Filtering

Hi, I recently obtained a vcf file with 90 samples from different generations. I realized that the file is extremely large. This is because I have a lot of possible sequencing errors, observing AF really closed to 0 or 1, and only in one of the samples out of 90. For this reason, I want to filter the whole file, selecting only the SNPs which at least in one sample has a value between 0.15 and 0.85. Do you know how can I do it?

I tried with this code "bcftools view -i 'sum(AF >= 0.15 && AF <= 0.85) > 0'" but it didn't gave me the expected result. Thank you!

vcf filtering af
• 1,479 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by Marc Ibanez • Marc Ibanez

AF is usually an INFO field, not a FORMAT/genotype field. What is the definition of AF in your header file and please, give us an example of such genotypes.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Hi, the definition of AF in my header is this one:

##FORMAT=<ID=AF,Number=A,Type=Float,Description="Observed ALT frequencies"> . And here I add an example of the genotype: GT:PL:DP:AD:SCR:QS:SAD:AF:XF    ./.:.:48:19,21:12:540,524:40:0.525:0.511628     ./.:.:42:24,13:9:537,305:37:0.351351:0.35.

Thank you in advance

• 0 views • link
updated by GenoMax GenoMax • written by Marc Ibanez • Marc Ibanez

I don't understand what AF refers to in this scenario when the genotype is ./.

Are these diploid organisms?

• 0 views • link
written by Jeremy Leipzig Jeremy Leipzig

1 answer

using vcffilterjdk https://jvarkit.readthedocs.io/en/latest/VcfFilterJdk/

java -jar dist/jvarkit.jar vcffilterjdk -e 'return variant.getGenotypes().stream().anyMatch(G->{double v=G.getAttributeAsDouble("AF",-1);return v>=0.15 && v<=0.85;});'  in.vcf
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this