fasta文件id后面加东西 python
时间: 2024-01-07 13:02:37 浏览: 92
你可以使用 Python 的字符串操作来在 FASTA 文件中的 ID 后面添加东西。下面是一个示例代码,它可以在每个 ID 后面添加 "_new"。
```
import os
# 输入文件名和输出文件名
input_file = "example.fasta"
output_file = "example_new.fasta"
# 检查输入文件是否存在
if not os.path.exists(input_file):
print("Error: Input file does not exist!")
exit()
# 打开输入文件和输出文件
with open(input_file, "r") as f_in, open(output_file, "w") as f_out:
# 初始化变量
current_id = ""
current_seq = ""
# 读取输入文件的每一行
for line in f_in:
# 如果这一行是 ID 行,则保存上一个序列并更新当前 ID
if line.startswith(">"):
if current_id:
f_out.write(current_id + "_new\n")
f_out.write(current_seq + "\n")
current_id = line.strip()
current_seq = ""
# 如果这一行是序列行,则将序列加入当前序列中
else:
current_seq += line.strip()
# 保存最后一个序列
if current_id:
f_out.write(current_id + "_new\n")
f_out.write(current_seq + "\n")
print("Done!")
```
在这个代码中,我们首先检查输入文件是否存在,然后打开输入文件和输出文件。接下来,我们逐行读取输入文件,并检查每一行是否是 ID 行。如果是 ID 行,我们将上一个序列写入输出文件,并更新当前 ID。如果这一行是序列行,则将序列加入当前序列中。在处理完所有行之后,我们将最后一个序列写入输出文件中,并关闭输入文件和输出文件。
在这个示例代码中,我们在每个 ID 后面添加了 "_new"。你可以根据需要修改这个代码,在 ID 后面添加其他字符串。
阅读全文