Posted by Ron E on Sep 30In 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(),...