将gbk格式文件转换成fasta格式文件
flag=0fasta = open ('/mnt/f/biopython/AY810830.fasta','w')with open ('/mnt/f/biopython/AY810830.gbk','r') as f: for line in f: if line[0:9] =='ACCESSION':
fasta.writelines('>'+line.split()[1]+'n') elif line[0:6] =='ORIGIN':
flag = 1
elif flag == 1:
s = line.split() if s != []:
print(s)
seq=''.join(s[1:])
fasta.writelines(seq.upper()+'n')
fasta.close()将多个fasta序列中提取含有'Homo sapiens'的序列
fasta_file=open('multi.fasta','r')
out_file=open('human.fasta','w')
seq = ''for line in fasta_file: if line[0]=='>' and seq =='':
header = line elif line[0] != '>':
seq=seq+line elif line[0] =='>' and seq !='': if "Homo sapiens" in header:
out_file.write(header+seq)
seq = ''
header = lineif "Homo sapiens" in header :
out_file.write(header+seq)
out_file.close()
作者：tianzhanlan
链接：https://www.jianshu.com/p/492a6e4a1b19
