美文网首页生物信息学
bcftools根据样品拆分vcf文件

bcftools根据样品拆分vcf文件

作者: 17号小行星 | 来源:发表于2022-01-24 18:16 被阅读0次
  • 文件1:需要提出的样品的id,文件名DB_sample.list
SYH10
SYH11
SYH12

文件2:含有不同样本的变异位点的vcf文件


image.png
  • 利用bcftools提取DB_sample.list中的样品的vcf文件
bcftools view  -S DB_sample.list samples.gatk.con.snp.12.vcf -O v -o samples.gatk.con.DB.snp.vcf

结果


image.png
  • 计算Pi值
vcftools --vcf samples.gatk.con.DB.snp.vcf  --window-pi 1000  --out samples.gatk.con.DB.snp.vcf

相关文章

网友评论

    本文标题:bcftools根据样品拆分vcf文件

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