用dict函数编写python脚本,读取多个蛋白序列fasta格式的序列文件,并输出无回车和空行的fasta文件

你可以使用Python的`csv`模块配合`dict`和文件操作来实现这个功能。首先，你需要遍历指定目录下的所有 `.fasta` 文件，然后逐行解析每个文件中的内容，并将它们添加到一个字典中，键为条目标题（通常是序列ID），值为序列。最后，将处理后的数据写入一个新的`.fasta` 文件，去除回车和空行。

下面是一个简单的示例：

import csv

def read_fasta_files(input_dir, output_file):
    seq_dict = {}
    
    # 遍历输入目录中的所有 .fasta 文件
    for filename in sorted(os.listdir(input_dir)):
        if filename.endswith(".fasta"):
            with open(os.path.join(input_dir, filename), 'r') as file:
                reader = csv.reader(file, delimiter='\n', quotechar=None)  # 使用'\n'作为分隔符

                # 翻页读取，跳过头部
                next(reader)
                
                # 读取每一条序列
                for i, row in enumerate(reader):
                    if i == 0:  # 标题行
                        title = row[0]
                    else:  # 序列数据
                        seq = ''.join(row).strip()  # 去除首尾空白并连接成字符串
                        seq_dict[title] = seq

    # 写入新文件
    with open(output_file, 'w') as out_file:
        for title, seq in seq_dict.items():
            out_file.write(f">{title}\n{seq}\n")

# 调用函数，例如
read_fasta_files('input_folder', 'output.fasta')

在这个脚本中，我们假设输入文件夹(`input_folder`)下有多个`.fasta` 文件，且文件内顺序是按标题排序的。脚本会创建一个新的名为`output.fasta` 的文件，其中包含合并和格式化的序列数据。