ls /mnt
ecoli-v41.1.ebwt ecoli-v41.fa qc-reads.map raw-reads.map.count ecoli-v41.2.ebwt ecoli-v41.rev.1.ebwt qc-reads.map.count ecoli-v41.3.ebwt ecoli-v41.rev.2.ebwt raw-reads.fq ecoli-v41.4.ebwt qc-reads.fq raw-reads.map
cd /mnt
/mnt
qc_counts = numpy.loadtxt('qc-reads.map.count')
raw_counts = numpy.loadtxt('raw-reads.map.count')
plot(qc_counts[:,0], qc_counts[:,1], label='qc')
plot(raw_counts[:, 0], raw_counts[:, 1], label='raw')
axis(ymin=0, ymax=80000, xmin=0, xmax=160)
xlabel('position in read')
ylabel('number of reads with mismatches at that position')
legend()
<matplotlib.legend.Legend at 0x7fe5c81910d0>