prographs 0.12.0 → 0.13.0

package/dist/{chunk-3BYXOISI.cjs → chunk-Q3HB3IX6.cjs}

@@ -22,7 +22,7 @@ var _chunkEEF5WRRZcjs = require('./chunk-EEF5WRRZ.cjs');
 var _chunkPY4R3IFFcjs = require('./chunk-PY4R3IFF.cjs');
 
 
-var _chunkUMEKIIHJcjs = require('./chunk-UMEKIIHJ.cjs');
+var _chunkLNWTC72Ycjs = require('./chunk-LNWTC72Y.cjs');
 
 // src/charts/lifescience/DoseResponseFit.ts
 var _d3 = require('d3');
@@ -215,7 +215,7 @@ var STRAY = { compound: "Compound B", dose: 10, replicate: 1, shift: 45 };
 var DoseResponseSamples = class {
   /** `{ compound, dose, response }`, three replicates per dose. */
   static assay(seed) {
-    const { n } = _chunkUMEKIIHJcjs.SampleData.rng(seed);
+    const { n } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
     return COMPOUNDS.flatMap(({ compound, bottom, top, ec50, hill }) => DOSES.flatMap((dose) => Array.from({ length: REPLICATES }, (_, replicate) => {
       const model = bottom + (top - bottom) * (dose ** hill / (dose ** hill + ec50 ** hill));
       const stray = compound === STRAY.compound && dose === STRAY.dose && replicate === STRAY.replicate ? STRAY.shift : 0;
@@ -508,7 +508,7 @@ var MaPlotSamples = class _MaPlotSamples {
   }
   /** `{ gene, red, green, baseMean, a, m, pvalue }`, with no p-value for the faintest genes, as after independent filtering. */
   static expression(seed, count = 1200) {
-    const z = _chunkUMEKIIHJcjs.SampleData.distribution(seed, count * DRAWS, 1).map((row) => (row.value - 40) / 8);
+    const z = _chunkLNWTC72Ycjs.SampleData.distribution(seed, count * DRAWS, 1).map((row) => (row.value - 40) / 8);
     return Array.from({ length: count }, (_, i) => {
       const [level, noise, pick, size] = z.slice(i * DRAWS, (i + 1) * DRAWS);
       const a = Math.max(3, Math.min(15, 8.5 + 2.2 * level));
@@ -904,7 +904,7 @@ var QUIET = 0.1;
 var OncoPrintSamples = class {
   /** `{ sample, gene, alteration, source }`, long form. A row with no gene marks a profiled sample with nothing to show. */
   static cohort(seed, samples = 64) {
-    const draws = _chunkUMEKIIHJcjs.SampleData.distribution(seed, samples * 20, 1).map((row) => _chunkAL4XQZH6cjs.StatMath.normalCdf((row.value - 40) / 8));
+    const draws = _chunkLNWTC72Ycjs.SampleData.distribution(seed, samples * 20, 1).map((row) => _chunkAL4XQZH6cjs.StatMath.normalCdf((row.value - 40) / 8));
     let next = 0;
     const u = () => draws[next++ % draws.length];
     const genes = [TP53, MYC, ...PATHWAYS.flat()];

package/dist/{chunk-RXCOVONX.cjs → chunk-RBMUOIXE.cjs}

@@ -1,4 +1,4 @@
-"use strict";Object.defineProperty(exports, "__esModule", {value: true}); function _nullishCoalesce(lhs, rhsFn) { if (lhs != null) { return lhs; } else { return rhsFn(); } } function _optionalChain(ops) { let lastAccessLHS = undefined; let value = ops[0]; let i = 1; while (i < ops.length) { const op = ops[i]; const fn = ops[i + 1]; i += 2; if ((op === 'optionalAccess' || op === 'optionalCall') && value == null) { return undefined; } if (op === 'access' || op === 'optionalAccess') { lastAccessLHS = value; value = fn(value); } else if (op === 'call' || op === 'optionalCall') { value = fn((...args) => value.call(lastAccessLHS, ...args)); lastAccessLHS = undefined; } } return value; } var _class; var _class2; var _class3; var _class4; var _class5; var _class6; var _class7; var _class8; var _class9; var _class10;
+"use strict";Object.defineProperty(exports, "__esModule", {value: true}); function _nullishCoalesce(lhs, rhsFn) { if (lhs != null) { return lhs; } else { return rhsFn(); } } function _optionalChain(ops) { let lastAccessLHS = undefined; let value = ops[0]; let i = 1; while (i < ops.length) { const op = ops[i]; const fn = ops[i + 1]; i += 2; if ((op === 'optionalAccess' || op === 'optionalCall') && value == null) { return undefined; } if (op === 'access' || op === 'optionalAccess') { lastAccessLHS = value; value = fn(value); } else if (op === 'call' || op === 'optionalCall') { value = fn((...args) => value.call(lastAccessLHS, ...args)); lastAccessLHS = undefined; } } return value; } var _class; var _class2; var _class3; var _class4; var _class5; var _class6; var _class7; var _class8; var _class9; var _class10; var _class11;
 
 
 var _chunkTOX7G74Vcjs = require('./chunk-TOX7G74V.cjs');
@@ -23,7 +23,7 @@ var _chunkXX2YTCORcjs = require('./chunk-XX2YTCOR.cjs');
 var _chunkEEF5WRRZcjs = require('./chunk-EEF5WRRZ.cjs');
 
 
-var _chunkUMEKIIHJcjs = require('./chunk-UMEKIIHJ.cjs');
+var _chunkLNWTC72Ycjs = require('./chunk-LNWTC72Y.cjs');
 
 // src/charts/biologycharts/biologyEnums.ts
 var HelicalWheelHelix = /* @__PURE__ */ ((HelicalWheelHelix2) => {
@@ -105,6 +105,24 @@ var SurvivorshipAge = /* @__PURE__ */ ((SurvivorshipAge2) => {
   SurvivorshipAge2["Relative"] = "relative-age";
   return SurvivorshipAge2;
 })(SurvivorshipAge || {});
+var PunnettCross = /* @__PURE__ */ ((PunnettCross2) => {
+  PunnettCross2["Monohybrid"] = "monohybrid";
+  PunnettCross2["Dihybrid"] = "dihybrid";
+  PunnettCross2["Trihybrid"] = "trihybrid";
+  return PunnettCross2;
+})(PunnettCross || {});
+var PunnettDominance = /* @__PURE__ */ ((PunnettDominance2) => {
+  PunnettDominance2["Complete"] = "complete";
+  PunnettDominance2["Incomplete"] = "incomplete";
+  PunnettDominance2["Codominant"] = "codominant";
+  return PunnettDominance2;
+})(PunnettDominance || {});
+var PunnettRatios = /* @__PURE__ */ ((PunnettRatios2) => {
+  PunnettRatios2["List"] = "ratio-list";
+  PunnettRatios2["Bars"] = "ratio-bars";
+  PunnettRatios2["None"] = "no-ratios";
+  return PunnettRatios2;
+})(PunnettRatios || {});
 
 // src/charts/biologycharts/EpidemicCurve.ts
 var _d3 = require('d3');
@@ -164,22 +182,24 @@ var WEB = [
 ];
 var COHORTS = [["Mountain ungulate", 13, 4.2, 9.2], ["Songbird", 8, 1, 1.9], ["Perennial herb", 10, 0.35, 0.12]];
 var REASONS = ["Wrong population", "Wrong outcome", "Not a comparative study"];
+var PEA_GENES = [["Seed shape", "R", "Round", "Wrinkled"], ["Seed colour", "Y", "Yellow", "Green"], ["Flower colour", "P", "Purple", "White"]];
+var PEA_CROSSES = [[[1, 1], [1, 1], [1, 1]], [[1, 0], [1, 0], [1, 1]], [[1, 1], [2, 0], [1, 0]]];
 var BiologySamples = class {
   /** `{ position, residue }`: an 18-residue amphipathic helix, hydrophobic along one face. */
   static peptide(seed) {
-    const { u } = _chunkUMEKIIHJcjs.SampleData.rng(seed), face = u() * 360;
+    const { u } = _chunkLNWTC72Ycjs.SampleData.rng(seed), face = u() * 360;
     return residues(Array.from({ length: 18 }, (_, i) => pick(u, Math.cos((i * 100 + face) * Math.PI / 180) > 0.15 ? "LLLAAIVFWM" : "KKEERSQDNT")).join(""));
   }
   /** `{ position, residue }`: a protein of four membrane-spanning stretches between polar loops. */
   static membraneProtein(seed) {
-    const { u } = _chunkUMEKIIHJcjs.SampleData.rng(seed);
+    const { u } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
     let sequence = `M${letters(u, 16 + Math.floor(u() * 10), LOOP)}`;
     for (let i = 0; i < 4; i++) sequence += letters(u, 20 + Math.floor(u() * 4), SPAN) + letters(u, 14 + Math.floor(u() * 16), LOOP);
     return residues(sequence);
   }
   /** `{ gene, metric, inSet }`: 300 genes with a ranking metric, and a gene set enriched near the top of the list. */
   static rankedGenes(seed) {
-    const { u, n } = _chunkUMEKIIHJcjs.SampleData.rng(seed), count = 300;
+    const { u, n } = _chunkLNWTC72Ycjs.SampleData.rng(seed), count = 300;
     const metrics = Array.from({ length: count }, (_, i) => {
       const t = 1 - 2 * i / (count - 1);
       return Math.sign(t) * Math.abs(t) ** 1.7 * 2.8 + n() * 0.05;
@@ -188,12 +208,12 @@ var BiologySamples = class {
   }
   /** `{ site, species, count }`: three forest plots, from many even species to a few dominant ones. */
   static communities(seed) {
-    const { n } = _chunkUMEKIIHJcjs.SampleData.rng(seed);
+    const { n } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
     return COMMUNITIES.flatMap(([site, richness, spread]) => Array.from({ length: richness }, (_, i) => ({ site, species: `Species ${i + 1}`, count: Math.max(1, Math.round(Math.exp(2.3 + n() * spread))) })));
   }
   /** A lake food web: links run from prey to consumer, `value` is biomass on nodes and the share of the diet on links. */
   static foodWeb(seed) {
-    const { u } = _chunkUMEKIIHJcjs.SampleData.rng(seed);
+    const { u } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
     const nodes = WEB.map(([id, group, biomass]) => ({ id, group, value: +(biomass * (0.75 + u() * 0.5)).toFixed(1) }));
     const links = WEB.flatMap(([target, , , prey]) => {
       const shares = prey.map(() => 0.3 + u());
@@ -204,7 +224,7 @@ var BiologySamples = class {
   }
   /** `{ name, sequence }`: eight aligned protein fragments with a conserved core, gaps written as `-`. */
   static alignment(seed) {
-    const { u } = _chunkUMEKIIHJcjs.SampleData.rng(seed), length = 30;
+    const { u } = _chunkLNWTC72Ycjs.SampleData.rng(seed), length = 30;
     const ancestor = letters(u, length, AMINO), gapAt = 4 + Math.floor(u() * 4);
     return TAXA.map((name, taxon) => {
       const rate = 0.04 + taxon * 0.045;
@@ -218,7 +238,7 @@ var BiologySamples = class {
   }
   /** A line list, `{ onset, classification }` per case, plus `{ onset, intervention }` rows for what was done and when. */
   static outbreak(seed) {
-    const { u, n } = _chunkUMEKIIHJcjs.SampleData.rng(seed), start = Date.UTC(2024, 2, 4 + seed % 5);
+    const { u, n } = _chunkLNWTC72Ycjs.SampleData.rng(seed), start = Date.UTC(2024, 2, 4 + seed % 5);
     const date = (day) => new Date(start + day * DAY).toISOString().slice(0, 10);
     const cases = Array.from({ length: 84 }, (_, i) => {
       const day = Math.min(33, Math.floor(Math.exp(Math.log(8) + n() * 0.4)) + (i % 4 === 0 ? 9 + Math.floor(u() * 6) : 0));
@@ -229,7 +249,7 @@ var BiologySamples = class {
   }
   /** `{ stage, item, n }`: the counts of a systematic review, a blank `item` being a box's own total. */
   static prisma(seed) {
-    const { u } = _chunkUMEKIIHJcjs.SampleData.rng(seed), part = (total, low, high) => Math.round(total * (low + u() * (high - low)));
+    const { u } = _chunkLNWTC72Ycjs.SampleData.rng(seed), part = (total, low, high) => Math.round(total * (low + u() * (high - low)));
     const databases = 1400 + part(1200, 0, 1), registers = 90 + part(160, 0, 1), identified = databases + registers;
     const duplicates = part(identified, 0.18, 0.26), automated = part(identified, 0.02, 0.05), other = part(identified, 4e-3, 0.012);
     const screened = identified - duplicates - automated - other, excluded = part(screened, 0.74, 0.84), sought = screened - excluded;
@@ -254,7 +274,7 @@ var BiologySamples = class {
   }
   /** `{ population, age, survivors }`: three cohorts of 1,000 followed by age class, one of each survivorship type. */
   static lifeTables(seed) {
-    const { u } = _chunkUMEKIIHJcjs.SampleData.rng(seed);
+    const { u } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
     return COHORTS.flatMap(([population, oldest, shape, scale]) => {
       let alive = 1e3;
       return Array.from({ length: oldest + 1 }, (_, age) => {
@@ -263,6 +283,26 @@ var BiologySamples = class {
       });
     });
   }
+  /**
+   * `{ gene, parent1, parent2, dominance, dominantTrait, recessiveTrait, mixedTrait }`: three pea genes and the genotype of each parent at them.
+   * Odd seeds cross two plants heterozygous at every gene; other seeds change the parents, and make flower colour incompletely dominant.
+   */
+  static cross(seed) {
+    const parents = PEA_CROSSES[seed % 2 ? 0 : 1 + Math.floor(seed / 2) % 2], blended = seed % 2 === 0;
+    const genotype = (letter, dominant) => letter.repeat(dominant) + letter.toLowerCase().repeat(2 - dominant);
+    return PEA_GENES.map(([gene, letter, dominantTrait, recessiveTrait], i) => {
+      const incomplete = blended && i === 2;
+      return {
+        gene,
+        parent1: genotype(letter, parents[i][0]),
+        parent2: genotype(letter, parents[i][1]),
+        dominance: incomplete ? "incomplete" : "complete",
+        dominantTrait,
+        recessiveTrait,
+        mixedTrait: incomplete ? "Pink" : null
+      };
+    });
+  }
 };
 
 // src/charts/biologycharts/EpidemicCurve.ts
@@ -1037,9 +1077,219 @@ var PrismaFlowDiagram = (_class7 = class _PrismaFlowDiagram extends _chunkEEF5WR
   }
 }, _class7.__initStatic7(), _class7);
 
-// src/charts/biologycharts/RankAbundancePlot.ts
-var RankAbundancePlot = (_class8 = class _RankAbundancePlot extends _chunkUXFWDEP6cjs.StatisticalChart {
+// src/charts/biologycharts/PunnettSquare.ts
+var GENES = { ["monohybrid" /* Monohybrid */]: 1, ["dihybrid" /* Dihybrid */]: 2, ["trihybrid" /* Trihybrid */]: 3 };
+var CHAR2 = 6.2;
+var CAPTION = 15;
+var HEAD3 = 20;
+var GAP = 16;
+var LINE = 14;
+var BAR = 20;
+var gcd = (a, b) => b ? gcd(b, a % b) : a;
+var PunnettSquare = (_class8 = class _PunnettSquare extends _chunkEEF5WRRZcjs.Chart2D {
   static __initStatic8() {this.meta = {
+    id: "punnett-square",
+    name: "Punnett square",
+    family: "matrix" /* Matrix */,
+    dimension: 2 /* Two */,
+    shape: "table" /* Table */,
+    description: "Crosses two parents from their genotypes: gametes along the top and side, offspring genotypes in the cells coloured by phenotype, and the genotype and phenotype ratios, for one gene or several with dominance set per gene.",
+    tags: ["genetics", "mendel", "monohybrid", "dihybrid", "genotype", "phenotype", "inheritance", "biology"],
+    sample: (seed) => BiologySamples.cross(seed),
+    variants: {
+      cross: ["dihybrid" /* Dihybrid */, "monohybrid" /* Monohybrid */],
+      ratios: Object.values(PunnettRatios)
+    }
+  }}
+  /** True when the ratio list is drawn, which names every phenotype beside its colour. */
+  #keyed = false;
+  defaults() {
+    return {
+      gene: "gene",
+      parent1: "parent1",
+      parent2: "parent2",
+      dominance: "dominance",
+      dominantTrait: "dominantTrait",
+      recessiveTrait: "recessiveTrait",
+      mixedTrait: "mixedTrait",
+      parents: ["Parent 1", "Parent 2"],
+      cross: "dihybrid" /* Dihybrid */,
+      ratios: "ratio-list" /* List */
+    };
+  }
+  plotMargin() {
+    return { top: 8, right: 12, bottom: 8, left: 12 };
+  }
+  get legendPosition() {
+    return this.#keyed ? "none" /* None */ : super.legendPosition;
+  }
+  /** The two alleles of a genotype: `Rr` gives R and r, and `IA/i` or `IA i` gives alleles of more than one letter. Anything else gives none. */
+  static alleles(genotype) {
+    const parts = /[\s/|,]/.test(genotype) ? genotype.split(/[\s/|,]+/).filter(Boolean) : [...genotype];
+    return parts.length === 2 ? parts : [];
+  }
+  /** Counts as the lowest whole-number ratio: 12, 4 gives `3 : 1`. */
+  static ratio(counts) {
+    const shared = counts.reduce(gcd, 0) || 1;
+    return counts.map((count) => count / shared).join(" : ");
+  }
+  #genes() {
+    const { data, gene, parent1, parent2, dominance, dominantTrait, recessiveTrait, mixedTrait, cross } = this.options, seen = /* @__PURE__ */ new Set();
+    const upper = (allele) => allele === allele.toUpperCase() && allele !== allele.toLowerCase();
+    return (Array.isArray(data) ? data : []).flatMap((row, i) => {
+      if (row === null || typeof row !== "object") return [];
+      const name = BiologyRows.text(row[gene], `Gene ${i + 1}`);
+      const parents = [_PunnettSquare.alleles(BiologyRows.text(row[parent1])), _PunnettSquare.alleles(BiologyRows.text(row[parent2]))];
+      if (seen.has(name) || !parents[0].length || !parents[1].length) return [];
+      seen.add(name);
+      const written = [...new Set(parents.flat())], lead = _nullishCoalesce(written.find(upper), () => ( written[0])), rest = written.filter((allele) => allele !== lead);
+      const said = BiologyRows.text(row[dominance]).toLowerCase();
+      const kind = said.startsWith("inc") ? "incomplete" /* Incomplete */ : said.startsWith("cod") ? "codominant" /* Codominant */ : "complete" /* Complete */;
+      const weak = _nullishCoalesce(rest[0], () => ( lead)), complete = kind === "complete" /* Complete */;
+      return [{
+        dominance: kind,
+        alleles: [lead, ...rest],
+        parents,
+        traits: [BiologyRows.text(row[dominantTrait], complete ? `${lead}_` : `${lead}${lead}`), BiologyRows.text(row[recessiveTrait], `${weak}${weak}`), BiologyRows.text(row[mixedTrait], `${lead}${weak}`)]
+      }];
+    }).slice(0, _nullishCoalesce(GENES[cross], () => ( 2)));
+  }
+  /** What an offspring carrying `pair` shows at a gene. Only complete dominance hides the second allele of a heterozygote. */
+  static #trait({ alleles, dominance, traits }, pair) {
+    const carried = pair.filter((allele) => allele === alleles[0]).length;
+    if (carried === 2 || carried === 1 && dominance === "complete" /* Complete */) return traits[0];
+    return carried ? traits[2] : traits[1];
+  }
+  /** The cross worked out: every gamete of each parent, the offspring of each pairing, and the tallies behind the ratios. */
+  get square() {
+    const genes = this.#genes(), single = genes.every((gene) => gene.alleles.every((allele) => allele.length === 1));
+    const gametes = (parent) => genes.length ? genes.reduce((all, gene) => all.flatMap((gamete) => gene.parents[parent].map((allele) => [...gamete, allele])), [[]]) : [];
+    const pair = (gene, a, b) => [a, b].sort((p, q) => gene.alleles.indexOf(p) - gene.alleles.indexOf(q));
+    const written = (pairs) => pairs.map((alleles) => alleles.join(single ? "" : "/")).join(single ? "" : " ");
+    const down2 = gametes(0), across = gametes(1);
+    const cells = down2.flatMap((a, row) => across.map((b, column) => {
+      const pairs = genes.map((gene, g) => pair(gene, a[g], b[g]));
+      return { row, column, genotype: written(pairs), phenotype: genes.map((gene, g) => _PunnettSquare.#trait(gene, pairs[g])).join(", ") };
+    }));
+    const tally = (name) => {
+      const classes = /* @__PURE__ */ new Map();
+      for (const cell of cells) (_nullishCoalesce(classes.get(name(cell)), () => ( classes.set(name(cell), { name: name(cell), count: 0, phenotype: cell.phenotype }).get(name(cell))))).count++;
+      return [...classes.values()];
+    };
+    const phenotypes = tally((cell) => cell.phenotype), rank = new Map(phenotypes.map((entry, i) => [entry.name, i]));
+    const text = (gamete) => gamete.join(single ? "" : " ");
+    return {
+      side: down2.map(text),
+      top: across.map(text),
+      cells,
+      phenotypes,
+      parents: [written(genes.map((gene) => pair(gene, gene.parents[0][0], gene.parents[0][1]))), written(genes.map((gene) => pair(gene, gene.parents[1][0], gene.parents[1][1])))],
+      genotypes: tally((cell) => cell.genotype).sort((a, b) => rank.get(a.phenotype) - rank.get(b.phenotype))
+    };
+  }
+  draw(ctx) {
+    const { ratios, parents: names } = this.options, { width, height } = ctx.plot, { palette, theme } = ctx;
+    const { side, top, parents, cells, genotypes, phenotypes } = this.square, n = side.length, total = cells.length;
+    for (const entry of phenotypes) palette.color(entry.name);
+    const count = new Map(genotypes.map((entry) => [entry.name, entry.count]));
+    const letters2 = Math.max(1, ...side.map((gamete) => gamete.length), ...top.map((gamete) => gamete.length));
+    const head = Math.max(18, Math.min(13, 60 / letters2) * 0.62 * letters2 + 10), sideBand = CAPTION + head, topBand = CAPTION + HEAD3;
+    const panel = ratios !== "no-ratios" /* None */ && n ? Math.max(118, Math.min(230, width * 0.38)) : 0;
+    const beside = panel > 0 && width - sideBand - panel - GAP >= 80 ? panel + GAP : 0;
+    const size = Math.max(10, Math.min(height - topBand, width - sideBand - beside)), cell = size / (n || 1), pad = cell >= 12 ? 2 : 1;
+    const gx = sideBand + (width - sideBand - size - beside) / 2, gy = topBand + (height - topBand - size) / 2;
+    this.#keyed = beside > 0 && ratios === "ratio-list" /* List */;
+    const key = (d) => `${d.row}/${d.column}`, cx = (d) => gx + (d.column + 0.5) * cell, cy = (d) => gy + (d.row + 0.5) * cell;
+    const rects = ctx.marks("cells", "rect", cells, key).attr("fill", (d) => palette.color(d.phenotype)).attr("rx", Math.min(4, cell / 8));
+    ctx.set(rects, { x: (d) => cx(d) - (cell - pad) / 2, y: (d) => cy(d) - (cell - pad) / 2, width: Math.max(0.5, cell - pad), height: Math.max(0.5, cell - pad), opacity: 1 }, { opacity: 0 }, 240 / (total || 1));
+    ctx.hover(rects, (d) => ({
+      title: d.genotype,
+      color: palette.color(d.phenotype),
+      rows: [["Phenotype", d.phenotype], ["Gametes", `${side[d.row]} and ${top[d.column]}`], ["Share of offspring", `${count.get(d.genotype)} in ${total}`]]
+    }));
+    const longest = Math.max(1, ...cells.map((d) => d.genotype.length)), font = Math.min(14, (cell - 6) / (longest * 0.62));
+    const lettered = font >= 6.5 ? cells : [], captioned = cell >= 46 && phenotypes.every((entry) => _chunkEEF5WRRZcjs.TextFit.width(entry.name, 5.6) <= cell - 8);
+    const onCell = (name, data, size2, shift, words) => {
+      const marks = ctx.text(ctx.marks(name, "text", data, key), size2, "ink").attr("text-anchor", "middle").attr("dy", "0.35em").attr("fill", (d) => palette.inkOn(palette.color(d.phenotype))).style("pointer-events", "none").text(words);
+      ctx.set(marks, { x: cx, y: (d) => cy(d) + shift });
+      return marks;
+    };
+    onCell("cell-genotypes", lettered, font, captioned ? -7 : 0, (d) => d.genotype).attr("font-weight", 600);
+    onCell("cell-phenotypes", captioned ? lettered : [], 10, 9, (d) => d.phenotype);
+    const gameteFont = Math.max(7, Math.min(13, (Math.min(cell, head) - 3) / (letters2 * 0.62)));
+    const headers = [
+      ...top.map((text, i) => ({ key: `top-${i}`, text, x: gx + (i + 0.5) * cell, y: gy - 9, anchor: "middle" })),
+      ...side.map((text, i) => ({ key: `side-${i}`, text, x: gx - 7, y: gy + (i + 0.5) * cell, anchor: "end" }))
+    ];
+    const gametes = ctx.text(ctx.marks("axis-gametes", "text", cell >= gameteFont * 0.9 ? headers : [], (d) => d.key), gameteFont, "ink").attr("font-weight", 600).attr("dy", "0.35em").attr("text-anchor", (d) => d.anchor).text((d) => d.text);
+    ctx.set(gametes, { x: (d) => d.x, y: (d) => d.y });
+    const captions = n ? [
+      { key: "top", text: `${names[1]}: ${parents[1]}`, room: Math.min(width, size + 2 * Math.min(gx, width - gx - size)), transform: `translate(${(gx + size / 2).toFixed(1)},${(gy - HEAD3 - 5).toFixed(1)})` },
+      { key: "side", text: `${names[0]}: ${parents[0]}`, room: height, transform: `translate(${(gx - head - 5).toFixed(1)},${(gy + size / 2).toFixed(1)}) rotate(-90)` }
+    ] : [];
+    const titled = ctx.text(ctx.marks("axis-parents", "text", captions, (d) => d.key)).attr("text-anchor", "middle").text((d) => _chunkEEF5WRRZcjs.TextFit.fit(d.text, d.room, CHAR2));
+    ctx.set(titled, { transform: (d) => d.transform });
+    const notes = n ? [] : [{ key: "empty", text: "No genotypes to cross", x: 4, y: 12, muted: true }];
+    const chips = [];
+    const px = gx + size + GAP;
+    let y = gy - topBand + 10;
+    const heading = (title, tallies) => {
+      const summary = _PunnettSquare.ratio(tallies.map((entry) => entry.count)), lead = _chunkEEF5WRRZcjs.TextFit.width(title, 6.9) + 8, inline = lead + _chunkEEF5WRRZcjs.TextFit.width(summary, CHAR2) <= panel;
+      notes.push({ key: title, text: title, x: px, y, bold: true });
+      const lines = inline ? [summary] : summary.split(" : ").reduce((rows, part) => {
+        const last = rows.length - 1, joined = `${rows[last]} : ${part}`;
+        if (rows.length && (_chunkEEF5WRRZcjs.TextFit.width(joined, CHAR2) + 8 <= panel || rows.length === 3)) rows[last] = joined;
+        else rows.push(part);
+        return rows;
+      }, []);
+      lines.forEach((text, i) => {
+        if (!inline) y += LINE;
+        notes.push({ key: `${title}-ratio-${i}`, text: _chunkEEF5WRRZcjs.TextFit.fit(i < lines.length - 1 ? `${text} :` : text, inline ? panel - lead : panel, CHAR2), x: inline ? px + lead : px, y });
+      });
+    };
+    const sections = beside ? [["Phenotypes", phenotypes], ["Genotypes", genotypes]] : [];
+    sections.forEach(([title, tallies], s) => {
+      heading(title, tallies);
+      if (ratios === "ratio-bars" /* Bars */) {
+        const inner = (entry) => Math.max(0.5, panel * entry.count / (total || 1) - 2), fits = (entry, text) => _chunkEEF5WRRZcjs.TextFit.width(text, 5.8) + 4 <= inner(entry);
+        const named = s > 0 && tallies.every((entry) => fits(entry, entry.name));
+        let x = px;
+        for (const entry of tallies) {
+          const label = named ? entry.name : fits(entry, `${entry.count}`) ? `${entry.count}` : void 0;
+          chips.push({ key: `${title}/${entry.name}`, x, y: y + 7, width: inner(entry), height: BAR, fill: palette.color(entry.phenotype), tally: entry, label });
+          x += inner(entry) + 2;
+        }
+        y += 7 + BAR + 24;
+        return;
+      }
+      y += LINE + 2;
+      const labels = tallies.map((entry) => `${entry.count} ${entry.name}`), widest = Math.max(...labels.map((label) => _chunkEEF5WRRZcjs.TextFit.width(label, CHAR2)));
+      const columns = Math.max(1, Math.min(tallies.length, Math.floor((panel + 8) / (widest + 22)))), step = panel / columns;
+      const kept = s ? 0 : Math.min(2, Math.ceil(genotypes.length / 2)) + 3;
+      const slots = Math.max(1, Math.floor((height - y) / LINE + 1) - kept) * columns;
+      const shown = tallies.length <= slots ? tallies.length : slots - 1;
+      tallies.slice(0, shown).forEach((entry, i) => {
+        const at = px + i % columns * step, row = y + Math.floor(i / columns) * LINE;
+        chips.push({ key: `${title}/${entry.name}`, x: at, y: row - 9, width: 10, height: 10, fill: palette.color(entry.phenotype), tally: entry });
+        notes.push({ key: `${title}/${entry.name}`, text: _chunkEEF5WRRZcjs.TextFit.fit(labels[i], step - 16, CHAR2), x: at + 14, y: row });
+      });
+      if (shown < tallies.length) notes.push({ key: `${title}-more`, text: `+${tallies.length - shown} more`, x: px + shown % columns * step, y: y + Math.floor(shown / columns) * LINE, muted: true });
+      y += Math.ceil((shown + (shown < tallies.length ? 1 : 0)) / columns) * LINE + 8;
+    });
+    const swatches = ctx.marks("ratio-marks", "rect", chips, (d) => d.key).attr("fill", (d) => d.fill).attr("rx", 3);
+    ctx.set(swatches, { x: (d) => d.x, y: (d) => d.y, width: (d) => d.width, height: (d) => d.height });
+    ctx.hover(swatches, (d) => ({ title: d.tally.name, color: d.fill, rows: [["Offspring", `${d.tally.count} in ${total}`], ...d.tally.name === d.tally.phenotype ? [] : [["Phenotype", d.tally.phenotype]]] }));
+    const inside = ctx.text(ctx.marks("ratio-labels", "text", chips.filter((chip) => chip.label), (d) => d.key), 10, "ink").attr("text-anchor", "middle").attr("dy", "0.35em").attr("fill", (d) => palette.inkOn(d.fill)).style("pointer-events", "none").text((d) => d.label);
+    ctx.set(inside, { x: (d) => d.x + d.width / 2, y: (d) => d.y + d.height / 2 });
+    const written = ctx.marks("axis-ratios", "text", notes, (d) => d.key).attr("font-weight", (d) => d.bold ? 600 : 400).text((d) => d.text);
+    ctx.text(written).attr("fill", (d) => d.bold ? theme.ink : d.muted ? theme.inkMuted : theme.inkSecondary);
+    ctx.set(written, { x: (d) => d.x, y: (d) => d.y });
+  }
+}, _class8.__initStatic8(), _class8);
+
+// src/charts/biologycharts/RankAbundancePlot.ts
+var RankAbundancePlot = (_class9 = class _RankAbundancePlot extends _chunkUXFWDEP6cjs.StatisticalChart {
+  static __initStatic9() {this.meta = {
     id: "rank-abundance-plot",
     name: "Rank-abundance plot",
     family: "distribution" /* Distribution */,
@@ -1154,19 +1404,19 @@ var RankAbundancePlot = (_class8 = class _RankAbundancePlot extends _chunkUXFWDE
       { r: 3.5, key: (d) => `${d.community.site}/${d.item.species}` }
     );
   }
-}, _class8.__initStatic8(), _class8);
+}, _class9.__initStatic9(), _class9);
 
 // src/charts/biologycharts/SequenceAlignmentView.ts
-var GAP = "-";
-var CHAR2 = 6.2;
+var GAP2 = "-";
+var CHAR3 = 6.2;
 var TICKS = 13;
 var TRACK = 22;
 var SPACE = 4;
 var KEY = 18;
 var KEY_STEPS = 5;
 var KEY_TEXT = "low to high conservation";
-var SequenceAlignmentView = (_class9 = class _SequenceAlignmentView extends _chunkEEF5WRRZcjs.Chart2D {
-  static __initStatic9() {this.meta = {
+var SequenceAlignmentView = (_class10 = class _SequenceAlignmentView extends _chunkEEF5WRRZcjs.Chart2D {
+  static __initStatic10() {this.meta = {
     id: "sequence-alignment",
     name: "Sequence alignment view",
     family: "matrix" /* Matrix */,
@@ -1200,19 +1450,19 @@ var SequenceAlignmentView = (_class9 = class _SequenceAlignmentView extends _chu
     return Array.from({ length }, (_, column) => {
       const counts = /* @__PURE__ */ new Map();
       for (const letters2 of sequences) {
-        const letter = _nullishCoalesce(letters2[column], () => ( GAP));
-        if (letter !== GAP) counts.set(letter, (_nullishCoalesce(counts.get(letter), () => ( 0))) + 1);
+        const letter = _nullishCoalesce(letters2[column], () => ( GAP2));
+        if (letter !== GAP2) counts.set(letter, (_nullishCoalesce(counts.get(letter), () => ( 0))) + 1);
       }
       const ranked = [...counts].sort((a, b) => b[1] - a[1] || a[0].localeCompare(b[0])), filled = ranked.reduce((sum2, [, count]) => sum2 + count, 0);
       const entropy = -ranked.reduce((sum2, [, count]) => sum2 + count / filled * Math.log2(count / filled), 0);
-      return { consensus: _nullishCoalesce(_optionalChain([ranked, 'access', _17 => _17[0], 'optionalAccess', _18 => _18[0]]), () => ( GAP)), gaps: sequences.length - filled, conservation: filled ? Math.max(0, 1 - entropy / most) * (filled / sequences.length) : 0 };
+      return { consensus: _nullishCoalesce(_optionalChain([ranked, 'access', _17 => _17[0], 'optionalAccess', _18 => _18[0]]), () => ( GAP2)), gaps: sequences.length - filled, conservation: filled ? Math.max(0, 1 - entropy / most) * (filled / sequences.length) : 0 };
     });
   }
   #read() {
     const { data, name, sequence } = this.options;
     return (Array.isArray(data) ? data : []).flatMap((row, i) => {
       if (row === null || typeof row !== "object") return [];
-      const letters2 = BiologyRows.text(row[sequence]).toUpperCase().replace(/\./g, GAP).replace(/[^A-Z-]/g, "");
+      const letters2 = BiologyRows.text(row[sequence]).toUpperCase().replace(/\./g, GAP2).replace(/[^A-Z-]/g, "");
       return letters2 ? [{ name: BiologyRows.text(row[name], `Sequence ${i + 1}`), letters: letters2 }] : [];
     });
   }
@@ -1221,7 +1471,7 @@ var SequenceAlignmentView = (_class9 = class _SequenceAlignmentView extends _chu
     const rows = this.#read(), columns = _SequenceAlignmentView.columns(rows.map((row) => row.letters)), nucleic = _SequenceAlignmentView.nucleic(rows.map((row) => row.letters));
     const byConservation = colorBy === "by-conservation" /* Conservation */, differences = colorBy === "differences" /* Differences */;
     const labels = [...rows.map((row) => row.name), ...consensus ? ["Consensus"] : [], ...track ? ["Conservation"] : []];
-    const gutter = rows.length ? Math.min(96, width * 0.27, Math.max(...labels.map((label) => _chunkEEF5WRRZcjs.TextFit.width(label, CHAR2))) + 10) : 0;
+    const gutter = rows.length ? Math.min(96, width * 0.27, Math.max(...labels.map((label) => _chunkEEF5WRRZcjs.TextFit.width(label, CHAR3))) + 10) : 0;
     const cellW = Math.min(22, (width - gutter) / (columns.length || 1));
     const fixed = TICKS + (track ? TRACK + SPACE : 0) + (consensus ? SPACE : 0) + (byConservation ? KEY : 0);
     const cellH = Math.max(3, Math.min(20, (height - fixed) / (rows.length + (consensus ? 1 : 0) || 1)));
@@ -1230,14 +1480,14 @@ var SequenceAlignmentView = (_class9 = class _SequenceAlignmentView extends _chu
     const columnX = (column) => gutter + column * cellW;
     const kind = (letter) => nucleic ? letter : ResidueChemistry.classOf(letter);
     if (!byConservation) {
-      const present = new Set(rows.flatMap((row) => [...row.letters].filter((letter) => letter !== GAP).map(kind)));
+      const present = new Set(rows.flatMap((row) => [...row.letters].filter((letter) => letter !== GAP2).map(kind)));
       for (const name of nucleic ? [...present].sort() : ResidueChemistry.classes) if (present.has(name)) palette.color(name);
     }
     const fill = (letter, column) => byConservation ? palette.sequential(columns[column].conservation) : palette.color(kind(letter));
     const mark = (owner, letter, column, y) => ({ key: `${owner}/${column}`, owner, column, letter, x: columnX(column), y, fill: fill(letter, column) });
-    const body = rows.flatMap((row, r) => [...row.letters].flatMap((letter, column) => letter === GAP ? [] : [mark(`${r}:${row.name}`, letter, column, top + r * cellH)]));
+    const body = rows.flatMap((row, r) => [...row.letters].flatMap((letter, column) => letter === GAP2 ? [] : [mark(`${r}:${row.name}`, letter, column, top + r * cellH)]));
     const matches = differences ? body.filter((cell) => cell.letter === columns[cell.column].consensus) : [];
-    const summary = consensus ? columns.flatMap((column, i) => column.consensus === GAP ? [] : [mark("Consensus", column.consensus, i, consensusY)]) : [];
+    const summary = consensus ? columns.flatMap((column, i) => column.consensus === GAP2 ? [] : [mark("Consensus", column.consensus, i, consensusY)]) : [];
     const cells = [...differences ? body.filter((cell) => cell.letter !== columns[cell.column].consensus) : body, ...summary];
     const tell = (cell) => {
       const column = columns[cell.column];
@@ -1264,17 +1514,17 @@ var SequenceAlignmentView = (_class9 = class _SequenceAlignmentView extends _chu
       ...cellH >= 9 ? rows.map((row, r) => ({ key: `row-${r}`, text: row.name, y: top + (r + 0.5) * cellH, bold: false })) : [],
       ...consensus ? [{ key: "consensus", text: "Consensus", y: consensusY + cellH / 2, bold: true }] : []
     ];
-    const names = ctx.text(ctx.marks("axis-names", "text", named, (d) => d.key)).attr("text-anchor", "end").attr("dy", "0.32em").attr("font-weight", (d) => d.bold ? 600 : 400).text((d) => _chunkEEF5WRRZcjs.TextFit.fit(d.text, gutter - 6, CHAR2));
+    const names = ctx.text(ctx.marks("axis-names", "text", named, (d) => d.key)).attr("text-anchor", "end").attr("dy", "0.32em").attr("font-weight", (d) => d.bold ? 600 : 400).text((d) => _chunkEEF5WRRZcjs.TextFit.fit(d.text, gutter - 6, CHAR3));
     ctx.set(names, { x: gutter - 6, y: (d) => d.y });
     const keyY = consensusY + (consensus ? cellH : -SPACE) + 7;
     const steps = byConservation && rows.length ? Array.from({ length: KEY_STEPS }, (_, i) => i) : [];
     const swatches = ctx.marks("conservation-key", "rect", steps, (d) => d).attr("fill", (d) => palette.sequential(d / (KEY_STEPS - 1))).attr("rx", 2);
     ctx.set(swatches, { x: (d) => gutter + d * 14, y: keyY, width: 12, height: 8 });
-    const fits = gutter + KEY_STEPS * 14 + 4 + _chunkEEF5WRRZcjs.TextFit.width(KEY_TEXT, CHAR2) <= width;
+    const fits = gutter + KEY_STEPS * 14 + 4 + _chunkEEF5WRRZcjs.TextFit.width(KEY_TEXT, CHAR3) <= width;
     const caption = ctx.text(ctx.marks("conservation-key-text", "text", steps.length && fits ? [KEY_TEXT] : []), 10, "muted").text((d) => d);
     ctx.set(caption, { x: gutter + KEY_STEPS * 14 + 4, y: keyY + 8 });
   }
-}, _class9.__initStatic9(), _class9);
+}, _class10.__initStatic10(), _class10);
 
 // src/charts/biologycharts/SurvivorshipCurve.ts
 var RADIX = 1e3;
@@ -1284,8 +1534,8 @@ var LABELS = {
   ["linear-scale" /* Linear */]: [[0.62, 8, -4, "start"], [0.2, 8, -8, "start"], [0.02, 8, 3, "start"]]
 };
 var SAMPLES = 60;
-var SurvivorshipCurve = (_class10 = class _SurvivorshipCurve extends _chunkUXFWDEP6cjs.StatisticalChart {
-  static __initStatic10() {this.meta = {
+var SurvivorshipCurve = (_class11 = class _SurvivorshipCurve extends _chunkUXFWDEP6cjs.StatisticalChart {
+  static __initStatic11() {this.meta = {
     id: "survivorship-curve",
     name: "Survivorship curve",
     family: "scientific" /* Scientific */,
@@ -1378,7 +1628,7 @@ var SurvivorshipCurve = (_class10 = class _SurvivorshipCurve extends _chunkUXFWD
       return { key: name, text: name, x: fx(at(life)) + dx, y: fy(_SurvivorshipCurve.reference(shape, life)) + dy, anchor, tone: "muted" };
     }));
   }
-}, _class10.__initStatic10(), _class10);
+}, _class11.__initStatic11(), _class11);
 
 // src/charts/biologycharts/index.ts
 var biologychartsCatalog = [
@@ -1390,7 +1640,8 @@ var biologychartsCatalog = [
   SurvivorshipCurve,
   FoodWebChart,
   EpidemicCurve,
-  PrismaFlowDiagram
+  PrismaFlowDiagram,
+  PunnettSquare
 ];
 
 
@@ -1421,4 +1672,8 @@ var biologychartsCatalog = [
 
 
 
-exports.HelicalWheelHelix = HelicalWheelHelix; exports.HelicalWheelLayout = HelicalWheelLayout; exports.HydropathyScale = HydropathyScale; exports.HydropathyForm = HydropathyForm; exports.GseaWeighting = GseaWeighting; exports.RankAbundanceMeasure = RankAbundanceMeasure; exports.RankAbundanceModel = RankAbundanceModel; exports.RankAbundanceRank = RankAbundanceRank; exports.FoodWebLevel = FoodWebLevel; exports.SequenceAlignmentColor = SequenceAlignmentColor; exports.EpidemicCurveBin = EpidemicCurveBin; exports.EpidemicCurveForm = EpidemicCurveForm; exports.PrismaDetail = PrismaDetail; exports.SurvivorshipScale = SurvivorshipScale; exports.SurvivorshipAge = SurvivorshipAge; exports.EpidemicCurve = EpidemicCurve; exports.FoodWebChart = FoodWebChart; exports.GseaEnrichmentPlot = GseaEnrichmentPlot; exports.ResidueClass = ResidueClass; exports.HelicalWheelPlot = HelicalWheelPlot; exports.HydropathyPlot = HydropathyPlot; exports.PrismaStage = PrismaStage; exports.PrismaFlowDiagram = PrismaFlowDiagram; exports.RankAbundancePlot = RankAbundancePlot; exports.SequenceAlignmentView = SequenceAlignmentView; exports.SurvivorshipCurve = SurvivorshipCurve; exports.biologychartsCatalog = biologychartsCatalog;
+
+
+
+
+exports.HelicalWheelHelix = HelicalWheelHelix; exports.HelicalWheelLayout = HelicalWheelLayout; exports.HydropathyScale = HydropathyScale; exports.HydropathyForm = HydropathyForm; exports.GseaWeighting = GseaWeighting; exports.RankAbundanceMeasure = RankAbundanceMeasure; exports.RankAbundanceModel = RankAbundanceModel; exports.RankAbundanceRank = RankAbundanceRank; exports.FoodWebLevel = FoodWebLevel; exports.SequenceAlignmentColor = SequenceAlignmentColor; exports.EpidemicCurveBin = EpidemicCurveBin; exports.EpidemicCurveForm = EpidemicCurveForm; exports.PrismaDetail = PrismaDetail; exports.SurvivorshipScale = SurvivorshipScale; exports.SurvivorshipAge = SurvivorshipAge; exports.PunnettCross = PunnettCross; exports.PunnettDominance = PunnettDominance; exports.PunnettRatios = PunnettRatios; exports.EpidemicCurve = EpidemicCurve; exports.FoodWebChart = FoodWebChart; exports.GseaEnrichmentPlot = GseaEnrichmentPlot; exports.ResidueClass = ResidueClass; exports.HelicalWheelPlot = HelicalWheelPlot; exports.HydropathyPlot = HydropathyPlot; exports.PrismaStage = PrismaStage; exports.PrismaFlowDiagram = PrismaFlowDiagram; exports.PunnettSquare = PunnettSquare; exports.RankAbundancePlot = RankAbundancePlot; exports.SequenceAlignmentView = SequenceAlignmentView; exports.SurvivorshipCurve = SurvivorshipCurve; exports.biologychartsCatalog = biologychartsCatalog;

package/dist/{chunk-MSBNHCSE.cjs → chunk-RJHGCUUD.cjs}

@@ -18,7 +18,7 @@ var _chunkEEF5WRRZcjs = require('./chunk-EEF5WRRZ.cjs');
 var _chunkPY4R3IFFcjs = require('./chunk-PY4R3IFF.cjs');
 
 
-var _chunkUMEKIIHJcjs = require('./chunk-UMEKIIHJ.cjs');
+var _chunkLNWTC72Ycjs = require('./chunk-LNWTC72Y.cjs');
 
 // src/charts/trackscharts/tracksEnums.ts
 var GeneTrackLayout = /* @__PURE__ */ ((GeneTrackLayout2) => {
@@ -259,7 +259,7 @@ var PLASMID_LENGTH = 5370;
 var TracksSamples = class {
   /** Invented feature rows: `{ gene, transcript, feature, start, end, strand }`, for two neighbouring genes on opposite strands. */
   static geneModels(seed) {
-    const { u } = _chunkUMEKIIHJcjs.SampleData.rng(seed);
+    const { u } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
     const exons = (count, from) => {
       let at = from;
       return Array.from({ length: count }, () => {
@@ -287,7 +287,7 @@ var TracksSamples = class {
   }
   /** Invented rows for one protein: mutations with `{ position, change, type, count, cohort }`, domains with `{ domain, start, end }`, and `length` on all. */
   static mutations(seed) {
-    const { u } = _chunkUMEKIIHJcjs.SampleData.rng(seed);
+    const { u } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
     const pick = (text5) => text5[Math.floor(u() * text5.length)];
     const kind = () => {
       let roll = u();
@@ -326,7 +326,7 @@ var TracksSamples = class {
   }
   /** Two invented DNA sequences as `{ name, sequence }`. The second differs from the first by substitutions, an insertion, a deletion, an inversion and a tandem duplication. */
   static sequences(seed) {
-    const { u } = _chunkUMEKIIHJcjs.SampleData.rng(seed);
+    const { u } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
     const base = () => BASES[Math.floor(u() * 4)];
     const random = (count) => Array.from({ length: count }, base).join("");
     const worn = (text5) => [...text5].map((letter) => u() < 0.04 ? base() : letter).join("");
@@ -345,7 +345,7 @@ var TracksSamples = class {
   }
   /** Invented map of a cloning vector: `{ plasmid, name, type, start, end, strand, length }`. Rows with no `end` are restriction sites. */
   static plasmid(seed) {
-    const { u } = _chunkUMEKIIHJcjs.SampleData.rng(seed);
+    const { u } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
     const nudge = (value) => value + Math.round(u() * 24);
     const feature = (name, type, start, end, strand) => ({ plasmid: "pDEMO-12", name, type, start: nudge(start), end: nudge(end), strand, length: PLASMID_LENGTH });
     const site = (name, at) => ({ plasmid: "pDEMO-12", name, type: "Restriction site", start: nudge(at), end: null, strand: null, length: PLASMID_LENGTH });
@@ -371,7 +371,7 @@ var TracksSamples = class {
   }
   /** Invented RNA-seq depth: `{ position, depth, sample, strand }` every 20 bp over 12 kb, with most reads on exons of a plus-strand gene. */
   static coverage(seed) {
-    const { u, n } = _chunkUMEKIIHJcjs.SampleData.rng(seed);
+    const { u, n } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
     const gain = (position) => {
       const exon = COVERAGE_EXONS.find(([start, end]) => position >= start && position < end);
       if (!exon) return 0.02;