美文网首页GWAS生物信息学
使用bcftools 将vcf文件拆分成单个vcf

使用bcftools 将vcf文件拆分成单个vcf

作者: 花生学生信 | 来源:发表于2021-07-12 10:45 被阅读0次

使用bcftools query -l test.vcf > 453.list 将vcf文件中的样品名提取出来

使用sh脚本
# 将453份数据写入变量$sample 

for sample in $(cat 453.list);

#创建每个样本的list文件

do touch list/$sample.list;

#将样本写入每个文件

#echo $sample >> list/$sample.list;

#使用bcftools提取

do bcftools view -S list/$sample.list /public/home/lianglunping/work/SV/453.vcf>/public/home/lianglunping/work/SV/453/$sample.vcf;

done

相关文章

网友评论

    本文标题:使用bcftools 将vcf文件拆分成单个vcf

    本文链接:https://www.haomeiwen.com/subject/wrnkpltx.html