Python 脚本：快速氨基酸序列突变

以下是一个用于快速针对氨基酸序列进行突变的 Python 脚本：

import csv

def mutate_sequence(sequence, mutations):
    mutated_sequence = list(sequence)
    for mutation in mutations:
        position = int(mutation[1:]) - 1
        mutated_sequence[position] = mutation[0]
    return ''.join(mutated_sequence)

def update_mut_csv(input_csv, output_csv):
    with open(input_csv, 'r') as file:
        reader = csv.reader(file)
        mutations = next(reader)
        sequence = next(reader)

        mutated_sequence = mutate_sequence(sequence[0], mutations)

        with open(output_csv, 'w', newline='') as outfile:
            writer = csv.writer(outfile)
            writer.writerow(mutations)
            writer.writerow(sequence)
            writer.writerow([mutated_sequence])

# 示例用法
input_csv = 'mut.csv'
output_csv = 'new_mut.csv'
update_mut_csv(input_csv, output_csv)

使用时，将原始的 mut.csv 文件放在与脚本相同的目录下，并确保 mut.csv 文件的格式为：

mutation1,mutation2,...
sequence

例如，'L2P' 表示将第二个氨基酸由 L 替换为 P。

然后运行脚本，它将生成一个新的 mut.csv 文件（new_mut.csv），其中 C3 单元格将包含突变后的序列。