This is a test version of Biostars. For the public version, visit https://www.biostars.org.
extracting snps from multiple vcf file but present in proportion of samples

Hi,

I have vcf files of 25 samples (all of them prepared using freebayes with same reference). I want to extract SNPs which are present in at least 80% of samples (i.e. present in any 20 samples). Kindly help me with it.

I have tried "bcftools isec". It gives me output of those SNPs which are present in at least 20 samples (what I want). But whichever sample was inputted first in file list will be used as a reference. Because of these, only SNPs which are present in my first sample along with any other 19 samples are outputted (This is what I don't want). It should output SNPs present in any 20 samples.

Hope I have explained my problem clearly.

Ankit.

vcf snp variant
• 1,963 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by ankit hinsu • ankit hinsu

1 answer

using vcffilterjdk : http://lindenb.github.io/jvarkit/VcfFilterJdk.html

java -jar jvarkit-git/dist/vcffilterjdk.jar -e 'return variant.getGenotypes().stream().filter(G->!(G.isNoCall() || G.isHomRef())).count()>=20;' input.vcf
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Thanks for reply...

I am guessing I need to merge all vcf file and then use this...

• 0 views • link
written by ankit hinsu • ankit hinsu

Log in to answer this question.

More posts like this