Samtools v1.22.1 における過剰なヒストグラム... ノート

Samtools v1.22.1 における過剰なヒストグラムビン数の不適切な処理がスタックオーバーフローを引き起こす

Ron Eさんが9月30日に投稿した内容です。samtools coverageサブコマンドの-w / --n-binsオプションを使用すると、カバレッジヒストグラムで生成する「ビン」の数を指定できます。コードは次のように計算します。stats[tid].bin_width = (stats[tid].end - stats[tid].beg) / n_bins;ビンの数(n_bins)が領域長(end - beg)に対して非常に大きい場合、この整数除算はゼロを生成したり、後続の算術演算で予期しない動作を引き起こしたりする可能性があります。後でprint_hist()内で…