biopipen 0.27.9__py3-none-any.whl → 0.28.1__py3-none-any.whl

biopipen/__init__.py

@@ -1 +1 @@
-__version__ = "0.27.9"
+__version__ = "0.28.1"

biopipen/ns/cellranger.py

@@ -7,12 +7,15 @@ class CellRangerCount(Proc):
     """Run cellranger count
 
     to count gene expression and/or feature barcode reads
+    requires cellranger v7+.
 
     Input:
         fastqs: The input fastq files
             Either a list of fastq files or a directory containing fastq files
             If a directory is provided, it should be passed as a list with one
             element.
+        id: The id defining output directory. If not provided, it is inferred
+            from the fastq files.
 
     Output:
         outdir: The output directory
@@ -23,21 +26,28 @@ class CellRangerCount(Proc):
         ref: Path of folder containing 10x-compatible transcriptome reference
         tmpdir: Path to temporary directory, used to save the soft-lined fastq files
             to pass to cellranger
-        include_introns: Set to false to exclude intronic reads in count.
+        include_introns (flag): Set to false to exclude intronic reads in count.
+        create_bam (flag): Enable or disable BAM file generation.
+            This is required by cellrange v8+. When using cellrange v8-, it will be
+            transformed to `--no-bam`.
         <more>: Other environment variables required by `cellranger count`
             See `cellranger count --help` for more details or
             https://www.10xgenomics.com/support/software/cell-ranger/advanced/cr-command-line-arguments#count
     """  # noqa: E501
-    input = "fastqs:files"
+    input = "fastqs:files, id"
     output = """outdir:dir:
         {%- set fastqs = in.fastqs -%}
         {%- if len(fastqs) == 1 and isdir(fastqs[0]) -%}
             {%- set fastqs = fastqs[0] | glob: "*.fastq.gz" -%}
         {%- endif -%}
-        {%- set sample = commonprefix(*fastqs) |
-            regex_replace: "_L\\d+_?$", "" |
-            regex_replace: "_S\\d+$", "" -%}
-        {{- sample -}}
+        {%- if in.id -%}
+            {{in.id}}
+        {%- else -%}
+            {%- set id = commonprefix(*fastqs) |
+                regex_replace: "_L\\d+(:?_.*)?$", "" |
+                regex_replace: "_S\\d+$", "" -%}
+            {{- id -}}
+        {%- endif -%}
     """
     lang = config.lang.python
     envs = {
@@ -45,7 +55,8 @@ class CellRangerCount(Proc):
         "cellranger": config.exe.cellranger,
         "ref": config.ref.ref_cellranger_gex,
         "tmpdir": config.path.tmpdir,
-        "include_introns": "true",
+        "include_introns": True,
+        "create_bam": False,
     }
     script = "file://../scripts/cellranger/CellRangerCount.py"
     plugin_opts = {
@@ -57,13 +68,16 @@ class CellRangerCount(Proc):
 class CellRangerVdj(Proc):
     """Run cellranger vdj
 
-    to perform sequence assembly and paired clonotype calling
+    to perform sequence assembly and paired clonotype calling.
+    requires cellranger v7+.
 
     Input:
         fastqs: The input fastq files
             Either a list of fastq files or a directory containing fastq files
             If a directory is provided, it should be passed as a list with one
             element.
+        id: The id determining the output directory. If not provided, it is inferred
+            from the fastq files.
 
     Output:
         outdir: The output directory
@@ -78,16 +92,20 @@ class CellRangerVdj(Proc):
             See `cellranger vdj --help` for more details or
             https://www.10xgenomics.com/support/software/cell-ranger/advanced/cr-command-line-arguments#vdj
     """  # noqa: E501
-    input = "fastqs:files"
+    input = "fastqs:files, id"
     output = """outdir:dir:
         {%- set fastqs = in.fastqs -%}
         {%- if len(fastqs) == 1 and isdir(fastqs[0]) -%}
             {%- set fastqs = fastqs[0] | glob: "*.fastq.gz" -%}
         {%- endif -%}
-        {%- set sample = commonprefix(*fastqs) |
-            regex_replace: "_L\\d+_?$", "" |
-            regex_replace: "_S\\d+$", "" -%}
-        {{- sample -}}
+        {%- if in.id -%}
+            {{in.id}}
+        {%- else -%}
+            {%- set id = commonprefix(*fastqs) |
+                regex_replace: "_L\\d+(:?_.*)?$", "" |
+                regex_replace: "_S\\d+$", "" -%}
+            {{- id -}}
+        {%- endif -%}
     """
     lang = config.lang.python
     envs = {
@@ -129,4 +147,5 @@ class CellRangerSummary(Proc):
     envs = {"group": None}
     plugin_opts = {
         "report": "file://../reports/cellranger/CellRangerSummary.svelte",
+        "report_paging": 8,
     }

biopipen/ns/cellranger_pipeline.py

@@ -7,6 +7,7 @@ from __future__ import annotations
 from typing import TYPE_CHECKING
 
 from diot import Diot
+from pipen.utils import mark, is_loading_pipeline
 from pipen_args.procgroup import ProcGroup
 
 if TYPE_CHECKING:
@@ -17,28 +18,43 @@ class CellRangerCountPipeline(ProcGroup):
     """The cellranger count pipeline
 
     Run cellranger count for multiple samples and summarize the metrics.
+
+    Args:
+        input (type=list): The list of lists of fastq files.
+            or the list of comma-separated string of fastq files.
+        ids (type=list): The list of ids for the samples.
     """
-    DEFAULTS = Diot(input=None)
+    DEFAULTS = Diot(input=None, ids=None)
 
     def post_init(self):
         """Check if the input is a list of fastq files"""
-        if (
+        if not is_loading_pipeline() and (
             not isinstance(self.opts.input, (list, tuple))
             or len(self.opts.input) == 0
-            or not isinstance(self.opts.input[0], (list, tuple))
         ):
             raise TypeError(
                 "The input of `CellRangerCountPipeline` should be a list of lists of "
                 "fastq files."
             )
 
+        if isinstance(self.opts.input, (list, tuple)):
+            self.opts.input = [
+                [y.strip() for y in x.split(",")]
+                if isinstance(x, str)
+                else x
+                for x in self.opts.input
+            ]
+
     @ProcGroup.add_proc
     def p_cellranger_count(self) -> Proc:
         """Build CellRangerCount process"""
         from .cellranger import CellRangerCount as _CellRangerCount
 
         class CellRangerCount(_CellRangerCount):
-            input_data = self.opts.input
+            if self.opts.ids:
+                input_data = list(zip(self.opts.input, self.opts.ids))
+            else:
+                input_data = self.opts.input
 
         return CellRangerCount
 
@@ -58,28 +74,43 @@ class CellRangerVdjPipeline(ProcGroup):
     """The cellranger vdj pipeline
 
     Run cellranger vdj for multiple samples and summarize the metrics.
+
+    Args:
+        input (type=list): The list of lists of fastq files.
+            or the list of comma-separated string of fastq files.
+        ids (type=list): The list of ids for the samples.
     """
-    DEFAULTS = Diot(input=None)
+    DEFAULTS = Diot(input=None, ids=None)
 
     def post_init(self):
         """Check if the input is a list of fastq files"""
-        if (
+        if not is_loading_pipeline() and (
             not isinstance(self.opts.input, (list, tuple))
             or len(self.opts.input) == 0
-            or not isinstance(self.opts.input[0], (list, tuple))
         ):
             raise TypeError(
                 "The input of `CellRangerVdjPipeline` should be a list of lists of "
                 "fastq files."
             )
 
+        if isinstance(self.opts.input, (list, tuple)):
+            self.opts.input = [
+                [y.strip() for y in x.split(",")]
+                if isinstance(x, str)
+                else x
+                for x in self.opts.input
+            ]
+
     @ProcGroup.add_proc
     def p_cellranger_vdj(self) -> Proc:
         """Build CellRangerVdj process"""
         from .cellranger import CellRangerVdj as _CellRangerVdj
 
         class CellRangerVdj(_CellRangerVdj):
-            input_data = self.opts.input
+            if self.opts.ids:
+                input_data = list(zip(self.opts.input, self.opts.ids))
+            else:
+                input_data = self.opts.input
 
         return CellRangerVdj
 

biopipen/ns/scrna.py

@@ -1689,11 +1689,11 @@ class CellTypeAnnotation(Proc):
             ///
 
         sccatch_args (ns): The arguments for `scCATCH::findmarkergene()` if `tool` is `sccatch`.
-            - species (choice): The specie of cells.
-                - Human: Human cells.
-                - Mouse: Mouse cells.
+            - species: The specie of cells.
             - cancer: If the sample is from cancer tissue, then the cancer type may be defined.
             - tissue: Tissue origin of cells must be defined.
+            - marker: The marker genes for cell type identification.
+            - if_use_custom_marker (flag): Whether to use custom marker genes. If `True`, no `species`, `cancer`, and `tissue` are needed.
             - <more>: Other arguments for [`scCATCH::findmarkergene()`](https://rdrr.io/cran/scCATCH/man/findmarkergene.html).
                 You can pass an RDS file to `sccatch_args.marker` to work as custom marker. If so,
                 `if_use_custom_marker` will be set to `TRUE` automatically.
@@ -1743,7 +1743,13 @@ class CellTypeAnnotation(Proc):
         "sctype_tissue": None,
         "sctype_db": config.ref.sctype_db,
         "cell_types": [],
-        "sccatch_args": {},
+        "sccatch_args": {
+            "species": None,
+            "cancer": "Normal",
+            "tissue": None,
+            "marker": None,
+            "if_use_custom_marker": False,
+        },
         "hitype_tissue": None,
         "hitype_db": None,
         "celltypist_args": {

biopipen/scripts/cellranger/CellRangerCount.py

@@ -5,23 +5,17 @@ from biopipen.utils.misc import run_command
 
 fastqs = {{in.fastqs | repr}}  # pyright: ignore  # noqa
 outdir = {{out.outdir | quote}}  # pyright: ignore
+id = {{out.outdir | basename | quote}}  # pyright: ignore
 
 cellranger = {{envs.cellranger | quote}}  # pyright: ignore
 tmpdir = Path({{envs.tmpdir | quote}})  # pyright: ignore
 ref = {{envs.ref | quote}}  # pyright: ignore
 ncores = {{envs.ncores | int}}  # pyright: ignore
+include_introns = {{envs.include_introns | repr}}
+create_bam = {{envs.create_bam | repr}}
 
-{% if "id" in envs -%}
-id = {{envs.id | quote}}  # pyright: ignore
-{%- else -%}
-id = {{out.outdir | basename | quote}}  # pyright: ignore
-{%- endif %}
-
-{% if "sample" in envs -%}
-sample = {{envs.sample | quote}}  # pyright: ignore
-{%- else -%}
-sample = {{out.outdir | basename | quote}}  # pyright: ignore
-{%- endif %}
+include_introns = str(include_introns).lower()
+create_bam = str(create_bam).lower()
 
 # create a temporary unique directory to store the soft-linked fastq files
 fastqdir = tmpdir / f"cellranger_count_{uuid.uuid4()}"
@@ -34,25 +28,35 @@ for fastq in fastqs:
     fastq = Path(fastq)
     (fastqdir / fastq.name).symlink_to(fastq)
 
-other_args = {{envs | dict_to_cli_args: dashify=True, exclude=['cellranger', 'transcriptome', 'ref', 'tmpdir', 'id', 'sample', 'ncores']}}  # pyright: ignore
+other_args = {{envs | dict_to_cli_args: dashify=True, exclude=['no_bam', 'create_bam', 'include_introns', 'cellranger', 'transcriptome', 'ref', 'tmpdir', 'id', 'ncores']}}  # pyright: ignore
 
 command = [
     cellranger,
     "count",
     "--id",
     id,
-    "--sample",
-    sample,
     "--fastqs",
     fastqdir,
     "--transcriptome",
-    ref,
+    Path(ref).resolve(),
     "--localcores",
     ncores,
     "--disable-ui",
+    "--include-introns",
+    include_introns,
     *other_args,
 ]
 
+# check cellranger version
+#   cellranger cellranger-7.2.0
+version = run_command([cellranger, "--version"], stdout = "RETURN")
+version = version.replace("cellranger", "").replace("-", "").strip()
+version = list(map(int, version.split(".")))
+if version[0] >= 8:
+    command += ["--create-bam", create_bam]
+elif not create_bam:
+    command += ["--no-bam"]
+
 run_command(command, fg=True, cwd=str(Path(outdir).parent))
 
 web_summary_html = Path(outdir) / "outs" / "web_summary.html"
@@ -68,7 +72,7 @@ print("# Modify web_summary.html to move javascript to a separate file")
 try:
     web_summary_js = Path(outdir) / "outs" / "web_summary.js"
     web_summary_content = web_summary_html.read_text()
-    regex = re.compile(r"<script>(?=/\*! For license)(.+)</script>", re.DOTALL)
+    regex = re.compile(r"<script>(.+)</script>", re.DOTALL)
     web_summary_html.write_text(regex.sub(
         '<script src="web_summary.js"></script>',
         web_summary_content,

biopipen/scripts/cellranger/CellRangerVdj.py

@@ -5,24 +5,13 @@ from biopipen.utils.misc import run_command
 
 fastqs = {{in.fastqs | repr}}  # pyright: ignore  # noqa
 outdir = {{out.outdir | quote}}  # pyright: ignore
+id = {{out.outdir | basename | quote}}  # pyright: ignore
 
 cellranger = {{envs.cellranger | quote}}  # pyright: ignore
 tmpdir = Path({{envs.tmpdir | quote}})  # pyright: ignore
 ref = {{envs.ref | quote}}  # pyright: ignore
 ncores = {{envs.ncores | int}}  # pyright: ignore
 
-{% if "id" in envs -%}
-id = {{envs.id | quote}}  # pyright: ignore
-{%- else -%}
-id = {{out.outdir | basename | quote}}  # pyright: ignore
-{%- endif %}
-
-{% if "sample" in envs -%}
-sample = {{envs.sample | quote}}  # pyright: ignore
-{%- else -%}
-sample = {{out.outdir | basename | quote}}  # pyright: ignore
-{%- endif %}
-
 # create a temporary unique directory to store the soft-linked fastq files
 fastqdir = tmpdir / f"cellranger_count_{uuid.uuid4()}"
 fastqdir.mkdir(parents=True, exist_ok=True)
@@ -34,19 +23,17 @@ for fastq in fastqs:
     fastq = Path(fastq)
     (fastqdir / fastq.name).symlink_to(fastq)
 
-other_args = {{envs | dict_to_cli_args: dashify=True, exclude=['cellranger', 'reference', 'ref', 'tmpdir', 'id', 'sample', 'ncores']}}  # pyright: ignore
+other_args = {{envs | dict_to_cli_args: dashify=True, exclude=['cellranger', 'reference', 'ref', 'tmpdir', 'id', 'ncores']}}  # pyright: ignore
 
 command = [
     cellranger,
     "vdj",
     "--id",
     id,
-    "--sample",
-    sample,
     "--fastqs",
     fastqdir,
     "--reference",
-    ref,
+    Path(ref).resolve(),
     "--localcores",
     ncores,
     "--disable-ui",
@@ -68,7 +55,7 @@ print("# Modify web_summary.html to move javascript to a separate file")
 try:
     web_summary_js = Path(outdir) / "outs" / "web_summary.js"
     web_summary_content = web_summary_html.read_text()
-    regex = re.compile(r"<script>(?=/\*! For license)(.+)</script>", re.DOTALL)
+    regex = re.compile(r"<script>(.+)</script>", re.DOTALL)
     web_summary_html.write_text(regex.sub(
         '<script src="web_summary.js"></script>',
         web_summary_content,

biopipen/scripts/scrna/CellTypeAnnotation-sccatch.R

@@ -7,19 +7,12 @@ outfile = {{out.outfile | r}}
 sccatch_args = {{envs.sccatch_args | r}}
 newcol = {{envs.newcol | r}}
 
-if (is.null(sccatch_args$tissue)) { stop("`envs.sccatch_args.tissue` origin of cells must be defined.") }
-if (is.null(sccatch_args$species)) {
-    sccatch_args$species = "Human"
-}
 if (!is.null(sccatch_args$marker)) {
     cellmatch = readRDS(sccatch_args$marker)
     sccatch_args$if_use_custom_marker = TRUE
 }
 sccatch_args$marker = cellmatch
 
-if (is.null(sccatch_args$cancer)) {
-    sccatch_args$cancer = "Normal"
-}
 if (is.integer(sccatch_args$use_method)) {
     sccatch_args$use_method = as.character(sccatch_args$use_method)
 }
@@ -30,6 +23,8 @@ obj = createscCATCH(data = GetAssayData(sobj), cluster = as.character(Idents(sob
 sccatch_args$object = obj
 
 obj = do_call(findmarkergene, sccatch_args)
+obj = findcelltype(object = obj)
+
 write.table(
     obj@celltype,
     file = file.path(dirname(outfile), "cluster2celltype.tsv"),

biopipen/scripts/scrna/CellsDistribution.R

@@ -360,6 +360,7 @@ do_case <- function(name, case) {
     }
 
     log_info("  Merging and saving pie charts ...")
+    devpars = case$devpars
     # assemble and save pie chart plots
     res <- devpars$res %||% 100
     #                         legend, cells_by names
@@ -405,6 +406,7 @@ do_case <- function(name, case) {
     }
 
     col_fun <- colorRamp2(c(0, max(hmdata, na.rm = T)), c("lightyellow", "purple"))
+    hm_devpars <- case$hm_devpars
     hm_res <- hm_devpars$res %||% 100
     hm_width <- hm_devpars$width %||% (600 + 15 * length(unique(meta$seurat_clusters)) + extra_width)
     hm_height <- hm_devpars$height %||% (450 + 15 * cells_rows + extra_height)

{biopipen-0.27.9.dist-info → biopipen-0.28.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: biopipen
-Version: 0.27.9
+Version: 0.28.1
 Summary: Bioinformatics processes/pipelines that can be run from `pipen run`
 License: MIT
 Author: pwwang

{biopipen-0.27.9.dist-info → biopipen-0.28.1.dist-info}/RECORD

@@ -1,4 +1,4 @@
-biopipen/__init__.py,sha256=sJWdNHkiTx9zQyc6_YJknVzen19cpBrSswdfVQaZ7S8,23
+biopipen/__init__.py,sha256=ZRQKbgDaGz_yuLk-cUKuk6ZBKCSRKZC8nQd041NRNXk,23
 biopipen/core/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 biopipen/core/config.py,sha256=edK5xnDhM8j27srDzsxubi934NMrglLoKrdcC8qsEPk,1069
 biopipen/core/config.toml,sha256=20RCI30Peee1EQdfb_UbV3Hf74XUPndJnYZlUThytsw,1781
@@ -10,8 +10,8 @@ biopipen/ns/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 biopipen/ns/bam.py,sha256=5AsYrB0mtr_mH6mCL6gjJ5rC4NywpjFkpFjUrBGp7Fk,9301
 biopipen/ns/bcftools.py,sha256=puCDfIL-1z6cz2y1Rlz-ESNIr8xJgeIjEQ440qicCvM,3467
 biopipen/ns/bed.py,sha256=UN38qUChDeE-ipuSBY8RVLwvJqM2wxSRmlhOiDo4JG0,5395
-biopipen/ns/cellranger.py,sha256=2FwE8YNsGEXn-vS0AmO7m5xmlTgNQjL1lhWUeJ6XNJA,4796
-biopipen/ns/cellranger_pipeline.py,sha256=D6gvIeasHjDCdro7f4wjomxRYTtsJT77Ld47XzlU3Fw,2891
+biopipen/ns/cellranger.py,sha256=yPBoNzVSY74J7uyVucaob5lqZKKru5-hYSM4f4Nr2OY,5553
+biopipen/ns/cellranger_pipeline.py,sha256=lyRUdJEeQaTWGUz5UbrSdHMH_cgzXAXLlERmanppEBs,4037
 biopipen/ns/cnv.py,sha256=vq6dZfEOyuVuqg3nP6FQtNmQ-JocpBJMX9IYlZ0OPD0,6803
 biopipen/ns/cnvkit.py,sha256=5mA2Q8-YDs4g1HoxtpB_NWnyZYwEThNr3s3wlubLQrQ,31130
 biopipen/ns/cnvkit_pipeline.py,sha256=2fJLn70L2jJ81ZMNdnU84Sf3HoKA2CSnHuDzLGR8jmw,36854
@@ -21,7 +21,7 @@ biopipen/ns/gsea.py,sha256=EsNRAPYsagaV2KYgr4Jv0KCnZGqayM209v4yOGGTIOI,7423
 biopipen/ns/misc.py,sha256=fzn0pXvdghMkQhu-e3MMapPNMyO6IAJbtTzVU3GbFa0,3246
 biopipen/ns/plot.py,sha256=fzJAKKl4a_tsVkLREGCQTFVHP049m33LdWgeYRb6v7M,5483
 biopipen/ns/rnaseq.py,sha256=bKAa6friFWof4yDTWZQahm1MS-lrdetO1GqDKdfxXYc,7708
-biopipen/ns/scrna.py,sha256=KL5Eu0mnIITLLSHAIz_sgr4ssmEU6AuBDXwedqYU7BI,105633
+biopipen/ns/scrna.py,sha256=OCEuE-TcqR7ADEEsgA6qL9rDK5R75YWRmmmvaNIFtZU,105925
 biopipen/ns/scrna_metabolic_landscape.py,sha256=6AhaynGG3lNRi96N2tReVT46BJMuEwooSSd2irBoN80,28347
 biopipen/ns/snp.py,sha256=EQ2FS0trQ7YThPmBVTpS66lc2OSfgQ6lCh6WnyP-C2g,5499
 biopipen/ns/stats.py,sha256=yJ6C1CXF84T7DDs9mgufqUOr89Rl6kybE5ji8Vnx6cw,13693
@@ -81,9 +81,9 @@ biopipen/scripts/bed/Bed2Vcf.py,sha256=u0mp_2Y4UtEA839zq9UENesH6Gyiwd4sZQW9wFnBV
 biopipen/scripts/bed/BedConsensus.py,sha256=gfAxuIalvCEpS0tiOyAJGPYGgHN0L-hm0K37Iteh5yw,2386
 biopipen/scripts/bed/BedLiftOver.sh,sha256=Y4gBsz9w4zhE29UmWojO6F4PXMMMWC1uCzjrxa19eOs,256
 biopipen/scripts/bed/BedtoolsMerge.py,sha256=TjKO5MpUzDj931bQAWku2660MVSiZzdMHt_v2Xbt0IE,355
-biopipen/scripts/cellranger/CellRangerCount.py,sha256=ZDcry8suLhulXiTsl01LGKmSJkewJ-TgHazLtfsBr6U,2516
+biopipen/scripts/cellranger/CellRangerCount.py,sha256=Mi6PiufFPp-cb4WRgZhniXi2ja5VUusRDW9RdAqFyX0,2831
 biopipen/scripts/cellranger/CellRangerSummary.R,sha256=8EqP4j074RLi5Adw47xyqQwkDEcvSS6Yfk0UgKKYOuo,11045
-biopipen/scripts/cellranger/CellRangerVdj.py,sha256=-QbhPKqFBZ15Es6NJaU7Lwf1KQW_3Lyv0aISh-Urk2M,2504
+biopipen/scripts/cellranger/CellRangerVdj.py,sha256=ipGkPzkiVy5IpPrYKUK4zgEIvz-aMJOKbMEN-bMStFE,2200
 biopipen/scripts/cnv/AneuploidyScore.R,sha256=liAN8u8_lj8voJ01oBW9Dw09yi388KF5f_gwPOv0wdE,8437
 biopipen/scripts/cnv/AneuploidyScoreSummary.R,sha256=9Zni5zqYfzevs5XSAt3fqD9WZ_RWr_ByUnXReKLLWoY,12337
 biopipen/scripts/cnv/TMADScore.R,sha256=uCLHQR6sMt-4uVUAEJlJxYXlai9ZE5J7xBl0sl-EkjU,1065
@@ -121,10 +121,10 @@ biopipen/scripts/scrna/AnnData2Seurat.R,sha256=16ASr8B7B7tj7hZpwL31v6GASfh8wcBwI
 biopipen/scripts/scrna/CellTypeAnnotation-celltypist.R,sha256=Hzqt-I41ThiSGPzqrg5ma4IQbCDnM1ztfTscPqthdfM,9009
 biopipen/scripts/scrna/CellTypeAnnotation-direct.R,sha256=ST7hJo6IY3-lfj1ymFCxxTH4cAsu8CVTcOKpELBeyuo,1855
 biopipen/scripts/scrna/CellTypeAnnotation-hitype.R,sha256=6_DBAlLKcHqaMyWGZWvTd4gFfHymfz9s2XLja8aj1qA,1869
-biopipen/scripts/scrna/CellTypeAnnotation-sccatch.R,sha256=1ejye0hs-EOwzzdP9gFWSLPcF6dOAA6VmNKXEjmS11E,1654
+biopipen/scripts/scrna/CellTypeAnnotation-sccatch.R,sha256=HO1ta0VWdBY1SuUNDqeb_Rp3UrbcqJgGpglegvMUGpw,1435
 biopipen/scripts/scrna/CellTypeAnnotation-sctype.R,sha256=1nKRtzhVoJ9y0yMg1sgI6u7czsrk2cN0FvNUCZo8l-o,3878
 biopipen/scripts/scrna/CellTypeAnnotation.R,sha256=OwLM_G4D7TG4HhIJjQxgIQM92X86lsWp9MVyXTTkLSc,618
-biopipen/scripts/scrna/CellsDistribution.R,sha256=isDr5-EWvOeWwVZdjOSsdX3QUpEaDvQFulIYawqFaQc,18854
+biopipen/scripts/scrna/CellsDistribution.R,sha256=lCobUWmNv97PF0HwwRbn9CGAJ1HGRuzlJe6CvlUijf0,18915
 biopipen/scripts/scrna/DimPlots.R,sha256=-mXOTMnpPxvR30XLjwcohFfFx7xTqWKKiICwJiD6yEo,1554
 biopipen/scripts/scrna/ExprImputation-alra.R,sha256=w3W1txJcdWg52-SETY2Z0lO7maDNfiMjBYIGN587YW0,843
 biopipen/scripts/scrna/ExprImputation-rmagic.R,sha256=jYIfqZpnvjKJkvItLnemPVtUApHBYQi1_L8rHVbEe1M,735
@@ -240,7 +240,7 @@ biopipen/utils/reference.py,sha256=6bPSwQa-GiDfr7xLR9a5T64Ey40y24yn3QfQ5wDFZkU,4
 biopipen/utils/rnaseq.R,sha256=Ro2B2dG-Z2oVaT5tkwp9RHBz4dp_RF-JcizlM5GYXFs,1298
 biopipen/utils/single_cell.R,sha256=pJjYP8bIZpNAtTQ32rOXhZxaM1Y-6D-xUcK3pql9tbk,4316
 biopipen/utils/vcf.py,sha256=ajXs0M_QghEctlvUlSRjWQIABVF02wPdYd-0LP4mIsU,9377
-biopipen-0.27.9.dist-info/METADATA,sha256=-Uvdu_dbpiM0EPBLSFLhhpySPAPxGtA6_FjRlMo0bME,882
-biopipen-0.27.9.dist-info/WHEEL,sha256=sP946D7jFCHeNz5Iq4fL4Lu-PrWrFsgfLXbbkciIZwg,88
-biopipen-0.27.9.dist-info/entry_points.txt,sha256=wu70aoBcv1UahVbB_5237MY-9M9_mzqmWjDD-oi3yz0,621
-biopipen-0.27.9.dist-info/RECORD,,
+biopipen-0.28.1.dist-info/METADATA,sha256=UPLio_HdqVCyeNe2VCcUTZXz1CGakrQFH7_Dbb9qnHQ,882
+biopipen-0.28.1.dist-info/WHEEL,sha256=sP946D7jFCHeNz5Iq4fL4Lu-PrWrFsgfLXbbkciIZwg,88
+biopipen-0.28.1.dist-info/entry_points.txt,sha256=wu70aoBcv1UahVbB_5237MY-9M9_mzqmWjDD-oi3yz0,621
+biopipen-0.28.1.dist-info/RECORD,,