Python新手使用def定义GC含量计算函数无预期结果问题求助

• 核心问题：你仅完成了GC计算函数的定义，没有主动调用该函数传入读取的fasta内容执行计算，自然不会输出预期结果。
• 潜在隐藏问题：如果读取的是标准FASTA格式文件，文件内容包含>开头的序列注释行、换行符等非目标碱基字符，直接统计会导致碱基计数、序列长度计算错误，和你注释中直接用纯序列测试的结果不一致。

修正后可运行代码 #

# 定义GC计算函数
def GC(fasta): 
    gc = sum(fasta.count(x) for x in["G","C"])
    return gc * 100.0 / len(fasta)

# 读取文件
with open ('/Users/files') as f: 
    fasta = f.read()
    print("读取的原始文件内容：\n", fasta) 

# 调用函数计算并打印结果
result = GC(fasta)
print("GC含量结果：", result)

# 如果是标准FASTA文件，需要先清洗数据再计算，示例如下：
# 清洗逻辑：过滤掉>开头的注释行，去掉所有空白字符，仅保留碱基序列
# cleaned_seq = ''.join([line.strip() for line in fasta.splitlines() if not line.startswith('>')])
# cleaned_result = GC(cleaned_seq)
# print("清洗序列后的GC含量结果：", cleaned_result)

验证说明 #

你注释中的计算逻辑和封装的函数计算逻辑完全等价，只要传入的fasta参数一致，得到的结果必然和你预期的27.75相同，按上述代码补充函数调用步骤即可得到正确结果。

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