1、 bedtools makewindows -g genome.txt -w 10000000 -s 1000000 > windows.bed
samtools depth xx.sort.rm.bam > depth (1)
利用脚本统计滑动窗口的平均depth
利用R作图.
12、bedtools coverage -a gene.bed -b bam -d -counts -sorted > coverage 计算区间的平均测序深度
samtools bedcov gene.bed sample.bam 计算bed区间所有位点测序深度的加和 & 手动算 得到bed文件中的平均测序深度!
与 samtools depth xx.sort.rm.bam > depth (1) 利用脚本统计滑动窗口的平均depth 方法所得结果相同!!!
网友评论