@jbrowse/bandage-core 10.2.7 → 10.3.1

package/dist/chunks/{chunk-RPRHR56C.js → chunk-BS6HPYCZ.js}

@@ -1,13 +1,13 @@
 import {
   polylineSlice
-} from "./chunk-N73TI3VG.js";
+} from "./chunk-BZFOFC4D.js";
 import {
   svgPath
-} from "./chunk-7SRWJLSC.js";
+} from "./chunk-UBQSRCWR.js";
 import {
   firstNodeAtOrAfter,
   isAnchored
-} from "./chunk-FYQWPEEY.js";
+} from "./chunk-RJWPJCIB.js";
 
 // src/genes/genePins.ts
 function genePins(graph, genes, positions) {

package/dist/chunks/chunk-BSOZW3FC.js

@@ -0,0 +1,30 @@
+// src/util/groups.ts
+function pushTo(map, key, value) {
+  const list = map.get(key);
+  if (list) {
+    list.push(value);
+  } else {
+    map.set(key, [value]);
+  }
+}
+function addTo(map, key, value) {
+  const set = map.get(key);
+  if (set) {
+    set.add(value);
+  } else {
+    map.set(key, /* @__PURE__ */ new Set([value]));
+  }
+}
+function groupBy(items, keyOf) {
+  const groups = /* @__PURE__ */ new Map();
+  for (const item of items) {
+    pushTo(groups, keyOf(item), item);
+  }
+  return groups;
+}
+
+export {
+  pushTo,
+  addTo,
+  groupBy
+};

package/dist/chunks/{chunk-XXT2JCAC.js → chunk-BX4AUMZE.js}

@@ -1,5 +1,9 @@
 // src/tubeMap/mismatches.ts
-import { forward, getXCoordinateOfBaseWithinNode } from "@jbrowse/tubemap-core";
+import {
+  forward,
+  getXCoordinateOfBaseWithinNode,
+  mirroredMismatch
+} from "@jbrowse/tubemap-core";
 function readMismatches(read, layout) {
   const marks = [];
   const entries = read.sequenceNew ?? [];
@@ -7,22 +11,16 @@ function readMismatches(read, layout) {
   entries.forEach((entry, i) => {
     const nodeIndex = layout.nodeMap.get(forward(entry.nodeName));
     const node = nodeIndex === void 0 ? void 0 : layout.nodes[nodeIndex];
-    const y = visits[i]?.y;
-    if (node && y !== void 0) {
-      const base = { readId: read.id, y, height: read.width };
+    const visit = visits[i];
+    if (node && visit?.y !== void 0) {
+      const base = { readId: read.id, y: visit.y, height: read.width };
       const at = (pos) => getXCoordinateOfBaseWithinNode(node, pos);
-      for (const mm of entry.mismatches) {
+      for (const edit of entry.mismatches) {
+        const mm = visit.isForward ? edit : mirroredMismatch(edit, node.sequenceLength);
         const x0 = at(mm.pos);
         if (x0 !== null) {
           if (mm.type === "insertion") {
-            const softClip = i === 0 && mm.pos === read.firstNodeOffset || i === entries.length - 1 && mm.pos === read.finalNodeCoverLength;
-            marks.push({
-              ...base,
-              kind: "insertion",
-              x: x0,
-              seq: mm.seq,
-              softClip
-            });
+            marks.push({ ...base, kind: "insertion", x: x0, seq: mm.seq });
           } else if (mm.type === "deletion" && mm.length !== void 0) {
             const x1 = at(mm.pos + mm.length);
             if (x1 !== null) {

package/dist/chunks/{chunk-N73TI3VG.js → chunk-BZFOFC4D.js}

@@ -1,10 +1,13 @@
 import {
   ownStrand,
   strandSides
-} from "./chunk-7SRWJLSC.js";
+} from "./chunk-UBQSRCWR.js";
+import {
+  pushTo
+} from "./chunk-BSOZW3FC.js";
 import {
   isBackbone
-} from "./chunk-FYQWPEEY.js";
+} from "./chunk-RJWPJCIB.js";
 
 // src/layout/mergeRuns.ts
 function flip(strand) {
@@ -14,13 +17,7 @@ function mergeRuns(graph, alone = /* @__PURE__ */ new Set()) {
   const sides = graph.edges.map(strandSides);
   const at = /* @__PURE__ */ new Map();
   const attach = (id, side, ei) => {
-    const key = `${id}|${side}`;
-    const list = at.get(key);
-    if (list) {
-      list.push(ei);
-    } else {
-      at.set(key, [ei]);
-    }
+    pushTo(at, `${id}|${side}`, ei);
   };
   graph.edges.forEach((edge, ei) => {
     attach(edge.from, sides[ei].from, ei);

package/dist/chunks/{chunk-4ETPXDO6.js → chunk-CGFPFBT2.js}

@@ -1,7 +1,7 @@
 import {
   formatBp,
   groupDigits
-} from "./chunk-PTABFYQ7.js";
+} from "./chunk-PFRQBCGC.js";
 
 // src/walkKey.ts
 function rangeText(contig, start, end) {

package/dist/chunks/{chunk-24WITXQY.js → chunk-CGMSKIEO.js}

@@ -3,27 +3,20 @@ import {
   edgeJoinOf,
   routedEdgeCurves,
   strandSides
-} from "./chunk-7SRWJLSC.js";
+} from "./chunk-UBQSRCWR.js";
+import {
+  groupBy
+} from "./chunk-BSOZW3FC.js";
 import {
   firstNodeAtOrAfter,
   isBackbone
-} from "./chunk-FYQWPEEY.js";
+} from "./chunk-RJWPJCIB.js";
 
 // src/deletionEdges.ts
 var MIN_DELETION_BP = 1;
 function bypassedPoints(nodePositions, bypassed) {
   return bypassed.flatMap((id) => nodePositions[id] ?? []);
 }
-function deletionArcCurves(nodePositions, deletion, axis, route, stranded = false) {
-  return edgeCurvesOf(
-    nodePositions,
-    deletion,
-    axis,
-    deletion.bypassed,
-    route,
-    stranded
-  );
-}
 function edgeCurvesOf(nodePositions, edge, axis, bypassed, route, stranded) {
   const from = nodePositions[edge.from];
   const to = nodePositions[edge.to];
@@ -44,17 +37,10 @@ function edgeCurvesOf(nodePositions, edge, axis, bypassed, route, stranded) {
   );
 }
 function backboneBySequence(graph) {
-  const bySequence = /* @__PURE__ */ new Map();
-  for (const node of graph.nodes) {
-    if (isBackbone(node)) {
-      const list = bySequence.get(node.stable.refName);
-      if (list) {
-        list.push(node);
-      } else {
-        bySequence.set(node.stable.refName, [node]);
-      }
-    }
-  }
+  const bySequence = groupBy(
+    graph.nodes.filter(isBackbone),
+    (node) => node.stable.refName
+  );
   for (const list of bySequence.values()) {
     list.sort((a, b) => a.stable.start - b.stable.start);
   }
@@ -167,8 +153,6 @@ function deletionDrawing(graph, deletions, show = true) {
 }
 
 export {
-  bypassedPoints,
-  deletionArcCurves,
   edgeCurvesOf,
   deletionEdges,
   closableDeletions,

package/dist/chunks/{chunk-VV6KGJ7M.js → chunk-CHNESV7L.js}

@@ -3,40 +3,37 @@ import {
 } from "./chunk-JMABOQWL.js";
 import {
   layoutModeByValue
-} from "./chunk-TKIJ5535.js";
+} from "./chunk-AZ4EEF6Q.js";
 import {
   parseRegion
-} from "./chunk-VQWHMWLF.js";
+} from "./chunk-M5NPD55W.js";
 import {
   version
-} from "./chunk-CG35S2S2.js";
+} from "./chunk-AQCOIE5R.js";
 import {
   axisScaleOf
-} from "./chunk-A7CFWFQS.js";
-import {
-  walkKey
-} from "./chunk-4ETPXDO6.js";
+} from "./chunk-BKEZLETL.js";
 import {
   figureScene
-} from "./chunk-23YFB5BB.js";
+} from "./chunk-QH2MIXBA.js";
 import {
   abgrAlpha,
   abgrBlue,
   abgrGreen,
   abgrRed
 } from "./chunk-G2FD7S3C.js";
+import {
+  FIGURE_DEFAULTS
+} from "./chunk-D7WIYPG3.js";
 import {
   engineSettingsOf,
   fitTransform,
   forceLayout,
   loadGraph
-} from "./chunk-QXUEF3AJ.js";
+} from "./chunk-KMBE2Q42.js";
 import {
   isBackbone
-} from "./chunk-FYQWPEEY.js";
-import {
-  facetCells
-} from "./chunk-XV4P5KEQ.js";
+} from "./chunk-RJWPJCIB.js";
 
 // src/figureData.ts
 function columns(keys) {
@@ -59,19 +56,18 @@ function figureData(graph, layout, o = {}) {
   if (layout.tubeMap) {
     throw new Error("figure data has no tables for a tube map yet");
   }
-  const width = o.width ?? 1200;
-  const height = o.height ?? 800;
+  const width = o.width ?? FIGURE_DEFAULTS.width;
+  const height = o.height ?? FIGURE_DEFAULTS.height;
   const pixelRows = layout.pixelRows ?? false;
   const {
     bounds,
     colorScheme,
     drawnDeletions,
     lift,
-    facet,
-    panels: facetted,
+    facets,
     labelOf,
     referenceName,
-    reference,
+    keyText,
     geometry
   } = figureScene(graph, layout, o);
   const fit = fitTransform(bounds, width, height, pixelRows);
@@ -80,20 +76,15 @@ function figureData(graph, layout, o = {}) {
   }
   const axis = axisScaleOf(fit.scale, pixelRows);
   const panels = [];
-  if (facetted) {
-    const place = facetCells(
-      facetted.map((p) => p.walks[0].name),
-      facet.field || "walk",
-      facet.domain
-    );
-    const ordered = facetted.map((panel, i) => ({ panel, cell: place.cells[i] })).sort((a, b) => a.cell - b.cell);
+  if (facets) {
+    const ordered = facets.panels.map((panel, i) => ({ panel, cell: facets.place.cells[i] })).sort((a, b) => a.cell - b.cell);
     for (const { panel } of ordered) {
       const walk = panel.walks[0];
-      const key = walkKey(walk, reference);
+      const key = keyText(walk);
       panels.push({
         highlight: panel,
         name: labelOf(walk.name),
-        key: key.delta + key.outside + key.reversed,
+        key: key.note,
         range: key.scale ?? ""
       });
     }

package/dist/chunks/{chunk-IGVV37VI.js → chunk-CKUVEHSS.js}

@@ -1,12 +1,12 @@
 import {
   buildNeighbors
-} from "./chunk-ES2N6V56.js";
+} from "./chunk-VX656OAZ.js";
 import {
   drawnNodeLength
-} from "./chunk-6JEWNXCQ.js";
+} from "./chunk-34BOZKIE.js";
 import {
   isBackbone
-} from "./chunk-FYQWPEEY.js";
+} from "./chunk-RJWPJCIB.js";
 
 // src/layout/referenceSeeds.ts
 var LANE_GAP = 30;

package/dist/chunks/{chunk-JBCZIBFC.js → chunk-CRCD3ODN.js}

@@ -51,7 +51,6 @@ function nodeInk(graph, nodeById, contigThickness, nodeWidth) {
 export {
   NODE_WIDTHS,
   meanDepth,
-  depthWidthFactor,
   nodeWidthPx,
   maxNodeWidthPx,
   nodeInk

package/dist/chunks/{chunk-AUFAPUO4.js → chunk-D7WIYPG3.js}

@@ -1,10 +1,12 @@
 import {
   ENGINE_DEFAULTS
-} from "./chunk-QXUEF3AJ.js";
+} from "./chunk-KMBE2Q42.js";
 
 // src/figureSettings.ts
 var FIGURE_DEFAULTS = {
   ...ENGINE_DEFAULTS,
+  width: 1200,
+  height: 800,
   colorScheme: "auto",
   nodeWidth: "depth",
   contigThickness: 6,

package/dist/chunks/{chunk-OIV2GBAJ.js → chunk-DBLV7PGG.js}

@@ -1,14 +1,14 @@
 import {
   PathVisits
-} from "./chunk-RAEABQNW.js";
+} from "./chunk-ZPZ7GNBT.js";
 import {
+  panSNFields,
   panSNSample
-} from "./chunk-BPLIXMH7.js";
+} from "./chunk-4XWZNPQ2.js";
 
 // src/gfa/graphTables.ts
-function walkFields(panSN) {
-  const [sample = "", haplotype = "", ...contig] = panSN.split("#");
-  return contig.length > 0 ? { sample, haplotype: +haplotype, contig: contig.join("#") } : void 0;
+function pLineName(name, start, end) {
+  return start === 0 ? name : `${name}:${start}-${end}`;
 }
 function namesByList(lists, names, listCount, pathNames) {
   const starts = new Int32Array(listCount + 1);
@@ -160,17 +160,27 @@ function graphFromTables({ nodes: n, links: l, walks: w }, name = "Imported GFA"
   const visitOrder = [];
   const blocks = [0];
   const pathSamples = [];
+  const edgeNameIndex = /* @__PURE__ */ new Map();
+  const edgeNames = [];
   for (const p of drawnWalks) {
-    const fields = walkFields(w.names[p]);
-    const pathName = fields ? `${fields.sample}#${fields.haplotype}#${fields.contig}` : w.names[p];
-    let named = nameIndex.get(pathName);
+    const walkName = w.names[p];
+    const fields = panSNFields(walkName);
+    const start = w.starts[p];
+    const visitName = fields ? `${fields.sample}#${fields.haplotype}#${fields.contig}` : walkName;
+    const pathName = fields ? visitName : pLineName(walkName, start, w.ends[p]);
+    let named = nameIndex.get(visitName);
     if (named === void 0) {
       named = pathNames.length;
-      nameIndex.set(pathName, named);
-      pathNames.push(pathName);
-      pathSamples.push(panSNSample(pathName));
+      nameIndex.set(visitName, named);
+      pathNames.push(visitName);
+      pathSamples.push(panSNSample(visitName));
+    }
+    let edgeNamed = edgeNameIndex.get(pathName);
+    if (edgeNamed === void 0) {
+      edgeNamed = edgeNames.length;
+      edgeNameIndex.set(pathName, edgeNamed);
+      edgeNames.push(pathName);
     }
-    const start = w.starts[p];
     const visitSample = pathSamples[named];
     const first = w.offsets[p];
     const nodeIds = new Array(w.offsets[p + 1] - first);
@@ -194,14 +204,14 @@ function graphFromTables({ nodes: n, links: l, walks: w }, name = "Imported GFA"
         if (list >= 0 && lastWalk[list] !== p) {
           lastWalk[list] = p;
           crossedList[crossings] = list;
-          crossedName[crossings++] = named;
+          crossedName[crossings++] = edgeNamed;
         }
       }
       prev = node;
     }
     paths.push({ name: pathName, nodeIds, start, ...fields });
     anchorPaths.push({
-      name: pathName,
+      name: visitName,
       sample: visitSample,
       start,
       length: pos - start
@@ -211,7 +221,7 @@ function graphFromTables({ nodes: n, links: l, walks: w }, name = "Imported GFA"
     crossedList.subarray(0, crossings),
     crossedName,
     lists,
-    pathNames
+    edgeNames
   );
   edges.forEach((edge, i) => {
     const names = listNames[edgeLists[i]];
@@ -254,7 +264,7 @@ function graphTablesGFA({ nodes, links, walks }) {
   }
   for (let p = 0; p < walks.names.length; p++) {
     const name = walks.names[p];
-    const fields = walkFields(name);
+    const fields = panSNFields(name);
     const start = walks.starts[p];
     const end = walks.ends[p];
     const steps = [];
@@ -262,7 +272,7 @@ function graphTablesGFA({ nodes, links, walks }) {
       steps.push([names[walks.steps[s]], walks.reversed[s]]);
     }
     lines.push(
-      fields ? `W	${fields.sample}	${fields.haplotype}	${fields.contig}	${start}	${end}	${steps.map(([id, rev]) => (rev ? "<" : ">") + id).join("")}` : `P	${start === 0 ? name : `${name}:${start}-${end}`}	${steps.map(([id, rev]) => id + (rev ? "-" : "+")).join(",")}	*`
+      fields ? `W	${fields.sample}	${fields.haplotype}	${fields.contig}	${start}	${end}	${steps.map(([id, rev]) => (rev ? "<" : ">") + id).join("")}` : `P	${pLineName(name, start, end)}	${steps.map(([id, rev]) => id + (rev ? "-" : "+")).join(",")}	*`
     );
   }
   return lines.join("\n");

package/dist/chunks/{chunk-VQC7KELB.js → chunk-DO62DG26.js}

@@ -1,11 +1,15 @@
 import {
   REFERENCE_RAMP_ALT_COLOR,
-  REFERENCE_RAMP_MAX_HUE,
   hslToRgb
-} from "./chunk-A4PFFC4Q.js";
+} from "./chunk-SWCULW6S.js";
 import {
   packAbgr
 } from "./chunk-G2FD7S3C.js";
+import {
+  REFERENCE_RAMP_LIGHTNESS,
+  REFERENCE_RAMP_SATURATION,
+  rampHueAt
+} from "./chunk-GQZVR5DW.js";
 
 // src/walkEncoding.ts
 var WALK_FIELDS = [
@@ -61,7 +65,11 @@ function schemeRgb(scheme, t) {
   const x = Math.max(0, Math.min(1, t));
   const family = WALK_SCHEMES.find((s) => s.value === scheme);
   if (!family || !("ends" in family)) {
-    return hslToRgb(x * REFERENCE_RAMP_MAX_HUE, 0.7, 0.5).map((c) => c * 255);
+    return hslToRgb(
+      rampHueAt(x),
+      REFERENCE_RAMP_SATURATION,
+      REFERENCE_RAMP_LIGHTNESS
+    ).map((c) => c * 255);
   }
   const from = oklab(family.ends[0]);
   const to = oklab(family.ends[1]);

package/dist/chunks/{chunk-JA2HJIYS.js → chunk-E4ZM4OOJ.js}

@@ -1,13 +1,13 @@
 import {
   bubbleSegmentIds,
   classifyBubble
-} from "./chunk-URZ3CH2R.js";
+} from "./chunk-A4JDEUAI.js";
 import {
   svgPath
-} from "./chunk-7SRWJLSC.js";
+} from "./chunk-UBQSRCWR.js";
 import {
   isBackbone
-} from "./chunk-FYQWPEEY.js";
+} from "./chunk-RJWPJCIB.js";
 
 // src/bubbles/bubbleHalos.ts
 var WHOLE_FRACTION = 0.9;
@@ -64,8 +64,6 @@ function bubbleHalos(graph, bubbles, positions, repeats = []) {
 }
 
 export {
-  WHOLE_FRACTION,
   coversGraph,
-  SMALL_VARIANT_BP,
   bubbleHalos
 };

package/dist/chunks/{chunk-QLWI2PZ3.js → chunk-EJRJAAN7.js}

@@ -46,18 +46,22 @@ var Canvas2DRenderer = class {
     this.highlightFactor = factor;
   }
   updateTransform(transform) {
+    this.transform = transform;
+  }
+  fitted() {
+    const t = this.transform;
     const k = this.fit;
-    this.transform = k === 1 ? transform : {
-      ...transform,
-      scaleX: transform.scaleX * k,
-      scaleY: transform.scaleY * k,
-      translateX: transform.translateX * k,
-      translateY: transform.translateY * k,
-      dpr: transform.dpr * k
+    return !t || k === 1 ? t : {
+      ...t,
+      scaleX: t.scaleX * k,
+      scaleY: t.scaleY * k,
+      translateX: t.translateX * k,
+      translateY: t.translateY * k,
+      dpr: t.dpr * k
     };
   }
   render(clearColor) {
-    const t = this.transform;
+    const t = this.fitted();
     if (!t) {
       return;
     }
@@ -84,7 +88,7 @@ var Canvas2DRenderer = class {
       weight: e.thickness,
       last: inRun(edgeRun, i)
     }))) {
-      this.strokeCurves(items, color);
+      this.strokeCurves(items, color, t);
     }
     const lifted = this.liftedStrokes(batch);
     for (const { color, items } of groupByPaint(batch.nodeStrokes, (s, i) => {
@@ -95,7 +99,7 @@ var Canvas2DRenderer = class {
         last: f !== void 0
       };
     })) {
-      this.strokeNodes(items, color);
+      this.strokeNodes(items, color, t);
     }
     const arrowRun = this.edgeRun(batch.arrowRuns);
     for (const { color, items } of groupByPaint(batch.arrows, (a, i) => ({
@@ -103,7 +107,7 @@ var Canvas2DRenderer = class {
       weight: 0,
       last: inRun(arrowRun, i)
     }))) {
-      this.fillArrows(items, color);
+      this.fillArrows(items, color, t);
     }
   }
   // Only what the highlights lift, over a transparent clear, for a hover layer
@@ -113,7 +117,8 @@ var Canvas2DRenderer = class {
     const ctx = this.ctx;
     ctx.clearRect(0, 0, ctx.canvas.width, ctx.canvas.height);
     const batch = this.batch;
-    if (!this.transform || !batch) {
+    const t = this.fitted();
+    if (!t || !batch) {
       return;
     }
     ctx.lineCap = "round";
@@ -125,7 +130,7 @@ var Canvas2DRenderer = class {
       weight: e.thickness,
       last: false
     }))) {
-      this.strokeCurves(items, color);
+      this.strokeCurves(items, color, t);
     }
     const lifted = [...this.liftedStrokes(batch)];
     for (const { color, items } of groupByPaint(lifted, ([i, f]) => ({
@@ -135,7 +140,8 @@ var Canvas2DRenderer = class {
     }))) {
       this.strokeNodes(
         items.map(([i]) => batch.nodeStrokes[i]),
-        color
+        color,
+        t
       );
     }
     const arrows = runItems(batch.arrows, this.edgeRun(batch.arrowRuns));
@@ -144,7 +150,7 @@ var Canvas2DRenderer = class {
       weight: 0,
       last: false
     }))) {
-      this.fillArrows(items, color);
+      this.fillArrows(items, color, t);
     }
   }
   edgeRun(runs) {
@@ -163,8 +169,7 @@ var Canvas2DRenderer = class {
     }
     return lifted;
   }
-  strokeCurves(items, color) {
-    const t = this.transform;
+  strokeCurves(items, color, t) {
     const ctx = this.ctx;
     ctx.strokeStyle = abgrToCssRgba(color);
     ctx.lineWidth = items[0].thickness * 2 * t.dpr;
@@ -188,8 +193,7 @@ var Canvas2DRenderer = class {
     }
     ctx.stroke();
   }
-  strokeNodes(items, color) {
-    const t = this.transform;
+  strokeNodes(items, color, t) {
     const ctx = this.ctx;
     ctx.strokeStyle = abgrToCssRgba(color);
     ctx.lineWidth = items[0].thickness * 2 * t.dpr;
@@ -204,8 +208,7 @@ var Canvas2DRenderer = class {
     }
     ctx.stroke();
   }
-  fillArrows(items, color) {
-    const t = this.transform;
+  fillArrows(items, color, t) {
     const ctx = this.ctx;
     ctx.fillStyle = abgrToCssRgba(color);
     ctx.beginPath();