import screed

filename = '../data/25k.fq.gz'

n = 0
for record in screed.open(filename):
    n += 1
    if n > 10:
        break
    print record.name
    print record.sequence
    print record.accuracy

import calc_gc

gc_list = []
for record in screed.open(filename):
    seq = record.sequence
    gc = calc_gc.calc_gc(seq)
    gc_list.append(gc)

hist(gc_list, bins=100, range=(0, 1))

import calc_gc

gc_list = []
for record in screed.open(filename):
    seq = record.sequence[:10]
    gc = calc_gc.calc_gc(seq)
    gc_list.append(gc)

hist(gc_list, bins=10, range=(0, 1))