Snakemake Wrappers在screen外正常运行但screen内报错求助

我刚帮同事排查过几乎一模一样的问题，你的核心症结就是Screen启动时没有正确继承conda环境的环境变量，导致流程调用了linuxbrew的Python而非你激活的conda Python，进而引发了依赖不匹配和pickle协议不兼容的错误。

问题拆解 #

• 为什么部分规则正常（比如fastqc、bwa/aln）？这类规则要么是调用独立的二进制文件（不依赖Python），要么其wrapper脚本对Python版本兼容性要求低；而cutadapt、bwa/samse的wrapper依赖特定版本的Python环境，一旦调用了错误的Python解释器，就会触发CalledProcessError或者pickle协议错误（比如你遇到的ValueError: unsupported pickle protocol: 3，大概率是因为linuxbrew的Python是2.x版本，无法解析Python3生成的pickle文件）。
• 为什么Screen内外PATH不一样？Screen默认启动时会读取系统的默认shell配置（比如/etc/profile、~/.bashrc），而不会自动继承当前终端已经激活的conda环境变量——除非你明确配置它这么做。

具体解决方案 #

1. 临时解决：启动Screen前强制传递conda环境变量

在激活conda环境后，用以下命令启动Screen，让它完全继承当前终端的环境：

screen -L -d -m -S snakemake_session bash -c "source ~/.bashrc; conda activate your_env_name; exec bash"

然后进入这个screen会话：

screen -r snakemake_session

之后再运行Snakemake，就能确保用的是conda环境里的Python了。

2. 永久解决：修改Screen配置，默认继承环境变量

编辑你的~/.screenrc文件（如果没有就新建一个），添加以下内容：

defshell -bash

这个配置会让Screen启动时使用登录shell，从而读取你当前用户的shell配置，包括conda的环境变量。保存后，下次启动Screen就会自动继承当前激活的conda环境了。

3. 流程层面加固：在Snakemake中指定Python路径

为了彻底避免环境变量的问题，你可以在Snakemake的规则里直接指定conda环境的Python解释器路径。比如在cutadapt的规则里：

rule cutadapt:
    input:
        # 你的输入文件
    output:
        # 你的输出文件
    params:
        # 你的参数
    wrapper:
        "0.74.0/bio/cutadapt"
    shell:
        "{conda_prefix}/bin/python {wrapper_script}"

或者更简单的方式，在Snakemake命令前显式指定Python路径：

{conda_prefix}/bin/python -m snakemake --cores 10

这里的{conda_prefix}可以通过echo $CONDA_PREFIX在激活环境后获取具体路径。

4. 验证环境是否正确

进入Screen后，先执行以下命令确认Python路径：

which python
python --version

确保输出的是conda环境下的Python路径和版本，而不是linuxbrew的。

额外提示 #

如果你的HPC集群有模块管理系统（比如module），记得在启动Screen前也要加载好必要的模块（比如bwa、samtools等），避免出现二进制文件找不到的问题。

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