# HG changeset patch
# User mvdbeek
# Date 1446219396 14400
# Node ID ed151db39c7e1d5d7fbbc3d2ffa639da9a0f2af9
# Parent  e041f9f3271f5208fb3c048346b4ccb5ced61260
planemo upload for repository https://github.com/mvdbeek/dapars commit b1b007c561ea6c9db145c88b6b128d66ecd05e24-dirty

diff -r e041f9f3271f -r ed151db39c7e dapars.py
--- a/dapars.py	Fri Oct 30 10:46:41 2015 -0400
+++ b/dapars.py	Fri Oct 30 11:36:36 2015 -0400
@@ -327,17 +327,17 @@
     ax2.set_ylabel("normalized nucleotide coverage")
     mse_control = [ condition[0] for condition in mse_list]
     mse_treatment = [ condition[1] for condition in mse_list]
-    minima_control = get_minima(np.array(mse_control))
-    minima_treatment = get_minima(np.array(mse_treatment))
+    #minima_control = get_minima(np.array(mse_control))
+    #minima_treatment = get_minima(np.array(mse_treatment))
     control = normalized_utr_coverage[:num_control]
     treatment = normalized_utr_coverage[num_control:]
     ax1.plot(mse_control, "b-")
     ax1.plot(mse_treatment, "r-")
     [ax2.plot(cov, "b-") for cov in control]
     [ax2.plot(cov, "r-") for cov in treatment]
-    [ax2.axvline(val, color="b", alpha=0.25) for val in minima_control]
+    #[ax2.axvline(val, color="b", alpha=0.25) for val in minima_control]
     ax2.axvline(mse_control.index(min(mse_control)), color="b", alpha=1)
-    [ax2.axvline(val, color="r", alpha=0.25) for val in minima_treatment]
+    #[ax2.axvline(val, color="r", alpha=0.25) for val in minima_treatment]
     ax2.axvline(mse_treatment.index(min(mse_treatment)), color="r", alpha=1)
     fig.add_subplot(ax1)
     fig.add_subplot(ax2)