2022-06-08修改vcf文件样本名称,提取指定vcf样本
2022-06-07 本文已影响0人
汪大山
1.修改vcf样本名
使用bcftools的 rehead指令
bcftools reheader -s sample.file old.vcf -o new.sample.vcf
-s 后面的为sample.file 指定要替换的样本名,内容如下
2.提取指定样本vcf文件
准备好样本ID文件,samplelist.txt,一个样本名一行,如下
bcftools view -S samplelist.txt old.vcf.gz -Ov > new.vcf