bcftools isec -p dir A.vcf.gz B.vcf.gz # 并集 bcftools isec -p dir -n=2 -w1 A.vcf.gz B.vcf.gz # 差集 bcftools isec -p dir -n-1 -c all A.vcf.gz B.vcf.gz dir:输出文件夹 -c:选择变异类型进行处理 -C:求第一个文件有其他文件没有的集合 -f:可筛选FILTER列的信息,如果想筛选P...
Pysam is a Python package for reading, manipulating, and writing genomics data such as SAM/BAM/CRAM and VCF/BCF files. It's a lightweight wrapper of the HTSlib API, the same one that powers samtools, bcftools, and tabix. - pysam/bcftools/vcfisec.c at 6b4
For example, when performing line intersections, the desire may be to consider as identical all sites with matching positions (*bcftools isec -c* 'all'), or only sites with matching variant type (*bcftools isec -c* 'snps'{nbsp} *-c* 'indels'), or only sites with all alleles identical...