alternative to BwaMemOptions

docs/api.rst

@@ -66,17 +66,19 @@ It is constructed directly and options set on the object:
    recs = aln.align(queries=["GATTACA"], opt=opt)
 
 
-The :meth:`pybwa.BwaMem.align` method accepts custom options provided as a :class:`~pybwa.BwaMemOptions` object.
-It is constructed via the :class:`~pybwa.BwaMemOptionsBuilder` class, to support scaling gap open and extend penalties
-when a custom match score, or the specification of presets (via `mode`).
+Similarly, the :meth:`pybwa.BwaMem.align` method accepts custom options provided as a :class:`~pybwa.BwaMemOptions` object.
+It is constructed directly and options set on the object:
 
 .. code-block:: python
 
-   builder = BwaMemOptionsBuilder()
-   builder.min_seed_len = 32
-   opt: BwaMemOptions = builder.build()
+   opt = BwaMemOptions()
+   opt.min_seed_len = 32
    recs = aln.align(queries=["GATTACA"], opt=opt)
 
+Note: the match score when set only scales other options (:code:`-TdBOELU`) unless the latter are specified
+**and** only after the :meth:`~pybwa.BwaMemOptions.finalize` has been called.
+The options become immutable unless :code:`copy=True` is passed to :meth:`~pybwa.BwaMemOptions.finalize`.
+
 The :class:`~pybwa.BwaIndex` object is useful when re-using the same index, such that it only needs to be loaded into memory
 once.
 Both constructors for the :class:`~pybwa.BwaAln` and :class:`~pybwa.BwaMem` objects accept an index.
@@ -109,9 +111,6 @@ Bwa Aln
 Bwa Mem
 =======
 
-.. autoclass:: pybwa.BwaMemOptionsBuilder
-   :members:
-
 .. autoclass:: pybwa.BwaMemOptions
    :members:
 

pybwa/libbwamem.pyi

@@ -18,42 +18,136 @@ class BwaMemMode(enum.Enum):
 
 class BwaMemOptions:
     def __init__(self, finalize: bool = False) -> None: ...
+    _finalized: bool
     _ignore_alt: bool
-    min_seed_len: int
-    mode: BwaMemMode
-    band_width: int
-    match_score: int
-    mismatch_penalty: int
-    minimum_score: int
-    unpaired_penalty: int
-    n_threads: int
-    skip_pairing: bool
-    output_all_for_fragments: bool
-    interleaved_paired_end: bool
-    short_split_as_secondary: bool
-    skip_mate_rescue: bool
-    soft_clip_supplementary: bool
-    with_xr_tag: bool
-    query_coord_as_primary: bool
-    keep_mapq_for_supplementary: bool
-    with_xb_tag: bool
-    max_occurrences: int
-    off_diagonal_x_dropoff: float
-    ignore_alternate_contigs: bool
-    internal_seed_split_factor: float
-    drop_chain_fraction: float
-    max_mate_rescue_rounds: int
-    min_seeded_bases_in_chain: int
-    seed_occurrence_in_3rd_round: int
-    xa_max_hits: int | tuple[int, int]
-    xa_drop_ratio: float
-    gap_open_penalty: int | tuple[int, int]
-    gap_extension_penalty: int | tuple[int, int]
-    clipping_penalty: int | tuple[int, int]
-
-class BwaMemOptionsBuilder(BwaMemOptions):
-    def __init__(self, options: BwaMemOptions | None = None) -> None: ...
-    def build(self) -> BwaMemOptions: ...
+    _mode: BwaMemMode | None
+    @property
+    def finalized(self) -> bool: ...
+    @property
+    def min_seed_len(self) -> int: ...
+    @min_seed_len.setter
+    def min_seed_len(self, value: int) -> None: ...
+    @property
+    def mode(self) -> BwaMemMode: ...
+    @mode.setter
+    def mode(self, value: BwaMemMode) -> None: ...
+    @property
+    def band_width(self) -> int: ...
+    @band_width.setter
+    def band_width(self, value: int) -> None: ...
+    @property
+    def match_score(self) -> int: ...
+    @match_score.setter
+    def match_score(self, value: int) -> None: ...
+    @property
+    def mismatch_penalty(self) -> int: ...
+    @mismatch_penalty.setter
+    def mismatch_penalty(self, value: int) -> None: ...
+    @property
+    def minimum_score(self) -> int: ...
+    @minimum_score.setter
+    def minimum_score(self, value: int) -> None: ...
+    @property
+    def unpaired_penalty(self) -> int: ...
+    @unpaired_penalty.setter
+    def unpaired_penalty(self, value: int) -> None: ...
+    @property
+    def n_threads(self) -> int: ...
+    @n_threads.setter
+    def n_threads(self, value: int) -> None: ...
+    @property
+    def skip_pairing(self) -> bool: ...
+    @skip_pairing.setter
+    def skip_pairing(self, value: bool) -> None: ...
+    @property
+    def output_all_for_fragments(self) -> bool: ...
+    @output_all_for_fragments.setter
+    def output_all_for_fragments(self, value: bool) -> None: ...
+    @property
+    def interleaved_paired_end(self) -> bool: ...
+    @interleaved_paired_end.setter
+    def interleaved_paired_end(self, value: bool) -> None: ...
+    @property
+    def short_split_as_secondary(self) -> bool: ...
+    @short_split_as_secondary.setter
+    def short_split_as_secondary(self, value: bool) -> None: ...
+    @property
+    def skip_mate_rescue(self) -> bool: ...
+    @skip_mate_rescue.setter
+    def skip_mate_rescue(self, value: bool) -> None: ...
+    @property
+    def soft_clip_supplementary(self) -> bool: ...
+    @soft_clip_supplementary.setter
+    def soft_clip_supplementary(self, value: bool) -> None: ...
+    @property
+    def with_xr_tag(self) -> bool: ...
+    @with_xr_tag.setter
+    def with_xr_tag(self, value: bool) -> None: ...
+    @property
+    def query_coord_as_primary(self) -> bool: ...
+    @query_coord_as_primary.setter
+    def query_coord_as_primary(self, value: bool) -> None: ...
+    @property
+    def keep_mapq_for_supplementary(self) -> bool: ...
+    @keep_mapq_for_supplementary.setter
+    def keep_mapq_for_supplementary(self, value: bool) -> None: ...
+    @property
+    def with_xb_tag(self) -> bool: ...
+    @with_xb_tag.setter
+    def with_xb_tag(self, value: bool) -> None: ...
+    @property
+    def max_occurrences(self) -> int: ...
+    @max_occurrences.setter
+    def max_occurrences(self, value: int) -> None: ...
+    @property
+    def off_diagonal_x_dropoff(self) -> float: ...
+    @off_diagonal_x_dropoff.setter
+    def off_diagonal_x_dropoff(self, value: float) -> None: ...
+    @property
+    def ignore_alternate_contigs(self) -> bool: ...
+    @ignore_alternate_contigs.setter
+    def ignore_alternate_contigs(self, value: bool) -> None: ...
+    @property
+    def internal_seed_split_factor(self) -> float: ...
+    @internal_seed_split_factor.setter
+    def internal_seed_split_factor(self, value: float) -> None: ...
+    @property
+    def drop_chain_fraction(self) -> float: ...
+    @drop_chain_fraction.setter
+    def drop_chain_fraction(self, value: float) -> None: ...
+    @property
+    def max_mate_rescue_rounds(self) -> int: ...
+    @max_mate_rescue_rounds.setter
+    def max_mate_rescue_rounds(self, value: int) -> None: ...
+    @property
+    def min_seeded_bases_in_chain(self) -> int: ...
+    @min_seeded_bases_in_chain.setter
+    def min_seeded_bases_in_chain(self, value: int) -> None: ...
+    @property
+    def seed_occurrence_in_3rd_round(self) -> int: ...
+    @seed_occurrence_in_3rd_round.setter
+    def seed_occurrence_in_3rd_round(self, value: int) -> None: ...
+    @property
+    def xa_max_hits(self) -> int | tuple[int, int]: ...
+    @xa_max_hits.setter
+    def xa_max_hits(self, value: int | tuple[int, int]) -> None: ...
+    @property
+    def xa_drop_ratio(self) -> float: ...
+    @xa_drop_ratio.setter
+    def xa_drop_ratio(self, value: float) -> None: ...
+    @property
+    def gap_open_penalty(self) -> int | tuple[int, int]: ...
+    @gap_open_penalty.setter
+    def gap_open_penalty(self, value: int | tuple[int, int]) -> None: ...
+    @property
+    def gap_extension_penalty(self) -> int | tuple[int, int]: ...
+    @gap_extension_penalty.setter
+    def gap_extension_penalty(self, value: int | tuple[int, int]) -> None: ...
+    @property
+    def clipping_penalty(self) -> int | tuple[int, int]: ...
+    @clipping_penalty.setter
+    def clipping_penalty(self, value: int | tuple[int, int]) -> None: ...
+    def finalize(self, copy: bool = False) -> BwaMemOptions: ...
 
 class BwaMem:
     _index: BwaIndex