基因组一过 10 Gb，常规变异检测的默认设置就开始崩，麻烦的地方不在比对慢，而在几个 32 位的坐标上限。下面这批数据约 23 Gb、十二条染色体，每条接近 2 Gb，比人类大七倍。会崩的地方有这么几处： bwa 的索引有 4 Gb 上限，整个基因组建不了索引。 .tbi 索引的坐标是 32 位，单条染色体超过 2,147,483,647 bp 就无法建索引，会直接报坐标溢出。 BAM 和...