Add -Y to minimap2 by default

Expose `--minimap2_args` on the command line.
Sentieon · Sep 17, 2024 · ee26dea · ee26dea
1 parent 9e49edb
commit ee26dea
Show file tree

Hide file tree

Showing 2 changed files with 13 additions and 0 deletions.
diff --git a/sentieon_cli/command_strings.py b/sentieon_cli/command_strings.py
@@ -260,6 +260,7 @@ def cmd_samtools_fastq_minimap2(
     sample_name: str,
     input_ref: Optional[pathlib.Path] = None,
     fastq_taglist: str = "*",
+    minimap2_args: str = "-Y",
     util_sort_args: str = "--cram_write_options version=3.0,compressor=rans",
 ) -> str:
     """Re-align an input BAM/CRAM/uBAM/uCRAM file with minimap2"""
@@ -288,6 +289,7 @@ def cmd_samtools_fastq_minimap2(
         "-t",
         str(cores),
         "-a",
+        minimap2_args,
         "-x",
         f"{model_bundle}/minimap2.model",
         str(reference),
@@ -434,6 +436,7 @@ def cmd_fastq_minimap2(
     model_bundle: pathlib.Path,
     cores: int,
     unzip: str = "gzip",
+    minimap2_args: str = "-Y",
     util_sort_args: str = "--cram_write_options version=3.0,compressor=rans",
 ) -> str:
     """Align an input fastq file with minimap2"""
@@ -449,6 +452,7 @@ def cmd_fastq_minimap2(
         "-t",
         str(cores),
         "-a",
+        minimap2_args,
         "-x",
         f"{model_bundle}/minimap2.model",
         "-R",

diff --git a/sentieon_cli/dnascope_longread.py b/sentieon_cli/dnascope_longread.py
@@ -72,6 +72,7 @@ def align_inputs(
     skip_version_check: bool = False,
     bam_format: bool = False,
     fastq_taglist: str = "*",
+    minimap2_args: str = "-Y",
     util_sort_args: str = "--cram_write_options version=3.0,compressor=rans",
     input_ref: Optional[pathlib.Path] = None,
     **_kwargs: Any,
@@ -107,6 +108,7 @@ def align_inputs(
                 sample_name,
                 input_ref,
                 fastq_taglist,
+                minimap2_args,
                 util_sort_args,
             )
         )
@@ -124,6 +126,7 @@ def align_fastq(
     readgroups: Optional[List[str]] = None,
     skip_version_check: bool = False,
     bam_format: bool = False,
+    minimap2_args: str = "-Y",
     util_sort_args: str = "--cram_write_options version=3.0,compressor=rans",
     **_kwargs: Any,
 ) -> List[pathlib.Path]:
@@ -166,6 +169,7 @@ def align_fastq(
                 model_bundle,
                 cores,
                 unzip,
+                minimap2_args,
                 util_sort_args,
             )
         )
@@ -726,6 +730,10 @@ def mosdepth(
     help="Use the BAM format instead of CRAM for output aligned files",
     action="store_true",
 )
+@arg(
+    "--minimap2_args",
+    help="Extra arguments for sentieon minimap2",
+)
 @arg(
     "--util_sort_args",
     help="Extra arguments for sentieon util sort",
@@ -766,6 +774,7 @@ def dnascope_longread(
     input_ref: Optional[pathlib.Path] = None,
     fastq_taglist: str = "*",  # pylint: disable=W0613
     bam_format: bool = False,  # pylint: disable=W0613
+    minimap2_args: str = "-Y",  # pylint: disable=W0613
     util_sort_args: str = (
         "--cram_write_options version=3.0,compressor=rans"
     ),  # pylint: disable=W0613