I have several VCF files generated from WES. Some share a condition, and others could be considered controls. I want to know which SNPs are shared between the files with the condition, and not in the controls. By "shared" I mean for a given SNP, I want to know whether they share the same genotype at that location. Many of the tools that I have tried, like bcftools isec seem to intersect by position, not by the actual genotype. Other tools like bcftools stats seem to only compare 2 files at once. Is there any way to do this for multiple files at once?
•
written
by
ccagg
1 answer
merge your VCF files with bcftools merge and the select the variants with vcffilterjdk http://lindenb.github.io/jvarkit/VcfFilterJdk.html
e.g:
java -jar vcffilterjdk.jar -e 'return variant.getGenotype("sample1").sameGenotype(variant.getGenotype("sample2")) && variant.getGenotype("sample2").sameGenotype(variant.getGenotype("sample3"));'
Log in to answer this question.
-
intersect variant from vcf files (2500 individual)
written by dragon940401 •
-
Finding Shared Patterns Between Patients with Commonality
written by fouerghi20 •
-
add gene names to 'isec' output files of bcftools'
written by minoo •
-
Tool: all2vcf: a tool to convert non-standard variant outputs (mummer. bcftools isec) to VCF
written by Matteo Schiavinato
-
bcftools isec not producing sites.txt
written by steve
-
Fasted method to intersect VCF files
written by nailu90 •
-
Using BCFtools to intersect multiple VCF files
written by Hadi M •
-
bcftools isec tries to create a directory called README.txt inside target directory
written by douwevanderleest •
-
comparing 3 VCF files for concordance and VENN plotting
written by drowl1 •
-
How bcftools isec works ?
written by sacha
