bedtools 软件可以统计每个窗口内SNP变异位点数量,具体命令如下:
bedtools coverage -a window.bed -b raw.filtered.snp.vcf.gz -counts >result.txt
其中bed文件格式如下:
结果文件如下:
最后一列为每个窗口内SNP数量。
348 篇文章
如果觉得我的文章对您有用,请随意打赏。你的支持将鼓励我继续创作!