def extract_reads_to_file(filename, start, stop): _, tmpfile = mapreads.get_temp_filename('readchunk.fa') fp = open(tmpfile, 'w') for record in retrieve_records(filename, start, stop): fp.write('>%s\n%s\n' % (record.name, record.sequence)) fp.close() return tmpfile
def extract_reads_to_file(filename, start, stop): _, tmpfile = mapreads.get_temp_filename('readchunk.fa') fp = open(tmpfile, 'w') for data in retrieve_bytes(filename, start, stop): fp.write(data) fp.close() return tmpfile