第一题

# 窗口拆分 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
posted on 2022-09-07 15:54  kid;)  阅读(112)  评论(0)    收藏  举报