Example #1
0
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
Example #2
0
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