python处理fasta文件，ID和序列放在一行

#!/usr/bin/python#-*- coding:utf-8 -*-\"处理fasta⽂件，将ID号和序列放在⼀⾏\"import sys

with open(sys.argv[1]) as f: fw=open('out.fasta', 'w') line=f.read()

line=line.replace('\\n', '').replace('>', '\\n>') for aa in line: fw.write(aa) fw.close()\"\"\"

>chr1|hos107.1#gene1

ACACTCCCGGGCCCCCCCCCCCCACCTTTCAAAAAAAAAAAAAAAAATTTTCCCCCCAAAGGGG>chr1|hos107.2#gene2

ACACTCCCGGGCCCCCCCCCCCCACCTTTCAAAAAAAAAAAAAAAAATTTTC

>chr1|hos107.4#gene3

ACACTCCCGGGCCCCCCCCCCCCACCTTTCAAAAAAAAAAAAAAAAATTTTC

>chr1|hos107.5#gene4

ACACTCCCGGGCCCCCCCCCCCCACCTTTCAAAAAAAAAAAAAAAAATTTTC\"\"\"\"\"\"

>chr1|hos107.1#gene1ACACTCCCGGGCCCCCCCCCCCCACCTTTCAAAAAAAAAAAAAAAAATTTTCCCCCCAAAGGGG>chr1|hos107.2#gene2ACACTCCCGGGCCCCCCCCCCCCACCTTTCAAAAAAAAAAAAAAAAATTTTC>chr1|hos107.4#gene3ACACTCCCGGGCCCCCCCCCCCCACCTTTCAAAAAAAAAAAAAAAAATTTTC>chr1|hos107.5#gene4ACACTCCCGGGCCCCCCCCCCCCACCTTTCAAAAAAAAAAAAAAAAATTTTC\"\"\"

#提取⽬标序列

f=open('./out.fasta', 'r')fw=open('target.fasta', 'w') for line in f.readlines():

if line.startswith('>chr1|hos107.1'): fw.write(line)f.close()fw.close()

\"\"\"可以从上述处理好的单⾏⽂件out.fasta中提取指定⽬标ID的⽂件，并将其写⼊到target.fasta⽂件中\"\"\"

#整体思路：

#先统⼀fasta⽂件格式从test.fasta----out.fasta#取出⽬标ID序列：out.fasta----target.fasta