vcftools过滤指定缺失率的变异位点
vcftools过滤指定缺失率的变异位点
vcf 文件中很多snp在某些样品中是缺失的，也就是基因型为 "./." 。如果缺失率较高，这种snp位点在很多分析中是不能用的，需要去掉。这时候就可以使用vcftools进行过滤。用到的选项为--max-missing 。
具体用法如下：
下图为原始vcf文件。
运行以下命令：
vcftools --vcf snp.vcf --recode --recode-INFO-all --stdout --max-missing 1 > snp.new.vcf
--max-missing 后跟的值为 0-1 ，1代表不允许缺失，0代表允许全部缺失
最后，得到的vcf文件中snp位点都是没有缺失的。
此外，我们在网易云课堂上有各种教学视频，有兴趣可以了解一下：
1. 文章越来越难发？是你没发现新思路，基因家族分析发2-4分文章简单快速，学习链接：基因家族分析实操课程
2. 转录组数据理解不深入？图表看不懂？点击链接学习深入解读数据结果文件，学习链接：转录组（有参）结果解读；转录组（无参）结果解读
3. 转录组数据深入挖掘技能-WGCNA，提升你的文章档次，学习链接：WGCNA-加权基因共表达网络分析
4. 转录组数据怎么挖掘？学习链接：转录组标准分析后的数据挖掘
6. 更多学习内容：linux、perl、R语言画图，更多免费课程请点击以下链接：
- 发表于 2019-01-04 16:05
- 阅读 ( 9948 )
- 分类：软件工具
