写一个严谨的python脚本，不依赖模块提取cds.fa文件中每个基因的最长转录本输出

以下是一个可能的python脚本:

# 打开cds.fa文件
with open("cds.fa") as file:
    genome = file.read()

# 将基因组按照每个基因分割为列表
genes = genome.split(">")

# 创建一个空字典来存储每个基因的最长转录本
longest_transcripts = {}

# 循环每个基因，提取其转录本
for gene in genes[1:]:  # 跳过第一个空字符串
    gene_lines = gene.split("\n")
    gene_name = gene_lines[0]
    seq = "".join(gene_lines[1:]).upper()  # 合并行并转换为大写

    # 初始化一个字典来存储每个转录本的长度
    transcript_lengths = {}

    # 循环每个可能转录本的开头（注意此处使用range(len(seq) - 2)，因为每个开头必须有一个起始密码子）
    for i in range(len(seq) - 2):
        # 如果开头是一个起始密码子，从这个位置开始查找最近的终止密码子
        if seq[i:i+3] == "ATG":
            for j in range(i+3, len(seq), 3):
                if seq[j:j+3] in ["TAA", "TAG", "TGA"]:
                    transcript = seq[i:j+3]
                    length = len(transcript)
                    transcript_lengths[transcript] = length
                    break  # 终止密码子已找到，不必再搜索

    # 使用max函数找到最长的转录本并添加到longest_transcripts字典中
    longest_transcript = max(transcript_lengths, key=transcript_lengths.get)
    longest_transcripts[gene_name] = longest_transcript

# 输出每个基因的最长转录本
for gene_name, longest_transcript in longest_transcripts.items():
    print(">" + gene_name + " longest transcript")
    print(longest_transcript)

这个脚本会读取名为"cds.fa"的文件，循环每个基因并提取其转录本，然后找到每个基因的最长转录本。最后，它会将每个基因的最长转录本输出到屏幕上。此脚本不依赖于第三方模块。