Snakemake BaseRecalibrator wrapper配置java_opts运行失败排查

你当前的配置存在3个直接导致运行失败、耗时过长的问题：

• 参数传递完全失效：0.74.0版本的GATK BaseRecalibrator wrapper约定Java启动参数必须写在params.java字段下，你写的params.java_opts不会被wrapper读取，配置的4G内存、GC线程参数根本没有传给Java进程，等于没有配置。
• 内存配置逻辑错误：你向集群调度申请了32G内存，但仅给Java堆分配4G，不仅28G申请到的内存完全闲置，4G堆内存也完全不足以支撑全基因组级别BAM的碱基质量重校正计算，进程会因为堆内存溢出被操作系统强制杀死。这类强制终止发生时进程来不及将错误信息刷入磁盘，就会出现你遇到的规则日志为空的情况，可以直接查集群调度的作业审计日志，如果退出码是137，就能确认是OOM强制终止。
• 多线程配置未生效：该版本wrapper不会自动将规则的threads: 8参数传递给GATK的原生计算线程参数，BaseRecalibrator默认单线程运行核心的PairHMM计算步骤，这是规则耗时过长的核心原因。你配置的-XX:ParallelGCThreads=8仅控制垃圾回收的线程数，和业务计算线程无关，配到8个反而会挤占计算资源。

rule recalibrate_base_qualities:
    (unchanged)
    params:
        (unchanged)
        # Java启动参数放在params.java下，堆内存配28G，留4G给堆外内存、元空间等系统开销，不超过32G的总申请值
        java="-Xmx28G -XX:ParallelGCThreads=2",
        # 显式传递GATK原生计算线程数，和申请的8线程匹配
        extra="--native-pair-hmm-threads 8"
    threads: 8
    resources:
        mem_mb=32768,
        time="12:00:00",
    wrapper:
        "0.74.0/bio/gatk/baserecalibrator"

补充说明：如果是外显子、靶向测序这类小体量数据，可以按比例降低内存申请值，全基因组样本BQSR步骤Java堆内存建议不低于24G。

内容的提问来源于stack exchange，提问作者gernophil