pen-stack 6.10.3__tar.gz → 6.10.4__tar.gz

{pen_stack-6.10.3 → pen_stack-6.10.4}/CHANGELOG.md

@@ -3,6 +3,30 @@
 All notable changes to PEN-STACK are documented here. This file follows
 [Keep a Changelog](https://keepachangelog.com/) and the program's phase structure.
 
+## [6.10.4] - 2026-06-20 - Chromatin incremental-value test (annotation, not a re-ranker) — chromatin work complete
+
+**PATCH — the final chromatin question, answered.** v6.10.3 validated accessibility as a moderate standalone
+predictor; v6.10.4 tests whether it adds **incremental** value *on top of* the CRISOT sequence score — i.e. whether
+it should re-rank.
+
+### Added
+- **Incremental-value analysis** (`scripts/offtarget_chromatin_incremental.py`; result in
+  `benchmarks/offtarget/chromatin_incremental.json`). On GUIDE-seq (HEK293T-matched), per off-target CRISOT +
+  accessibility: (A) a conditional logistic regression `active ~ z(CRISOT) + z(accessibility)` and (B) a
+  leave-one-guide-out held-out AUPRC of CRISOT-only vs a CRISOT+accessibility combiner — tested at two candidate
+  imbalances (1:16 and a realistic 1:123). **Result:** accessibility carries a **small, real conditional signal**
+  (coef ≈ 0.35, bootstrap CI excludes 0 at both imbalances) but adds **NO held-out ranking improvement** over CRISOT
+  (AUPRC gap CI includes 0 at both: −0.0025 [−0.011, +0.005] and +0.0027 [−0.014, +0.021]).
+
+### Decision
+- **Chromatin is a validated ANNOTATION, NOT a re-ranker.** The CRISOT sequence score already captures the
+  practically-relevant nomination ranking; a CRISOT+accessibility combiner does not improve held-out AUPRC, so it is
+  **not** wired into the numeric risk score (`CHROMATIN_VALIDATION.changes_numeric_risk_score = False`). The fitted
+  combiner coefficients are recorded for transparency but intentionally not applied. Ledger + docs updated.
+- **This completes the chromatin context work** (v6.10.1 wired → v6.10.2 cross-cell weak → v6.10.3 cell-type-matched
+  VALIDATED moderate → v6.10.4 incremental tested, annotation-only). Nothing about chromatin is now open or
+  undisclosed.
+
 ## [6.10.3] - 2026-06-20 - Chromatin context: cell-type-matched validation (VALIDATED, moderate)
 
 **PATCH — the definitive chromatin test.** v6.10.2's controlled experiment was ambiguous (GUIDE-seq positive, TTISS

{pen_stack-6.10.3 → pen_stack-6.10.4}/CITATION.cff

@@ -1,7 +1,7 @@
 cff-version: 1.2.0
 message: "If you use PEN-STACK, please cite it as below."
 title: "PEN-STACK: open infrastructure for genome writing"
-version: 6.10.3
+version: 6.10.4
 date-released: 2026-06-20
 authors:
   - family-names: "Mahaboob Ali"

{pen_stack-6.10.3 → pen_stack-6.10.4}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pen-stack
-Version: 6.10.3
+Version: 6.10.4
 Summary: Open infrastructure for genome writing: the Writable Genome atlas, the Writer Atlas, and the Write Planner.
 Author-email: Anees Ahmed Mahaboob Ali <ahmedaneesm@gmail.com>
 License: MIT
@@ -90,7 +90,7 @@ every design against rule-grounded mechanism, reports calibrated confidence, cit
 [![codecov](https://codecov.io/gh/ahmedanees-m/pen-stack/branch/main/graph/badge.svg)](https://codecov.io/gh/ahmedanees-m/pen-stack)
 [![License: MIT](https://img.shields.io/badge/License-MIT-informational.svg)](LICENSE)
 [![Python 3.11+](https://img.shields.io/badge/python-3.11%2B-blue.svg)](https://www.python.org/)
-[![Version](https://img.shields.io/badge/version-6.10.3-blue.svg)](CHANGELOG.md)
+[![Version](https://img.shields.io/badge/version-6.10.4-blue.svg)](CHANGELOG.md)
 [![Status](https://img.shields.io/badge/status-1.0%20First%20Stable-success.svg)](docs/STABILITY.md)
 [![Tests](https://img.shields.io/badge/tests-378%20passing-success.svg)](tests/)
 [![Lint: ruff](https://img.shields.io/badge/lint-ruff-purple.svg)](https://github.com/astral-sh/ruff)

{pen_stack-6.10.3 → pen_stack-6.10.4}/README.md

@@ -15,7 +15,7 @@ every design against rule-grounded mechanism, reports calibrated confidence, cit
 [![codecov](https://codecov.io/gh/ahmedanees-m/pen-stack/branch/main/graph/badge.svg)](https://codecov.io/gh/ahmedanees-m/pen-stack)
 [![License: MIT](https://img.shields.io/badge/License-MIT-informational.svg)](LICENSE)
 [![Python 3.11+](https://img.shields.io/badge/python-3.11%2B-blue.svg)](https://www.python.org/)
-[![Version](https://img.shields.io/badge/version-6.10.3-blue.svg)](CHANGELOG.md)
+[![Version](https://img.shields.io/badge/version-6.10.4-blue.svg)](CHANGELOG.md)
 [![Status](https://img.shields.io/badge/status-1.0%20First%20Stable-success.svg)](docs/STABILITY.md)
 [![Tests](https://img.shields.io/badge/tests-378%20passing-success.svg)](tests/)
 [![Lint: ruff](https://img.shields.io/badge/lint-ruff-purple.svg)](https://github.com/astral-sh/ruff)

{pen_stack-6.10.3 → pen_stack-6.10.4}/docs/DEVIATIONS_AND_DISCLOSURES.md

@@ -16,7 +16,7 @@ applicable (e.g. single contributor).
 |---|---|---|
 | Nuclease predictor named in plan (CCLMoff / CRISMER) vs shipped (CRISOT) | **DISCLOSED** | We use **CRISOT-Score** (Chen 2023, CC-BY-NC, CPU), not the plan's CCLMoff/CRISMER. Justified: **CRISMER ships no license** (cannot be redistributed/wrapped); **CCLMoff** is a GPU RNA-language-model stack AND was trained on these very assays (evaluating it on them would be **leakage**), whereas **CRISOT-Score is MD-physics, assay-agnostic → a leakage-clean held-out evaluation**. CRISOT is the more rigorous and license-appropriate choice. A genuine substitution, disclosed. |
 | Assay coverage (plan listed GUIDE/CIRCLE/CHANGE/SITE/SURRO/TTISS/Digenome-seq) | **FIXED (v6.10.1)** | v6.10.0 used GUIDE-seq + CIRCLE-seq only. v6.10.1 adds **CHANGE-seq + SITE-seq** (independent broad guide panels) — CRISOT beats homology on **all four** (per-guide bootstrap CI excludes 0). SURRO-seq is targeted/biased (kept as an orthogonal reference, not genome-wide truth); TTISS/Digenome-seq remain **DISCLOSED** as not-yet-folded-in. |
-| "Chromatin-aware via the Stage B feature store" | **VALIDATED (moderate, cell-type-matched; v6.10.3)** | v6.10.0 overclaimed "chromatin-aware" with only a caller hook → v6.10.1 added a real `locus_accessibility()` (verified on the VM) → v6.10.2 ran a controlled test with the cross-cell K562 proxy and got a *weak/inconsistent* result (GUIDE-seq 0.58, TTISS 0.346, in-vitro null) → **v6.10.3 settled it with a CELL-TYPE-MATCHED track** (ENCODE HEK293T DNase, ENCFF529BOG; matches GUIDE-seq/TTISS). **Cell-type matching lifts the canonical WT-Cas9 cell-based assay GUIDE-seq from AUROC 0.58 → 0.671 (CI [0.642, 0.701])**, the in-vitro control stays null (0.494). **VERDICT: VALIDATED (moderate, cell-type-matched) for WT-Cas9 cell-based off-target activity** — the cross-cell K562 proxy was dampening a real effect. Honest caveats: the effect is **moderate** (the sequence/CRISOT score still dominates); it does **not** transfer to TTISS (0.383, an expected outlier — a Cas9-*variant* specificity assay); and chromatin is surfaced as a validated **annotation** that does **not yet change the numeric risk score** (a calibrated CRISOT+accessibility combination is the remaining refinement, deferred). Full result: `benchmarks/offtarget/chromatin_validation.json`; reproducible scripts `scripts/offtarget_chromatin_{validation,matched}.py`. |
+| "Chromatin-aware via the Stage B feature store" | **VALIDATED (moderate, cell-type-matched; v6.10.3)** | v6.10.0 overclaimed "chromatin-aware" with only a caller hook → v6.10.1 added a real `locus_accessibility()` (verified on the VM) → v6.10.2 ran a controlled test with the cross-cell K562 proxy and got a *weak/inconsistent* result (GUIDE-seq 0.58, TTISS 0.346, in-vitro null) → **v6.10.3 settled it with a CELL-TYPE-MATCHED track** (ENCODE HEK293T DNase, ENCFF529BOG; matches GUIDE-seq/TTISS). **Cell-type matching lifts the canonical WT-Cas9 cell-based assay GUIDE-seq from AUROC 0.58 → 0.671 (CI [0.642, 0.701])**, the in-vitro control stays null (0.494). **VERDICT: VALIDATED (moderate, cell-type-matched) for WT-Cas9 cell-based off-target activity** — the cross-cell K562 proxy was dampening a real effect. Honest caveats: the effect is **moderate** (the sequence/CRISOT score dominates); it does **not** transfer to TTISS (0.383, an expected outlier — a Cas9-*variant* specificity assay). **v6.10.4 closed the last piece** — does accessibility add *incremental* value over CRISOT? On GUIDE-seq (HEK293T-matched), accessibility carries a **small real conditional signal** (logreg coef ~0.35, CI excludes 0 at two imbalances) but adds **NO held-out ranking improvement** over CRISOT (AUPRC gap CI includes 0 at both). **Decision: chromatin is a validated ANNOTATION, NOT a re-ranker** — it does **not** change the numeric risk score (CRISOT captures the ranking); the fitted combiner is recorded but intentionally not applied. Full results: `benchmarks/offtarget/chromatin_{validation,incremental}.json`; reproducible scripts `scripts/offtarget_chromatin_{validation,matched,incremental}.py`. **This is the final, complete state — nothing about chromatin is now open or undisclosed.** |
 | Held-out **guide** AND **locus** splits | **DISCLOSED** | Held-out-**guide** is implemented; a genomic-coordinate **locus** split is **not possible** — the harmonized assay data ships sequences, not coordinates. Instead we add **cross-assay generalization** (the assay-agnostic CRISOT-Score evaluated on GUIDE/CIRCLE canonical guides + CHANGE/SITE independent panels). This is the leakage-clean substitute and is stated in `benchmarks/offtarget/split.json`. |
 | Learned **integrase** off-target scorer (plan: "learned on HIDE/Cryptic-seq") | **DISCLOSED** | Shipped as a documented **pseudo-attB cryptic scan** on the real Bxb1 attB core, not a learned model — Cryptic-seq/HIDE-seq are recent preprints and IntQuery has **no public weights**, so a real learned integrase predictor is not groundable. Flagged extrapolative; IntQuery cited as paper-only. |
 | Bench "joins the Challenge" | **FIXED (v6.10.1)** | An `offtarget` nomination task is added to the Genome-Writing Challenge (`benchmarks/genome_writing_challenge/harness.py`) — non-circular (label = wet-lab Active call), data-gated on the fixture. |

{pen_stack-6.10.3 → pen_stack-6.10.4}/docs/cards/offtarget_data.md

@@ -61,10 +61,16 @@ AUROC of accessibility for active-vs-inactive off-targets per assay, with in-vit
 DNase track (`ENCFF529BOG`, matching GUIDE-seq's HEK293 and TTISS's HEK293T): **cell-type matching lifts the
 canonical WT-Cas9 cell-based assay from 0.58 → 0.671** (the cross-cell proxy was dampening a real effect), with the
 in-vitro control still null. **Verdict: VALIDATED (moderate, cell-type-matched)** for WT-Cas9 cell-based off-target
-activity. Honest caveats: the effect is **moderate** (the sequence/CRISOT score still dominates nomination); it does
-**not** transfer to TTISS (a Cas9-variant specificity assay, the expected outlier); and chromatin is surfaced as a
-validated **annotation** that does **not yet change the numeric risk score** (a calibrated CRISOT+accessibility
-combination is the remaining refinement). Reproducible: `scripts/offtarget_chromatin_{validation,matched}.py`.
+activity. Honest caveats: the effect is **moderate** (the sequence/CRISOT score dominates nomination); it does
+**not** transfer to TTISS (a Cas9-variant specificity assay, the expected outlier).
+
+**Incremental value over CRISOT (v6.10.4, `chromatin_incremental.json`):** on GUIDE-seq (HEK293T-matched),
+accessibility carries a **small real conditional signal** (logistic-regression coefficient ~0.35, bootstrap CI
+excludes 0 at both 1:16 and 1:123 candidate imbalance) but adds **NO held-out ranking improvement** over CRISOT
+(leave-one-guide-out AUPRC gap CI includes 0 at both). **Decision: chromatin is a validated ANNOTATION, NOT a
+re-ranker** — it does **not** change the numeric risk score (CRISOT already captures the practically-relevant ranking
+signal); the fitted CRISOT+accessibility combiner is recorded but intentionally not applied. Reproducible:
+`scripts/offtarget_chromatin_{validation,matched,incremental}.py`.
 
 ## Risk calibration (grounded)
 The nomination risk band IS the empirical fraction of candidates at *k* mismatches that were validated-active

{pen_stack-6.10.3 → pen_stack-6.10.4}/pen_stack/__init__.py

@@ -1,2 +1,2 @@
 """PEN-STACK v3.0 - open infrastructure for genome writing."""
-__version__ = "6.10.3"
+__version__ = "6.10.4"

{pen_stack-6.10.3 → pen_stack-6.10.4}/pen_stack/wgenome/offtarget_data.py

@@ -75,12 +75,24 @@ CHROMATIN_VALIDATION = {
         "siteseq": {"modality": "in_vitro_control", "auroc": 0.494, "k562_cross_cell": 0.469, "ci95": [0.475, 0.514]},
     },
     "matched_track": "ENCODE HEK293T DNase-seq (ENCFF529BOG, GRCh38)",
+    # v6.10.4 — does accessibility add INCREMENTAL value over the CRISOT sequence score? Tested at two imbalances:
+    # a small REAL conditional signal (logreg acc coef ~0.35, CI excludes 0) but NO held-out ranking improvement
+    # (CRISOT+acc AUPRC gap CI includes 0). DECISION: annotation, NOT a re-ranker. Full result:
+    # benchmarks/offtarget/chromatin_incremental.json.
+    "incremental_over_crisot": {
+        "conditional_acc_coef": 0.351, "conditional_acc_coef_ci95": [0.2385, 0.5584], "adds_conditional_signal": True,
+        "heldout_auprc_gap": 0.0027, "heldout_auprc_gap_ci95": [-0.014, 0.0214], "improves_ranking": False,
+        "decision": "annotation only — accessibility carries a small real conditional signal but does NOT improve "
+                    "held-out nomination ranking over the CRISOT sequence score (at realistic imbalance); a "
+                    "re-ranking combiner is NOT wired in (no demonstrated benefit).",
+    },
+    "changes_numeric_risk_score": False,
     "note": "cell-type-matched accessibility predicts WT-Cas9 cell-based off-target activity (GUIDE-seq AUROC "
             "0.58 cross-cell -> 0.671 matched; in-vitro control null -> method sound). MODERATE effect (the "
             "sequence/CRISOT score dominates nomination) and CELL-TYPE-SPECIFIC. Does NOT transfer to the "
-            "Cas9-VARIANT assay TTISS (0.383, the expected outlier). Surfaced as a validated qualitative annotation; "
-            "it does NOT yet change the numeric risk score (a calibrated CRISOT+accessibility combination is the "
-            "remaining refinement).",
+            "Cas9-VARIANT assay TTISS (0.383, the expected outlier). v6.10.4 tested the incremental value over "
+            "CRISOT: a small real conditional signal but NO held-out ranking improvement -> chromatin is surfaced "
+            "as a VALIDATED ANNOTATION and does NOT change the numeric risk score (CRISOT captures the ranking).",
 }
 
 # ---- Off-Target-Bench headline (REAL full-data result; per-guide AUPRC, held-out-guide bootstrap CI) -----

{pen_stack-6.10.3 → pen_stack-6.10.4}/pen_stack/wgenome/offtarget_predict.py

@@ -103,10 +103,11 @@ def _chromatin_modifier(accessibility: float | None = None, locus_acc: dict | No
     Stage B accessibility (`locus_acc`) when available, else a caller-supplied 0..1 scalar; abstains (None) when
     neither is given. It is QUALITATIVE and does NOT change the numeric risk score — a controlled test on our data
     found the signal weak/inconsistent (see `offtarget_data.CHROMATIN_VALIDATION`)."""
-    val = {"validated": True, "effect_size": "moderate",
-           "validation": "validated (moderate, cell-type-matched): GUIDE-seq off-target AUROC 0.58 cross-cell -> "
-           "0.671 with matched HEK293T DNase (CI [0.642,0.701]); in-vitro control null. Annotation only — does not "
-           "yet change the numeric risk score (sequence/CRISOT dominates).",
+    val = {"validated": True, "effect_size": "moderate", "changes_risk_score": False,
+           "validation": "validated standalone (moderate, cell-type-matched: GUIDE-seq off-target AUROC 0.58 "
+           "cross-cell -> 0.671 matched HEK293T DNase). ANNOTATION ONLY — v6.10.4 tested the incremental value over "
+           "CRISOT: small real conditional signal but NO held-out ranking improvement, so it does NOT change the "
+           "numeric risk score (CRISOT sequence score captures the ranking).",
            "doi": "10.1038/s41587-020-0555-7"}
     if locus_acc is not None:
         raises = bool(locus_acc.get("accessible"))

{pen_stack-6.10.3 → pen_stack-6.10.4}/pen_stack.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pen-stack
-Version: 6.10.3
+Version: 6.10.4
 Summary: Open infrastructure for genome writing: the Writable Genome atlas, the Writer Atlas, and the Write Planner.
 Author-email: Anees Ahmed Mahaboob Ali <ahmedaneesm@gmail.com>
 License: MIT
@@ -90,7 +90,7 @@ every design against rule-grounded mechanism, reports calibrated confidence, cit
 [![codecov](https://codecov.io/gh/ahmedanees-m/pen-stack/branch/main/graph/badge.svg)](https://codecov.io/gh/ahmedanees-m/pen-stack)
 [![License: MIT](https://img.shields.io/badge/License-MIT-informational.svg)](LICENSE)
 [![Python 3.11+](https://img.shields.io/badge/python-3.11%2B-blue.svg)](https://www.python.org/)
-[![Version](https://img.shields.io/badge/version-6.10.3-blue.svg)](CHANGELOG.md)
+[![Version](https://img.shields.io/badge/version-6.10.4-blue.svg)](CHANGELOG.md)
 [![Status](https://img.shields.io/badge/status-1.0%20First%20Stable-success.svg)](docs/STABILITY.md)
 [![Tests](https://img.shields.io/badge/tests-378%20passing-success.svg)](tests/)
 [![Lint: ruff](https://img.shields.io/badge/lint-ruff-purple.svg)](https://github.com/astral-sh/ruff)

{pen_stack-6.10.3 → pen_stack-6.10.4}/pen_stack.egg-info/SOURCES.txt

@@ -480,6 +480,7 @@ prereg/ws_writer.yaml
 prereg/ws_wv.yaml
 scripts/calibrate_immune_axes.py
 scripts/fetch_licensed_sources.py
+scripts/offtarget_chromatin_incremental.py
 scripts/offtarget_chromatin_matched.py
 scripts/offtarget_chromatin_validation.py
 scripts/p1_build_atlas.py

{pen_stack-6.10.3 → pen_stack-6.10.4}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "pen-stack"
-version = "6.10.3"
+version = "6.10.4"
 description = "Open infrastructure for genome writing: the Writable Genome atlas, the Writer Atlas, and the Write Planner."
 readme = "README.md"
 requires-python = ">=3.11"

pen_stack-6.10.4/scripts/offtarget_chromatin_incremental.py

@@ -0,0 +1,148 @@
+"""Reproducible analysis (v6.10.4): does chromatin accessibility add INCREMENTAL value over the CRISOT sequence
+score for nominating off-targets? On GUIDE-seq (cell-based) with cell-type-matched HEK293T DNase accessibility:
+
+  (A) conditional logistic regression  active ~ z(CRISOT) + z(accessibility)  -> bootstrap-over-guides 95% CI on
+      the accessibility coefficient (does accessibility carry signal CONDITIONAL on CRISOT?).
+  (B) leave-one-guide-out held-out AUPRC of CRISOT-only vs a CRISOT+accessibility logistic combiner, bootstrap
+      95% CI on the per-guide gap (does the combiner IMPROVE held-out nomination ranking?).
+
+RESULT (committed in benchmarks/offtarget/chromatin_incremental.json): accessibility carries a SMALL, REAL
+conditional signal (coef ~0.35, CI excludes 0 at both 1:16 and 1:123 imbalance) but adds NO held-out ranking
+improvement over CRISOT (AUPRC gap CI includes 0 at both). DECISION: chromatin is a validated ANNOTATION, NOT a
+re-ranker — the CRISOT sequence score already captures the practically-relevant ranking signal.
+
+Needs: the CRISOT repo (CRISOT-Score), hg38 GRCh38.fa, the HEK293T DNase bigWig, the GUIDE-seq CSV, plus
+xgboost / pyBigWig / scikit-learn. Run on the VM. (Inactive cap controls the candidate imbalance.)
+"""
+from __future__ import annotations
+
+import json
+import os
+import subprocess
+import sys
+
+import numpy as np
+import pandas as pd
+import pyBigWig
+from sklearn.linear_model import LogisticRegression
+from sklearn.metrics import average_precision_score
+
+sys.path.insert(0, os.environ.get("CRISOT", "/crisot"))
+from crisot_modules import CRISOT  # noqa: E402
+from utils import load_pkl  # noqa: E402
+
+FA = os.environ.get("FA", "/ref/GRCh38.fa")
+BW = os.environ.get("BW", "/ref/hek293t_dnase.bigWig")
+DS = os.environ.get("DS", "/d")
+INACT_CAP = int(os.environ.get("INACT_CAP", "4000"))   # 4000/guide ~ realistic imbalance
+COMP = str.maketrans("ACGT", "TGCA")
+VALID = {f"chr{i}" for i in list(range(1, 23)) + ["X", "Y"]}
+
+
+def rc(s: str) -> str:
+    return s.translate(COMP)[::-1]
+
+
+def main() -> dict:
+    rng = np.random.RandomState(7)
+    pr, _ab, _bins, _w = load_pkl(os.environ.get("CRISOT", "/crisot") + "/models/crisot_score_param.pkl")
+    model = CRISOT(param=pr, ref_genome=os.environ.get("CRISOT", "/crisot") + "/script/hg38.na")
+    df = pd.read_csv(f"{DS}/guideseq.csv")
+    df["On20"] = df["On"].str[:20]
+    g = df.groupby("On20")["Active"].sum()
+    df = df[df["On20"].isin(list(g[g >= 5].index))].copy()
+    parts = []
+    for _gd, sub in df.groupby("On20"):
+        inact = sub[sub.Active == 0]
+        parts.append(pd.concat([sub[sub.Active == 1], inact.sample(n=min(INACT_CAP, len(inact)), random_state=rng)]))
+    d = pd.concat(parts).reset_index(drop=True)
+    d["crisot"] = model.score(data_df=d, On="On", Off="Off")
+    d["off23"] = d["Off"].str[:23]
+    pat = {}
+    for off in d["off23"].unique():
+        if len(off) == 23 and set(off) <= set("ACGT"):
+            pat[off] = off
+            pat[rc(off)] = off
+    open("/tmp/pats.txt", "w").write("\n".join(pat) + "\n")
+    hits: dict = {}
+    cur, buf = None, []
+
+    def flush(chrom, seq):
+        c = chrom if chrom and chrom.startswith("chr") else f"chr{chrom}"
+        if not seq or c not in VALID:
+            return
+        open("/tmp/chr.txt", "w").write(seq)
+        p = subprocess.run(["grep", "-boFf", "/tmp/pats.txt", "/tmp/chr.txt"], capture_output=True, text=True)
+        for ln in p.stdout.splitlines():
+            ob, _, m = ln.partition(":")
+            if pat.get(m):
+                hits.setdefault(pat[m], []).append((c, int(ob)))
+    with open(FA) as fh:
+        for line in fh:
+            if line.startswith(">"):
+                if cur:
+                    flush(cur, "".join(buf))
+                cur, buf = line[1:].split()[0], []
+            else:
+                buf.append(line.strip().upper())
+        if cur:
+            flush(cur, "".join(buf))
+    bw = pyBigWig.open(BW)
+    bwc = set(bw.chroms().keys())
+    cache: dict = {}
+
+    def acc(off):
+        vals = []
+        for c, p in hits.get(off, []):
+            bc = c if c in bwc else (c[3:] if c[3:] in bwc else None)
+            if bc is None:
+                continue
+            b0 = (p // 1000) * 1000
+            if (bc, b0) not in cache:
+                try:
+                    v = bw.stats(bc, b0, min(b0 + 1000, bw.chroms(bc)), type="mean")[0]
+                    cache[(bc, b0)] = float(v) if v is not None else None
+                except Exception:  # noqa: BLE001
+                    cache[(bc, b0)] = None
+            if cache[(bc, b0)] is not None:
+                vals.append(cache[(bc, b0)])
+        return max(vals) if vals else None
+    d["acc"] = d["off23"].map(acc)
+    d = d.dropna(subset=["acc"]).reset_index(drop=True)
+    d["zc"] = (d.crisot - d.crisot.mean()) / d.crisot.std()
+    d["za"] = (d.acc - d.acc.mean()) / d.acc.std()
+    guides = sorted(d.On20.unique())
+
+    def fit(dd):
+        lr = LogisticRegression(max_iter=2000, class_weight="balanced").fit(dd[["zc", "za"]], dd.Active)
+        return lr.coef_[0], lr.intercept_[0]
+    (cc, ca), _b = fit(d)
+    accb = []
+    for _ in range(1000):
+        gs = rng.choice(guides, len(guides), replace=True)
+        try:
+            accb.append(fit(pd.concat([d[d.On20 == x] for x in gs]))[0][1])
+        except Exception:  # noqa: BLE001
+            pass
+    acc_ci = [round(float(np.percentile(accb, 2.5)), 4), round(float(np.percentile(accb, 97.5)), 4)]
+    cris_ap, comb_ap = [], []
+    for held in guides:
+        tr, te = d[d.On20 != held], d[d.On20 == held]
+        if te.Active.sum() == 0 or len(te) < 5:
+            continue
+        lr = LogisticRegression(max_iter=2000, class_weight="balanced").fit(tr[["zc", "za"]], tr.Active)
+        comb_ap.append(average_precision_score(te.Active, lr.predict_proba(te[["zc", "za"]])[:, 1]))
+        cris_ap.append(average_precision_score(te.Active, te.zc.values))
+    gap = np.array(comb_ap) - np.array(cris_ap)
+    gb = [float(np.mean(gap[rng.randint(0, len(gap), len(gap))])) for _ in range(1000)]
+    out = {"n_offtargets": int(len(d)), "n_actives": int(d.Active.sum()), "n_guides": len(guides),
+           "conditional_acc_coef": round(float(ca), 4), "acc_coef_ci95": acc_ci,
+           "mean_crisot_only_auprc": round(float(np.mean(cris_ap)), 4),
+           "mean_crisot_plus_acc_auprc": round(float(np.mean(comb_ap)), 4),
+           "auprc_gap_ci95": [round(float(np.percentile(gb, 2.5)), 4), round(float(np.percentile(gb, 97.5)), 4)]}
+    print(json.dumps(out, indent=2))
+    return out
+
+
+if __name__ == "__main__":
+    main()