trackplot 0.5.2__tar.gz → 0.5.4__tar.gz

{trackplot-0.5.2 → trackplot-0.5.4}/PKG-INFO

@@ -1,34 +1,30 @@
-Metadata-Version: 2.1
+Metadata-Version: 2.4
 Name: trackplot
-Version: 0.5.2
+Version: 0.5.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/
-License: BSD-3
-Author: ygidtu
-Author-email: ygidtu@gmail.com
-Requires-Python: >=3.10,<3.13
-Classifier: License :: Other/Proprietary License
-Classifier: Programming Language :: Python :: 3
-Classifier: Programming Language :: Python :: 3.10
-Classifier: Programming Language :: Python :: 3.11
-Classifier: Programming Language :: Python :: 3.12
-Requires-Dist: adjusttext (>=0.7.3,<0.8.0)
-Requires-Dist: cairocffi (>=1.4.0,<2.0.0)
-Requires-Dist: click (>=8.1.3,<9.0.0)
-Requires-Dist: click-option-group (>=0.5.5,<0.6.0)
-Requires-Dist: filetype (>=1.2.0,<2.0.0)
-Requires-Dist: flask (>=2.3.2,<3.0.0)
-Requires-Dist: hicmatrix (>=17,<18)
-Requires-Dist: loguru (>=0.6.0,<0.7.0)
-Requires-Dist: matplotlib (>=3.6.3,<4.0.0)
-Requires-Dist: numpy (>=1.24.1,<2.0.0)
-Requires-Dist: pandas (>=1.5.3,<2.0.0)
-Requires-Dist: pybigwig (>=0.3.18,<0.4.0)
-Requires-Dist: pysam (>=0.21.0,<0.22.0)
-Requires-Dist: requests (>=2.28.2,<3.0.0)
-Requires-Dist: scipy (>=1.10.0,<2.0.0)
-Requires-Dist: seaborn (>=0.12.2,<0.13.0)
-Requires-Dist: xmltodict (>=0.13.0,<0.14.0)
+Author-email: ygidtu <ygidtu@gmail.com>
+License-Expression: BSD-3-Clause
+Requires-Python: >=3.11
 Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: adjusttext>=1.3.0
+Requires-Dist: cairocffi>=1.7.1
+Requires-Dist: click>=8.2.1
+Requires-Dist: click-option-group>=0.5.7
+Requires-Dist: filetype>=1.2.0
+Requires-Dist: flask>=3.1.1
+Requires-Dist: hicmatrix>=17.2
+Requires-Dist: loguru>=0.7.3
+Requires-Dist: matplotlib>=3.10.3
+Requires-Dist: numpy>=2.3.0
+Requires-Dist: pandas>=2.3.0
+Requires-Dist: pybigwig>=0.3.24
+Requires-Dist: pysam>=0.23.3
+Requires-Dist: requests>=2.32.4
+Requires-Dist: scipy>=1.15.3
+Requires-Dist: seaborn>=0.13.2
+Requires-Dist: xmltodict>=0.14.2
+Dynamic: license-file
 
 # trackplot
 
@@ -298,4 +294,3 @@ contact [Yiming Zhang](https://github.com/ygidtu) or
 If you use the tool in your publication, please cite by
 
 [Zhang Y, Zhou R, Liu L, et al. Trackplot: A flexible toolkit for combinatorial analysis of genomic data[J]. PLoS computational biology, 2023, 19(9): e1011477.](https://journals.plos.org/ploscompbiol/article?id=10.1371/journal.pcbi.1011477)
-

trackplot-0.5.4/pyproject.toml

@@ -0,0 +1,48 @@
+[project]
+name = "trackplot"
+version = "0.5.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" }
+]
+license = "BSD-3-Clause"
+readme = "README.md"
+requires-python = ">=3.11"
+dependencies = [
+    "adjusttext>=1.3.0",
+    "cairocffi>=1.7.1",
+    "click>=8.2.1",
+    "click-option-group>=0.5.7",
+    "filetype>=1.2.0",
+    "flask>=3.1.1",
+    "hicmatrix>=17.2",
+    "loguru>=0.7.3",
+    "matplotlib>=3.10.3",
+    "numpy>=2.3.0",
+    "pandas>=2.3.0",
+    "pybigwig>=0.3.24",
+    "pysam>=0.23.3",
+    "requests>=2.32.4",
+    "scipy>=1.15.3",
+    "seaborn>=0.13.2",
+    "xmltodict>=0.14.2",
+]
+
+
+[build-system]
+requires = ["setuptools>=65.0", "wheel"]
+build-backend = "setuptools.build_meta"
+
+
+[tool.setuptools.packages.find]
+# 搜索目前目录下的trackplot和ui及其子文件夹
+where = ["."]
+include = ["trackplot*", "ui*"]  # 通配符匹配所有子包
+
+
+[tool.setuptools.package-data]
+# 不限定父目录(等号前为父目录)，包含其下匹配特定正则的文件
+"*" = [
+    "pyproject.toml", "*"
+]
+

trackplot-0.5.4/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+

{trackplot-0.5.2 → trackplot-0.5.4}/trackplot/base/Junction.py

@@ -134,6 +134,19 @@ class Junction(object):
 
         return self.start > other.end
 
+    def str(self, with_strand: bool = True) -> str:
+        u"""convert junctions to string, with or without strand
+        
+        Keyword arguments:
+        :param with_strand: as name says
+        Return: chr:strat-end:strand or chr:start-end
+        """
+        if with_strand:
+            return str(self)
+        return f"{self.chromosome}:{self.start}-{self.end}"
+        
+        
+
 
 if __name__ == '__main__':
     pass

{trackplot-0.5.2 → trackplot-0.5.4}/trackplot/file/Annotation.py

@@ -2,8 +2,10 @@
 # -*- coding:utf-8 -*-
 u"""
 Created by ygidtu@gmail.com at 2020.05.07
-
 This scripts contains the class handle the reference file
+
+Modified by AD 2025/01/12 to include only specified transcripts, otherwise their exon coordinates still occupy the alignments.
+
 """
 import glob
 import gzip
@@ -349,16 +351,15 @@ class Annotation(File):
 
         return output_gtf
 
-    def __load_gtf__(self):
-
+    def __load_gtf__(self, transcripts_to_show: list[str]|None = None):
         u"""
         Load transcripts inside of region from gtf file
         :param region: target region
         :return: list of Transcript
-        """
+        """ 
+        # AD - passing transcripts_to_show
         transcripts = {}
         exons = {}
-
         for rec in Reader.read_gtf(self.path, self.region):
             start = max(rec.start, self.region.start)
             end = min(rec.end, self.region.end)
@@ -369,6 +370,13 @@ class Annotation(File):
                 break
 
             if re.search(r"(rna|transcript|cds)", rec.feature, re.I):
+
+                if transcripts_to_show:
+                    _name = rec.transcript_name if "transcript_name" in rec.attributes else rec.transcript_id
+                    if _name not in transcripts_to_show:
+                        logger.info(f"Skipping transcript {_name}")
+                        continue
+
                 if rec.transcript_id not in transcripts.keys():
                     transcripts[rec.transcript_id] = Transcript(
                         chromosome=rec.contig,
@@ -380,7 +388,8 @@ class Annotation(File):
                         gene=rec.gene_name if "gene_name" in rec.attributes else "",
                         transcript=rec.transcript_name if "transcript_name" in rec.attributes else "",
                         exons=[]
-                    )
+                 )
+
             elif re.search(r"(exon)", rec.feature, re.I):
                 if rec.transcript_id not in exons.keys():
                     exons[rec.transcript_id] = []
@@ -581,9 +590,16 @@ class Annotation(File):
         assert isinstance(region, GenomicLoci), "region should be a GenomicLoci object"
         if transcripts is None:
             transcripts = []
+        elif isinstance(transcripts, str):
+            # AD when no transcripts are specified, transcripts is an empty string, not None
+            if transcripts == "":
+                transcripts = []
+            else:
+                transcripts = transcripts.split(",")
         self.region = region
+        # AD - pass specified transcripts early
         if self.category == "gtf":
-            self.__load_gtf__()
+            self.__load_gtf__(transcripts)
             if self.add_local_domain:
                 self.__load_local_domain__(region)
 

{trackplot-0.5.2 → trackplot-0.5.4}/trackplot/file/Bam.py

@@ -20,7 +20,7 @@ from trackplot.base.Junction import Junction
 from trackplot.base.ReadDepth import ReadDepth
 from trackplot.base.Readder import Reader
 from trackplot.conf.config import NORMALIZATION
-from trackplot.file.File import SingleCell
+from trackplot.file.File import SingleCell  
 
 
 class Bam(SingleCell):
@@ -275,14 +275,32 @@ class Bam(SingleCell):
                     site_minus[start - region.start] += 1
 
             for k, v in spanned_junctions.items():
-                if included_junctions and str(k) not in included_junctions:
+                kept = v >= threshold
+                
+                # if the number of junctiosn is lower than threshold, then skip
+                if not kept:
                     continue
-
-                if v >= threshold:
+                
+                # if included_junctions is provided, then skip all junctions by default
+                if included_junctions:
+                    kept = False
+                
+                # check whether junctions should be kept
+                if k.str(with_strand = True) in included_junctions:
+                    logger.debug(f"{str(k)} is included")
+                    kept = True
+                elif k.str(with_strand = False) in included_junctions:
+                    logger.debug(f"{str(k)} is included, but strand is ignored")
+                    kept = True
+
+                if not kept:
+                    logger.debug(f"{str(k)} is not included")
+                else:
                     if k.strand == "+":
                         spanned_junctions_plus[k] = 1 + spanned_junctions_plus.get(k, v)
                     elif k.strand == "-":
                         spanned_junctions_minus[k] = -1 + spanned_junctions_minus.get(k, v)
+                        
         except IOError as err:
             logger.error('There is no .bam file at {0}'.format(self.path))
             logger.error(err)

{trackplot-0.5.2 → trackplot-0.5.4}/trackplot/file/ReadSegments.py

@@ -382,7 +382,9 @@ class ReadSegment(File):
                     current_ignore_num = min(
                         [self.deletion_ignore, read.query_alignment_length * self.del_ratio_ignore])
                 else:
-                    current_ignore_num = np.Inf
+                    current_ignore_num = np.inf
+                    # AD - AttributeError: `np.Inf` was removed in the NumPy 2.0 release. Use `np.inf` instead.
+
 
                 exon_bound = []
                 intron_bound = []

{trackplot-0.5.2 → trackplot-0.5.4}/trackplot/plot.py

@@ -30,7 +30,7 @@ logging.getLogger('matplotlib.font_manager').setLevel(logging.ERROR)
 faulthandler.enable()
 
 
-__version__ = "0.5.2"
+__version__ = "0.5.4"
 __author__ = "ygidtu & Ran Zhou"
 __email__ = "ygidtu@gmail.com"
 
@@ -1437,7 +1437,7 @@ class Plot(object):
 
         if output:
             logger.info(f"saving fig into {output}")
-            fig.savefig(output, transparent=True, bbox_inches='tight')
+            fig.savefig(output, transparent=False, bbox_inches='tight') # AD
         elif return_image:
             output = io.BytesIO()
             if return_image == "png":

{trackplot-0.5.2 → trackplot-0.5.4}/trackplot/plot_func.py

@@ -2,11 +2,13 @@
 # -*- coding:utf-8 -*-
 u"""
 This script contains the functions to draw different images
+Modified by AD 2025/01/13
 """
-import gzip
+import itertools # AD - for cumulative sum of graph_coords
+#import gzip
 import math
 from copy import deepcopy
-from decimal import Decimal
+#from decimal import Decimal
 from typing import Dict, List, Optional, Union, Tuple, Set
 
 import matplotlib as mpl
@@ -25,6 +27,7 @@ from matplotlib.textpath import TextPath
 from matplotlib.transforms import Affine2D
 from scipy.cluster.hierarchy import dendrogram, linkage
 from scipy.stats import gaussian_kde, zscore
+from scipy.stats import beta # AD
 
 from trackplot.anno.theme import Theme
 from trackplot.base.GenomicLoci import GenomicLoci
@@ -100,7 +103,6 @@ def cubic_bezier(pts, t):
     p3 = np.array(p3)
     return p0 * (1 - t) ** 3 + 3 * t * p1 * (1 - t) ** 2 + 3 * t ** 2 * (1 - t) * p2 + t ** 3 * p3
 
-
 def __merge_exons__(exons: List[List[int]]):
     u"""
     merge the overlap exons into one
@@ -124,7 +126,7 @@ def __merge_exons__(exons: List[List[int]]):
 
 
 def init_graph_coords(region: GenomicLoci, exons: Optional[List[List[int]]] = None, exon_scale=1,
-                      intron_scale=.5) -> np.array:
+                      intron_scale=0.5) -> np.array:
     u"""
     init the default
     :param region: the plot region
@@ -135,33 +137,59 @@ def init_graph_coords(region: GenomicLoci, exons: Optional[List[List[int]]] = No
     graph_coords = np.zeros(len(region), dtype=int)
 
     if exons:
-        # if there have exons, init graph_coords by exon and intron scales
-        for i in range(0, exons[0][0] - region.start):
-            graph_coords[i] = (i - 0) * intron_scale
-
-        exons = __merge_exons__(exons)
-        for i in range(0, len(exons)):
-            exon = exons[i]
-
-            if i > 0:
-                intron = [exons[i - 1][1], exons[i][0]]
-
-                for j in range(intron[0], intron[1]):
+        # AD - if user specifies an intron scale in line with Trackplot, addition is transparent
+        if intron_scale <= 1:
+            # if there have exons, init graph_coords by exon and intron scales
+            for i in range(0, exons[0][0] - region.start):
+                graph_coords[i] = (i - 0) * intron_scale
+            exons = __merge_exons__(exons)
+            for i in range(0, len(exons)):
+                exon = exons[i]
+                if i > 0:
+                    intron = [exons[i - 1][1], exons[i][0]]
+
+                    for j in range(intron[0], intron[1]):
+                        if j >= region.start:
+                            graph_coords[j - region.start] = graph_coords[intron[0] - region.start - 1] + (
+                                    j - intron[0] + 1) * intron_scale
+                for j in range(exon[0], exon[1] + 1):
                     if j >= region.start:
-                        graph_coords[j - region.start] = graph_coords[intron[0] - region.start - 1] + (
-                                j - intron[0] + 1) * intron_scale
-
-            for j in range(exon[0], exon[1] + 1):
-                if j >= region.start:
-                    graph_coords[j - region.start] = graph_coords[exon[0] - region.start - 1] + (
-                            j - exon[0] + 1) * exon_scale
-
-        intron = [exons[-1][-1], region.end]
-        for i in range(intron[0], intron[1]):
-            if i >= region.start:
-                graph_coords[i - region.start] = graph_coords[intron[0] - region.start - 1] + (
-                        i - intron[0] + 1) * intron_scale
-
+                        graph_coords[j - region.start] = graph_coords[exon[0] - region.start - 1] + (
+                                j - exon[0] + 1) * exon_scale
+            intron = [exons[-1][-1], region.end]
+            for i in range(intron[0], intron[1]):
+                if i >= region.start:
+                    graph_coords[i - region.start] = graph_coords[intron[0] - region.start - 1] + (
+                            i - intron[0] + 1) * intron_scale
+        else:
+            exons = __merge_exons__(exons)
+            while exons[0][1] < region.start:
+                exons = exons[1:]
+            while exons[-1][0] > region.end:
+                exons = exons[:-1]
+            exons[0][0] = max(exons[0][0], region.start)
+            exons[-1][1] = min(exons[-1][1], region.end)
+            last_interval = region.end - exons[-1][1]
+            for i, e in enumerate(exons):
+                exons[i][0] -= region.start
+                exons[i][1] -= region.start
+
+            steps = [float(exon_scale)] * len(region)
+            if exons[0][0] > 0:
+                step = intron_scale / exons[0][0]
+                for i in range(exons[0][0]):
+                    steps[i] = step
+            for e in range(1, len(exons)):
+                interval = exons[e][0] - exons[e-1][1] - 2
+                step = intron_scale / interval
+                for i in range(exons[e-1][1], exons[e][0]):
+                    steps[i] = step
+                interval = exons[e][0] - exons[e-1][1] - 2
+            if last_interval := len(region) - exons[-1][1] - 1:
+                step = intron_scale / last_interval
+                for i in range(exons[1][1] +1, len(region)):
+                    steps[i] = step
+            graph_coords = list(map(int, itertools.accumulate(steps)))
     else:
         # if there is not any exons, just init graph_coords by region
         for i, j in enumerate(range(region.start, region.end + 1)):
@@ -175,7 +203,6 @@ def init_graph_coords(region: GenomicLoci, exons: Optional[List[List[int]]] = No
                 graph_coords[i] = max(graph_coords) + 1
     return graph_coords
 
-
 def set_x_ticks(
         ax: mpl.axes.Axes,
         region: GenomicLoci,
@@ -284,6 +311,8 @@ def set_y_ticks(
         universal_y_ticks = sorted(universal_y_ticks)
 
         for lab in universal_y_ticks:
+            curr_y_tick_labels.append(f"{int(lab)}")
+            """
             if y_axis_skip_zero and lab == 0:
                 # Exclude label for 0
                 curr_y_tick_labels.append("")
@@ -291,6 +320,7 @@ def set_y_ticks(
                 curr_y_tick_labels.append(f"{Decimal(lab):.2E}")
             else:
                 curr_y_tick_labels.append(f"{lab:.1f}" if lab % 1 != 0 else f"{int(lab)}")
+            """
         u"""
         @2019.01.04
         If there is no bam file, draw a blank y-axis 
@@ -306,6 +336,7 @@ def set_y_ticks(
     """
 
     if show_y_label:
+        
         def __dynamic_distance__(distance_between_label_axis: float, label: str, scale: int = 100) -> float:
             if distance_between_label_axis != 0:
                 return -distance_between_label_axis
@@ -314,9 +345,7 @@ def set_y_ticks(
 
         curr_y_tick_labels = sorted(curr_y_tick_labels, key=lambda x: len(x), reverse=True)
         ax.text(
-            x=__dynamic_distance__(
-                distance_between_label_axis, curr_y_tick_labels[0] if curr_y_tick_labels else ""
-            ) * max(graph_coords),
+            x=__dynamic_distance__(distance_between_label_axis, curr_y_tick_labels[0] if curr_y_tick_labels else "") * max(graph_coords),
             y=(max_used_y_val + min_used_y_val) / 2,
             s=label, fontsize=font_size, ha="right"
         )
@@ -629,7 +658,7 @@ def plot_density(
         color="blue",
         font_size: int = 8,
         show_junction_number: bool = True,
-        junction_number_font_size: int = 5,
+        junction_number_font_size: int = 12,
         n_y_ticks: int = 4,
         distance_between_label_axis: float = .1,
         show_y_label: bool = True,
@@ -664,6 +693,8 @@ def plot_density(
     :param kwargs:
     :return:
     """
+    # max_used_y_val is None
+    # distance_between_label_axis: 0
     if obj:
         assert obj.region is not None, "please load data first"
         region = obj.region
@@ -684,7 +715,17 @@ def plot_density(
     except AttributeError:
         # depth do not have junctions
         jxns = {}
+    
+    # AD - convert junction counts to log scale early if requested
+    if obj.log_trans and obj.log_trans.isdigit() and int(obj.log_trans) > 0: # in ["2", "10"]:
+
+        y_label += f" (log{obj.log_trans})"
+        denominator = np.log(int(obj.log_trans))
+        for k, v in jxns.items():
+            jxns[k] = np.log1p(v)  / denominator
 
+    min_used_y_val = 0 # AD
+    fixed_min_used_y = True #  AD
     fixed_max_used_y, fixed_min_used_y = max_used_y_val is not None, min_used_y_val is not None
     if max_used_y_val is None:
         if isinstance(data, dict):
@@ -720,7 +761,7 @@ def plot_density(
         if data.strand_aware:
             max_used_y_val = max(abs(min_used_y_val), max_used_y_val)
             min_used_y_val = -max(abs(min_used_y_val), max_used_y_val) if data.minus is not None else 0
-
+            
     if jxns:
         # sort the junctions by intron length for better plotting look
         jxns_sorted_list = sorted(jxns.keys(), key=lambda x: (x.end - x.start, x.start, x.end), reverse=True)
@@ -732,10 +773,8 @@ def plot_density(
             min_junction_count = min(jxns.values())
         junction_count_gap = max_junction_count - min_junction_count
 
-        recorded_pts = set()
         jxn_numbers = []
         for jxn_idx, jxn in enumerate(jxns_sorted_list):
-            logger.debug(f"junctions of {y_label}: {jxn} - {round(jxns[jxn], 2)}")
             leftss, rightss = jxn.start, jxn.end
 
             # junction must at least have one anchor located in plotted region
@@ -753,15 +792,18 @@ def plot_density(
             add two new variables to make it clear which one is index, which one is genomic site 
             """
             ss1, ss2 = graph_coords[ss1_idx], graph_coords[ss2_idx]
+            # AD = keep junction arcs on top
+            jxn_on_top = True
 
-            # draw junction on bottom
+            # draw junction on bottom 
             if kwargs.get("density_by_strand"):
                 jxn_on_top = jxn.strand == "+"
             else:
                 jxn_on_top = jxn_idx % 2 == 0
+                #jxn_on_top = True # AD - keep all junctions on same strand
                 if abs(min_used_y_val) < max_used_y_val:
                     min_used_y_val = -max_used_y_val
-
+            
             if fill_step == "post":
                 ss1_idx = max(0, ss1_idx - 1)
             elif fill_step == "pre":
@@ -785,27 +827,22 @@ def plot_density(
                     -right_dens - current_height,
                     -right_dens if not ss2_modified else -right_dens - current_height
                 ]
-
-            if sum(pts) > 0:
-                pts_ = "_".join([str(x) for x in pts])
-                while pts_ in recorded_pts:
-                    pts[1], pts[2] = pts[1] * 1.1, pts[2] * 1.1
-                    pts_ = "_".join([str(x) for x in pts])
-                recorded_pts.add(pts_)
-
+  
             """
             @2018.12.26
             scale the junctions line width
             """
             if junction_count_gap > 0:
-                line_width = (jxns[jxn] - min_junction_count) / junction_count_gap
+                #line_width = (jxns[jxn] - min_junction_count) / junction_count_gap
+                line_width = .5 + 1.5 * (jxns[jxn] - min_junction_count) / junction_count_gap
             else:
                 line_width = 0
+            #line_width = max(.5,np.log())
 
             pts = [(ss1, pts[0]), (ss1, pts[1]), (ss2, pts[2]), (ss2, pts[3])]
             ax.add_patch(PathPatch(Path(pts, [Path.MOVETO, Path.CURVE4, Path.CURVE4, Path.CURVE4]),
                                    ec=color, lw=line_width + 0.2, fc='none'))
-
+            
             if show_junction_number:
                 midpt = cubic_bezier(pts, .5)
 
@@ -821,10 +858,6 @@ def plot_density(
                 t.set_bbox(dict(alpha=0))
                 jxn_numbers.append(t)
 
-        try:
-            adjust_text(jxn_numbers, force_text=(0.2, 0.2), arrowprops=dict(arrowstyle="-", color='black', lw=1))
-        except IndexError as err:
-            logger.debug(err)
 
     if obj and obj.title:
         ax.text(max(graph_coords) - len(obj.title), max_used_y_val, obj.title, color=color, fontsize=font_size)
@@ -849,13 +882,13 @@ def plot_density(
     set_y_ticks(
         ax, label=y_label, theme=theme,
         graph_coords=graph_coords,
-        max_used_y_val=max_used_y_val * 1.1,  # keep a buffer at top and bottom of junctions
-        min_used_y_val=min_used_y_val * 1.1,
+        max_used_y_val=max_used_y_val, # * 1.1,  # keep a buffer at top and bottom of junctions
+        min_used_y_val=min_used_y_val, # * 1.1, # AD
         n_y_ticks=n_y_ticks,
         distance_between_label_axis=distance_between_label_axis,
         font_size=font_size,
         show_y_label=show_y_label,
-        y_axis_skip_zero=False if not isinstance(data, dict) and data.strand_aware else True
+        y_axis_skip_zero=False #if not isinstance(data, dict) and data.strand_aware else True
     )