!python sandbox/bin-reads-by-abundance.py data/stamps-reads.fa.gz 

!python scripts/abundance-dist-single.py -s -k 20 -x 1e8 stamps-reads.fa.gz.bink.2 bink2.hist
!python scripts/abundance-dist-single.py -s -k 20 -x 1e8 stamps-reads.fa.gz.bink.1 bink1.hist

raw = numpy.loadtxt('raw-reads.hist')
bink = numpy.loadtxt('bink-reads.hist')
bink1 = numpy.loadtxt('bink1.hist')
bink2 = numpy.loadtxt('bink2.hist')

plot(raw[:,0], raw[:,1])
plot(bink[:,0], bink[:,1])
axis(ymax=200)

plot(raw[:,0], raw[:,1])
plot(bink[:,0], bink[:,1])
axis(ymax=200, xmax=200)



plot(bink1[:,0], bink1[:,1])
plot(bink2[:,0], bink2[:,1])
axis(ymax=800)

!python scripts/load-into-counting.py -x 1e7 -k 20 bink1.kh stamps-reads.fa.gz.bink.1
!python scripts/abundance-dist.py -s bink1.kh a.fa bink1.x.a.hist
!python scripts/abundance-dist.py -s bink1.kh b.fa bink1.x.b.hist

!python scripts/load-into-counting.py -x 1e7 -k 20 bink2.kh stamps-reads.fa.gz.bink.2
!python scripts/abundance-dist.py -s bink2.kh a.fa bink2.x.a.hist
!python scripts/abundance-dist.py -s bink2.kh b.fa bink2.x.b.hist

bink1xa = numpy.loadtxt('bink1.x.a.hist')
bink1xb = numpy.loadtxt('bink1.x.b.hist')
bink2xa = numpy.loadtxt('bink2.x.a.hist')
bink2xb = numpy.loadtxt('bink2.x.b.hist')


plot(bink1xa[:,0], bink1xa[:,1], label='bink1 x a')
plot(bink1xb[:,0], bink1xb[:,1], label='bink1 x b')
plot(bink2xa[:,0], bink2xa[:,1], label='bink2 x a')
plot(bink2xb[:,0], bink2xb[:,1], label='bink2 x b')
legend()



`