【发布时间】:2021-04-07 19:35:35
【问题描述】:
我有一个 fasta 文件,我想为每个序列制作单独的文件。每个序列的文件名应该是唯一的,所以我使用了序列的描述。它给出了一个错误:output=open(header,"w") 错误说:OSError: [Errno 22] Invalid argument: '>g1.fasta'
代码:
in_file = open("test.fasta","r")
lines=in_file.read().strip()
lines=lines.replace("\r\n","\n")
in_file.close()
sequences=lines.split("\n>")
in_file.close()
print(sequences)
for sequence in sequences:
elements=sequence.split("\n")
header = elements[0]
seq = "".join(elements[1:])
name=header+".fasta"
output=open(name,"w")
output.write(seq)
output.close()
【问题讨论】:
-
请提供一个test.fasta文件的例子
标签: python python-2.7 loops for-loop