第一种
with open(seq,‘r’) as f:
line = f.read()
line = line.replace(’\n’,’’).replace(’>’,’\n>’)
第二种在这里插入代码片
def read_fastq(): #读取序列
seqs = []
with open(r’seq’,‘r’) as f:
seq_ls = []
for i in f:
i = i.strip()
if not i:
continue
if i[0] == ‘>’:
if seq_ls :
seqs.append([line,’’.join(seq_ls)])
line = i
else:
seq_ls.append(i)
else:
seqs.append([line,’’.join(seq_ls)])
return seqs
结果