如何统计字符串列表各位置G/C字符的占比？现有代码待修正

需求：对字符串列表中的每个字符串逐字符比对，判断字符是否为G或C，最终计算列表中每个字符位置上匹配该指定字符的百分比。

用户初始代码：

def GC_content_pos(reads_list):
    for read in reads_list:
        for position in range(len(read)):
            if read[position] == "G" or read[position] == "C":
                #do something

实现步骤与完整代码 #

1. 初始化计数数组：先确认所有输入字符串长度一致（这是按位置统计的前提），创建一个和字符串长度相同的列表，初始值全为0，用于记录每个位置的GC字符数量。
2. 遍历统计GC数量：逐个遍历列表中的字符串，再遍历每个字符位置，若字符是G或C，对应位置的计数加1。
3. 计算百分比：用每个位置的GC计数除以总字符串数量，再乘以100得到该位置的GC百分比。

完整代码：

def GC_content_pos(reads_list):
    # 处理空输入的情况
    if not reads_list:
        return []
    
    # 假设所有字符串长度一致，取第一个字符串的长度初始化计数数组
    pos_counts = [0] * len(reads_list[0])
    total_reads = len(reads_list)
    
    for read in reads_list:
        # 校验当前字符串长度是否和预期一致，避免报错
        if len(read) != len(pos_counts):
            raise ValueError("所有字符串长度必须一致")
        
        for position in range(len(read)):
            if read[position] in ("G", "C"):
                pos_counts[position] += 1
    
    # 计算每个位置的GC百分比，可按需调整小数位数
    gc_percentages = [round((count / total_reads) * 100, 2) for count in pos_counts]
    return gc_percentages

代码说明 #

• 增加空输入判断，避免后续逻辑报错；
• 加入字符串长度校验，确保输入符合统计要求，否则抛出明确错误；
• 用in ("G", "C")替代连续或判断，代码更简洁易读；
• 通过列表推导式快速计算百分比，用round()控制小数位数，适配不同精度需求。

内容的提问来源于stack exchange，提问作者C. van Veen