genome-wide coverage metric for WGS data?
Hi guys

Thanks a lot for developing mosdepth - I have a potentially silly question. I've started analysing some pacbio WGS data and was hoping mosdepth would be helpful for producing a genome wide average coverage figure. Essentially I'm looking for the best approximation to the 'we got 15X coverage for this sample' in order to decide with my biology-based colleagues if we need to sequence more or not.

Is this something I can do with mosdepth? The closest to this in my understanding would be to look at the cumulative plot produced by your python script and then set a threshold somewhere. Is there a better way? Thanks!