Open aprilwzh opened 5 years ago
`fr = open('fastq','r') fw = open('81.out','w') fr = fr.readlines() n = len(fr) i = 0 while i < n : m = i%4 if m == 0: name = fr[i] fw.write(name) fw.write('\n') i = i +1 elif m == 1: seq = fr[i] fw.write(seq) fw.write('\n') i = i +1 else: i = i +1
else: fw.close()`
fastq格式为四行,前两列有用,所以就把文件转为列表,然后统计行数,然后进行除四取余输出。
`fr = open('fastq','r') fw = open('81.out','w') fr = fr.readlines() n = len(fr) i = 0 while i < n : m = i%4 if m == 0: name = fr[i] fw.write(name) fw.write('\n') i = i +1 elif m == 1: seq = fr[i] fw.write(seq) fw.write('\n') i = i +1 else: i = i +1
else: fw.close()`