修复DNA_Sequence空列表赋值错误及Python序列翻译脚本问题

这是此前帖子《可通过1152种不同密码子表将1000条DNA序列翻译为蛋白质的Python脚本无法运行》的续篇，目前我在PyCharm中编辑的脚本内容如下：

import yaml
file = open('/home/ihor/Документи/біоінформатика/DNA_Sequence72 - Копія.py', 'r')
DNA = file.read()
DNA_Sequences = DNA.split(',')
Genetic_Codes = open('/home/ihor/Документи/біоінформатика/Genetic_Codes.py', 'r')
Genetic_Codes = Genet...

现在要解决DNA_Sequences被赋值为空列表的问题，咱们一步步排查修复：

先搞懂空列表出现的原因 #

大概率是这两种情况导致的：

• 读取的DNA文件内容为空，或者内容里没有逗号分隔符，使得split(',')返回空/全是空字符串的列表
• 文件路径写错、没有读取权限，导致DNA变量拿到的是空字符串

具体修复方案 #

1. 给文件读取加安全校验和规范处理

手动打开文件容易遗漏关闭操作，换成with语句自动管理文件生命周期，同时加入内容校验逻辑：

import yaml

try:
    # 读取DNA序列文件，用with自动处理关闭
    with open('/home/ihor/Документи/біоінформатика/DNA_Sequence72 - Копія.py', 'r') as dna_file:
        DNA_content = dna_file.read().strip()  # 去掉首尾空白，避免空行干扰
    
    # 检查文件内容是否为空
    if not DNA_content:
        raise ValueError("DNA文件内容是空的，请检查文件是否正确！")
    
    # 分割序列，同时过滤掉空字符串（比如文件末尾多了个逗号的情况）
    DNA_Sequences = [seq.strip() for seq in DNA_content.split(',') if seq.strip()]
    
    # 再次确认是否得到有效序列
    if not DNA_Sequences:
        raise ValueError("没提取到有效DNA序列！请确认文件里的序列是用逗号分隔的，且内容不为空。")

    # 读取密码子表文件（同样用with更安全）
    with open('/home/ihor/Документи/біоінформатика/Genetic_Codes.py', 'r') as gc_file:
        Genetic_Codes = gc_file.read()
        # 这里补充你后续解析密码子表、翻译序列的逻辑

# 捕获文件找不到的错误
except FileNotFoundError as e:
    print(f"出错啦：找不到指定文件 - {e.filename}")
# 捕获其他异常
except Exception as e:
    print(f"发生错误：{str(e)}")

2. 检查文件内容格式

打开你的DNA_Sequence72 - Копія.py文件，确认格式是否符合预期：

• 确认序列之间确实用逗号分隔（比如ATCG123,GGCC456,TTAA789）
• 如果序列是按换行分隔的，就把split(',')改成split('\n')
• 删掉文件末尾多余的逗号或空行，避免分割出无效的空字符串

3. 验证文件路径和权限

因为你的路径包含非英文字符，建议直接在PyCharm里拖拽文件到代码中，自动生成正确路径，避免手动输入出错。同时要确认当前用户有读取该文件的权限。

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