如何通过Python3执行复杂Shell命令？语法错误求助

首先，咱们先拆解你原代码里的几个问题：

• 你写的system.os(shell)是错误的，Python里调用系统命令的正确方式是os.system()，这是直接触发SyntaxError的核心原因。
• 另外，你的Shell命令字符串用双引号包裹，里面又嵌套了双引号，虽然Shell本身支持这种写法，但Python字符串里会引发语法冲突，要么转义内部双引号，要么换用单引号包裹整个命令。

下面给你几种可行的解决方案：

方案1：修正os.system的写法 #

如果还是想用os.system，可以把Shell命令用三个单引号包裹（避免转义麻烦），同时修正调用方式：

import os

# 用三单引号包裹Shell命令，规避双引号嵌套的语法问题
shell_cmd = '''ls *R1.fastq.gz | while read a; do b=${a%R1.fastq.gz}R2.fastq.gz; c=${a%R1.fastq.gz}R1.out.fq; d=${a%R1.fastq.gz}R2.out.fq; e=${a%R1.fastq.gz}.s.fq; echo "sickle pe -t sanger -f $a -r $b -o $c -p $d -s $e"; done > SICKLE.sh'''
os.system(shell_cmd)

方案2：使用更推荐的subprocess模块 #

subprocess是Python官方推荐的执行外部命令的模块，比os.system更灵活、安全，还能更好地处理输入输出。这里用subprocess.run的现代写法：

import subprocess

shell_cmd = '''ls *R1.fastq.gz | while read a; do b=${a%R1.fastq.gz}R2.fastq.gz; c=${a%R1.fastq.gz}R1.out.fq; d=${a%R1.fastq.gz}R2.out.fq; e=${a%R1.fastq.gz}.s.fq; echo "sickle pe -t sanger -f $a -r $b -o $c -p $d -s $e"; done > SICKLE.sh'''

# shell=True表示执行Shell命令，check=True会在命令执行失败时抛出异常，方便排查问题
subprocess.run(shell_cmd, shell=True, check=True)

方案3：更Pythonic的写法（推荐） #

如果你的文件名里可能包含空格或特殊字符，用Shell的ls遍历会有潜在问题。不如直接用Python自己的文件遍历功能，更安全可靠：

import glob

# 用glob匹配所有R1结尾的fastq.gz文件
for r1_file in glob.glob('*R1.fastq.gz'):
    # 提取文件名前缀，构造对应的R2和输出文件
    base_name = r1_file[:-len('R1.fastq.gz')]
    r2_file = f'{base_name}R2.fastq.gz'
    r1_out = f'{base_name}R1.out.fq'
    r2_out = f'{base_name}R2.out.fq'
    s_out = f'{base_name}.s.fq'
    
    # 构造sickle命令
    sickle_cmd = f'sickle pe -t sanger -f {r1_file} -r {r2_file} -o {r1_out} -p {r2_out} -s {s_out}'
    
    # 追加写入到SICKLE.sh
    with open('SICKLE.sh', 'a') as f:
        f.write(f'{sickle_cmd}\n')

这种方式完全脱离了Shell的依赖，能更好地处理各种特殊文件名，也更容易调试和维护。

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