dist_hist = get_distance_histogram(data_folder,
adaID,
fragment,
VERBOSE=VERBOSE)
except IOError:
continue
dist_hists.append((samplename_seq, fragment, dist_hist))
dist_hists.sort(key=itemgetter(1))
fig, ax = plt.subplots()
for i, (samplename_seq, fragment, h) in enumerate(dist_hists):
plot_distance_histogram(h,
ax=ax,
color=cm.jet(1.0 * i / len(dist_hists)),
label=', '.join([samplename_seq,
fragment]))
ax.set_title(samplename)
ax.legend(loc=1, fontsize=10)
if use_save:
foldername = sample.get_foldername() + 'figures/'
mkdirs(foldername)
fn = foldername + 'distance_to_consensus_seqsamples.png'
fig.savefig(fn)
plt.close(fig)
if not use_save:
plt.ion()
plt.show()
