trackplot 0.6.3__tar.gz → 0.6.4__tar.gz

{trackplot-0.6.3/trackplot.egg-info → trackplot-0.6.4}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: trackplot
-Version: 0.6.3
+Version: 0.6.4
 Summary: The trackplot is a tool for visualizing various next-generation sequencing (NGS) data, including DNA-seq, RNA-seq, single-cell RNA-seq and full-length sequencing datasets. https://sashimi.readthedocs.io/
 Author-email: ygidtu <ygidtu@gmail.com>
 License-Expression: BSD-3-Clause

{trackplot-0.6.3 → trackplot-0.6.4}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "trackplot"
-version = "0.6.3"
+version = "0.6.4"
 description = "The trackplot is a tool for visualizing various next-generation sequencing (NGS) data, including DNA-seq, RNA-seq, single-cell RNA-seq and full-length sequencing datasets. https://sashimi.readthedocs.io/"
 authors = [
     { name = "ygidtu", email = "ygidtu@gmail.com" }

{trackplot-0.6.3 → trackplot-0.6.4}/trackplot/cli.py

@@ -1214,6 +1214,30 @@ def _add_interval(p, f):
 @optgroup.option(
     "--title", type=click.STRING, default=None, help="Title", show_default=True
 )
+@optgroup.option(
+    "--no-title",
+    type=click.BOOL,
+    is_flag=True,
+    default=False,
+    help="Disable the plot title entirely",
+    show_default=True,
+)
+@optgroup.option(
+    "--no-region-text",
+    type=click.BOOL,
+    is_flag=True,
+    default=False,
+    help="Disable the --event region text on the x-axis ruler",
+    show_default=True,
+)
+@optgroup.option(
+    "--junctions-on-top",
+    type=click.BOOL,
+    is_flag=True,
+    default=False,
+    help="Force all junction arcs to be drawn above the density plot (only effective when --density-by-strand is not set)",
+    show_default=True,
+)
 @optgroup.option(
     "--font", type=click.STRING, default=None, help="Fonts", show_default=True
 )
@@ -1362,6 +1386,10 @@ def main(**kwargs):
         normalize_format=kwargs.get("normalize_format"),
         fill_step=kwargs.get("fill_step", "post"),
         smooth_bin=kwargs["smooth_bin"],
+        title=kwargs["title"],
+        no_title=kwargs["no_title"],
+        no_region_text=kwargs["no_region_text"],
+        junctions_on_top=kwargs["junctions_on_top"],
     )
 
 

{trackplot-0.6.3 → trackplot-0.6.4}/trackplot/plot/coord.py

@@ -134,14 +134,17 @@ def set_x_ticks(
         x_label = f"{x_label}, y axis is {log_trans} transformed"
 
     ax.hlines(y=0, xmin=0, xmax=max(graph_coords), color="black", lw=1)
-    ax.text(
-        x=graph_coords[len(graph_coords) // 2],
-        y=-2.8,
-        s=x_label,
-        fontsize=font_size,
-        ha="center",
-        va="top",
-    )
+
+    no_region_text = kwargs.pop("no_region_text", False)
+    if not no_region_text:
+        ax.text(
+            x=graph_coords[len(graph_coords) // 2],
+            y=-2.8,
+            s=x_label,
+            fontsize=font_size,
+            ha="center",
+            va="top",
+        )
 
     bk = 1
     if not sequence and nx_ticks > 1:

{trackplot-0.6.3 → trackplot-0.6.4}/trackplot/plot/core.py

@@ -986,6 +986,11 @@ class Plot(object):
 
         igv_height_scale = kwargs.get("igv_height_scale")
 
+        # Extract / pop custom params before they flow into downstream **kwargs
+        title = kwargs.pop("title", None)
+        no_title = kwargs.pop("no_title", False)
+        junctions_on_top = kwargs.pop("junctions_on_top", False)
+
         # ====== Phase 1: Load data ======
         self._load_plot_data(n_jobs, *args, **kwargs)
 
@@ -1048,8 +1053,8 @@ class Plot(object):
             else:
                 ax_var = plt.subplot(gs[curr_idx, 0])
 
-            if curr_idx == 0:
-                ax_var.set_title(str(self.region), loc="left")
+            if curr_idx == 0 and not no_title:
+                ax_var.set_title(title or str(self.region), loc="left")
 
             max_y_val_, min_y_val_ = self._resolve_plot_y_limits(
                 p,
@@ -1088,6 +1093,7 @@ class Plot(object):
                             distance_between_label_axis=distance_between_label_axis,
                             raster=raster,
                             fill_step=fill_step,
+                            junctions_on_top=junctions_on_top,
                             **temp_params,
                         )
                         curr_idx += 1
@@ -1103,6 +1109,7 @@ class Plot(object):
                         distance_between_label_axis=distance_between_label_axis,
                         raster=raster,
                         fill_step=fill_step,
+                        junctions_on_top=junctions_on_top,
                         **self.params.get(p, {}),
                     )
             elif p.type == "hic":
@@ -1123,6 +1130,7 @@ class Plot(object):
                     min_used_y_val=min_y_val_,
                     distance_between_label_axis=distance_between_label_axis,
                     raster=raster,
+                    junctions_on_top=junctions_on_top,
                     **self.params.get(p, {}),
                 )
                 curr_idx += 1

{trackplot-0.6.3 → trackplot-0.6.4}/trackplot/plot/limits.py

@@ -26,6 +26,7 @@ def _compute_y_limits(
     density_by_strand: bool = False,
     fill_step: str = "post",
     show_mean_jxn_number: bool = False,
+    junctions_on_top: bool = False,
 ) -> Tuple[float, float]:
     """
     Compute y-axis limits from ReadDepth data, including junction arc extents.
@@ -119,6 +120,8 @@ def _compute_y_limits(
 
             if density_by_strand:
                 jxn_on_top = jxn.strand == "+"
+            elif junctions_on_top:
+                jxn_on_top = True
             else:
                 jxn_on_top = jxn_idx % 2 == 0
                 if abs(min_used_y_val) < max_used_y_val:
@@ -163,7 +166,7 @@ def _compute_y_limits(
         if density_by_strand:
             max_used_y_val = max(abs(min_used_y_val), max_used_y_val)
             min_used_y_val = -max_used_y_val
-    elif not density_by_strand and not jxns:
+    elif not density_by_strand and (junctions_on_top or not jxns):
         min_used_y_val = 0
 
     # Small expansion so the top / bottom arcs have visual breathing room
@@ -207,4 +210,5 @@ def precompute_y_limits(
         density_by_strand=kwargs.get("density_by_strand", False),
         fill_step=fill_step,
         show_mean_jxn_number=False,
+        junctions_on_top=kwargs.get("junctions_on_top", False),
     )

{trackplot-0.6.3 → trackplot-0.6.4}/trackplot/plot/render.py

@@ -493,6 +493,8 @@ def plot_density(
     # Compute y limits using shared logic
     fixed_max_used_y = max_used_y_val is not None
 
+    junctions_on_top = kwargs.get("junctions_on_top", False)
+
     max_used_y_val, min_used_y_val = _compute_y_limits(
         data=data,
         region=region,
@@ -502,6 +504,7 @@ def plot_density(
         density_by_strand=kwargs.get("density_by_strand", False),
         fill_step=fill_step,
         show_mean_jxn_number=show_mean_jxn_number,
+        junctions_on_top=junctions_on_top,
     )
     
     # Draw fill
@@ -557,6 +560,8 @@ def plot_density(
 
             if kwargs.get("density_by_strand"):
                 jxn_on_top = jxn.strand == "+"
+            elif junctions_on_top:
+                jxn_on_top = True
             else:
                 jxn_on_top = jxn_idx % 2 == 0
                 if abs(min_used_y_val) < max_used_y_val:
@@ -659,7 +664,7 @@ def plot_density(
     ):
         max_used_y_val = max(abs(min_used_y_val), max_used_y_val)
         min_used_y_val = -max_used_y_val
-    elif not kwargs.get("density_by_strand") and not jxns:
+    elif not kwargs.get("density_by_strand") and (junctions_on_top or not jxns):
         min_used_y_val = 0
 
     set_y_ticks(

{trackplot-0.6.3 → trackplot-0.6.4/trackplot.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: trackplot
-Version: 0.6.3
+Version: 0.6.4
 Summary: The trackplot is a tool for visualizing various next-generation sequencing (NGS) data, including DNA-seq, RNA-seq, single-cell RNA-seq and full-length sequencing datasets. https://sashimi.readthedocs.io/
 Author-email: ygidtu <ygidtu@gmail.com>
 License-Expression: BSD-3-Clause