编写一个python程序，读取系列文件，然后统计各密码子的使用频率，并将频率表输出到结果文件（CSV文件），输出为每行一个密码子的使用频率

好的，以下是 Python 代码实现，假设密码子文件名为 "codons.txt"，结果文件名为 "frequency.csv"：

# 密码子与对应的氨基酸
codon_table = {
    "TTT": "F", "TTC": "F", "TTA": "L", "TTG": "L",
    "CTT": "L", "CTC": "L", "CTA": "L", "CTG": "L",
    "ATT": "I", "ATC": "I", "ATA": "I", "ATG": "M",
    "GTT": "V", "GTC": "V", "GTA": "V", "GTG": "V",
    "TCT": "S", "TCC": "S", "TCA": "S", "TCG": "S",
    "CCT": "P", "CCC": "P", "CCA": "P", "CCG": "P",
    "ACT": "T", "ACC": "T", "ACA": "T", "ACG": "T",
    "GCT": "A", "GCC": "A", "GCA": "A", "GCG": "A",
    "TAT": "Y", "TAC": "Y", "TAA": "*", "TAG": "*",
    "CAT": "H", "CAC": "H", "CAA": "Q", "CAG": "Q",
    "AAT": "N", "AAC": "N", "AAA": "K", "AAG": "K",
    "GAT": "D", "GAC": "D", "GAA": "E", "GAG": "E",
    "TGT": "C", "TGC": "C", "TGA": "*", "TGG": "W",
    "CGT": "R", "CGC": "R", "CGA": "R", "CGG": "R",
    "AGT": "S", "AGC": "S", "AGA": "R", "AGG": "R",
    "GGT": "G", "GGC": "G", "GGA": "G", "GGG": "G",
}

# 初始化密码子频率字典
codon_freq = {codon: 0 for codon in codon_table}

# 读取密码子文件，统计频率
with open("codons.txt", "r") as f:
    for line in f:
        for i in range(0, len(line), 3):
            codon = line[i:i+3]
            if codon in codon_table:
                codon_freq[codon] += 1

# 将频率表输出到结果文件
with open("frequency.csv", "w") as f:
    f.write("Codon,Frequency\n")
    for codon in codon_table:
        f.write(f"{codon},{codon_freq[codon]}\n")

代码中，我们首先定义了一个密码子与氨基酸的对应字典 `codon_table`，然后初始化密码子频率字典 `codon_freq`，将每个密码子的初始频率设为 0。接着，我们打开密码子文件，循环读取每行密码子，统计每个密码子出现的次数，将结果保存在 `codon_freq` 中。最后，我们打开结果文件，将密码子与频率写入文件中，每行一个。