I have 3 VCF files.
A1.vcf A2.vcf A3.vcf
I want to get the common SNPs that are present in all these three files. And output must be in vcf format.
Bioinformatics Stack Exchange is a question and answer site for researchers, developers, students, teachers, and end users interested in bioinformatics. It only takes a minute to sign up.Sign up to join this community
for files in A1.vcf A2.vcf A3.vcf do bgzip "$files" tabix "$files" done bcftools isec -n=3 -c all -o all_common_SNPs.vcf A1.vcf.gz A2.vcf.gz A3.vcf.gz ```