@jbrowse/bandage-core 10.2.7 → 10.3.0

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (190) hide show
  1. package/dist/alleleProjection/projectAlleles.js +8 -4
  2. package/dist/anchoredNodes.js +2 -1
  3. package/dist/bubbleSpreads.js +2 -2
  4. package/dist/bubbles/bubbleHalos.js +10 -11
  5. package/dist/bubbles/bubblesFromGraph.js +8 -3
  6. package/dist/bubbles/classifyBubble.js +9 -6
  7. package/dist/bubbles/popBubble.js +7 -5
  8. package/dist/bubbles/superbubbles.js +18 -14
  9. package/dist/chunks/{chunk-RAKOB6IV.js → chunk-234NFZS4.js} +8 -6
  10. package/dist/chunks/{chunk-A4PFFC4Q.js → chunk-2XEQR54X.js} +60 -102
  11. package/dist/chunks/{chunk-6JEWNXCQ.js → chunk-34BOZKIE.js} +0 -2
  12. package/dist/chunks/{chunk-RVS2QCOX.js → chunk-3GLG3ZLE.js} +2 -11
  13. package/dist/chunks/{chunk-PUG7LXDF.js → chunk-3PZSBNLI.js} +4 -1
  14. package/dist/chunks/{chunk-KBEP7VYK.js → chunk-3VCKBSML.js} +4 -4
  15. package/dist/chunks/{chunk-BPLIXMH7.js → chunk-4XWZNPQ2.js} +14 -0
  16. package/dist/chunks/{chunk-LKW2TP4N.js → chunk-5YILQTHE.js} +18 -16
  17. package/dist/chunks/{chunk-B6IW3IRE.js → chunk-65ILXEIN.js} +14 -15
  18. package/dist/chunks/{chunk-4DRQTBBE.js → chunk-6EOB6ANX.js} +1 -1
  19. package/dist/chunks/{chunk-QBHPHLAQ.js → chunk-7JGLJK44.js} +41 -36
  20. package/dist/chunks/{chunk-T2PCSAFY.js → chunk-7XR26JST.js} +2 -2
  21. package/dist/chunks/{chunk-URZ3CH2R.js → chunk-A4JDEUAI.js} +1 -1
  22. package/dist/chunks/{chunk-JLSJMCDJ.js → chunk-A54O52XZ.js} +3 -3
  23. package/dist/chunks/chunk-A55OZVYX.js +327 -0
  24. package/dist/chunks/{chunk-CWFNEKST.js → chunk-A7H753BA.js} +3 -3
  25. package/dist/chunks/{chunk-VV6KGJ7M.js → chunk-AMX4OQAW.js} +18 -27
  26. package/dist/chunks/{chunk-RPRHR56C.js → chunk-BS6HPYCZ.js} +3 -3
  27. package/dist/chunks/chunk-BSOZW3FC.js +30 -0
  28. package/dist/chunks/{chunk-CU6P5RUG.js → chunk-BVYR7IFN.js} +22 -16
  29. package/dist/chunks/{chunk-XXT2JCAC.js → chunk-BX4AUMZE.js} +11 -13
  30. package/dist/chunks/{chunk-N73TI3VG.js → chunk-BZFOFC4D.js} +6 -9
  31. package/dist/chunks/{chunk-4ETPXDO6.js → chunk-CGFPFBT2.js} +1 -1
  32. package/dist/chunks/{chunk-24WITXQY.js → chunk-CGMSKIEO.js} +9 -25
  33. package/dist/chunks/{chunk-IGVV37VI.js → chunk-CKUVEHSS.js} +3 -3
  34. package/dist/chunks/{chunk-JBCZIBFC.js → chunk-CRCD3ODN.js} +0 -1
  35. package/dist/chunks/{chunk-OIV2GBAJ.js → chunk-DBLV7PGG.js} +27 -17
  36. package/dist/chunks/{chunk-JA2HJIYS.js → chunk-E4ZM4OOJ.js} +3 -5
  37. package/dist/chunks/{chunk-ABUIN5AK.js → chunk-EHXBZ5JZ.js} +29 -24
  38. package/dist/chunks/{chunk-QLWI2PZ3.js → chunk-EJRJAAN7.js} +24 -21
  39. package/dist/chunks/{chunk-Q7LYHGGQ.js → chunk-GIHTHYXY.js} +5 -4
  40. package/dist/chunks/{chunk-L3O7AWLI.js → chunk-GJLZ2MKE.js} +2 -2
  41. package/dist/chunks/{chunk-LOEA5IHZ.js → chunk-GJU6BJJA.js} +8 -8
  42. package/dist/chunks/{chunk-XL472HHD.js → chunk-GQZVR5DW.js} +12 -7
  43. package/dist/chunks/{chunk-UVBGRJQX.js → chunk-HKYXWZQS.js} +1 -1
  44. package/dist/chunks/chunk-HMXCN6PG.js +17 -0
  45. package/dist/chunks/{chunk-WSIFAZU7.js → chunk-HNWJ4B7S.js} +12 -13
  46. package/dist/chunks/chunk-HOAIOEI7.js +18 -0
  47. package/dist/chunks/{chunk-XMW266WS.js → chunk-HSVFCA2J.js} +18 -20
  48. package/dist/chunks/{chunk-3XSM4OZU.js → chunk-IYALCKNB.js} +5 -5
  49. package/dist/chunks/{chunk-QXUEF3AJ.js → chunk-KWHWNOTY.js} +22 -17
  50. package/dist/chunks/{chunk-M6LK4GII.js → chunk-LGJBB5BP.js} +8 -9
  51. package/dist/chunks/{chunk-VQWHMWLF.js → chunk-M5NPD55W.js} +18 -14
  52. package/dist/chunks/{chunk-TKIJ5535.js → chunk-MD5VPXXT.js} +6 -6
  53. package/dist/chunks/{chunk-HK5Y2YWV.js → chunk-MYDHU2UV.js} +6 -11
  54. package/dist/chunks/{chunk-TZFWUT2C.js → chunk-MYLAJDKB.js} +2 -2
  55. package/dist/chunks/{chunk-4HWM3FNK.js → chunk-OMUPUB4J.js} +39 -56
  56. package/dist/chunks/{chunk-PTABFYQ7.js → chunk-PFRQBCGC.js} +49 -60
  57. package/dist/chunks/{chunk-T6V4KI2L.js → chunk-PPDE4IVT.js} +1 -1
  58. package/dist/chunks/{chunk-KC7T5RWG.js → chunk-PSCSZSVK.js} +9 -16
  59. package/dist/chunks/{chunk-CG35S2S2.js → chunk-PTQPTDZH.js} +1 -1
  60. package/dist/chunks/{chunk-L2HDM2R3.js → chunk-QQCOVMMT.js} +7 -7
  61. package/dist/chunks/{chunk-CQHUZHET.js → chunk-RFG42VIO.js} +46 -30
  62. package/dist/chunks/{chunk-FYQWPEEY.js → chunk-RJWPJCIB.js} +5 -11
  63. package/dist/chunks/{chunk-AR4BUBMU.js → chunk-RWKP7UO4.js} +34 -35
  64. package/dist/chunks/{chunk-WBSWROE6.js → chunk-SHKGYGM4.js} +14 -18
  65. package/dist/chunks/{chunk-2CMYNW3E.js → chunk-SI3BB4DO.js} +1 -1
  66. package/dist/chunks/{chunk-DY3JCAMX.js → chunk-SOMZGDVS.js} +7 -16
  67. package/dist/chunks/{chunk-EFQLEGWM.js → chunk-SZHQK6GL.js} +7 -8
  68. package/dist/chunks/{chunk-5GMFEPWR.js → chunk-TCOLHD7N.js} +16 -28
  69. package/dist/chunks/{chunk-XV4P5KEQ.js → chunk-TGUFT6KK.js} +6 -3
  70. package/dist/chunks/{chunk-CM6LO6OC.js → chunk-UAFATXYA.js} +10 -14
  71. package/dist/chunks/{chunk-7SRWJLSC.js → chunk-UBQSRCWR.js} +11 -8
  72. package/dist/chunks/{chunk-RUCH6ZFL.js → chunk-V4YML3GR.js} +2 -1
  73. package/dist/chunks/{chunk-VQC7KELB.js → chunk-V7J6U4XC.js} +11 -3
  74. package/dist/chunks/{chunk-ES2N6V56.js → chunk-VX656OAZ.js} +6 -11
  75. package/dist/chunks/chunk-WA5WNP6V.js +37 -0
  76. package/dist/chunks/{chunk-EQ2NIFRH.js → chunk-WJ4TIAVA.js} +1 -1
  77. package/dist/chunks/{chunk-A7CFWFQS.js → chunk-WONB6LMC.js} +1 -1
  78. package/dist/chunks/{chunk-23YFB5BB.js → chunk-XUT4HXVE.js} +29 -13
  79. package/dist/chunks/{chunk-G45YY72K.js → chunk-Z6XZWPXK.js} +17 -25
  80. package/dist/chunks/{chunk-XI3ZUTWF.js → chunk-ZNJ4JAPT.js} +14 -20
  81. package/dist/chunks/{chunk-AUFAPUO4.js → chunk-ZPBEUOBM.js} +3 -1
  82. package/dist/chunks/{chunk-RAEABQNW.js → chunk-ZPZ7GNBT.js} +1 -0
  83. package/dist/cli/figure.js +59 -56
  84. package/dist/cli/fold.js +19 -17
  85. package/dist/deletionEdges.js +5 -7
  86. package/dist/facetGrid.js +2 -1
  87. package/dist/figure.js +45 -42
  88. package/dist/figureData.js +48 -44
  89. package/dist/figureScene.js +34 -29
  90. package/dist/figureSettings.js +19 -17
  91. package/dist/foldVariants.js +5 -3
  92. package/dist/gbzCut.js +12 -7
  93. package/dist/gbzJoin.js +2 -1
  94. package/dist/gbzWindow.js +11 -6
  95. package/dist/genes/genePins.js +6 -4
  96. package/dist/gfa/gfaConverter.js +7 -5
  97. package/dist/gfa/gfaTables.js +2 -1
  98. package/dist/gfa/graphTables.js +3 -3
  99. package/dist/graphLabels.js +8 -5
  100. package/dist/index.js +88 -85
  101. package/dist/labelLayout.js +12 -9
  102. package/dist/layout/anchoredLayout.js +8 -5
  103. package/dist/layout/drawnScale.js +1 -5
  104. package/dist/layout/forceCache.js +19 -17
  105. package/dist/layout/mergeRuns.js +5 -3
  106. package/dist/layout/orderedLayout.js +5 -2
  107. package/dist/layout/orientToReference.js +3 -2
  108. package/dist/layout/placeOffReference.js +7 -4
  109. package/dist/layout/referenceSeeds.js +6 -4
  110. package/dist/layout/sampleRowLayout.js +8 -5
  111. package/dist/layout/trimToWindow.js +6 -4
  112. package/dist/layout/tubeMapLayout.js +17 -19
  113. package/dist/layout/walkRowDraw.js +18 -17
  114. package/dist/layout/walkRowLayout.js +10 -6
  115. package/dist/layout/walkRows.js +9 -5
  116. package/dist/layout/walkStrip.js +21 -18
  117. package/dist/layoutModes.js +23 -20
  118. package/dist/nodeWidths.js +1 -3
  119. package/dist/pansn.js +5 -1
  120. package/dist/pathAnchoring.js +6 -4
  121. package/dist/pathColors.js +12 -9
  122. package/dist/pathVisits.js +3 -1
  123. package/dist/pipeline.js +18 -16
  124. package/dist/reference.js +5 -3
  125. package/dist/referenceRampCss.js +9 -17
  126. package/dist/referenceSpan.js +4 -2
  127. package/dist/referenceStrip.js +18 -14
  128. package/dist/renderer/Canvas2DRenderer.js +1 -1
  129. package/dist/renderer/GeometryBuilder.js +19 -14
  130. package/dist/renderer/svgCanvas.js +1 -1
  131. package/dist/tubeMap/axis.js +11 -6
  132. package/dist/tubeMap/bundle.js +16 -12
  133. package/dist/tubeMap/coarsen.js +20 -345
  134. package/dist/tubeMap/deviations.js +53 -23
  135. package/dist/tubeMap/draw.js +11 -8
  136. package/dist/tubeMap/frame.js +3 -2
  137. package/dist/tubeMap/genes.js +10 -7
  138. package/dist/tubeMap/mismatches.js +1 -1
  139. package/dist/tubeMap/nodeColors.js +18 -14
  140. package/dist/tubeMap/reads.js +1 -1
  141. package/dist/tubeMap/warp.js +2 -1
  142. package/dist/types/alleleProjection/projectAlleles.d.ts +1 -4
  143. package/dist/types/bubbles/bubbleHalos.d.ts +0 -2
  144. package/dist/types/bubbles/bubblesFromGraph.d.ts +6 -4
  145. package/dist/types/bubbles/classifyBubble.d.ts +2 -1
  146. package/dist/types/deletionEdges.d.ts +0 -8
  147. package/dist/types/figureScene.d.ts +17 -7
  148. package/dist/types/figureSettings.d.ts +2 -0
  149. package/dist/types/gbzCut.d.ts +4 -0
  150. package/dist/types/gbzWindow.d.ts +1 -0
  151. package/dist/types/gfa-core/gfaParser.d.ts +1 -1
  152. package/dist/types/graphLabels.d.ts +2 -1
  153. package/dist/types/layout/drawnScale.d.ts +0 -2
  154. package/dist/types/layout/referenceSeeds.d.ts +2 -1
  155. package/dist/types/layout/tubeMapLayout.d.ts +0 -7
  156. package/dist/types/layout/walkRowDraw.d.ts +0 -1
  157. package/dist/types/layout/walkRowLayout.d.ts +1 -0
  158. package/dist/types/layout/walkRows.d.ts +1 -0
  159. package/dist/types/nodeWidths.d.ts +0 -1
  160. package/dist/types/pansn.d.ts +9 -0
  161. package/dist/types/pathVisits.d.ts +6 -0
  162. package/dist/types/referenceRampCss.d.ts +4 -0
  163. package/dist/types/renderer/Canvas2DRenderer.d.ts +1 -0
  164. package/dist/types/renderer/GeometryBuilder.d.ts +3 -2
  165. package/dist/types/tubeMap/axis.d.ts +3 -1
  166. package/dist/types/tubeMap/coarsen.d.ts +3 -1
  167. package/dist/types/tubeMap/deviations.d.ts +2 -1
  168. package/dist/types/tubeMap/mismatches.d.ts +0 -1
  169. package/dist/types/tubeMap/warp.d.ts +2 -1
  170. package/dist/types/util/SpatialIndex.d.ts +1 -1
  171. package/dist/types/util/geometry.d.ts +4 -0
  172. package/dist/types/util/groups.d.ts +3 -0
  173. package/dist/types/util/hitDetection.d.ts +2 -2
  174. package/dist/types/util/partitionPoint.d.ts +1 -0
  175. package/dist/types/util/weakMemo.d.ts +1 -0
  176. package/dist/types/version.d.ts +1 -1
  177. package/dist/util/SpatialIndex.js +8 -5
  178. package/dist/util/edgeCurves.js +7 -4
  179. package/dist/util/geometry.js +3 -1
  180. package/dist/util/groups.js +10 -0
  181. package/dist/util/hitDetection.js +9 -6
  182. package/dist/util/partitionPoint.js +6 -0
  183. package/dist/util/weakMemo.js +6 -0
  184. package/dist/version.js +1 -1
  185. package/dist/viewport.js +19 -17
  186. package/dist/walkEncoding.js +18 -14
  187. package/dist/walkHighlight.js +19 -15
  188. package/dist/walkKey.js +9 -6
  189. package/package.json +2 -2
  190. package/dist/chunks/chunk-ZY2AFL4Y.js +0 -43
@@ -1,20 +1,24 @@
1
1
  import {
2
2
  bundleRoutes
3
- } from "../chunks/chunk-5GMFEPWR.js";
4
- import "../chunks/chunk-XI3ZUTWF.js";
5
- import "../chunks/chunk-CWFNEKST.js";
3
+ } from "../chunks/chunk-TCOLHD7N.js";
4
+ import "../chunks/chunk-ZNJ4JAPT.js";
5
+ import "../chunks/chunk-A7H753BA.js";
6
6
  import "../chunks/chunk-C44YRP5F.js";
7
- import "../chunks/chunk-3XSM4OZU.js";
8
- import "../chunks/chunk-M6LK4GII.js";
9
- import "../chunks/chunk-RAEABQNW.js";
10
- import "../chunks/chunk-BPLIXMH7.js";
11
- import "../chunks/chunk-PTABFYQ7.js";
7
+ import "../chunks/chunk-5YILQTHE.js";
8
+ import "../chunks/chunk-IYALCKNB.js";
9
+ import "../chunks/chunk-LGJBB5BP.js";
10
+ import "../chunks/chunk-ZPZ7GNBT.js";
11
+ import "../chunks/chunk-PFRQBCGC.js";
12
+ import "../chunks/chunk-HMXCN6PG.js";
12
13
  import "../chunks/chunk-AB35TJST.js";
13
- import "../chunks/chunk-N73TI3VG.js";
14
- import "../chunks/chunk-24WITXQY.js";
15
- import "../chunks/chunk-7SRWJLSC.js";
14
+ import "../chunks/chunk-BZFOFC4D.js";
15
+ import "../chunks/chunk-CGMSKIEO.js";
16
+ import "../chunks/chunk-UBQSRCWR.js";
16
17
  import "../chunks/chunk-CPMHA4W5.js";
17
- import "../chunks/chunk-FYQWPEEY.js";
18
+ import "../chunks/chunk-BSOZW3FC.js";
19
+ import "../chunks/chunk-RJWPJCIB.js";
20
+ import "../chunks/chunk-HOAIOEI7.js";
21
+ import "../chunks/chunk-4XWZNPQ2.js";
18
22
  export {
19
23
  bundleRoutes
20
24
  };
@@ -1,351 +1,26 @@
1
1
  import {
2
- pathSteps
3
- } from "../chunks/chunk-XI3ZUTWF.js";
4
- import "../chunks/chunk-CWFNEKST.js";
2
+ coarsenTubeMap,
3
+ lastAtMost
4
+ } from "../chunks/chunk-A55OZVYX.js";
5
+ import "../chunks/chunk-ZNJ4JAPT.js";
6
+ import "../chunks/chunk-A7H753BA.js";
5
7
  import "../chunks/chunk-C44YRP5F.js";
6
- import "../chunks/chunk-3XSM4OZU.js";
7
- import {
8
- pathOrigin
9
- } from "../chunks/chunk-M6LK4GII.js";
10
- import {
11
- PathVisitsBuilder
12
- } from "../chunks/chunk-RAEABQNW.js";
13
- import "../chunks/chunk-BPLIXMH7.js";
14
- import "../chunks/chunk-PTABFYQ7.js";
8
+ import "../chunks/chunk-5YILQTHE.js";
9
+ import "../chunks/chunk-IYALCKNB.js";
10
+ import "../chunks/chunk-LGJBB5BP.js";
11
+ import "../chunks/chunk-ZPZ7GNBT.js";
12
+ import "../chunks/chunk-PFRQBCGC.js";
13
+ import "../chunks/chunk-HMXCN6PG.js";
15
14
  import "../chunks/chunk-AB35TJST.js";
16
- import "../chunks/chunk-N73TI3VG.js";
17
- import "../chunks/chunk-24WITXQY.js";
18
- import "../chunks/chunk-7SRWJLSC.js";
15
+ import "../chunks/chunk-BZFOFC4D.js";
16
+ import "../chunks/chunk-CGMSKIEO.js";
17
+ import "../chunks/chunk-UBQSRCWR.js";
19
18
  import "../chunks/chunk-CPMHA4W5.js";
20
- import {
21
- isBackbone
22
- } from "../chunks/chunk-FYQWPEEY.js";
23
-
24
- // src/tubeMap/coarsen.ts
25
- function readWalk(path, steps) {
26
- const runs = [];
27
- const open = () => ({ a: -Infinity, b: Infinity, detour: [], bp: 0 });
28
- const gaps = [open()];
29
- for (const { node, strand } of steps) {
30
- const gap = gaps.at(-1);
31
- const last = runs.at(-1);
32
- if (isBackbone(node)) {
33
- const start = node.stable.start;
34
- const end = start + node.length;
35
- const inverted = strand !== (node.stable.strand ?? "+");
36
- const joins = last !== void 0 && gap.detour.length === 0;
37
- if (joins && !last.inverted && !inverted && last.end === start) {
38
- last.end = end;
39
- } else if (joins && last.inverted && inverted && last.start === end) {
40
- last.start = start;
41
- } else {
42
- runs.push({ start, end, inverted });
43
- gaps.push(open());
44
- }
45
- } else {
46
- gap.detour.push({ node, strand });
47
- gap.bp += node.length;
48
- }
49
- }
50
- runs.forEach((run, i) => {
51
- gaps[i].b = run.start;
52
- gaps[i + 1].a = run.end;
53
- });
54
- return { path, runs, gaps };
55
- }
56
- function isSmall(gap, sigma) {
57
- return Number.isFinite(gap.a) && Number.isFinite(gap.b) && gap.b >= gap.a && Math.max(gap.b - gap.a, gap.bp) < sigma;
58
- }
59
- function strictlyInside(cuts, a, b) {
60
- let lo = 0;
61
- let hi = cuts.length;
62
- while (lo < hi) {
63
- const mid = lo + hi >> 1;
64
- if (cuts[mid] <= a) {
65
- lo = mid + 1;
66
- } else {
67
- hi = mid;
68
- }
69
- }
70
- return lo < cuts.length && cuts[lo] < b;
71
- }
72
- function referenceGaps(backbone) {
73
- const cuts = /* @__PURE__ */ new Map();
74
- for (let i = 1; i < backbone.length; i++) {
75
- const prev = backbone[i - 1];
76
- const node = backbone[i];
77
- const end = prev.stable.start + prev.length;
78
- if (end !== node.stable.start) {
79
- cuts.set(end, "reference-gap");
80
- cuts.set(node.stable.start, "reference-gap");
81
- }
82
- }
83
- return cuts;
84
- }
85
- function findCuts(walks, sigma, cuts) {
86
- const cut = (bp, cause) => {
87
- if (Number.isFinite(bp) && !cuts.has(bp)) {
88
- cuts.set(bp, cause);
89
- }
90
- };
91
- for (const { runs, gaps } of walks) {
92
- const first = runs[0];
93
- const last = runs.at(-1);
94
- if (first && last) {
95
- cut(first.start, "walk-end");
96
- cut(last.end, "walk-end");
97
- }
98
- for (const run of runs) {
99
- if (run.inverted && run.end - run.start >= sigma) {
100
- cut(run.start, "inversion");
101
- cut(run.end, "inversion");
102
- }
103
- }
104
- for (const gap of gaps) {
105
- if (gap.b < gap.a) {
106
- cut(gap.a, "backward");
107
- cut(gap.b, "backward");
108
- } else if (!isSmall(gap, sigma)) {
109
- cut(gap.a, "variant");
110
- cut(gap.b, "variant");
111
- }
112
- }
113
- }
114
- const spans = walks.flatMap(({ runs, gaps }) => [
115
- ...gaps.filter((g) => isSmall(g, sigma)).map((g) => [g.a, g.b]),
116
- ...runs.filter((r) => r.inverted && r.end - r.start < sigma).map((r) => [r.start, r.end])
117
- ]);
118
- for (let changed = true; changed; ) {
119
- changed = false;
120
- const sorted = [...cuts.keys()].sort((x, y) => x - y);
121
- for (const [a, b] of spans) {
122
- if (strictlyInside(sorted, a, b)) {
123
- changed ||= !cuts.has(a) || !cuts.has(b);
124
- cut(a, "overlap");
125
- cut(b, "overlap");
126
- }
127
- }
128
- }
129
- return cuts;
130
- }
131
- function mergeReference(backbone, cuts) {
132
- const merged = [];
133
- let group = [];
134
- const close = () => {
135
- const first = group[0];
136
- const last = group.at(-1);
137
- if (first && last) {
138
- const bp0 = first.stable.start;
139
- const bp1 = last.stable.start + last.length;
140
- const name = `${bp0}-${bp1}`;
141
- merged.push({
142
- node: {
143
- id: `${name}+`,
144
- name,
145
- length: bp1 - bp0,
146
- depth: 0,
147
- stable: {
148
- refName: first.stable.refName,
149
- start: bp0,
150
- rank: 0,
151
- strand: "+"
152
- }
153
- },
154
- members: group
155
- });
156
- }
157
- group = [];
158
- };
159
- for (const node of backbone) {
160
- if (cuts.has(node.stable.start)) {
161
- close();
162
- }
163
- group.push(node);
164
- }
165
- close();
166
- return merged;
167
- }
168
- function coarsenTubeMap(graph, sigma) {
169
- const backbone = graph.nodes.filter(isBackbone).sort((a, b) => a.stable.start - b.stable.start);
170
- const paths = graph.paths ?? [];
171
- if (backbone.length === 0 || paths.length === 0) {
172
- return void 0;
173
- }
174
- const stepsOf = pathSteps(graph);
175
- const walks = paths.map((p, i) => readWalk(p, stepsOf[i]));
176
- const cuts = findCuts(walks, sigma, referenceGaps(backbone));
177
- const sortedCuts = [...cuts.keys()].sort((x, y) => x - y);
178
- const merged = mergeReference(backbone, cuts);
179
- const starts = merged.map((m) => m.node.stable.start);
180
- const mergedAt = (bp) => {
181
- let lo = 0;
182
- let hi = starts.length - 1;
183
- while (lo < hi) {
184
- const mid = lo + hi + 1 >> 1;
185
- if (starts[mid] <= bp) {
186
- lo = mid;
187
- } else {
188
- hi = mid - 1;
189
- }
190
- }
191
- return lo;
192
- };
193
- const members = /* @__PURE__ */ new Map();
194
- const coarseOf = /* @__PURE__ */ new Map();
195
- for (const m of merged) {
196
- members.set(
197
- m.node.id,
198
- m.members.map((n) => n.id)
199
- );
200
- for (const n of m.members) {
201
- coarseOf.set(n.id, m.node.id);
202
- }
203
- }
204
- const routes = /* @__PURE__ */ new Map();
205
- const routeFor = (gap) => {
206
- const key = `${gap.a} ${gap.b} ${gap.detour.map((s) => s.node.id + s.strand).join(",")}`;
207
- let node = routes.get(key);
208
- if (!node) {
209
- const name = `r${routes.size}`;
210
- node = {
211
- id: `${name}+`,
212
- name,
213
- length: gap.bp,
214
- depth: 0,
215
- stable: {
216
- refName: backbone[0].stable.refName,
217
- start: Number.isFinite(gap.a) ? gap.a : gap.b,
218
- rank: 1,
219
- strand: "+"
220
- }
221
- };
222
- routes.set(key, node);
223
- members.set(
224
- node.id,
225
- gap.detour.map((s) => s.node.id)
226
- );
227
- for (const s of gap.detour) {
228
- if (!coarseOf.has(s.node.id)) {
229
- coarseOf.set(s.node.id, node.id);
230
- }
231
- }
232
- }
233
- return node;
234
- };
235
- const deviations = /* @__PURE__ */ new Map();
236
- const coarsePaths = [];
237
- const pathVisits = new PathVisitsBuilder();
238
- const edges = /* @__PURE__ */ new Map();
239
- const depth = /* @__PURE__ */ new Map();
240
- for (const { path, runs, gaps } of walks) {
241
- const steps = [];
242
- const devs = [];
243
- let continuing = false;
244
- const append = (node, strand) => {
245
- const last = steps.at(-1);
246
- if (!(continuing && last?.node === node && last.strand === strand)) {
247
- steps.push({ node, strand });
248
- }
249
- continuing = true;
250
- };
251
- const structural = (gap) => {
252
- continuing = false;
253
- if (gap.detour.length > 0) {
254
- append(routeFor(gap), "+");
255
- continuing = false;
256
- }
257
- };
258
- const fold = (gap) => {
259
- if (gap.bp > 0 || gap.b > gap.a) {
260
- devs.push({ start: gap.a, end: gap.b, bp: gap.bp });
261
- }
262
- const home = merged[mergedAt(gap.b > gap.a ? gap.a : gap.a - 1)];
263
- for (const s of gap.detour) {
264
- if (!coarseOf.has(s.node.id)) {
265
- coarseOf.set(s.node.id, home.node.id);
266
- members.get(home.node.id).push(s.node.id);
267
- }
268
- }
269
- };
270
- runs.forEach((run, i) => {
271
- const gap = gaps[i];
272
- if (isSmall(gap, sigma) && !strictlyInside(sortedCuts, gap.a, gap.b)) {
273
- fold(gap);
274
- } else {
275
- structural(gap);
276
- }
277
- const i0 = mergedAt(run.start);
278
- const i1 = mergedAt(run.end - 1);
279
- const small = run.end - run.start < sigma;
280
- if (run.inverted && !small) {
281
- continuing = false;
282
- for (let j = i1; j >= i0; j--) {
283
- append(merged[j].node, "-");
284
- }
285
- continuing = false;
286
- } else {
287
- if (run.inverted) {
288
- devs.push({
289
- start: run.start,
290
- end: run.end,
291
- bp: run.end - run.start,
292
- inverted: true
293
- });
294
- }
295
- for (let j = i0; j <= i1; j++) {
296
- append(merged[j].node, "+");
297
- }
298
- }
299
- });
300
- structural(gaps.at(-1));
301
- const visitName = pathOrigin(path.name).name;
302
- let pos = 0;
303
- steps.forEach(({ node, strand }, i) => {
304
- depth.set(node.id, (depth.get(node.id) ?? 0) + 1);
305
- pathVisits.add(
306
- node.name,
307
- visitName,
308
- path.sample ?? visitName,
309
- pos,
310
- strand === "-"
311
- );
312
- pos += node.length;
313
- const next = steps[i + 1];
314
- if (next) {
315
- const key = `${node.id}${strand} ${next.node.id}${next.strand}`;
316
- const edge = edges.get(key);
317
- if (edge) {
318
- edge.pathIds.push(path.name);
319
- } else {
320
- edges.set(key, {
321
- from: node.id,
322
- to: next.node.id,
323
- fromStrand: strand,
324
- toStrand: next.strand,
325
- pathIds: [path.name]
326
- });
327
- }
328
- }
329
- });
330
- coarsePaths.push({ ...path, nodeIds: steps.map((s) => s.node.id) });
331
- deviations.set(path.name, devs);
332
- }
333
- const nodes = [...merged.map((m) => m.node), ...routes.values()].filter((n) => depth.has(n.id)).map((n) => ({ ...n, depth: depth.get(n.id) }));
334
- return {
335
- graph: {
336
- ...graph,
337
- nodes,
338
- edges: [...edges.values()],
339
- paths: coarsePaths,
340
- pathVisits: pathVisits.build(),
341
- reads: void 0
342
- },
343
- members,
344
- coarseOf,
345
- deviations,
346
- cuts
347
- };
348
- }
19
+ import "../chunks/chunk-BSOZW3FC.js";
20
+ import "../chunks/chunk-RJWPJCIB.js";
21
+ import "../chunks/chunk-HOAIOEI7.js";
22
+ import "../chunks/chunk-4XWZNPQ2.js";
349
23
  export {
350
- coarsenTubeMap
24
+ coarsenTubeMap,
25
+ lastAtMost
351
26
  };
@@ -1,36 +1,68 @@
1
1
  import {
2
+ lastAtMost
3
+ } from "../chunks/chunk-A55OZVYX.js";
4
+ import {
5
+ boxX,
2
6
  referenceBoxes
3
- } from "../chunks/chunk-XMW266WS.js";
4
- import "../chunks/chunk-PTABFYQ7.js";
7
+ } from "../chunks/chunk-HSVFCA2J.js";
8
+ import "../chunks/chunk-ZNJ4JAPT.js";
9
+ import "../chunks/chunk-A7H753BA.js";
10
+ import "../chunks/chunk-C44YRP5F.js";
11
+ import "../chunks/chunk-5YILQTHE.js";
12
+ import "../chunks/chunk-IYALCKNB.js";
13
+ import {
14
+ pathOrigin
15
+ } from "../chunks/chunk-LGJBB5BP.js";
16
+ import "../chunks/chunk-ZPZ7GNBT.js";
17
+ import "../chunks/chunk-PFRQBCGC.js";
18
+ import "../chunks/chunk-HMXCN6PG.js";
5
19
  import "../chunks/chunk-AB35TJST.js";
6
- import "../chunks/chunk-N73TI3VG.js";
7
- import "../chunks/chunk-24WITXQY.js";
8
- import "../chunks/chunk-7SRWJLSC.js";
20
+ import "../chunks/chunk-BZFOFC4D.js";
21
+ import "../chunks/chunk-CGMSKIEO.js";
22
+ import "../chunks/chunk-UBQSRCWR.js";
9
23
  import "../chunks/chunk-CPMHA4W5.js";
10
- import "../chunks/chunk-FYQWPEEY.js";
24
+ import {
25
+ pushTo
26
+ } from "../chunks/chunk-BSOZW3FC.js";
27
+ import "../chunks/chunk-RJWPJCIB.js";
28
+ import "../chunks/chunk-HOAIOEI7.js";
29
+ import "../chunks/chunk-4XWZNPQ2.js";
11
30
 
12
31
  // src/tubeMap/deviations.ts
13
32
  function boxAt(boxes, bp, visited) {
14
- let lo = 0;
15
- let hi = boxes.length - 1;
16
- while (lo < hi) {
17
- const mid = lo + hi + 1 >> 1;
18
- if (boxes[mid].bp0 <= bp) {
19
- lo = mid;
20
- } else {
21
- hi = mid - 1;
22
- }
23
- }
33
+ const lo = lastAtMost(boxes, (box) => box.bp0, bp);
24
34
  return [boxes[lo], boxes[lo - 1]].find(
25
35
  (box) => box && visited.has(box.index) && box.bp0 <= bp && bp <= box.bp1
26
36
  );
27
37
  }
38
+ function deviationsByPath(graph, deviations) {
39
+ const byWalk = /* @__PURE__ */ new Map();
40
+ for (const [name, devs] of deviations) {
41
+ for (const d of devs) {
42
+ pushTo(byWalk, pathOrigin(name).name, d);
43
+ }
44
+ }
45
+ const byPath = /* @__PURE__ */ new Map();
46
+ for (const path of graph.paths ?? []) {
47
+ if (path.members) {
48
+ const union = /* @__PURE__ */ new Map();
49
+ for (const d of path.members.flatMap((m) => byWalk.get(m) ?? [])) {
50
+ union.set(`${d.start}-${d.end}`, d);
51
+ }
52
+ byPath.set(path.name, [...union.values()]);
53
+ } else {
54
+ byPath.set(path.name, deviations.get(path.name) ?? []);
55
+ }
56
+ }
57
+ return byPath;
58
+ }
28
59
  function deviationMarks(drawing, deviations) {
29
- const { layout } = drawing;
60
+ const { graph, layout } = drawing;
30
61
  const boxes = [...referenceBoxes(drawing).values()].flat().sort((a, b) => a.bp0 - b.bp0);
62
+ const byPath = deviationsByPath(graph, deviations);
31
63
  const marks = [];
32
64
  for (const track of layout.tracks) {
33
- const devs = track.name ? deviations.get(track.name) : void 0;
65
+ const devs = track.name ? byPath.get(track.name) : void 0;
34
66
  if (devs?.length) {
35
67
  const visited = /* @__PURE__ */ new Map();
36
68
  for (const segment of track.path) {
@@ -41,12 +73,10 @@ function deviationMarks(drawing, deviations) {
41
73
  for (const d of devs) {
42
74
  const box = boxAt(boxes, d.start, visited);
43
75
  if (box) {
44
- const { bp0, bp1, x0, x1, index } = box;
45
- const at = (bp) => x0 + (Math.min(bp, bp1) - bp0) / (bp1 - bp0) * (x1 - x0);
46
76
  marks.push({
47
- x0: at(d.start),
48
- x1: at(d.end),
49
- y: visited.get(index),
77
+ x0: boxX(box, d.start),
78
+ x1: boxX(box, d.end),
79
+ y: visited.get(box.index),
50
80
  height: track.width
51
81
  });
52
82
  }
@@ -8,15 +8,18 @@ import {
8
8
  tubeMapMismatchAt,
9
9
  tubeMapPicture,
10
10
  tubeMapTrackAt
11
- } from "../chunks/chunk-B6IW3IRE.js";
12
- import "../chunks/chunk-XXT2JCAC.js";
13
- import "../chunks/chunk-HK5Y2YWV.js";
14
- import "../chunks/chunk-CWFNEKST.js";
11
+ } from "../chunks/chunk-65ILXEIN.js";
12
+ import "../chunks/chunk-BX4AUMZE.js";
13
+ import "../chunks/chunk-MYDHU2UV.js";
14
+ import "../chunks/chunk-A7H753BA.js";
15
15
  import "../chunks/chunk-C44YRP5F.js";
16
- import "../chunks/chunk-M6LK4GII.js";
17
- import "../chunks/chunk-RAEABQNW.js";
18
- import "../chunks/chunk-BPLIXMH7.js";
19
- import "../chunks/chunk-FYQWPEEY.js";
16
+ import "../chunks/chunk-5YILQTHE.js";
17
+ import "../chunks/chunk-LGJBB5BP.js";
18
+ import "../chunks/chunk-ZPZ7GNBT.js";
19
+ import "../chunks/chunk-BSOZW3FC.js";
20
+ import "../chunks/chunk-RJWPJCIB.js";
21
+ import "../chunks/chunk-HOAIOEI7.js";
22
+ import "../chunks/chunk-4XWZNPQ2.js";
20
23
  export {
21
24
  drawTubeMap,
22
25
  drawTubeMapHighlight,
@@ -1,8 +1,9 @@
1
1
  import {
2
2
  tubeMapFrame,
3
3
  tubeMapNodeAt
4
- } from "../chunks/chunk-2CMYNW3E.js";
5
- import "../chunks/chunk-HK5Y2YWV.js";
4
+ } from "../chunks/chunk-SI3BB4DO.js";
5
+ import "../chunks/chunk-MYDHU2UV.js";
6
+ import "../chunks/chunk-HOAIOEI7.js";
6
7
  export {
7
8
  tubeMapFrame,
8
9
  tubeMapNodeAt
@@ -3,15 +3,18 @@ import {
3
3
  drawTubeMapGenes,
4
4
  tubeMapGeneRows,
5
5
  tubeMapGenes
6
- } from "../chunks/chunk-4DRQTBBE.js";
7
- import "../chunks/chunk-XMW266WS.js";
8
- import "../chunks/chunk-PTABFYQ7.js";
6
+ } from "../chunks/chunk-6EOB6ANX.js";
7
+ import "../chunks/chunk-HSVFCA2J.js";
8
+ import "../chunks/chunk-PFRQBCGC.js";
9
+ import "../chunks/chunk-HMXCN6PG.js";
9
10
  import "../chunks/chunk-AB35TJST.js";
10
- import "../chunks/chunk-N73TI3VG.js";
11
- import "../chunks/chunk-24WITXQY.js";
12
- import "../chunks/chunk-7SRWJLSC.js";
11
+ import "../chunks/chunk-BZFOFC4D.js";
12
+ import "../chunks/chunk-CGMSKIEO.js";
13
+ import "../chunks/chunk-UBQSRCWR.js";
13
14
  import "../chunks/chunk-CPMHA4W5.js";
14
- import "../chunks/chunk-FYQWPEEY.js";
15
+ import "../chunks/chunk-BSOZW3FC.js";
16
+ import "../chunks/chunk-RJWPJCIB.js";
17
+ import "../chunks/chunk-HOAIOEI7.js";
15
18
  export {
16
19
  GENE_ROW_PX,
17
20
  drawTubeMapGenes,
@@ -1,6 +1,6 @@
1
1
  import {
2
2
  tubeMapMismatches
3
- } from "../chunks/chunk-XXT2JCAC.js";
3
+ } from "../chunks/chunk-BX4AUMZE.js";
4
4
  export {
5
5
  tubeMapMismatches
6
6
  };
@@ -1,22 +1,26 @@
1
1
  import {
2
2
  tubeMapNodeColors
3
- } from "../chunks/chunk-T2PCSAFY.js";
4
- import "../chunks/chunk-A4PFFC4Q.js";
3
+ } from "../chunks/chunk-7XR26JST.js";
4
+ import "../chunks/chunk-2XEQR54X.js";
5
5
  import "../chunks/chunk-G2FD7S3C.js";
6
- import "../chunks/chunk-ZY2AFL4Y.js";
7
- import "../chunks/chunk-3XSM4OZU.js";
8
- import "../chunks/chunk-JBCZIBFC.js";
9
- import "../chunks/chunk-M6LK4GII.js";
10
- import "../chunks/chunk-RAEABQNW.js";
11
- import "../chunks/chunk-ES2N6V56.js";
12
- import "../chunks/chunk-BPLIXMH7.js";
13
- import "../chunks/chunk-PTABFYQ7.js";
6
+ import "../chunks/chunk-GQZVR5DW.js";
7
+ import "../chunks/chunk-WA5WNP6V.js";
8
+ import "../chunks/chunk-IYALCKNB.js";
9
+ import "../chunks/chunk-CRCD3ODN.js";
10
+ import "../chunks/chunk-LGJBB5BP.js";
11
+ import "../chunks/chunk-VX656OAZ.js";
12
+ import "../chunks/chunk-ZPZ7GNBT.js";
13
+ import "../chunks/chunk-PFRQBCGC.js";
14
+ import "../chunks/chunk-HMXCN6PG.js";
14
15
  import "../chunks/chunk-AB35TJST.js";
15
- import "../chunks/chunk-N73TI3VG.js";
16
- import "../chunks/chunk-24WITXQY.js";
17
- import "../chunks/chunk-7SRWJLSC.js";
16
+ import "../chunks/chunk-BZFOFC4D.js";
17
+ import "../chunks/chunk-CGMSKIEO.js";
18
+ import "../chunks/chunk-UBQSRCWR.js";
18
19
  import "../chunks/chunk-CPMHA4W5.js";
19
- import "../chunks/chunk-FYQWPEEY.js";
20
+ import "../chunks/chunk-BSOZW3FC.js";
21
+ import "../chunks/chunk-RJWPJCIB.js";
22
+ import "../chunks/chunk-HOAIOEI7.js";
23
+ import "../chunks/chunk-4XWZNPQ2.js";
20
24
  export {
21
25
  tubeMapNodeColors
22
26
  };
@@ -5,7 +5,7 @@ import {
5
5
  readColor,
6
6
  stepName,
7
7
  tubeMapReads
8
- } from "../chunks/chunk-CWFNEKST.js";
8
+ } from "../chunks/chunk-A7H753BA.js";
9
9
  import "../chunks/chunk-C44YRP5F.js";
10
10
  export {
11
11
  FORWARD_READ_COLORS,
@@ -2,7 +2,8 @@ import {
2
2
  MAX_GAP_PX,
3
3
  referenceKnots,
4
4
  warpX
5
- } from "../chunks/chunk-HK5Y2YWV.js";
5
+ } from "../chunks/chunk-MYDHU2UV.js";
6
+ import "../chunks/chunk-HOAIOEI7.js";
6
7
  export {
7
8
  MAX_GAP_PX,
8
9
  referenceKnots,
@@ -1,4 +1,5 @@
1
1
  import type { Graph } from '../types.js';
2
+ export { parsePanSN } from '../pansn.js';
2
3
  export interface ProjectedAllele {
3
4
  start: number;
4
5
  refSpan: number;
@@ -6,8 +7,4 @@ export interface ProjectedAllele {
6
7
  nodeOffsets: number[];
7
8
  pathLength: number;
8
9
  }
9
- export declare function parsePanSN(refName: string): {
10
- sample: string;
11
- haplotype: number | undefined;
12
- };
13
10
  export declare function projectAlleles(graph: Graph): ProjectedAllele[];
@@ -12,7 +12,5 @@ export interface BubbleHalo {
12
12
  whole: boolean;
13
13
  tick: boolean;
14
14
  }
15
- export declare const WHOLE_FRACTION = 0.9;
16
15
  export declare function coversGraph(members: number, nodes: number): boolean;
17
- export declare const SMALL_VARIANT_BP = 50;
18
16
  export declare function bubbleHalos(graph: Graph, bubbles: MinigraphBubble[], positions: Record<string, NodeSegment[]>, repeats?: readonly RepeatSpan[]): BubbleHalo[];