phc-ingestion 0.8.33__py3-none-any.whl → 0.8.34__py3-none-any.whl

ingestion/nextgen/util/interpretation.py

@@ -2,6 +2,9 @@ from logging import Logger
 
 
 def map_interpretation(status: str, log: Logger):
+    """
+    Map interpretation for structural and copy number variants
+    """
     if status == "Pathogenic":
         return "Pathogenic"
     elif "VUS" in status:
@@ -9,3 +12,17 @@ def map_interpretation(status: str, log: Logger):
     else:
         log.error(f"Failed to resolve interpretation: {status}")
         return ""
+
+
+def map_vendsig(vendsig: str) -> str:
+    """
+    Map vendor significance for short variants
+    """
+    if vendsig in ["Pathogenic"]:
+        return "VENDSIG=Pathogenic"
+    elif vendsig in ["Likely Pathogenic", "LikelyPathogenic"]:
+        return "VENDSIG=Likely pathogenic"
+    elif vendsig in ["VUS"]:
+        return "VENDSIG=Uncertain significance"
+    else:
+        raise RuntimeError(f"Unable to map vendor significance: {vendsig}")

ingestion/nextgen/util/nextgen_specific_genes.py

@@ -14,7 +14,7 @@ nextgen_specific_genes_with_location: list[GeneWithLocation] = [
     {"gene": "CCND3", "chr": "chr6", "start": 41920534, "end": 42562008},
     {"gene": "MYC", "chr": "chr8", "start": 125309416, "end": 129673293},
     {"gene": "CCND1", "chr": "chr11", "start": 69090733, "end": 69656860},
-    {"gene": "IGH", "chr": "chr14", "start": 105516968, "end": 109902208},
+    {"gene": "IGH", "chr": "chr14", "start": 105325507, "end": 109902208},
     {"gene": "MAF", "chr": "chr16", "start": 78428398, "end": 79615096},
     {"gene": "MAFB", "chr": "chr20", "start": 39039005, "end": 40688948},
     {"gene": "IGL", "chr": "chr22", "start": 22012552, "end": 22965858},

ingestion/nextgen/util/process_manifest.py

@@ -65,7 +65,7 @@ def get_cell_purity(interpretation_lines: list):
         return float(00.00)
 
 
-def extract_patient_data(patient_info_lines: list):
+def extract_patient_data(patient_info_lines: list[str]):
     patient_data: dict = {}
     patient_data["patientInfo"] = {}
 

ingestion/nextgen/util/process_vcf.py

@@ -5,6 +5,7 @@ from logging import Logger
 from typing import Literal
 
 from ingestion.nextgen.util.alteration_table import AlterationTableRow, ShortVariantGene
+from ingestion.nextgen.util.interpretation import map_vendsig
 
 SequenceType = Literal["somatic", "germline"]
 
@@ -204,17 +205,6 @@ def process_vcf(
     return {"vcf_path_name": vcf_path, "vcf_line_count": line_count}
 
 
-def map_vendsig(vendsig: str) -> str:
-    if vendsig in ["Pathogenic"]:
-        return "VENDSIG=Pathogenic"
-    elif vendsig in ["Likely Pathogenic", "LikelyPathogenic"]:
-        return "VENDSIG=Likely pathogenic"
-    elif vendsig in ["VUS"]:
-        return "VENDSIG=Uncertain significance"
-    else:
-        raise RuntimeError(f"Unable to map vendor significance: {vendsig}")
-
-
 def add_vendsig_to_info(
     info: str,
     short_variant_table_rows: list[AlterationTableRow[ShortVariantGene]],

ingestion/vcf_standardization/Variant.py

@@ -18,7 +18,7 @@ class Variant:
         self.info = {x.split("=")[0]: x.split("=")[1] for x in fields[7].split(";") if "=" in x}
         self.frmt = fields[8].split(":")
         self.smpl = fields[9].split(":")
-        self.ad_af_dp = {"AD": False, "AF": False, "DP": False}
+        self.ad_af_dp: dict[str, bool | str] = {"AD": False, "AF": False, "DP": False}
 
     def standardize_allele_frequency(self, log):
         # Detect if allele frequency is present either in the INFO or FORMAT/SAMPLE fields
@@ -131,19 +131,19 @@ class Variant:
         )
         return updated_variant
 
+    @classmethod
+    def check_formatting(cls, var: str):
+        # Loose formatting check, return as Variant class object
+        split_var = var.split("\t")
+        if len(split_var) < 8 or not split_var[1].isdigit():
+            raise RuntimeError(f"Variant contains incorrect number, or invalid fields:  {var}")
 
-def check_formatting(var: str) -> Variant:
-    # Loose formatting check, return as Variant class object
-    split_var = var.split("\t")
-    if len(split_var) < 8 or not split_var[1].isdigit():
-        raise RuntimeError(f"Variant contains incorrect number, or invalid fields:  {var}")
+        if len(split_var) == 8:
+            split_var.append(".")  # Add placeholder for FORMAT
+            split_var.append(".")  # Add placeholder for SAMPLE
 
-    if len(split_var) == 8:
-        split_var.append(".")  # Add placeholder for FORMAT
-        split_var.append(".")  # Add placeholder for SAMPLE
+        elif len(split_var) == 9:
+            split_var.append(".")  # Add placeholder for SAMPLE
 
-    elif len(split_var) == 9:
-        split_var.append(".")  # Add placeholder for SAMPLE
-
-    working_variant = Variant(split_var)
-    return working_variant
+        working_variant = cls(split_var)
+        return working_variant

ingestion/vcf_standardization/standardize.py

@@ -9,15 +9,15 @@ from ingestion.vcf_standardization.util.read_write import (
     read_headers,
     read_variants,
 )
-from ingestion.vcf_standardization.Variant import check_formatting
+from ingestion.vcf_standardization.Variant import Variant
 
 
-def format_variant(variant: str, log: Logger, vendsig_dict: dict = None) -> Optional[str]:
+def format_variant(variant: str, log: Logger, vendsig_dict: dict | None = None) -> Optional[str]:
     # Ignore structural variants
     if "SVTYPE" in variant:
         return None
     # Working variant
-    wv = check_formatting(variant)
+    wv = Variant.check_formatting(variant)
 
     # Only process variants that aren't multiallelic
     if len(wv.alt.split(",")) == 1:
@@ -48,7 +48,7 @@ def standardize_vcf(
     out_path: str,
     case_id: str,
     log: Logger,
-    vendsig_dict: dict = None,
+    vendsig_dict: dict | None = None,
     compression: bool = False,
 ) -> Optional[int]:
     check_vcf(infile, log)

ingestion/vcf_standardization/util/read_write.py

@@ -1,10 +1,11 @@
 import gzip
+from logging import Logger
 import re
 import os
 from typing import Iterator, Optional
 
 
-def check_vcf(infile, log):
+def check_vcf(infile: str, log: Logger) -> None:
     log.info("Checking VCF file")
     # Check if file exists. Raise if it doesn't.
     if os.path.exists(infile) == False:
@@ -68,7 +69,7 @@ def write_vcf(
     compression: bool,
     line_count: int,
     log,
-):
+) -> int:
     log.info(f"Writing standardized VCF to {outfile}")
 
     with gzip.open(outfile, "wt") if compression else open(outfile, "w") as w:

{phc_ingestion-0.8.33.dist-info → phc_ingestion-0.8.34.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: phc-ingestion
-Version: 0.8.33
+Version: 0.8.34
 Summary: Functions for LifeOmic PHC genomic ingestions
 License: MIT
 Author-email: LifeOmic Development <development@lifeomic.com>

{phc_ingestion-0.8.33.dist-info → phc_ingestion-0.8.34.dist-info}/RECORD

@@ -31,14 +31,14 @@ ingestion/generic/utils.py,sha256=1MEIru7uq38IjUdL8lcHqDH0oTki9uWrz1f2e-pmRoU,28
 ingestion/nextgen/__init__.py,sha256=7LQ-h_Bvc5P1QcHMdzsqi1Qm4fTJn04-ozar2ty9wSc,59
 ingestion/nextgen/process.py,sha256=kDCnU685v7aqJ3i4HpFdb7HqgHRSBKqtYPpuyN7qWmM,3976
 ingestion/nextgen/util/alteration_table.py,sha256=OqstLK6cgoNvRWy8bW6_iABaAn5ggCi1xBM8GOU6wYQ,6060
-ingestion/nextgen/util/interpretation.py,sha256=ozuzb0vozff34zfP6AdOiUmI8Q77hI02jve_nCPZHfE,297
+ingestion/nextgen/util/interpretation.py,sha256=tQ3qzAUwBzwK1tQzH9ujZk_VCQ2wP8HzSZY0fImJ5Wo,818
 ingestion/nextgen/util/manifest_helpers.py,sha256=2xrpEtHbCb1Kea1wJeObkDfTiBklmffQt_o2hMgOSOE,1208
-ingestion/nextgen/util/nextgen_specific_genes.py,sha256=hgam7HVE324FwOf7G4Wk4cUArch9vHIjBZRUUyF3ukg,1206
+ingestion/nextgen/util/nextgen_specific_genes.py,sha256=1jFcqvtYAlJ7eBwOBm1UC2TzAbjHjdlvPBUzxr1G8dY,1206
 ingestion/nextgen/util/pre_filter_somatic_vcf.py,sha256=mIaUihmGLbS38D4Gy_Qtf1lFAfW0A-LgAgQmsrEiI-M,3529
 ingestion/nextgen/util/process_cnv.py,sha256=MIirc8e0k6lsaTZkRM3U3L3IvbrcHmKQ4xlIu585514,2430
-ingestion/nextgen/util/process_manifest.py,sha256=EGYaTcub4M08mFTAh4CNHPRkP8_a5r4jMJaExm9Nkko,8423
+ingestion/nextgen/util/process_manifest.py,sha256=EnV9I90vnanDvuoErbMfz6yAfjzM5LdhhUF4q5DJd8w,8428
 ingestion/nextgen/util/process_structural.py,sha256=FKjkK7BkIlocnLs8rFCjrMC39FCQnD0nQCeWvi7cRoA,7539
-ingestion/nextgen/util/process_vcf.py,sha256=SN0C13F45R_N5UaMaVSUDSCtIMmpHfaMTo7_5PkFkrM,8085
+ingestion/nextgen/util/process_vcf.py,sha256=ZZURSMnZhHDpFahzijZ4MvCfSWTPdIktzmnCKVVUbGs,7768
 ingestion/nextgen/util/types.py,sha256=SSzt5gv-kss1PR45eQUelypWrGI-dAfQMO3GSD-T-Wg,22
 ingestion/resources/GRCh37_map.csv.gz,sha256=JOEkjtbYrJpIdyoZdCvfJhvvz2dNfkSve7lXSXkCCD8,408290
 ingestion/resources/GRCh38_map.csv.gz,sha256=qriYO2_buCCb4T6WcuZ-pCwPxMsm0TL2OxAHvJ1cEfA,612373
@@ -47,13 +47,13 @@ ingestion/shared_util/coords_to_genes.py,sha256=vz9EfgFm3BS6pEPnslbEka8cJKlQZtHJ
 ingestion/shared_util/gene_to_coords.py,sha256=M-q5ateLSQ4fCF0uMk5TX2uBLRrcZzXqXEf05TPaLsU,876
 ingestion/shared_util/open_maybe_gzipped.py,sha256=FrOPJ4OgfpQGyT3f1Su1rFeuuYYu6QJ-nVIBIosbfhw,232
 ingestion/shared_util/types.py,sha256=u9AD2OrTQWMBtK_7VXHsD8Rv6HFs-7ZUItNl4KXdL7k,68
-ingestion/vcf_standardization/Variant.py,sha256=MBT8x25Ub1GRkTtnvPMV9SN7LUoF4Xkons-m3kYDkeY,5398
+ingestion/vcf_standardization/Variant.py,sha256=aoSqT4XAECxCF0JZgv1YRmxuw20WGeWkwFTw0x6FmWc,5475
 ingestion/vcf_standardization/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-ingestion/vcf_standardization/standardize.py,sha256=PaRqQRSrnI79WFgNvbvnXFcvXHYqshOOxCXup5eMZ_M,2289
+ingestion/vcf_standardization/standardize.py,sha256=zYzZxncq8USA1bUs26L-ByLPTnUlGyVvS3LJVGD19BM,2302
 ingestion/vcf_standardization/util/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 ingestion/vcf_standardization/util/af_helpers.py,sha256=dpTzoeIQVeBRt0ETF3a9rp5ojZqznHg4x_hCZ8OPcOg,1061
 ingestion/vcf_standardization/util/dp_helpers.py,sha256=Nq8oLOLObu4_pv16qwwgpALRlUoJVCULrd9cFOD-eoI,823
-ingestion/vcf_standardization/util/read_write.py,sha256=IQotJ27To1MoQcRstc5AbHZtUuJz5cqkkZiHsDNaBvI,2471
-phc_ingestion-0.8.33.dist-info/WHEEL,sha256=B19PGBCYhWaz2p_UjAoRVh767nYQfk14Sn4TpIZ-nfU,87
-phc_ingestion-0.8.33.dist-info/METADATA,sha256=CM6kTtndCIkqq55vXl_x2vehEZ_mL29fK7TblsCsz9E,552
-phc_ingestion-0.8.33.dist-info/RECORD,,
+ingestion/vcf_standardization/util/read_write.py,sha256=x3Pf6Dq8tmolblbCS5CrNmrcHS3FGfqBSFpFgvFGC4g,2526
+phc_ingestion-0.8.34.dist-info/WHEEL,sha256=B19PGBCYhWaz2p_UjAoRVh767nYQfk14Sn4TpIZ-nfU,87
+phc_ingestion-0.8.34.dist-info/METADATA,sha256=5afEPjKh3tGTH02AjkVcQ-AluTJ-JuEqSBiZkD7HHAU,552
+phc_ingestion-0.8.34.dist-info/RECORD,,