Python序列计数报错NameError：n_nt变量未定义问题排查

你的代码触发NameError的核心原因是：只有当序列中存在'N'时，n_nt变量才会被定义，如果序列里没有'N'，这个变量从未被赋值，打印时自然会报错。另外，原代码的循环逻辑完全冗余——每次遍历核苷酸都重复调用seq.count()，做了大量无意义的重复计算。

先初始化所有计数变量为0，确保即使某类核苷酸数量为0，变量也存在；直接用count()方法一次性获取各类核苷酸的数量，同时添加非法核苷酸的检查逻辑，最后统一打印结果。

修复后的代码：

import sys

seq = 'TGCCTTGGGCACCATGCAGTACCAAACGGAACGATAGTG'

# 初始化所有计数，即使N为0也会被赋值
a_nt = seq.count('A')
g_nt = seq.count('G')
c_nt = seq.count('C')
t_nt = seq.count('T')
n_nt = seq.count('N')

# 检查序列中是否存在非法核苷酸
for nucleotide in seq:
    if nucleotide not in 'AGCTN':
        sys.exit("no this code")

# 打印所有计数结果
print(a_nt, g_nt, c_nt, t_nt, n_nt)

• 直接计数：用seq.count()一次性获取各类核苷酸数量，避免冗余循环，代码更高效简洁。
• 提前初始化：所有计数变量一开始就被赋值，彻底解决未定义报错问题。
• 非法字符检查：单独遍历序列一次，发现非AGCTN字符立即退出，完全符合需求。

如果想通过循环统计来理解计数逻辑，也可以用字典管理计数（扩展性更好）：

import sys

seq = 'TGCCTTGGGCACCATGCAGTACCAAACGGAACGATAGTG'

# 初始化所有核苷酸计数为0
counts = {'A': 0, 'G': 0, 'C': 0, 'T': 0, 'N': 0}

for nucleotide in seq:
    if nucleotide in counts:
        counts[nucleotide] += 1
    else:
        sys.exit("no this code")

# 按要求顺序打印结果
print(counts['A'], counts['G'], counts['C'], counts['T'], counts['N'])

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