Discrepancies calculatin average coverage with Bedtools and Mosdepth

Question

[María José]

Hi,

I would like to calculate the average coverage of specific genomic regions. However I get different results with mosdepth and bedtools:

bedtools

bedtools coverage -a u1.narrowPeak -b u1_bwa.bam -mean > MeanCoverageu1.bedgraph

chrY 57216300 57217350 u1_1 4 4105 .       0 0 0 825 18933.2226562

mosdepth

mosdepth -t 8 -n --by u1.narrowPeak u1_mosdepth_output u1_bwa.bam

chrY 57216300 57217350  u1_1 1 5015.89

[Pierre Lindenbaum]

you asked a question a few hours ago, but we got no validation or feedback Error with bigWigAverageOverBed