Posted by Ron E on Sep 30
In the samtools coverage subcommand, the -w / –n-bins option allows the
user to specify how many “bins” to produce in the coverage histogram. The
code computes: stats[tid].bin_width = (stats[tid].end – stats[tid].beg) /
n_bins; When the number of bins (n_bins) is extremely large relative to the
region length (end – beg), this integer division can yield zero, or lead to
unexpected behavior in subsequent arithmetic. Later in print_hist(),…
– Read More – Full Disclosure
