Python用while与index计算等长序列一致性返回None排查

现有代码核心错误点 #

• 无返回值：Python函数未显式声明return时默认返回None，这是输出结果显示None%的直接原因
• while循环逻辑完全错误：初始值i=0、num_matches=0，循环条件i < num_matches从初始状态就不成立，循环体一次都不会执行；且循环内没有编写对应索引位置的字符比对、匹配计数更新的逻辑
• 异常分支无中断：序列长度不等时仅打印提示，没有终止函数执行，后续如果强行运行索引遍历会触发越界报错
• 缺失一致性计算逻辑：序列一致性的计算公式为(匹配位点数 / 序列总长度) * 100，原有代码没有实现该计算步骤

修正后可运行代码 #

严格满足提出的两个要求：先输出每条序列长度、使用while循环+索引遍历统计匹配位点

# 示例序列列表
sequences = ['AUUGCAUG', 'CGUGGCUA']

def compute_sequence_identity(sequence_a, sequence_b):
    # 计算并打印每条序列的长度
    length_a = len(sequence_a)
    length_b = len(sequence_b)
    print(f"序列1长度：{length_a}，序列2长度：{length_b}")

    # 长度不等直接返回
    if length_a != length_b:
        print("Sequences are not of equal length.")
        return None

    # 初始化索引、匹配计数器
    idx = 0
    match_count = 0
    # while循环遍历所有索引位置
    while idx < length_a:
        # 比对当前索引位置的字符是否一致
        if sequence_a[idx] == sequence_b[idx]:
            match_count += 1
        # 索引向后移动一位
        idx += 1

    # 计算一致性百分比，保留2位小数
    identity = round((match_count / length_a) * 100, 2)
    return identity

# 调用函数计算
sequence_identity = compute_sequence_identity(sequences[0], sequences[1])
print("Sequence identity: " + str(sequence_identity) + "%")

运行结果说明 #

针对给出的示例序列，代码运行后会先输出两条序列长度均为8，逐位比对后匹配位点数为4，最终输出序列一致性为50.0%。

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