Updating scripts with alignment_length as global (was undefined and t…

…hus taken by default) : now also for the plot title
bonsai-team · May 20, 2022 · f94fa77 · f94fa77
1 parent 5ed210c
commit f94fa77
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git 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
@@ -19,7 +19,7 @@ def seeds_coefficients(seed_pattern = "###-#--#-#--##-###", alignment_length = 6
     @return tuple(list<int>,list<int>,list<int>,list<float>)  : 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()