摘要: 1、bcftools提取指定区段的vcf文件 下载安装bcftools 见如下命令: 注意:输入的vcf以gz格式存在,不然会报错:Failed to open 1000Genomes.vcf: not compressed with bgzip 如何将vcf生成gz格式,见这篇文章bcftools 阅读全文
posted @ 2018-06-22 14:47 橙子牛奶糖 阅读(17190) 评论(0) 推荐(0) 编辑
摘要: 对于某些特殊vcf,想去除没有变异的基因型(主要形式为0|0或者0/0),则需要用到GATK的--excludeNonVariants参数,命令行如下: 参考链接: https://www.biostars.org/p/203809/ https://software.broadinstitute. 阅读全文
posted @ 2018-06-22 14:37 橙子牛奶糖 阅读(1374) 评论(0) 推荐(0) 编辑
摘要: 1、下载安装bcftools。 2、准备样本ID文件,这里命名为samplelistname.txt,一个样本一行,如下所示: sample1 sample2 sample3 3、输入命令: 参考链接: https://www.biostars.org/p/184950/ https://samto 阅读全文
posted @ 2018-06-22 10:21 橙子牛奶糖 阅读(9259) 评论(2) 推荐(0) 编辑