# 窗口拆分 test.bed 需要修改
bedtools makewindows -g test.bed -w 100000 -s 100000 >> 100k_window.bed
# tumor coverage count
bedtools coverage -a 100k_window.bed -b tumor.bam > tumor.counts.txt
# normal coverage count
bedtools coverage -a 100k_window.bed -b normal.bam > normal.counts.txt
原创文章，作者：ItWorker，如若转载，请注明出处：https://blog.ytso.com/288142.html