def DrawBaseStats(options, base_lens, base_freqs, log, min_mult):
hist_name1 = os.path.join(options.output_dir, "base_seq_lens.png")
len_hist_settings = drawing_utils.GetGraphicalSettings(
xlabel="Sequence length",
ylabel="# sequences",
output_filename=hist_name1)
drawing_utils.DrawHistogram(base_lens, len_hist_settings)
if os.path.exists(hist_name1):
log.info(
"* Histogram of distribution of base sequence lengths was written to "
+ hist_name1)
else:
log.info(
"ERROR: Histogram of distribution of base sequence lengths was not found"
)
sys.exit(1)
hist_name2 = os.path.join(options.output_dir, "base_seq_freqs.png")
freq_hist_settings = drawing_utils.GetGraphicalSettings(
xlabel="Base sequence frequency (>" + str(min_mult) + ")",
ylabel="# sequences",
output_filename=hist_name2)
drawing_utils.DrawHistogram(base_freqs, freq_hist_settings)
if os.path.exists(hist_name2):
log.info(
"* Histogram of distribution of base sequence frequencies was written to "
+ hist_name2)
else:
log.info(
"ERROR: Histogram of distribution of base sequence frequencies was not found"
)
sys.exit(1)
def DrawMutatedStats(options, mutated_freqs, mutation_pos, log, min_mult):
hist_name1 = os.path.join(options.output_dir, "mutated_seq_freqs.png")
freq_hist_settings = drawing_utils.GetGraphicalSettings(
xlabel="Mutated sequence frequency (>" + str(min_mult) + ")",
ylabel="# sequences",
output_filename=hist_name1)
drawing_utils.DrawHistogram(mutated_freqs, freq_hist_settings)
if os.path.exists(hist_name1):
log.info(
"* Histogram of distribution of mutated sequence frequencies was written to "
+ hist_name1)
else:
log.info(
"ERROR: Histogram of distribution of mutated sequence frequencies was not found"
)
sys.exit(1)
hist_name2 = os.path.join(options.output_dir, "shm_positions.png")
pos_hist_settings = drawing_utils.GetGraphicalSettings(
xlabel="Relative SHM position",
ylabel="# SHMs",
output_filename=hist_name2)
drawing_utils.DrawMutationHistogram(mutation_pos, pos_hist_settings,
options.chain_type)
if os.path.exists(hist_name2):
log.info(
"* Histogram of distribution of SHM positions was written to " +
hist_name2)
else:
log.info(
"ERROR: Histogram of distribution of SHM positions was not found")
sys.exit(1)