bioforge 3.0.0__tar.gz → 3.1.0__tar.gz

{bioforge-3.0.0 → bioforge-3.1.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: bioforge
-Version: 3.0.0
+Version: 3.1.0
 Summary: High-performance bioinformatics engine for Edge Computing — 5-bit encoding, vectorised NumPy core, optional C backend
 Author-email: Aarón Aranda Torrijos <noe9601@gmail.com>
 License: # PolyForm Noncommercial License 1.0.0

{bioforge-3.0.0 → bioforge-3.1.0}/bioforge/__init__.py

@@ -36,7 +36,7 @@ from .biocore import (
 from .genomemap import GenomeAligner, Mapping
 from .smart_translator import SmartTranslator
 
-__version__ = "3.0.0"
+__version__ = "3.1.0"
 __author__  = "Aarón Aranda Torrijos"
 
 __all__ = [

{bioforge-3.0.0 → bioforge-3.1.0}/bioforge/aligner.py

@@ -179,6 +179,7 @@ class SequenceAligner:
         seq_b: PackedSequence,
         mode:  Literal['global', 'semi-global'] = 'global',
         band:  int | None = None,
+        detect_mutations: bool = True,
     ) -> AlignmentResult:
         """
         Align seq_a (reference) against seq_b (query) — Needleman-Wunsch.
@@ -238,7 +239,8 @@ class SequenceAligner:
                 raise AlignmentError(f"band debe ser ≥ 1, se recibió {band}.")
             if C_AVAILABLE:
                 return cls._align_banded_c(
-                    codes_a, codes_b, m, n, band, mode, seq_a.seq_type
+                    codes_a, codes_b, m, n, band, mode, seq_a.seq_type,
+                    detect_mutations,
                 )
             # NumPy fallback: banded fill sobre matriz completa
             if m > cls._MAX_SAFE_LEN or n > cls._MAX_SAFE_LEN:
@@ -261,7 +263,8 @@ class SequenceAligner:
             )
 
         if C_AVAILABLE:
-            return cls._align_c(codes_a, codes_b, m, n, mode, seq_a.seq_type)
+            return cls._align_c(codes_a, codes_b, m, n, mode, seq_a.seq_type,
+                                detect_mutations)
 
         H = cls._fill_matrix(codes_a, codes_b, m, n, mode)
         return cls._traceback(H, codes_a, codes_b, m, n, mode, seq_a.seq_type)
@@ -361,6 +364,7 @@ class SequenceAligner:
         band: int,
         mode: str,
         seq_type: SeqType,
+        detect_mutations: bool = True,
     ) -> AlignmentResult:
         decode_bytes = (
             _NUC_DECODE_BYTES if seq_type == SeqType.NUCLEOTIDE else _AA_DECODE_BYTES
@@ -370,7 +374,7 @@ class SequenceAligner:
             int(cls.MATCH), int(cls.MISMATCH), int(cls.GAP),
             band, mode,
         )
-        mutations = cls._detect_mutations(aligned_a, aligned_b)
+        mutations = cls._detect_mutations(aligned_a, aligned_b) if detect_mutations else []
         aln_len   = n_matches + n_mismatches + n_gaps
         identity  = n_matches / aln_len if aln_len else 0.0
         return AlignmentResult(
@@ -389,6 +393,7 @@ class SequenceAligner:
         n: int,
         mode: str,
         seq_type: SeqType,
+        detect_mutations: bool = True,
     ) -> AlignmentResult:
         decode_bytes = (
             _NUC_DECODE_BYTES if seq_type == SeqType.NUCLEOTIDE else _AA_DECODE_BYTES
@@ -398,7 +403,7 @@ class SequenceAligner:
             int(cls.MATCH), int(cls.MISMATCH), int(cls.GAP),
             mode,
         )
-        mutations = cls._detect_mutations(aligned_a, aligned_b)
+        mutations = cls._detect_mutations(aligned_a, aligned_b) if detect_mutations else []
         aln_len   = n_matches + n_mismatches + n_gaps
         identity  = n_matches / aln_len if aln_len else 0.0
         return AlignmentResult(

{bioforge-3.0.0 → bioforge-3.1.0}/bioforge/genomemap.py

@@ -308,7 +308,8 @@ def _extend(ref: str, read: str, ch: Chain) -> Optional[Mapping]:
     band = min(256, abs(len(ref_sub) - len(read_sub)) + 32)
     try:
         res = SequenceAligner.align(_pack(ref_sub), _pack(read_sub),
-                                    mode="global", band=band)
+                                    mode="global", band=band,
+                                    detect_mutations=False)   # el mapeo no las usa
         cigar, n_match, block = _cigar(res.aligned_a, res.aligned_b)
         identity = res.identity
     except Exception:                       # noqa: BLE001 — extensión best-effort

{bioforge-3.0.0 → bioforge-3.1.0}/bioforge.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: bioforge
-Version: 3.0.0
+Version: 3.1.0
 Summary: High-performance bioinformatics engine for Edge Computing — 5-bit encoding, vectorised NumPy core, optional C backend
 Author-email: Aarón Aranda Torrijos <noe9601@gmail.com>
 License: # PolyForm Noncommercial License 1.0.0

{bioforge-3.0.0 → bioforge-3.1.0}/tests/test_aligner.py

@@ -44,6 +44,23 @@ def _prot(seq: str, header: str = "h") -> PackedSequence:
 # §1  PROPIEDADES MATEMÁTICAS
 # ══════════════════════════════════════════════════════════════════════════════
 
+def test_detect_mutations_flag():
+    """detect_mutations=False omite la lista de mutaciones sin tocar el resto."""
+    a = "ACGTACGTACGTACGTACGT"
+    b = "ACGTACGTTCGTACGTACGT"          # una sustitución
+    full = SequenceAligner.align(_nuc(a), _nuc(b))
+    fast = SequenceAligner.align(_nuc(a), _nuc(b), detect_mutations=False)
+    assert len(full.mutations) == 1
+    assert fast.mutations == []
+    # todo lo demás idéntico (identidad, score, matches)
+    assert fast.identity == full.identity
+    assert fast.score == full.score
+    assert fast.n_matches == full.n_matches
+    # y también en la ruta banded (la que usa el mapeador)
+    fastb = SequenceAligner.align(_nuc(a), _nuc(b), band=6, detect_mutations=False)
+    assert fastb.mutations == []
+
+
 def test_identical_sequences_perfect_score():
     """Dos secuencias idénticas → score máximo, identidad 1.0, 0 mutaciones."""
     seq = "ATGGTGCACCTGACTCCTGAGGAGAAGTCTGCC"