Python实现DNA序列中AGAT/AATG/TATC重复次数统计的代码求助

• 语法逻辑误区：判断条件lines in DNAsample=='AGAT'存在初学阶段常见的连续比较语法误解，Python会把这类连续判断解析为(lines in DNAsample) and (DNAsample == 'AGAT')——既没有判断文本内容中是否存在目标序列，也因为文件对象本身不可能等于字符串'AGAT'，所有判断分支永远不会触发。
• 流式文件读取的机制误解：Python的文件对象是单向流式结构，第一次逐行遍历完成后，文件指针会停在文件末尾，后续两次重复遍历DNAsample的循环根本读不到任何有效内容，对应的计数器永远停留在初始值0。
• 计数逻辑错误：针对TATC的计数分支，无论匹配条件是否满足，都写了tatc_Counter+=0，完全没有自增逻辑；同时原代码判断整行内容和目标序列完全相等，但实际DNA序列文件通常是一整段连续的长字符串，不会单独把4字符的目标序列作为整行存储，整行相等的判断完全不符合实际文件的存储场景。
• 冗余无效代码：所有else分支的+=0对计数器没有任何实际作用，属于完全多余的逻辑；原代码打开文件后没有自动关闭的机制，长期运行容易造成文件资源占用。

针对编程入门者的认知水平，以下代码做了逻辑简化，同时覆盖了序列重叠匹配的统计需求：

# 初始化三个目标序列的计数器
agat_Counter = 0
aatg_Counter = 0
tatc_Counter = 0

# 使用with语句打开文件，自动管理文件资源，无需手动调用close()
with open('DNA SEQUENCE FILE.txt', 'r', encoding='utf-8') as f:
    # 一次性读取全文件内容，去除首尾换行/空格，统一转大写避免大小写不匹配
    dna_str = f.read().strip().upper()

# 三个目标序列长度均为4，逐位置滑动窗口匹配
target_length = 4
# 遍历所有可能的起始位置，避免漏算
for start_idx in range(len(dna_str) - target_length + 1):
    # 截取当前位置长度为4的子串做匹配
    sub_seq = dna_str[start_idx:start_idx+target_length]
    if sub_seq == "AGAT":
        agat_Counter += 1
    elif sub_seq == "AATG":
        aatg_Counter += 1
    elif sub_seq == "TATC":
        tatc_Counter += 1

# 输出最终统计结果
print(f"AGAT 重复次数：{agat_Counter}")
print(f"AATG 重复次数：{aatg_Counter}")
print(f"TATC 重复次数：{tatc_Counter}")

实现说明：滑动窗口的匹配方式可以覆盖所有相邻位置的序列匹配，包括重叠出现的场景——例如序列中相邻位置出现重叠的目标片段时，也能被准确统计，不会出现漏算。

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