cool-seq-tool 0.14.1__tar.gz → 0.14.3__tar.gz

{cool_seq_tool-0.14.1 → cool_seq_tool-0.14.3}/.gitignore

@@ -135,3 +135,6 @@ docs/source/reference/api
 
 # Files created by doctests
 tpm3.fasta
+
+# VSCode
+.vscode

{cool_seq_tool-0.14.1 → cool_seq_tool-0.14.3}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: cool_seq_tool
-Version: 0.14.1
+Version: 0.14.3
 Summary: Common Operation on Lots of Sequences Tool
 Author: Kori Kuzma, James Stevenson, Katie Stahl, Alex Wagner
 License: MIT License
@@ -38,22 +38,20 @@ Classifier: Intended Audience :: Developers
 Classifier: Topic :: Scientific/Engineering :: Bio-Informatics
 Classifier: License :: OSI Approved :: MIT 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-Python: >=3.10
+Classifier: Programming Language :: Python :: 3.13
+Requires-Python: >=3.11
 Description-Content-Type: text/markdown
 License-File: LICENSE
 Requires-Dist: asyncpg
-Requires-Dist: aiofiles
 Requires-Dist: boto3
 Requires-Dist: agct>=0.1.0-dev1
 Requires-Dist: polars~=1.0
-Requires-Dist: hgvs
 Requires-Dist: biocommons.seqrepo
 Requires-Dist: pydantic<3.0,>=2.0
 Requires-Dist: ga4gh.vrs<3.0,>=2.1.3
-Requires-Dist: wags-tails~=0.3.2
+Requires-Dist: wags-tails~=0.4.0
 Requires-Dist: bioutils
 Provides-Extra: dev
 Requires-Dist: pre-commit>=4.2.0; extra == "dev"
@@ -64,7 +62,7 @@ Requires-Dist: ruff==0.12.1; extra == "dev"
 Provides-Extra: tests
 Requires-Dist: pytest; extra == "tests"
 Requires-Dist: pytest-cov; extra == "tests"
-Requires-Dist: pytest-asyncio==0.18.3; extra == "tests"
+Requires-Dist: pytest-asyncio; extra == "tests"
 Requires-Dist: mock; extra == "tests"
 Provides-Extra: docs
 Requires-Dist: sphinx==6.1.3; extra == "docs"

{cool_seq_tool-0.14.1 → cool_seq_tool-0.14.3}/docs/source/install.rst

@@ -6,7 +6,7 @@ Installation
 Prerequisites
 -------------
 
-* Python 3.10+
+* Python 3.11+
 * A UNIX-like environment (e.g. MacOS, WSL, Ubuntu)
 * A recent version of PostgreSQL (ideally at least 11+)
 

{cool_seq_tool-0.14.1 → cool_seq_tool-0.14.3}/docs/source/reference/index.rst

@@ -55,6 +55,7 @@ Data Mappers
 
    cool_seq_tool.mappers.alignment
    cool_seq_tool.mappers.exon_genomic_coords
+   cool_seq_tool.mappers.feature_overlap
    cool_seq_tool.mappers.liftover
    cool_seq_tool.mappers.mane_transcript
 

{cool_seq_tool-0.14.1 → cool_seq_tool-0.14.3}/pyproject.toml

@@ -16,24 +16,22 @@ classifiers = [
     "Topic :: Scientific/Engineering :: Bio-Informatics",
     "License :: OSI Approved :: MIT License",
     "Programming Language :: Python :: 3",
-    "Programming Language :: Python :: 3.10",
     "Programming Language :: Python :: 3.11",
     "Programming Language :: Python :: 3.12",
+    "Programming Language :: Python :: 3.13",
 ]
-requires-python = ">=3.10"
+requires-python = ">=3.11"
 description = "Common Operation on Lots of Sequences Tool"
 license = {file = "LICENSE"}
 dependencies = [
     "asyncpg",
-    "aiofiles",
     "boto3",
     "agct >= 0.1.0-dev1",
     "polars ~= 1.0",
-    "hgvs",
     "biocommons.seqrepo",
     "pydantic >=2.0,<3.0",
     "ga4gh.vrs >=2.1.3,<3.0",
-    "wags-tails ~= 0.3.2",
+    "wags-tails ~= 0.4.0",
     "bioutils",
 ]
 dynamic = ["version"]
@@ -46,7 +44,7 @@ dev = [
     "psycopg2-binary",
     "ruff==0.12.1",
 ]
-tests = ["pytest", "pytest-cov", "pytest-asyncio==0.18.3", "mock"]
+tests = ["pytest", "pytest-cov", "pytest-asyncio", "mock"]
 docs = [
     "sphinx==6.1.3",
     "sphinx-autodoc-typehints==1.22.0",
@@ -81,6 +79,9 @@ build-backend = "setuptools.build_meta"
 "cool_seq_tool.resources" = ["transcript_mapping.tsv"]
 
 [tool.pytest.ini_options]
+asyncio_mode = "auto"
+asyncio_default_fixture_loop_scope = "session"
+asyncio_default_test_loop_scope = "session"
 addopts = "--cov=src --cov-report term-missing"
 testpaths = ["tests"]
 

{cool_seq_tool-0.14.1 → cool_seq_tool-0.14.3}/src/cool_seq_tool/mappers/__init__.py

@@ -4,6 +4,13 @@ from .alignment import AlignmentMapper  # noqa: I001
 from .liftover import LiftOver
 from .mane_transcript import ManeTranscript
 from .exon_genomic_coords import ExonGenomicCoordsMapper
+from .feature_overlap import FeatureOverlap
 
 
-__all__ = ["AlignmentMapper", "ExonGenomicCoordsMapper", "LiftOver", "ManeTranscript"]
+__all__ = [
+    "AlignmentMapper",
+    "ExonGenomicCoordsMapper",
+    "FeatureOverlap",
+    "LiftOver",
+    "ManeTranscript",
+]

{cool_seq_tool-0.14.1 → cool_seq_tool-0.14.3}/src/cool_seq_tool/mappers/alignment.py

@@ -106,7 +106,7 @@ class AlignmentMapper:
         c_end_pos: int,
         cds_start: int | None = None,
         coordinate_type: CoordinateType = CoordinateType.RESIDUE,
-        target_genome_assembly: bool = Assembly.GRCH38,
+        target_genome_assembly: Assembly = Assembly.GRCH38,
     ) -> tuple[dict | None, str | None]:
         """Translate cDNA representation to genomic representation
 

cool_seq_tool-0.14.3/src/cool_seq_tool/mappers/feature_overlap.py

@@ -0,0 +1,252 @@
+"""Module for getting feature (gene/exon) overlap"""
+
+import re
+from pathlib import Path
+
+import polars as pl
+from ga4gh.core import ga4gh_identify
+from ga4gh.vrs.models import SequenceLocation, SequenceReference
+
+from cool_seq_tool.handlers import SeqRepoAccess
+from cool_seq_tool.resources.data_files import DataFile, get_data_file
+from cool_seq_tool.schemas import Assembly, CdsOverlap, CoordinateType
+
+# Pattern for chromosome
+CHR_PATTERN = r"X|Y|([1-9]|1[0-9]|2[0-2])"
+
+
+class FeatureOverlapError(Exception):
+    """Custom exception for the Feature Overlap class"""
+
+
+class FeatureOverlap:
+    """The class for getting feature overlap"""
+
+    def __init__(
+        self,
+        seqrepo_access: SeqRepoAccess,
+        mane_refseq_genomic_path: Path | None = None,
+        from_local: bool = False,
+    ) -> None:
+        """Initialize the FeatureOverlap class. Will load RefSeq data and store as df.
+
+        :param seqrepo_access: Client for accessing SeqRepo data
+        :param mane_refseq_genomic_path: Path to MANE RefSeq Genomic GFF data
+        :param from_local: if ``True``, don't check for or acquire latest version --
+            just provide most recent locally available file, if possible, and raise
+            error otherwise
+        """
+        if not mane_refseq_genomic_path:
+            mane_refseq_genomic_path = get_data_file(
+                DataFile.MANE_REFSEQ_GENOMIC, from_local
+            )
+        self.seqrepo_access = seqrepo_access
+        self.mane_refseq_genomic_path = mane_refseq_genomic_path
+        self.df = self._load_mane_refseq_gff_data()
+
+    def _load_mane_refseq_gff_data(self) -> pl.DataFrame:
+        """Load MANE RefSeq GFF data file into DataFrame.
+
+        :return: DataFrame containing MANE RefSeq Genomic GFF data for CDS. Columns
+            include `type`, `chromosome` (chromosome without 'chr' prefix), `cds_start`,
+            `cds_stop`, `info_name` (name of record), and `gene`. `cds_start` and
+            `cds_stop` use inter-residue coordinates.
+        """
+        df = pl.read_csv(
+            self.mane_refseq_genomic_path,
+            separator="\t",
+            has_header=False,
+            skip_rows=9,
+            columns=[0, 2, 3, 4, 8],
+        )
+        df.columns = ["chromosome", "type", "cds_start", "cds_stop", "info"]
+
+        # Restrict to only feature of interest: CDS (which has gene info)
+        df = df.filter(pl.col("type") == "CDS")
+
+        # Get name from the info field
+        # Get gene from the info field
+        # Get chromosome names without prefix and without suffix for alternate transcripts
+        # Convert start and stop to ints
+        # Convert to inter-residue coordinates
+        # Only return certain columns
+        return df.with_columns(
+            (pl.col("info").str.extract(r"Name=([^;]+)", 1).alias("info_name")),
+            (pl.col("info").str.extract(r"gene=([^;]+)", 1).alias("gene")),
+            (pl.col("chromosome").str.extract(r"^chr?([^_]+)", 1).alias("chromosome")),
+            (pl.col("cds_start").cast(pl.Int64) - 1).alias("cds_start"),
+            (pl.col("cds_stop").cast(pl.Int64).alias("cds_stop")),
+        ).select(
+            [
+                pl.col("type"),
+                pl.col("chromosome"),
+                pl.col("cds_start"),
+                pl.col("cds_stop"),
+                pl.col("info_name"),
+                pl.col("gene"),
+            ]
+        )
+
+    def _get_chr_from_alt_ac(self, identifier: str) -> str:
+        """Get chromosome given genomic identifier
+
+        :param identifier: Genomic identifier on GRCh38 assembly
+        :raises FeatureOverlapError: If unable to find associated GRCh38 chromosome
+        :return: Chromosome. 1..22, X, Y. No 'chr' prefix.
+        """
+        aliases, error_msg = self.seqrepo_access.translate_identifier(
+            identifier, Assembly.GRCH38.value
+        )
+
+        if error_msg:
+            raise FeatureOverlapError(str(error_msg))
+
+        if not aliases:
+            error_msg = (
+                f"Unable to find {Assembly.GRCH38.value} aliases for: {identifier}"
+            )
+            raise FeatureOverlapError(error_msg)
+
+        assembly_chr_pattern = (
+            rf"^{Assembly.GRCH38.value}:(?P<chromosome>{CHR_PATTERN})$"
+        )
+        for a in aliases:
+            chr_match = re.match(assembly_chr_pattern, a)
+            if chr_match:
+                break
+
+        if not chr_match:
+            error_msg = (
+                f"Unable to find {Assembly.GRCH38.value} chromosome for: {identifier}"
+            )
+            raise FeatureOverlapError(error_msg)
+
+        chr_groupdict = chr_match.groupdict()
+        return chr_groupdict["chromosome"]
+
+    def get_grch38_mane_gene_cds_overlap(
+        self,
+        start: int,
+        end: int,
+        chromosome: str | None = None,
+        identifier: str | None = None,
+        coordinate_type: CoordinateType = CoordinateType.RESIDUE,
+    ) -> dict[str, list[CdsOverlap]] | None:
+        """Given GRCh38 genomic data, find the overlapping MANE features (gene and cds).
+        The genomic data is specified as a sequence location by `chromosome`, `start`,
+        `end`. All CDS regions with which the input sequence location has nonzero base
+        pair overlap will be returned.
+
+        :param start: GRCh38 start position
+        :param end: GRCh38 end position
+        :param chromosome: Chromosome. 1..22, X, or Y. If not provided, must provide
+            `identifier`. If both `chromosome` and `identifier` are provided,
+            `chromosome` will be used.
+        :param identifier: Genomic identifier on GRCh38 assembly. If not provided, must
+            provide `chromosome`. If both `chromosome` and `identifier` are provided,
+            `chromosome` will be used.
+        :param coordinate_type: Coordinate type for ``start`` and ``end``
+        :raise FeatureOverlapError: If missing required fields or unable to find
+            associated ga4gh identifier
+        :return: MANE feature (gene/cds) overlap data represented as a dict. The
+            dictionary will be keyed by genes which overlap the input sequence location.
+            Each gene contains a list of the overlapping CDS regions with the beginning
+            and end of the input sequence location's overlap with each
+        """
+        ga4gh_seq_id = None
+        if chromosome:
+            if not re.match(f"^{CHR_PATTERN}$", chromosome):
+                error_msg = "`chromosome` must be 1, ..., 22, X, or Y"
+                raise FeatureOverlapError(error_msg)
+        else:
+            if identifier:
+                chromosome = self._get_chr_from_alt_ac(identifier)
+                if identifier.startswith("ga4gh:SQ."):
+                    ga4gh_seq_id = identifier
+            else:
+                error_msg = "Must provide either `chromosome` or `identifier`"
+                raise FeatureOverlapError(error_msg)
+
+        # Convert residue to inter-residue
+        if coordinate_type == CoordinateType.RESIDUE:
+            start -= 1
+
+        # Get feature dataframe (df uses inter-residue)
+        feature_df = self.df.filter(
+            (pl.col("chromosome") == chromosome)
+            & (pl.col("cds_start") <= end)
+            & (pl.col("cds_stop") >= start)
+        )
+
+        if feature_df.is_empty():
+            return None
+
+        # Add overlap columns
+        feature_df = feature_df.with_columns(
+            [
+                pl.when(pl.col("cds_start") < start)
+                .then(start)
+                .otherwise(pl.col("cds_start"))
+                .alias("overlap_start"),
+                pl.when(pl.col("cds_stop") > end)
+                .then(end)
+                .otherwise(pl.col("cds_stop"))
+                .alias("overlap_stop"),
+            ]
+        )
+
+        # Get ga4gh identifier for chromosome
+        if not ga4gh_seq_id:
+            grch38_chr = f"{Assembly.GRCH38.value}:{chromosome}"
+            ga4gh_aliases, error_msg = self.seqrepo_access.translate_identifier(
+                grch38_chr, "ga4gh"
+            )
+
+            # Errors should never happen but catching just in case
+            if error_msg:
+                raise FeatureOverlapError(str(error_msg))
+
+            if not ga4gh_aliases:
+                error_msg = f"Unable to find ga4gh identifier for: {grch38_chr}"
+                raise FeatureOverlapError(error_msg)
+
+            ga4gh_seq_id = ga4gh_aliases[0]
+
+        def _get_seq_loc(
+            start_pos: int, stop_pos: int, refget_ac: str
+        ) -> SequenceLocation:
+            """Get VRS Sequence Location
+
+            :param start_pos: Start position
+            :param stop_pos: Stop position
+            :param refget_ac: Refget Accession (SQ.)
+            :return: VRS Sequence Location
+            """
+            _sl = SequenceLocation(
+                sequenceReference=SequenceReference(
+                    refgetAccession=refget_ac,
+                ),
+                start=start_pos,
+                end=stop_pos,
+            )
+            ga4gh_identify(_sl)
+            return _sl
+
+        resp = {}
+        refget_ac = ga4gh_seq_id.split("ga4gh:")[-1]
+        for gene, group in feature_df.group_by("gene"):
+            gene = gene[0]
+            _gene_overlap_data = [
+                CdsOverlap(
+                    cds=_get_seq_loc(
+                        cds_row["cds_start"], cds_row["cds_stop"], refget_ac
+                    ),
+                    overlap=_get_seq_loc(
+                        cds_row["overlap_start"], cds_row["overlap_stop"], refget_ac
+                    ),
+                ).model_dump(by_alias=True, exclude_none=True)
+                for cds_row in group.iter_rows(named=True)
+            ]
+            resp[gene] = _gene_overlap_data
+
+        return resp

{cool_seq_tool-0.14.1 → cool_seq_tool-0.14.3}/src/cool_seq_tool/mappers/mane_transcript.py

@@ -55,7 +55,7 @@ class DataRepresentation(BaseModel):
     """Define object model for final output representation"""
 
     gene: str | None = None
-    refseq: str
+    refseq: str | None = None
     ensembl: str | None = None
     pos: tuple[int, int]
     strand: Strand
@@ -447,7 +447,7 @@ class ManeTranscript:
 
     async def _g_to_c(
         self,
-        g: dict,
+        g: GenomicTxMetadata,
         refseq_c_ac: str,
         status: TranscriptPriority,
         ensembl_c_ac: str | None = None,
@@ -590,16 +590,23 @@ class ManeTranscript:
         if mane_transcript:
             mane_start_pos = mane_transcript.pos[0]
             mane_end_pos = mane_transcript.pos[1]
-            if anno == AnnotationLayer.CDNA:
+            if anno == AnnotationLayer.CDNA and isinstance(
+                mane_transcript, CdnaRepresentation
+            ):
                 mane_cds = mane_transcript.coding_start_site
                 mane_start_pos += mane_cds
                 mane_end_pos += mane_cds
-            mane_ref, _ = self.seqrepo_access.get_reference_sequence(
-                mane_transcript.refseq,
-                start=mane_start_pos,
-                end=mane_end_pos if mane_start_pos != mane_end_pos else None,
-                coordinate_type=coordinate_type,
-            )
+
+            if mane_transcript.refseq:
+                mane_ref, _ = self.seqrepo_access.get_reference_sequence(
+                    mane_transcript.refseq,
+                    start=mane_start_pos,
+                    end=mane_end_pos if mane_start_pos != mane_end_pos else None,
+                    coordinate_type=coordinate_type,
+                )
+            else:
+                mane_ref = None
+
             if not mane_ref:
                 _logger.info("Unable to validate reference for MANE Transcript")
 
@@ -1330,7 +1337,7 @@ class ManeTranscript:
         gene: str | None = None,
         coordinate_type: CoordinateType = CoordinateType.RESIDUE,
         try_longest_compatible: bool = False,
-    ) -> dict | None:
+    ) -> ProteinAndCdnaRepresentation | None:
         """Given GRCh38 genomic representation, return protein representation.
 
         Will try MANE Select and then MANE Plus Clinical. If neither is found and

{cool_seq_tool-0.14.1 → cool_seq_tool-0.14.3}/src/cool_seq_tool/resources/data_files.py

@@ -6,7 +6,11 @@ from importlib import resources
 from os import environ
 from pathlib import Path
 
-from wags_tails import NcbiLrgRefSeqGeneData, NcbiManeSummaryData
+from wags_tails import (
+    NcbiLrgRefSeqGeneData,
+    NcbiManeRefSeqGenomicData,
+    NcbiManeSummaryData,
+)
 
 _logger = logging.getLogger(__name__)
 
@@ -16,6 +20,7 @@ class DataFile(str, Enum):
 
     TRANSCRIPT_MAPPINGS = "transcript_mappings"
     MANE_SUMMARY = "mane_summary"
+    MANE_REFSEQ_GENOMIC = "mane_refseq_genomic"
     LRG_REFSEQGENE = "lrg_refseqgene"
 
     def lower(self) -> str:
@@ -37,6 +42,12 @@ _resource_acquisition_params = {
             from_local=from_local
         )[0],
     ),
+    DataFile.MANE_REFSEQ_GENOMIC: (
+        "MANE_REFSEQ_GENOMIC_PATH",
+        lambda from_local: NcbiManeRefSeqGenomicData(silent=True).get_latest(
+            from_local=from_local
+        )[0],
+    ),
     DataFile.LRG_REFSEQGENE: (
         "LRG_REFSEQGENE_PATH",
         lambda from_local: NcbiLrgRefSeqGeneData(silent=True).get_latest(
@@ -53,6 +64,7 @@ def get_data_file(resource: DataFile, from_local: bool = False) -> Path:
 
     * ``Resource.TRANSCRIPT_MAPPINGS`` -> ``TRANSCRIPT_MAPPINGS_PATH``
     * ``Resource.MANE_SUMMARY`` -> ``MANE_SUMMARY_PATH``
+    * ``Resource.MANE_REFSEQ_GENOMIC`` -> ``MANE_REFSEQ_GENOMIC_PATH``
     * ``Resource.LRG_REFSEQGENE`` -> ``LRG_REFSEQGENE_PATH``
 
     Otherwise, this function falls back on default expected locations:

{cool_seq_tool-0.14.1 → cool_seq_tool-0.14.3}/src/cool_seq_tool/resources/status.py

@@ -24,6 +24,7 @@ ResourceStatus = namedtuple(
         DataFile.TRANSCRIPT_MAPPINGS.lower(),
         DataFile.MANE_SUMMARY.lower(),
         DataFile.LRG_REFSEQGENE.lower(),
+        DataFile.MANE_REFSEQ_GENOMIC.lower(),
         "liftover",
     ),
 )
@@ -37,6 +38,7 @@ async def check_status(
     sr: SeqRepo | None = None,
     chain_file_37_to_38: str | None = None,
     chain_file_38_to_37: str | None = None,
+    mane_refseq_genomic_path: str | None = None,
 ) -> ResourceStatus:
     """Perform basic status checks on availability of required data resources.
 
@@ -62,6 +64,7 @@ async def check_status(
         is used for ``agct``. If this is not provided, will check to see if
         ``LIFTOVER_CHAIN_38_TO_37`` env var is set. If neither is provided, will allow
         ``agct`` to download a chain file from UCSC
+    :param mane_refseq_genomic_path: Optional path to MANE RefSeq Genomic GFF data
     :return: boolean description of availability of each resource, given current
         environment configurations
     """
@@ -69,19 +72,21 @@ async def check_status(
         DataFile.TRANSCRIPT_MAPPINGS.lower(): transcript_file_path,
         DataFile.LRG_REFSEQGENE.lower(): lrg_refseqgene_path,
         DataFile.MANE_SUMMARY.lower(): mane_data_path,
+        DataFile.MANE_REFSEQ_GENOMIC.lower(): mane_refseq_genomic_path,
     }
 
     status = {
         DataFile.TRANSCRIPT_MAPPINGS.lower(): False,
         DataFile.LRG_REFSEQGENE.lower(): False,
         DataFile.MANE_SUMMARY.lower(): False,
+        DataFile.MANE_REFSEQ_GENOMIC.lower(): False,
         "liftover": False,
         "uta": False,
         "seqrepo": False,
     }
     for r in list(DataFile):
         name_lower = r.lower()
-        declared_path = file_path_params[name_lower]
+        declared_path = file_path_params.get(name_lower)
         if declared_path and declared_path.exists() and declared_path.is_file():
             status[name_lower] = True
             continue

{cool_seq_tool-0.14.1 → cool_seq_tool-0.14.3}/src/cool_seq_tool/schemas.py

@@ -4,6 +4,7 @@ import datetime
 from enum import Enum, IntEnum
 from typing import Literal
 
+from ga4gh.vrs.models import SequenceLocation
 from pydantic import (
     BaseModel,
     ConfigDict,
@@ -13,7 +14,7 @@ from pydantic import (
 
 from cool_seq_tool import __version__
 
-_now = str(datetime.datetime.now(tz=datetime.timezone.utc))
+_now = str(datetime.datetime.now(tz=datetime.UTC))
 
 
 class AnnotationLayer(str, Enum):
@@ -167,3 +168,37 @@ class ServiceMeta(BaseModelForbidExtra):
             }
         }
     )
+
+
+class CdsOverlap(BaseModelForbidExtra):
+    """Create model for representing CDS start/stop and Overlap start/stop"""
+
+    cds: SequenceLocation
+    overlap: SequenceLocation
+
+    model_config = ConfigDict(
+        json_schema_extra={
+            "example": {
+                "cds": {
+                    "id": "ga4gh:SL.fYRYzNIAoe6UQF9MT1XaYsFscoU68ZJv",
+                    "type": "SequenceLocation",
+                    "sequenceReference": {
+                        "refgetAccession": "SQ.F-LrLMe1SRpfUZHkQmvkVKFEGaoDeHul",
+                        "type": "SequenceReference",
+                    },
+                    "start": 140726493,
+                    "end": 140726516,
+                },
+                "overlap": {
+                    "id": "ga4gh:SL.fYRYzNIAoe6UQF9MT1XaYsFscoU68ZJv",
+                    "type": "SequenceLocation",
+                    "sequenceReference": {
+                        "refgetAccession": "SQ.F-LrLMe1SRpfUZHkQmvkVKFEGaoDeHul",
+                        "type": "SequenceReference",
+                    },
+                    "start": 140726493,
+                    "end": 140726516,
+                },
+            }
+        }
+    )

{cool_seq_tool-0.14.1 → cool_seq_tool-0.14.3}/src/cool_seq_tool/utils.py

@@ -47,7 +47,7 @@ def service_meta() -> ServiceMeta:
     """
     return ServiceMeta(
         version=__version__,
-        response_datetime=datetime.datetime.now(tz=datetime.timezone.utc),
+        response_datetime=datetime.datetime.now(tz=datetime.UTC),
     )
 
 

{cool_seq_tool-0.14.1 → cool_seq_tool-0.14.3}/src/cool_seq_tool.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: cool_seq_tool
-Version: 0.14.1
+Version: 0.14.3
 Summary: Common Operation on Lots of Sequences Tool
 Author: Kori Kuzma, James Stevenson, Katie Stahl, Alex Wagner
 License: MIT License
@@ -38,22 +38,20 @@ Classifier: Intended Audience :: Developers
 Classifier: Topic :: Scientific/Engineering :: Bio-Informatics
 Classifier: License :: OSI Approved :: MIT 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-Python: >=3.10
+Classifier: Programming Language :: Python :: 3.13
+Requires-Python: >=3.11
 Description-Content-Type: text/markdown
 License-File: LICENSE
 Requires-Dist: asyncpg
-Requires-Dist: aiofiles
 Requires-Dist: boto3
 Requires-Dist: agct>=0.1.0-dev1
 Requires-Dist: polars~=1.0
-Requires-Dist: hgvs
 Requires-Dist: biocommons.seqrepo
 Requires-Dist: pydantic<3.0,>=2.0
 Requires-Dist: ga4gh.vrs<3.0,>=2.1.3
-Requires-Dist: wags-tails~=0.3.2
+Requires-Dist: wags-tails~=0.4.0
 Requires-Dist: bioutils
 Provides-Extra: dev
 Requires-Dist: pre-commit>=4.2.0; extra == "dev"
@@ -64,7 +62,7 @@ Requires-Dist: ruff==0.12.1; extra == "dev"
 Provides-Extra: tests
 Requires-Dist: pytest; extra == "tests"
 Requires-Dist: pytest-cov; extra == "tests"
-Requires-Dist: pytest-asyncio==0.18.3; extra == "tests"
+Requires-Dist: pytest-asyncio; extra == "tests"
 Requires-Dist: mock; extra == "tests"
 Provides-Extra: docs
 Requires-Dist: sphinx==6.1.3; extra == "docs"

{cool_seq_tool-0.14.1 → cool_seq_tool-0.14.3}/src/cool_seq_tool.egg-info/SOURCES.txt

@@ -39,6 +39,7 @@ src/cool_seq_tool/handlers/seqrepo_access.py
 src/cool_seq_tool/mappers/__init__.py
 src/cool_seq_tool/mappers/alignment.py
 src/cool_seq_tool/mappers/exon_genomic_coords.py
+src/cool_seq_tool/mappers/feature_overlap.py
 src/cool_seq_tool/mappers/liftover.py
 src/cool_seq_tool/mappers/mane_transcript.py
 src/cool_seq_tool/resources/__init__.py
@@ -51,6 +52,7 @@ src/cool_seq_tool/sources/transcript_mappings.py
 src/cool_seq_tool/sources/uta_database.py
 tests/conftest.py
 tests/test_utils.py
+tests/handlers/test_feature_overlap.py
 tests/handlers/test_seqrepo_access.py
 tests/mappers/test_alignment.py
 tests/mappers/test_exon_genomic_coords.py