5. bcftools isec 求变异的交集、并集 # 求交集 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:求第一个文件有其他文件...
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
n=2 -w1 A.vcf.gz B.vcf.gz --- Extract records private to A or B comparing by position only --- bcftools isec -p dir -n-1 -c all A.vcf.gz B.vcf.gz --- Print a list of records which are present in A and B but not in C and D --- bcftools isec -n~1100 -c all A...