读取fasta文件 with open('your_fasta_file.fasta', 'r') as f: lines = f.readlines() # 将序列ID和序列分别存储到两个列表中 ids = [line.strip()[1:] for line in lines[::2]] sequences = [line.strip() for line in lines[1::2]] # 将序列ID和序列合并为一行 fasta = [f'>{id} {sequence}' for id, sequence in zip(ids, sequences)] # 将fasta文件写入到文件中 with open('output.fasta', 'w') as f: f.write('\n'.join(fasta))
时间: 2024-04-01 17:37:00 浏览: 97
通过指定目录,读取该目录下所有子文件夹中的内容 输入根目录和指定类别的名字,即可读取所有文件
这个问题是一个关于读取fasta文件并将其写入到输出文件的代码段,你想知道这段代码的作用是什么。实际上,这段代码可以帮助你将一个fasta格式的文件读取到内存中,并将其重新格式化为一个新的fasta文件输出。在这个代码段中,我们首先打开一个fasta文件并读取其中的所有行,然后将序列ID和序列分别存储到两个列表中。接下来,我们使用zip函数将这两个列表合并为一个fasta格式的条目,这个条目包含一个序列ID和一个序列。最后,我们将这些fasta格式的条目写入到一个新的文件中。
阅读全文