vcf文件提取SNP位点

测序得到VCF文件后,有时需要提取部分SNP位点,这里介绍的工具是VCFtools,安装以及基础功能见下贴:
vcftools安装及基础用法 - 简书 (jianshu.com)

1.准备SNP位点list

首先将需要提取的SNP位点的ID整理在一个文本文件中,一个ID一行。这里要注意ID和vcf文件中的ID要一致,Chr_1;chr_1还是只有数字编号1:

$ cd /your/path/vcftools
$ vim snp.id
Chr1__272173254
Chr10__132755179
Chr11__187947859
Chr12__233234905
Chr14__135335406

2.拆分SNP

 $ vcftools --vcf root.id.vcf --snps snp.id --recode --out sig.snp
  • --vcf:VCF文件全称
  • --snps:SNP的id列表
  • --out:输出vcf文件前缀

引用转载请注明出处,如有错误敬请指出。

©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容