From f94fa771c6a97090594757420ec3b34303f37567 Mon Sep 17 00:00:00 2001 From: noe Date: Fri, 20 May 2022 17:56:25 +0200 Subject: [PATCH] Updating scripts with alignment_length as global (was undefined and thus taken by default) : now also for the plot title --- plot_mow_seeds.py | 2 +- plot_spaced_seeds.py | 4 ++-- 2 files changed, 3 insertions(+), 3 deletions(-) diff --git a/plot_mow_seeds.py b/plot_mow_seeds.py index 79ece48..ecfca61 100755 --- a/plot_mow_seeds.py +++ b/plot_mow_seeds.py @@ -189,7 +189,7 @@ def plot_for_seed(plt = None, seed_pattern = "NNNnNnnNnNnnNNnNNN", color_plot = plt.ylim([-0.01,+1.01]) plt.ylabel("1st hit frequency/probability") -plt.title("1st hit frequency/probability (alignment length = 64)") +plt.title("1st hit frequency/probability (alignment length = "+str(alignment_length)+")") plt.legend(loc="upper left") plt.show() diff --git a/plot_spaced_seeds.py b/plot_spaced_seeds.py index 73e0375..b831581 100755 --- a/plot_spaced_seeds.py +++ b/plot_spaced_seeds.py @@ -19,7 +19,7 @@ def seeds_coefficients(seed_pattern = "###-#--#-#--##-###", alignment_length = 6 @return tuple(list,list,list,list) : the tuple with four lists """ list_tuples = () - err = os.system('./src/iedera -spaced -l '+str(alignment_length)+' -s 11,24 -m "'+seed_pattern+'" -p -o iedera_output.txt') + err = os.system('./src/iedera -spaced -l '+str(alignment_length)+' -s 1,64 -m "'+seed_pattern+'" -p -o iedera_output.txt') if err != 0: sys.exit("./src/iedera not running (please try to run \"./configure; make\" before)") @@ -159,7 +159,7 @@ def plot_for_seed(plt = None, seed_pattern = "###-#--#-#--##-###", color_plot = plt.ylim([-0.01,+1.01]) plt.ylabel("1st hit frequency/probability") -plt.title("1st hit frequency/probability (alignment length = 64)") +plt.title("1st hit frequency/probability (alignment length = "+str(alignment_length)+")") plt.legend(loc="upper left") plt.show()