草庐IT

bcftools

全部标签

取两个VCF的差集

参考:https://www.biostars.org/p/113509/#113509https://samtools.github.io/bcftools/bcftools.html#isecbcftoolsisec[option]A.vcf.gzB.vcf.gz[…]创建A和B的交集和补集,将输出保存在dir/*bcftoolsisec-pdirA.vcf.gzB.vcf.gz过滤存在于A(需要INFO/MAF>=0.01)和B(需要INFO/dbSNP)但不在C的位点,并创建交集,仅包括过滤后出现在至少两个文件中的位点bcftoolsisec-e'MAF精确匹配等位基因,从A中提取A和