Entering edit mode
6.0 years ago
mostafarafiepour
▴
180
Hi All Dear,
I have 5 VCF files, each of this files is for a population. So, I want to create a Venn Diagram based on SNPs.
I have been able to find this script by searching.
vcf-compare sample1.flt.vcf.gz sample2.flt.vcf.gz > venn.out
So, How do I edit this script for my case?
Best Regard
From each VCF subset,
"chrom, pos, alt, ref"
columns (bcftools), read into R (or Python). Paste those columns together, then use one of the many ready packages for Venn diagrams.There is already a relevant thread in Biostars How To Create A Venn Diagram Of Sample1.Vcf, Sample2.Vcf And Sample1&2(Union).Vcf Files . Also, use upSet plots while comparing more than 3 samples together.