mosdepth icon indicating copy to clipboard operation
mosdepth copied to clipboard

Feature Request: Add threshold for Whole Genome

Open harish0201 opened this issue 4 years ago • 1 comments

Hi!,

Currently the threshold feature can output the chromosome/scaffold wise threshold outcome for every/selected chromosome.

Currently it looks like this:

#chrom  start   end     region  1X      10X     20X     30X
1       0       197608386       unknown 195634794       183535477       51086214        2596469
2       0       149682049       unknown 148581500       139670841       37081610        1436371

I believe that we can add a line for the whole genome as well like this:

#chrom  start   end     region  1X      10X     20X     30X
Total  0       sum(1+2)         unknown value_1X         value_10X       value_20X       value_30X
1       0       197608386       unknown 195634794       183535477       51086214        2596469
2       0       149682049       unknown 148581500       139670841       37081610        1436371

harish0201 avatar Feb 12 '20 06:02 harish0201

I agree this is useful. But this should be simple to calculate externally.

brentp avatar Feb 12 '20 16:02 brentp