bcftools 为 vcf 文件建索引及合并 vcf 文件

本文转载自查看原文 2018-01-16 09:28 3076

1. bgzip 压缩 vcf 文件为 gz 文件

bgzip -c T8_B8_TN_20171226190133_TN_haplotyper.vcf >T8_B8_TN_20171226190133_TN_haplotyper.vcf.gz
bgzip -c T9_B9_TN_20171226195650_TN_haplotyper.vcf >T9_B9_TN_20171226195650_TN_haplotyper.vcf.gz

2. bcftools 为 gz 文件建索引

bcftools index -t T8_B8_TN_20171226190133_TN_haplotyper.vcf.gz
bcftools index -t T9_B9_TN_20171226195650_TN_haplotyper.vcf.gz

3. bcftools 合并 vcf 文件

bcftools merge -m snps -f PASS,. --force-samples T8_B8_TN_20171226190133_TN_haplotyper.vcf.gz T9_B9_TN_20171226195650_TN_haplotyper.vcf.gz >merge.vcf

免责声明！

本站转载的文章为个人学习借鉴使用，本站对版权不负任何法律责任。如果侵犯了您的隐私权益，请联系本站邮箱yoyou2525@163.com删除。

猜您在找 bcftools合并vcf文件【Bcftools】合并不同sample的vcf文件，通过bcftools 用bcftools将多个vcf文件合并成一个vcf文件或将多个vcf和合并成的vcf文件拆分成单个样本的vcf文件【BCFTOOLS】按样本拆分VCF文件 bcftools 提取vcf（snp/indel）文件子集使用bcftools查看vcf文件所含位点数量信息使用bcftools提取指定样本的vcf文件（extract specified samples in vcf format） VCF文件导入导出 tabix 操作VCF文件 bcftools或vcftools提取指定区段或者多个指定位置的vcf文件（extract specified position ）