У меня есть этот код, который читает из файла и записывает в другой файл, но файл записи пуст.
f=open(FASTQFILE, "r")
fo=open(FASTAFILE, "w")
while True:
header=f.read()
if not header:break
header=header.strip()
seq=f.readline()
seq=seq.strip()
ignore=f.readline()
qualities=f.readline()
result=filter_quality(qualities.strip())
if (result >MINSIZE):
header=header[0].replace("@", ">") + header[1:]
fo.write(header+ "\n")
fo.write(seq[:result] + "\n")
f.close()
fo.close()
в основном эта функция превращает fastq в файл fasta