mosdepth计算depth

##make pos.bed 

bedtools makewindows -g CS_genome.txt -w 5000000 -s 2500000 > CS_win.bed

CS_genome.txt

CS_win.bed

###calculate the depth of region

 mosdepth -b CS_win.bed H3_output sample1.Aligned.sortedByCoord.out.bam &

###result

 Per-base output will go to sample-output.per-base.bed.gz, the mean for each region will go to sample-output.regions.bed.gz; each of those will be written along with a CSI index that can be used for tabix queries. The distribution of depths will go to sample-output.mosdepth.dist.txt

###**regions.bed.gz