Does anyone have a link to a python script for calculating population private SNPs? Many articles indicate that ad-hoc scripts were used to obtain private SNPs but the code or the algorithm summary is always unavailable. I know about
bcftools isec and
vcftools isec. I'd like to have a look at the Python implementations. I have vcf_files for every population that I'm interested in.