1.制作保留区间bed文件
![](https://img.haomeiwen.com/i15061072/9e9f434763d0f2e7.png)
2.利用bedtools根据bed文件提取落在区间内的SNPs。
bedtools intersect -a myfile.vcf.gz -b mutil-region.bed -header > output.vcf
其他资料也提供了一些其他提取方法(没仔细研究):
https://www.cnblogs.com/miyuanbiotech/p/14530891.html
1.制作保留区间bed文件
2.利用bedtools根据bed文件提取落在区间内的SNPs。
bedtools intersect -a myfile.vcf.gz -b mutil-region.bed -header > output.vcf
其他资料也提供了一些其他提取方法(没仔细研究):
https://www.cnblogs.com/miyuanbiotech/p/14530891.html
本文标题:2023-07-05 根据多个指定位置提取vcf子集(实现过滤点
本文链接:https://www.haomeiwen.com/subject/qdkhudtx.html
网友评论