--- 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)