jbrowse-plugin-graphgenomeviewer 7.1.1 → 7.1.2

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 (230) hide show
  1. package/dist/chunks/AddTrackWorkflow-RPH4ZFSX.js +2 -0
  2. package/dist/chunks/{AddTrackWorkflow-Q6RY23LF.js.map → AddTrackWorkflow-RPH4ZFSX.js.map} +2 -2
  3. package/dist/chunks/{ChooseSamplesDialog-BNVYBGJQ.js → ChooseSamplesDialog-JRN4YLBY.js} +2 -2
  4. package/dist/chunks/{ChooseWalksDialog-EN3C6AQA.js → ChooseWalksDialog-UXLFZTZJ.js} +2 -2
  5. package/dist/chunks/GbzBaseSyntenyAdapter-LRECKL74.js +2 -0
  6. package/dist/chunks/GbzBaseSyntenyAdapter-LRECKL74.js.map +7 -0
  7. package/dist/chunks/{GraphGenomeView-BG7UHKCT.js → GraphGenomeView-NJ5XCREB.js} +2 -2
  8. package/dist/chunks/GraphTrackSettingsDialog-47PQ6LM7.js +6 -0
  9. package/dist/chunks/GraphTrackSettingsDialog-47PQ6LM7.js.map +7 -0
  10. package/dist/chunks/{HighlightColorDialog-YT55H62Y.js → HighlightColorDialog-UK5VDDV6.js} +2 -2
  11. package/dist/chunks/{LinearGraphDisplay-SPILNH4J.js → LinearGraphDisplay-TXTKAACH.js} +2 -2
  12. package/dist/chunks/MinigraphBubbleAdapter-HW7XQMAI.js +2 -0
  13. package/dist/chunks/MinigraphBubbleAdapter-HW7XQMAI.js.map +7 -0
  14. package/dist/chunks/RgfaTabixAdapter-PSKYWD2X.js +3 -0
  15. package/dist/chunks/RgfaTabixAdapter-PSKYWD2X.js.map +7 -0
  16. package/dist/chunks/{TrackWarningsDialog-QT4YEXGD.js → TrackWarningsDialog-TRYBOIHX.js} +2 -2
  17. package/dist/chunks/WalkTabixSyntenyAdapter-U53JLYVU.js +2 -0
  18. package/dist/chunks/WalkTabixSyntenyAdapter-U53JLYVU.js.map +7 -0
  19. package/dist/chunks/chunk-2MCOY3YG.js +3 -0
  20. package/dist/chunks/chunk-2MCOY3YG.js.map +7 -0
  21. package/dist/chunks/{chunk-S3QFRJR3.js → chunk-5O6YRJHL.js} +2 -2
  22. package/dist/chunks/{chunk-S3QFRJR3.js.map → chunk-5O6YRJHL.js.map} +1 -1
  23. package/dist/chunks/chunk-BUY2PV7D.js +2 -0
  24. package/dist/chunks/chunk-BUY2PV7D.js.map +7 -0
  25. package/dist/chunks/chunk-DVCZ3BOC.js +2 -0
  26. package/dist/chunks/chunk-DVCZ3BOC.js.map +7 -0
  27. package/dist/chunks/{chunk-WI7LRX6I.js → chunk-EGDLEQ47.js} +3 -3
  28. package/dist/chunks/{chunk-WI7LRX6I.js.map → chunk-EGDLEQ47.js.map} +4 -4
  29. package/dist/chunks/chunk-FHZXW7IG.js +2 -0
  30. package/dist/chunks/chunk-FHZXW7IG.js.map +7 -0
  31. package/dist/chunks/{chunk-SIGBYUXP.js → chunk-FITPB5X6.js} +2 -2
  32. package/dist/chunks/{chunk-SIGBYUXP.js.map → chunk-FITPB5X6.js.map} +1 -1
  33. package/dist/chunks/chunk-FVK2THE4.js +2 -0
  34. package/dist/chunks/chunk-FVK2THE4.js.map +7 -0
  35. package/dist/chunks/{chunk-AHFQHYJI.js → chunk-IBWORBGY.js} +2 -2
  36. package/dist/chunks/chunk-MND7GAOG.js +2 -0
  37. package/dist/chunks/chunk-MND7GAOG.js.map +7 -0
  38. package/dist/chunks/chunk-MURN3JPK.js +2 -0
  39. package/dist/chunks/chunk-MURN3JPK.js.map +7 -0
  40. package/dist/chunks/{chunk-JT3RXYIL.js → chunk-N6PQZSAX.js} +1 -1
  41. package/dist/chunks/{chunk-JT3RXYIL.js.map → chunk-N6PQZSAX.js.map} +1 -1
  42. package/dist/chunks/chunk-NWJYLP6Z.js +4 -0
  43. package/dist/chunks/chunk-NWJYLP6Z.js.map +7 -0
  44. package/dist/chunks/{chunk-FAVU7473.js → chunk-S75W27NT.js} +2 -2
  45. package/dist/chunks/{chunk-FAVU7473.js.map → chunk-S75W27NT.js.map} +1 -1
  46. package/dist/chunks/{chunk-77GS3NVS.js → chunk-TUKJRVBS.js} +1 -1
  47. package/dist/chunks/{chunk-77GS3NVS.js.map → chunk-TUKJRVBS.js.map} +1 -1
  48. package/dist/chunks/chunk-UGCNPRKC.js +2 -0
  49. package/dist/chunks/{chunk-E3WYW7Z6.js.map → chunk-UGCNPRKC.js.map} +2 -2
  50. package/dist/chunks/{chunk-GUNY323E.js → chunk-VJ3SO4CH.js} +1 -1
  51. package/dist/chunks/{chunk-GUNY323E.js.map → chunk-VJ3SO4CH.js.map} +1 -1
  52. package/dist/chunks/{chunk-ABNOMXT3.js → chunk-W6EQLQ6H.js} +1 -1
  53. package/dist/chunks/{chunk-ABNOMXT3.js.map → chunk-W6EQLQ6H.js.map} +1 -1
  54. package/dist/chunks/chunk-WAMBI7OX.js +2 -0
  55. package/dist/chunks/chunk-WAMBI7OX.js.map +7 -0
  56. package/dist/chunks/chunk-WCGVUZWP.js +2 -0
  57. package/dist/chunks/chunk-WCGVUZWP.js.map +7 -0
  58. package/dist/chunks/chunk-WEGWTC52.js +2 -0
  59. package/dist/chunks/chunk-WEGWTC52.js.map +7 -0
  60. package/dist/chunks/chunk-YU7T5WOQ.js +2 -0
  61. package/dist/chunks/{chunk-6OPB7FOJ.js.map → chunk-YU7T5WOQ.js.map} +2 -2
  62. package/dist/chunks/{chunk-IIGYYZI5.js → chunk-ZMOMFU2B.js} +1 -1
  63. package/dist/chunks/{chunk-IIGYYZI5.js.map → chunk-ZMOMFU2B.js.map} +1 -1
  64. package/dist/chunks/renderSvg-KSQI5DH6.js +2 -0
  65. package/dist/chunks/renderSvg-KSQI5DH6.js.map +7 -0
  66. package/dist/jbrowse-plugin-graphgenomeviewer.esm.js +4 -4
  67. package/dist/jbrowse-plugin-graphgenomeviewer.esm.js.map +4 -4
  68. package/package.json +4 -5
  69. package/src/GbzBaseSyntenyAdapter/GbzBaseSyntenyAdapter.test.ts +12 -10
  70. package/src/GbzBaseSyntenyAdapter/GbzBaseSyntenyAdapter.ts +31 -39
  71. package/src/GbzBaseSyntenyAdapter/configSchema.ts +5 -29
  72. package/src/GetGraphHaplotypes.ts +6 -15
  73. package/src/GetGraphReads.ts +7 -14
  74. package/src/GetSubgraph.ts +11 -20
  75. package/src/GraphAddTrackWorkflow/AddTrackWorkflow.test.tsx +73 -0
  76. package/src/GraphAddTrackWorkflow/AddTrackWorkflow.tsx +1 -0
  77. package/src/GraphAddTrackWorkflow/buildTrackConfig.test.ts +3 -2
  78. package/src/GraphAddTrackWorkflow/buildTrackConfig.ts +2 -11
  79. package/src/GraphAddTrackWorkflow/graphManifest.ts +15 -3
  80. package/src/GraphGenomeView/components/BubbleHalos.tsx +5 -24
  81. package/src/GraphGenomeView/components/ElTree.tsx +6 -1
  82. package/src/GraphGenomeView/components/FacetPanels.test.tsx +65 -0
  83. package/src/GraphGenomeView/components/FacetPanels.tsx +17 -12
  84. package/src/GraphGenomeView/components/GenePins.tsx +5 -22
  85. package/src/GraphGenomeView/components/GraphCanvas.tsx +40 -234
  86. package/src/GraphGenomeView/components/GraphSettingsDialog.tsx +102 -108
  87. package/src/GraphGenomeView/components/GraphToolbar.tsx +5 -6
  88. package/src/GraphGenomeView/components/HoverLayer.tsx +3 -7
  89. package/src/GraphGenomeView/components/HoverTooltip.test.tsx +39 -0
  90. package/src/GraphGenomeView/components/HoverTooltip.tsx +122 -0
  91. package/src/GraphGenomeView/components/LabelChip.tsx +3 -5
  92. package/src/GraphGenomeView/components/LabelLayer.tsx +4 -17
  93. package/src/GraphGenomeView/components/LabelledSelect.tsx +1 -1
  94. package/src/GraphGenomeView/components/Legend.tsx +14 -23
  95. package/src/GraphGenomeView/components/PaneContextMenu.test.tsx +48 -0
  96. package/src/GraphGenomeView/components/PaneContextMenu.tsx +144 -0
  97. package/src/GraphGenomeView/components/PaneLayer.tsx +91 -0
  98. package/src/GraphGenomeView/components/PaneSvg.tsx +32 -50
  99. package/src/GraphGenomeView/components/ReferenceStripOverlay.test.tsx +48 -0
  100. package/src/GraphGenomeView/components/ReferenceStripOverlay.tsx +13 -42
  101. package/src/GraphGenomeView/components/SettingsMenu.test.tsx +29 -0
  102. package/src/GraphGenomeView/components/SettingsMenu.tsx +56 -65
  103. package/src/GraphGenomeView/components/TubeMapOverlay.tsx +15 -46
  104. package/src/GraphGenomeView/components/WalkRowsOverlay.tsx +58 -132
  105. package/src/GraphGenomeView/components/WalkStrip.tsx +17 -20
  106. package/src/GraphGenomeView/components/WalkStripLocator.tsx +8 -15
  107. package/src/GraphGenomeView/components/usePaneGestures.test.tsx +35 -0
  108. package/src/GraphGenomeView/components/usePaneGestures.ts +88 -91
  109. package/src/GraphGenomeView/genes/geneFeatures.ts +27 -45
  110. package/src/GraphGenomeView/graphLayers.ts +4 -4
  111. package/src/GraphGenomeView/host.ts +1 -1
  112. package/src/GraphGenomeView/legendSections.test.ts +23 -1
  113. package/src/GraphGenomeView/legendSections.ts +68 -74
  114. package/src/GraphGenomeView/model.test.ts +373 -39
  115. package/src/GraphGenomeView/model.ts +67 -11
  116. package/src/GraphGenomeView/pane/fitViews.ts +57 -514
  117. package/src/GraphGenomeView/pane/graphViews.ts +219 -160
  118. package/src/GraphGenomeView/pane/hostViews.ts +2 -10
  119. package/src/GraphGenomeView/pane/keepInView.ts +5 -4
  120. package/src/GraphGenomeView/pane/launchActions.ts +117 -139
  121. package/src/GraphGenomeView/pane/launchViews.ts +8 -27
  122. package/src/GraphGenomeView/pane/loadActions.ts +208 -187
  123. package/src/GraphGenomeView/pane/morph.ts +3 -9
  124. package/src/GraphGenomeView/pane/paneBase.ts +22 -11
  125. package/src/GraphGenomeView/pane/renderingActions.ts +65 -11
  126. package/src/GraphGenomeView/pane/settingActions.ts +17 -42
  127. package/src/GraphGenomeView/pane/tubeMapViews.ts +445 -0
  128. package/src/GraphGenomeView/pane/walkRowViews.ts +29 -37
  129. package/src/GraphGenomeView/renderPipeline.test.ts +1 -1
  130. package/src/GraphGenomeView/repeats/repeatFeatures.test.ts +23 -0
  131. package/src/GraphGenomeView/repeats/repeatFeatures.ts +16 -23
  132. package/src/GraphGenomeView/repeats/walkCalls.test.ts +18 -0
  133. package/src/GraphGenomeView/repeats/walkCalls.ts +46 -38
  134. package/src/GraphGenomeView/subgraphLoad.test.ts +67 -3
  135. package/src/GraphGenomeView/trackFeatures.ts +24 -0
  136. package/src/GraphGenomeView/tubeMapPanels.ts +28 -36
  137. package/src/GraphGenomeView/viewModel.ts +8 -56
  138. package/src/GraphGenomeView/walkRowGroups.ts +104 -30
  139. package/src/GraphGenomeView/walkRowSpans.test.ts +6 -0
  140. package/src/GraphGenomeView/walkRowSpans.ts +22 -31
  141. package/src/GraphGenomeView/walkRowsCut.ts +8 -17
  142. package/src/LinearGraphDisplay/LinearGraphDisplay.test.ts +173 -26
  143. package/src/LinearGraphDisplay/components/SubgraphContextSelect.test.tsx +1 -1
  144. package/src/LinearGraphDisplay/components/SubgraphContextSelect.tsx +3 -3
  145. package/src/LinearGraphDisplay/components/SubgraphHaplotypesField.test.tsx +3 -3
  146. package/src/LinearGraphDisplay/components/SubgraphHaplotypesField.tsx +7 -10
  147. package/src/LinearGraphDisplay/configSchema.ts +0 -3
  148. package/src/LinearGraphDisplay/denseWindow.ts +8 -11
  149. package/src/LinearGraphDisplay/model.ts +89 -82
  150. package/src/MinigraphBubbleAdapter/MinigraphBubbleAdapter.test.ts +21 -0
  151. package/src/MinigraphBubbleAdapter/MinigraphBubbleAdapter.ts +6 -14
  152. package/src/MinigraphBubbleAdapter/configSchema.ts +3 -9
  153. package/src/MinigraphBubbleAdapter/index.ts +1 -2
  154. package/src/RgfaTabixAdapter/RgfaTabixAdapter.test.ts +45 -1
  155. package/src/RgfaTabixAdapter/RgfaTabixAdapter.ts +111 -276
  156. package/src/RgfaTabixAdapter/configSchema.ts +15 -49
  157. package/src/RgfaTabixAdapter/index.ts +1 -2
  158. package/src/RgfaTabixAdapter/rgfaBed.ts +25 -15
  159. package/src/RgfaTabixAdapter/rgfaCut.ts +139 -0
  160. package/src/RgfaTabixAdapter/test_data/walks2_twopart.GRCh38.walks.bed.gz +0 -0
  161. package/src/RgfaTabixAdapter/test_data/walks2_twopart.GRCh38.walks.bed.gz.tbi +0 -0
  162. package/src/RgfaTabixAdapter/walkFiles.test.ts +33 -0
  163. package/src/RgfaTabixAdapter/walkReader.ts +103 -27
  164. package/src/RgfaTabixAdapter/walkRowRuns.test.ts +53 -0
  165. package/src/RgfaTabixAdapter/walkRowRuns.ts +42 -24
  166. package/src/WalkTabixSyntenyAdapter/WalkTabixSyntenyAdapter.test.ts +37 -6
  167. package/src/WalkTabixSyntenyAdapter/WalkTabixSyntenyAdapter.ts +23 -45
  168. package/src/WalkTabixSyntenyAdapter/configSchema.ts +11 -24
  169. package/src/WalkTabixSyntenyAdapter/index.ts +1 -2
  170. package/src/WalkTabixSyntenyAdapter/walkLanes.ts +35 -43
  171. package/src/adapterImplementing.ts +27 -0
  172. package/src/graphTrackDefaults/index.ts +14 -20
  173. package/src/hoverSync/lgvHover.ts +2 -4
  174. package/src/launchFromGraph/contributors.ts +1 -5
  175. package/src/launchFromGraph/graphMenuItems.ts +3 -4
  176. package/src/launchFromGraph/launchFromGraph.test.ts +0 -1
  177. package/src/launchFromGraph/launchFromGraph.ts +8 -13
  178. package/src/panSNTabix.test.ts +1 -26
  179. package/src/panSNTabix.ts +54 -51
  180. package/src/synteny/lanePairs.ts +8 -0
  181. package/src/synteny/matchFeature.ts +52 -0
  182. package/src/tabixRange.ts +13 -10
  183. package/src/version.ts +1 -1
  184. package/dist/chunks/AddTrackWorkflow-Q6RY23LF.js +0 -2
  185. package/dist/chunks/GbzBaseSyntenyAdapter-A56NDFW7.js +0 -2
  186. package/dist/chunks/GbzBaseSyntenyAdapter-A56NDFW7.js.map +0 -7
  187. package/dist/chunks/GraphTrackSettingsDialog-YOJF7KCT.js +0 -6
  188. package/dist/chunks/GraphTrackSettingsDialog-YOJF7KCT.js.map +0 -7
  189. package/dist/chunks/MinigraphBubbleAdapter-2XGZPDI4.js +0 -2
  190. package/dist/chunks/MinigraphBubbleAdapter-2XGZPDI4.js.map +0 -7
  191. package/dist/chunks/RgfaTabixAdapter-ZTFVNM4G.js +0 -3
  192. package/dist/chunks/RgfaTabixAdapter-ZTFVNM4G.js.map +0 -7
  193. package/dist/chunks/WalkTabixSyntenyAdapter-I6E6HAN4.js +0 -2
  194. package/dist/chunks/WalkTabixSyntenyAdapter-I6E6HAN4.js.map +0 -7
  195. package/dist/chunks/chunk-2QQDYK5V.js +0 -2
  196. package/dist/chunks/chunk-2QQDYK5V.js.map +0 -7
  197. package/dist/chunks/chunk-4IJKHORO.js +0 -2
  198. package/dist/chunks/chunk-4IJKHORO.js.map +0 -7
  199. package/dist/chunks/chunk-6CAAGNOP.js +0 -4
  200. package/dist/chunks/chunk-6CAAGNOP.js.map +0 -7
  201. package/dist/chunks/chunk-6OPB7FOJ.js +0 -2
  202. package/dist/chunks/chunk-CPCTMFXV.js +0 -2
  203. package/dist/chunks/chunk-CPCTMFXV.js.map +0 -7
  204. package/dist/chunks/chunk-E3WYW7Z6.js +0 -2
  205. package/dist/chunks/chunk-EHCF6AXH.js +0 -2
  206. package/dist/chunks/chunk-EHCF6AXH.js.map +0 -7
  207. package/dist/chunks/chunk-FBGWMI6M.js +0 -2
  208. package/dist/chunks/chunk-FBGWMI6M.js.map +0 -7
  209. package/dist/chunks/chunk-G63ZHV6F.js +0 -3
  210. package/dist/chunks/chunk-G63ZHV6F.js.map +0 -7
  211. package/dist/chunks/chunk-LUTMKRXY.js +0 -2
  212. package/dist/chunks/chunk-LUTMKRXY.js.map +0 -7
  213. package/dist/chunks/chunk-NKBE5KYL.js +0 -2
  214. package/dist/chunks/chunk-NKBE5KYL.js.map +0 -7
  215. package/dist/chunks/chunk-TXI43PEN.js +0 -2
  216. package/dist/chunks/chunk-TXI43PEN.js.map +0 -7
  217. package/dist/chunks/chunk-WRH2EJKX.js +0 -2
  218. package/dist/chunks/chunk-WRH2EJKX.js.map +0 -7
  219. package/dist/chunks/chunk-YAADZ2AI.js +0 -2
  220. package/dist/chunks/chunk-YAADZ2AI.js.map +0 -7
  221. package/dist/chunks/renderSvg-B3WCFN65.js +0 -2
  222. package/dist/chunks/renderSvg-B3WCFN65.js.map +0 -7
  223. package/src/GraphGenomeView/components/WalkRowContextMenu.tsx +0 -105
  224. /package/dist/chunks/{ChooseSamplesDialog-BNVYBGJQ.js.map → ChooseSamplesDialog-JRN4YLBY.js.map} +0 -0
  225. /package/dist/chunks/{ChooseWalksDialog-EN3C6AQA.js.map → ChooseWalksDialog-UXLFZTZJ.js.map} +0 -0
  226. /package/dist/chunks/{GraphGenomeView-BG7UHKCT.js.map → GraphGenomeView-NJ5XCREB.js.map} +0 -0
  227. /package/dist/chunks/{HighlightColorDialog-YT55H62Y.js.map → HighlightColorDialog-UK5VDDV6.js.map} +0 -0
  228. /package/dist/chunks/{LinearGraphDisplay-SPILNH4J.js.map → LinearGraphDisplay-TXTKAACH.js.map} +0 -0
  229. /package/dist/chunks/{TrackWarningsDialog-QT4YEXGD.js.map → TrackWarningsDialog-TRYBOIHX.js.map} +0 -0
  230. /package/dist/chunks/{chunk-AHFQHYJI.js.map → chunk-IBWORBGY.js.map} +0 -0
@@ -0,0 +1,2 @@
1
+ import{b as j,c as O,d as z,e as J,g as D}from"./chunk-DVCZ3BOC.js";import"./chunk-5OI6LO4I.js";import"./chunk-7WIBGKBU.js";import"./chunk-OGQ6OQDN.js";import{d as R,e as v}from"./chunk-NDRLCGPC.js";import{b as G,h as S,j as A,m as w}from"./chunk-TENPGIPX.js";import{h as r}from"./chunk-YU7T5WOQ.js";import{A as _,B as W,C as H,F as b,a as x,f as M,g as p,o as B,s as I,t as P,u as N}from"./chunk-UGCNPRKC.js";import{b as t,c as g,d as E}from"./chunk-JMBDTNUR.js";var s=JBrowseExports["@jbrowse/core/util/tss-react"],$=o=>s!=null&&o in Object(s)?s[o]:Object.assign(function(){throw new Error("This JBrowse does not provide '"+o+"' from '@jbrowse/core/util/tss-react' to plugins: GraphGenomeView was built for a different JBrowse version, so update the plugin or JBrowse")},{jbrowseHostMissing:o});var V=$("makeStyles"),oe=s&&s.__esModule?s.default:s;var ee=V()(o=>({paper:{margin:o.spacing(),padding:o.spacing()},field:{marginTop:o.spacing(2)},submit:{marginTop:25,marginBottom:100,display:"block"}})),te=["RgfaTabixAdapter","GbzBaseSyntenyAdapter","MinigraphBubbleAdapter"],ae=E(function({model:a}){let{classes:l}=ee(),f=A(a),d=R(a),[n,X]=r("RgfaTabixAdapter"),[i,q]=r(),[y,k]=r(),[C,K]=r(),m=n==="GbzBaseSyntenyAdapter",[F,Q]=r(""),[h,U]=r("Pangenome graph"),[L,u]=r(),[Y,T]=r(!1);async function Z(){if(!i||!a.assembly)return;u(void 0),T(!0);let e=n==="RgfaTabixAdapter"?await j(i,S(a).pluginManager):void 0;if(T(!1),!(!v(a)||!a.assembly))try{let c=h.trim();G({model:a,session:f,conf:D({choice:n,loc:i,indexLoc:y,readsLoc:m?C:void 0,manifest:e,assembly:a.assembly,sample:F,trackId:w({name:c}),name:c})})}catch(c){u(c)}}return g(_,{className:l.paper,children:[L?t(M,{error:L}):null,g(I,{children:[t(N,{children:"File type"}),t(H,{value:n,onChange:e=>{X(e.target.value),k(void 0),u(void 0)},children:te.map(e=>t(P,{value:e,control:t(W,{}),label:O[e]},e))})]}),t(p,{location:i,name:z[n],rootModel:d,setLocation:q}),t(p,{location:y,name:J[n],rootModel:d,setLocation:k}),m?t(p,{location:C,name:"Path to reads aligned to the graph as GAF (optional; a .gz is read by the .tbi beside it)",rootModel:d,setLocation:K}):null,t(b,{className:l.field,value:F,onChange:e=>{Q(e.target.value)},label:"Sample name in the graph",slotProps:{htmlInput:{"data-testid":"graph-sample-input"}},helperText:m?"Optional. The graph's sample for this assembly, e.g. GRCh38 where the assembly is hg38; leave blank when they share a name or the graph has one reference sample":"Optional. The PanSN prefix the graph gives this assembly, e.g. GRCh38 for HPRC's GRCh38#0#chr1; leave blank when the graph's stable names are bare",placeholder:"GRCh38",fullWidth:!0}),t(b,{className:l.field,value:h,helperText:"Track name",slotProps:{htmlInput:{"data-testid":"graph-track-name-input"}},onChange:e=>{U(e.target.value)}}),t(x,{session:f,helperText:"Select assembly to add track to",selected:a.assembly,onChange:e=>{a.setAssembly(e)},fullWidth:!0}),t(B,{variant:"contained",className:l.submit,disabled:!i||!h.trim()||!a.assembly||Y,onClick:()=>{Z()},children:"Submit"})]})}),ue=ae;export{ue as default};
2
+ //# sourceMappingURL=AddTrackWorkflow-RPH4ZFSX.js.map
@@ -1,7 +1,7 @@
1
1
  {
2
2
  "version": 3,
3
3
  "sources": ["host:@jbrowse/core/util/tss-react", "../../src/GraphAddTrackWorkflow/AddTrackWorkflow.tsx"],
4
- "sourcesContent": ["const host = JBrowseExports[\"@jbrowse/core/util/tss-react\"]\nconst read = name =>\n host != null && name in Object(host)\n ? host[name]\n : Object.assign(\n function () {\n throw new Error(\n \"This JBrowse does not provide '\" + name + \"' from '\" + \"@jbrowse/core/util/tss-react\" + \"' to plugins: GraphGenomeView was built for a different JBrowse version, so update the plugin or JBrowse\",\n )\n },\n { jbrowseHostMissing: name },\n )\nexport const cx = /* @__PURE__ */ read(\"cx\")\nexport const keyframes = /* @__PURE__ */ read(\"keyframes\")\nexport const makeStyles = /* @__PURE__ */ read(\"makeStyles\")\nexport default host && host.__esModule ? host.default : host", "import { useState } from 'react'\n\nimport { AssemblySelector, ErrorMessage, FileSelector } from '@jbrowse/core/ui'\nimport {\n addTrackFromWidget,\n getEnv,\n getSession,\n makeTrackId,\n} from '@jbrowse/core/util'\nimport { makeStyles } from '@jbrowse/core/util/tss-react'\nimport { getRoot, isAlive } from '@jbrowse/mobx-state-tree'\nimport {\n Button,\n FormControl,\n FormControlLabel,\n FormLabel,\n Paper,\n Radio,\n RadioGroup,\n TextField,\n} from '@mui/material'\nimport { observer } from 'mobx-react'\n\nimport {\n GRAPH_FILE_FIELDS,\n GRAPH_FILE_LABELS,\n GRAPH_INDEX_FIELDS,\n buildTrackConfig,\n} from './buildTrackConfig'\nimport { readGraphManifest } from './graphManifest'\n\nimport type { GraphFileChoice } from './buildTrackConfig'\nimport type {\n AbstractRootModel,\n AddTrackWorkflowModel,\n FileLocation,\n} from '@jbrowse/core/util'\n\nconst useStyles = makeStyles()(theme => ({\n paper: {\n margin: theme.spacing(),\n padding: theme.spacing(),\n },\n field: {\n marginTop: theme.spacing(2),\n },\n submit: {\n marginTop: 25,\n marginBottom: 100,\n display: 'block',\n },\n}))\n\nconst CHOICES: GraphFileChoice[] = [\n 'RgfaTabixAdapter',\n 'GbzBaseSyntenyAdapter',\n 'MinigraphBubbleAdapter',\n]\n\nconst GraphAddTrackWidget = observer(function GraphAddTrackWidget({\n model,\n}: {\n model: AddTrackWorkflowModel\n}) {\n const { classes } = useStyles()\n const session = getSession(model)\n const rootModel = getRoot<AbstractRootModel>(model)\n const [choice, setChoice] = useState<GraphFileChoice>('RgfaTabixAdapter')\n const [loc, setLoc] = useState<FileLocation>()\n const [indexLoc, setIndexLoc] = useState<FileLocation>()\n const [readsLoc, setReadsLoc] = useState<FileLocation>()\n const gbz = choice === 'GbzBaseSyntenyAdapter'\n const [sample, setSample] = useState('')\n const [trackName, setTrackName] = useState('Pangenome graph')\n const [error, setError] = useState<unknown>()\n const [submitting, setSubmitting] = useState(false)\n\n async function handleSubmit() {\n if (!loc || !model.assembly) {\n return\n }\n setError(undefined)\n setSubmitting(true)\n const manifest =\n choice === 'RgfaTabixAdapter'\n ? await readGraphManifest(loc, getEnv(model).pluginManager)\n : undefined\n setSubmitting(false)\n if (!isAlive(model) || !model.assembly) {\n return\n }\n try {\n const name = trackName.trim()\n addTrackFromWidget({\n model,\n session,\n conf: buildTrackConfig({\n choice,\n loc,\n indexLoc,\n readsLoc: gbz ? readsLoc : undefined,\n manifest,\n assembly: model.assembly,\n sample,\n trackId: makeTrackId({ name }),\n name,\n }),\n })\n } catch (e) {\n setError(e)\n }\n }\n\n return (\n <Paper className={classes.paper}>\n {error ? <ErrorMessage error={error} /> : null}\n <FormControl>\n <FormLabel>File type</FormLabel>\n <RadioGroup\n value={choice}\n onChange={event => {\n setChoice(event.target.value as GraphFileChoice)\n setError(undefined)\n }}\n >\n {CHOICES.map(option => (\n <FormControlLabel\n key={option}\n value={option}\n control={<Radio />}\n label={GRAPH_FILE_LABELS[option]}\n />\n ))}\n </RadioGroup>\n </FormControl>\n <FileSelector\n location={loc}\n name={GRAPH_FILE_FIELDS[choice]}\n rootModel={rootModel}\n setLocation={setLoc}\n />\n <FileSelector\n location={indexLoc}\n name={GRAPH_INDEX_FIELDS[choice]}\n rootModel={rootModel}\n setLocation={setIndexLoc}\n />\n {gbz ? (\n <FileSelector\n location={readsLoc}\n name=\"Path to reads aligned to the graph as GAF (optional; a .gz is read by the .tbi beside it)\"\n rootModel={rootModel}\n setLocation={setReadsLoc}\n />\n ) : null}\n <TextField\n className={classes.field}\n value={sample}\n onChange={event => {\n setSample(event.target.value)\n }}\n label=\"Sample name in the graph\"\n slotProps={{ htmlInput: { 'data-testid': 'graph-sample-input' } }}\n helperText={\n gbz\n ? \"Optional. The graph's sample for this assembly, e.g. GRCh38 where the assembly is hg38; leave blank when they share a name or the graph has one reference sample\"\n : \"Optional. The PanSN prefix the graph gives this assembly, e.g. GRCh38 for HPRC's GRCh38#0#chr1; leave blank when the graph's stable names are bare\"\n }\n placeholder=\"GRCh38\"\n fullWidth\n />\n <TextField\n className={classes.field}\n value={trackName}\n helperText=\"Track name\"\n slotProps={{ htmlInput: { 'data-testid': 'graph-track-name-input' } }}\n onChange={event => {\n setTrackName(event.target.value)\n }}\n />\n <AssemblySelector\n session={session}\n helperText=\"Select assembly to add track to\"\n selected={model.assembly}\n onChange={arg => {\n model.setAssembly(arg)\n }}\n fullWidth\n />\n <Button\n variant=\"contained\"\n className={classes.submit}\n disabled={!loc || !trackName.trim() || !model.assembly || submitting}\n onClick={() => {\n void handleSubmit()\n }}\n >\n Submit\n </Button>\n </Paper>\n )\n})\n\nexport default GraphAddTrackWidget\n"],
5
- "mappings": "gdAAA,IAAMA,EAAO,eAAe,8BAA8B,EACpDC,EAAOC,GACXF,GAAQ,MAAQE,KAAQ,OAAOF,CAAI,EAC/BA,EAAKE,CAAI,EACT,OAAO,OACL,UAAY,CACV,MAAM,IAAI,MACR,kCAAoCA,EAAO,8IAC7C,CACF,EACA,CAAE,mBAAoBA,CAAK,CAC7B,EAGC,IAAMC,EAA6BC,EAAK,YAAY,EACpDC,GAAQC,GAAQA,EAAK,WAAaA,EAAK,QAAUA,ECuBxD,IAAMC,GAAYC,EAAW,EAAEC,IAAU,CACvC,MAAO,CACL,OAAQA,EAAM,QAAQ,EACtB,QAASA,EAAM,QAAQ,CACzB,EACA,MAAO,CACL,UAAWA,EAAM,QAAQ,CAAC,CAC5B,EACA,OAAQ,CACN,UAAW,GACX,aAAc,IACd,QAAS,OACX,CACF,EAAE,EAEIC,GAA6B,CACjC,mBACA,wBACA,wBACF,EAEMC,GAAsBC,EAAS,SAA6B,CAChE,MAAAC,CACF,EAEG,CACD,GAAM,CAAE,QAAAC,CAAQ,EAAIP,GAAU,EACxBQ,EAAUC,EAAWH,CAAK,EAC1BI,EAAYC,EAA2BL,CAAK,EAC5C,CAACM,EAAQC,CAAS,EAAIC,EAA0B,kBAAkB,EAClE,CAACC,EAAKC,CAAM,EAAIF,EAAuB,EACvC,CAACG,EAAUC,CAAW,EAAIJ,EAAuB,EACjD,CAACK,EAAUC,CAAW,EAAIN,EAAuB,EACjDO,EAAMT,IAAW,wBACjB,CAACU,EAAQC,CAAS,EAAIT,EAAS,EAAE,EACjC,CAACU,EAAWC,CAAY,EAAIX,EAAS,iBAAiB,EACtD,CAACY,EAAOC,CAAQ,EAAIb,EAAkB,EACtC,CAACc,EAAYC,CAAa,EAAIf,EAAS,EAAK,EAElD,eAAegB,GAAe,CAC5B,GAAI,CAACf,GAAO,CAACT,EAAM,SACjB,OAEFqB,EAAS,MAAS,EAClBE,EAAc,EAAI,EAClB,IAAME,EACJnB,IAAW,mBACP,MAAMoB,EAAkBjB,EAAKkB,EAAO3B,CAAK,EAAE,aAAa,EACxD,OAEN,GADAuB,EAAc,EAAK,EACf,GAACK,EAAQ5B,CAAK,GAAK,CAACA,EAAM,UAG9B,GAAI,CACF,IAAM6B,EAAOX,EAAU,KAAK,EAC5BY,EAAmB,CACjB,MAAA9B,EACA,QAAAE,EACA,KAAM6B,EAAiB,CACrB,OAAAzB,EACA,IAAAG,EACA,SAAAE,EACA,SAAUI,EAAMF,EAAW,OAC3B,SAAAY,EACA,SAAUzB,EAAM,SAChB,OAAAgB,EACA,QAASgB,EAAY,CAAE,KAAAH,CAAK,CAAC,EAC7B,KAAAA,CACF,CAAC,CACH,CAAC,CACH,OAASI,EAAG,CACVZ,EAASY,CAAC,CACZ,CACF,CAEA,OACEC,EAACC,EAAA,CAAM,UAAWlC,EAAQ,MACvB,UAAAmB,EAAQgB,EAACC,EAAA,CAAa,MAAOjB,EAAO,EAAK,KAC1Cc,EAACI,EAAA,CACC,UAAAF,EAACG,EAAA,CAAU,qBAAS,EACpBH,EAACI,EAAA,CACC,MAAOlC,EACP,SAAUmC,GAAS,CACjBlC,EAAUkC,EAAM,OAAO,KAAwB,EAC/CpB,EAAS,MAAS,CACpB,EAEC,SAAAxB,GAAQ,IAAI6C,GACXN,EAACO,EAAA,CAEC,MAAOD,EACP,QAASN,EAACQ,EAAA,EAAM,EAChB,MAAOC,EAAkBH,CAAM,GAH1BA,CAIP,CACD,EACH,GACF,EACAN,EAACU,EAAA,CACC,SAAUrC,EACV,KAAMsC,EAAkBzC,CAAM,EAC9B,UAAWF,EACX,YAAaM,EACf,EACA0B,EAACU,EAAA,CACC,SAAUnC,EACV,KAAMqC,EAAmB1C,CAAM,EAC/B,UAAWF,EACX,YAAaQ,EACf,EACCG,EACCqB,EAACU,EAAA,CACC,SAAUjC,EACV,KAAK,4FACL,UAAWT,EACX,YAAaU,EACf,EACE,KACJsB,EAACa,EAAA,CACC,UAAWhD,EAAQ,MACnB,MAAOe,EACP,SAAUyB,GAAS,CACjBxB,EAAUwB,EAAM,OAAO,KAAK,CAC9B,EACA,MAAM,2BACN,UAAW,CAAE,UAAW,CAAE,cAAe,oBAAqB,CAAE,EAChE,WACE1B,EACI,mKACA,qJAEN,YAAY,SACZ,UAAS,GACX,EACAqB,EAACa,EAAA,CACC,UAAWhD,EAAQ,MACnB,MAAOiB,EACP,WAAW,aACX,UAAW,CAAE,UAAW,CAAE,cAAe,wBAAyB,CAAE,EACpE,SAAUuB,GAAS,CACjBtB,EAAasB,EAAM,OAAO,KAAK,CACjC,EACF,EACAL,EAACc,EAAA,CACC,QAAShD,EACT,WAAW,kCACX,SAAUF,EAAM,SAChB,SAAUmD,GAAO,CACfnD,EAAM,YAAYmD,CAAG,CACvB,EACA,UAAS,GACX,EACAf,EAACgB,EAAA,CACC,QAAQ,YACR,UAAWnD,EAAQ,OACnB,SAAU,CAACQ,GAAO,CAACS,EAAU,KAAK,GAAK,CAAClB,EAAM,UAAYsB,EAC1D,QAAS,IAAM,CACRE,EAAa,CACpB,EACD,kBAED,GACF,CAEJ,CAAC,EAEM6B,GAAQvD",
4
+ "sourcesContent": ["const host = JBrowseExports[\"@jbrowse/core/util/tss-react\"]\nconst read = name =>\n host != null && name in Object(host)\n ? host[name]\n : Object.assign(\n function () {\n throw new Error(\n \"This JBrowse does not provide '\" + name + \"' from '\" + \"@jbrowse/core/util/tss-react\" + \"' to plugins: GraphGenomeView was built for a different JBrowse version, so update the plugin or JBrowse\",\n )\n },\n { jbrowseHostMissing: name },\n )\nexport const cx = /* @__PURE__ */ read(\"cx\")\nexport const keyframes = /* @__PURE__ */ read(\"keyframes\")\nexport const makeStyles = /* @__PURE__ */ read(\"makeStyles\")\nexport default host && host.__esModule ? host.default : host", "import { useState } from 'react'\n\nimport { AssemblySelector, ErrorMessage, FileSelector } from '@jbrowse/core/ui'\nimport {\n addTrackFromWidget,\n getEnv,\n getSession,\n makeTrackId,\n} from '@jbrowse/core/util'\nimport { makeStyles } from '@jbrowse/core/util/tss-react'\nimport { getRoot, isAlive } from '@jbrowse/mobx-state-tree'\nimport {\n Button,\n FormControl,\n FormControlLabel,\n FormLabel,\n Paper,\n Radio,\n RadioGroup,\n TextField,\n} from '@mui/material'\nimport { observer } from 'mobx-react'\n\nimport {\n GRAPH_FILE_FIELDS,\n GRAPH_FILE_LABELS,\n GRAPH_INDEX_FIELDS,\n buildTrackConfig,\n} from './buildTrackConfig'\nimport { readGraphManifest } from './graphManifest'\n\nimport type { GraphFileChoice } from './buildTrackConfig'\nimport type {\n AbstractRootModel,\n AddTrackWorkflowModel,\n FileLocation,\n} from '@jbrowse/core/util'\n\nconst useStyles = makeStyles()(theme => ({\n paper: {\n margin: theme.spacing(),\n padding: theme.spacing(),\n },\n field: {\n marginTop: theme.spacing(2),\n },\n submit: {\n marginTop: 25,\n marginBottom: 100,\n display: 'block',\n },\n}))\n\nconst CHOICES: GraphFileChoice[] = [\n 'RgfaTabixAdapter',\n 'GbzBaseSyntenyAdapter',\n 'MinigraphBubbleAdapter',\n]\n\nconst GraphAddTrackWidget = observer(function GraphAddTrackWidget({\n model,\n}: {\n model: AddTrackWorkflowModel\n}) {\n const { classes } = useStyles()\n const session = getSession(model)\n const rootModel = getRoot<AbstractRootModel>(model)\n const [choice, setChoice] = useState<GraphFileChoice>('RgfaTabixAdapter')\n const [loc, setLoc] = useState<FileLocation>()\n const [indexLoc, setIndexLoc] = useState<FileLocation>()\n const [readsLoc, setReadsLoc] = useState<FileLocation>()\n const gbz = choice === 'GbzBaseSyntenyAdapter'\n const [sample, setSample] = useState('')\n const [trackName, setTrackName] = useState('Pangenome graph')\n const [error, setError] = useState<unknown>()\n const [submitting, setSubmitting] = useState(false)\n\n async function handleSubmit() {\n if (!loc || !model.assembly) {\n return\n }\n setError(undefined)\n setSubmitting(true)\n const manifest =\n choice === 'RgfaTabixAdapter'\n ? await readGraphManifest(loc, getEnv(model).pluginManager)\n : undefined\n setSubmitting(false)\n if (!isAlive(model) || !model.assembly) {\n return\n }\n try {\n const name = trackName.trim()\n addTrackFromWidget({\n model,\n session,\n conf: buildTrackConfig({\n choice,\n loc,\n indexLoc,\n readsLoc: gbz ? readsLoc : undefined,\n manifest,\n assembly: model.assembly,\n sample,\n trackId: makeTrackId({ name }),\n name,\n }),\n })\n } catch (e) {\n setError(e)\n }\n }\n\n return (\n <Paper className={classes.paper}>\n {error ? <ErrorMessage error={error} /> : null}\n <FormControl>\n <FormLabel>File type</FormLabel>\n <RadioGroup\n value={choice}\n onChange={event => {\n setChoice(event.target.value as GraphFileChoice)\n setIndexLoc(undefined)\n setError(undefined)\n }}\n >\n {CHOICES.map(option => (\n <FormControlLabel\n key={option}\n value={option}\n control={<Radio />}\n label={GRAPH_FILE_LABELS[option]}\n />\n ))}\n </RadioGroup>\n </FormControl>\n <FileSelector\n location={loc}\n name={GRAPH_FILE_FIELDS[choice]}\n rootModel={rootModel}\n setLocation={setLoc}\n />\n <FileSelector\n location={indexLoc}\n name={GRAPH_INDEX_FIELDS[choice]}\n rootModel={rootModel}\n setLocation={setIndexLoc}\n />\n {gbz ? (\n <FileSelector\n location={readsLoc}\n name=\"Path to reads aligned to the graph as GAF (optional; a .gz is read by the .tbi beside it)\"\n rootModel={rootModel}\n setLocation={setReadsLoc}\n />\n ) : null}\n <TextField\n className={classes.field}\n value={sample}\n onChange={event => {\n setSample(event.target.value)\n }}\n label=\"Sample name in the graph\"\n slotProps={{ htmlInput: { 'data-testid': 'graph-sample-input' } }}\n helperText={\n gbz\n ? \"Optional. The graph's sample for this assembly, e.g. GRCh38 where the assembly is hg38; leave blank when they share a name or the graph has one reference sample\"\n : \"Optional. The PanSN prefix the graph gives this assembly, e.g. GRCh38 for HPRC's GRCh38#0#chr1; leave blank when the graph's stable names are bare\"\n }\n placeholder=\"GRCh38\"\n fullWidth\n />\n <TextField\n className={classes.field}\n value={trackName}\n helperText=\"Track name\"\n slotProps={{ htmlInput: { 'data-testid': 'graph-track-name-input' } }}\n onChange={event => {\n setTrackName(event.target.value)\n }}\n />\n <AssemblySelector\n session={session}\n helperText=\"Select assembly to add track to\"\n selected={model.assembly}\n onChange={arg => {\n model.setAssembly(arg)\n }}\n fullWidth\n />\n <Button\n variant=\"contained\"\n className={classes.submit}\n disabled={!loc || !trackName.trim() || !model.assembly || submitting}\n onClick={() => {\n void handleSubmit()\n }}\n >\n Submit\n </Button>\n </Paper>\n )\n})\n\nexport default GraphAddTrackWidget\n"],
5
+ "mappings": "gdAAA,IAAMA,EAAO,eAAe,8BAA8B,EACpDC,EAAOC,GACXF,GAAQ,MAAQE,KAAQ,OAAOF,CAAI,EAC/BA,EAAKE,CAAI,EACT,OAAO,OACL,UAAY,CACV,MAAM,IAAI,MACR,kCAAoCA,EAAO,8IAC7C,CACF,EACA,CAAE,mBAAoBA,CAAK,CAC7B,EAGC,IAAMC,EAA6BC,EAAK,YAAY,EACpDC,GAAQC,GAAQA,EAAK,WAAaA,EAAK,QAAUA,ECuBxD,IAAMC,GAAYC,EAAW,EAAEC,IAAU,CACvC,MAAO,CACL,OAAQA,EAAM,QAAQ,EACtB,QAASA,EAAM,QAAQ,CACzB,EACA,MAAO,CACL,UAAWA,EAAM,QAAQ,CAAC,CAC5B,EACA,OAAQ,CACN,UAAW,GACX,aAAc,IACd,QAAS,OACX,CACF,EAAE,EAEIC,GAA6B,CACjC,mBACA,wBACA,wBACF,EAEMC,GAAsBC,EAAS,SAA6B,CAChE,MAAAC,CACF,EAEG,CACD,GAAM,CAAE,QAAAC,CAAQ,EAAIP,GAAU,EACxBQ,EAAUC,EAAWH,CAAK,EAC1BI,EAAYC,EAA2BL,CAAK,EAC5C,CAACM,EAAQC,CAAS,EAAIC,EAA0B,kBAAkB,EAClE,CAACC,EAAKC,CAAM,EAAIF,EAAuB,EACvC,CAACG,EAAUC,CAAW,EAAIJ,EAAuB,EACjD,CAACK,EAAUC,CAAW,EAAIN,EAAuB,EACjDO,EAAMT,IAAW,wBACjB,CAACU,EAAQC,CAAS,EAAIT,EAAS,EAAE,EACjC,CAACU,EAAWC,CAAY,EAAIX,EAAS,iBAAiB,EACtD,CAACY,EAAOC,CAAQ,EAAIb,EAAkB,EACtC,CAACc,EAAYC,CAAa,EAAIf,EAAS,EAAK,EAElD,eAAegB,GAAe,CAC5B,GAAI,CAACf,GAAO,CAACT,EAAM,SACjB,OAEFqB,EAAS,MAAS,EAClBE,EAAc,EAAI,EAClB,IAAME,EACJnB,IAAW,mBACP,MAAMoB,EAAkBjB,EAAKkB,EAAO3B,CAAK,EAAE,aAAa,EACxD,OAEN,GADAuB,EAAc,EAAK,EACf,GAACK,EAAQ5B,CAAK,GAAK,CAACA,EAAM,UAG9B,GAAI,CACF,IAAM6B,EAAOX,EAAU,KAAK,EAC5BY,EAAmB,CACjB,MAAA9B,EACA,QAAAE,EACA,KAAM6B,EAAiB,CACrB,OAAAzB,EACA,IAAAG,EACA,SAAAE,EACA,SAAUI,EAAMF,EAAW,OAC3B,SAAAY,EACA,SAAUzB,EAAM,SAChB,OAAAgB,EACA,QAASgB,EAAY,CAAE,KAAAH,CAAK,CAAC,EAC7B,KAAAA,CACF,CAAC,CACH,CAAC,CACH,OAASI,EAAG,CACVZ,EAASY,CAAC,CACZ,CACF,CAEA,OACEC,EAACC,EAAA,CAAM,UAAWlC,EAAQ,MACvB,UAAAmB,EAAQgB,EAACC,EAAA,CAAa,MAAOjB,EAAO,EAAK,KAC1Cc,EAACI,EAAA,CACC,UAAAF,EAACG,EAAA,CAAU,qBAAS,EACpBH,EAACI,EAAA,CACC,MAAOlC,EACP,SAAUmC,GAAS,CACjBlC,EAAUkC,EAAM,OAAO,KAAwB,EAC/C7B,EAAY,MAAS,EACrBS,EAAS,MAAS,CACpB,EAEC,SAAAxB,GAAQ,IAAI6C,GACXN,EAACO,EAAA,CAEC,MAAOD,EACP,QAASN,EAACQ,EAAA,EAAM,EAChB,MAAOC,EAAkBH,CAAM,GAH1BA,CAIP,CACD,EACH,GACF,EACAN,EAACU,EAAA,CACC,SAAUrC,EACV,KAAMsC,EAAkBzC,CAAM,EAC9B,UAAWF,EACX,YAAaM,EACf,EACA0B,EAACU,EAAA,CACC,SAAUnC,EACV,KAAMqC,EAAmB1C,CAAM,EAC/B,UAAWF,EACX,YAAaQ,EACf,EACCG,EACCqB,EAACU,EAAA,CACC,SAAUjC,EACV,KAAK,4FACL,UAAWT,EACX,YAAaU,EACf,EACE,KACJsB,EAACa,EAAA,CACC,UAAWhD,EAAQ,MACnB,MAAOe,EACP,SAAUyB,GAAS,CACjBxB,EAAUwB,EAAM,OAAO,KAAK,CAC9B,EACA,MAAM,2BACN,UAAW,CAAE,UAAW,CAAE,cAAe,oBAAqB,CAAE,EAChE,WACE1B,EACI,mKACA,qJAEN,YAAY,SACZ,UAAS,GACX,EACAqB,EAACa,EAAA,CACC,UAAWhD,EAAQ,MACnB,MAAOiB,EACP,WAAW,aACX,UAAW,CAAE,UAAW,CAAE,cAAe,wBAAyB,CAAE,EACpE,SAAUuB,GAAS,CACjBtB,EAAasB,EAAM,OAAO,KAAK,CACjC,EACF,EACAL,EAACc,EAAA,CACC,QAAShD,EACT,WAAW,kCACX,SAAUF,EAAM,SAChB,SAAUmD,GAAO,CACfnD,EAAM,YAAYmD,CAAG,CACvB,EACA,UAAS,GACX,EACAf,EAACgB,EAAA,CACC,QAAQ,YACR,UAAWnD,EAAQ,OACnB,SAAU,CAACQ,GAAO,CAACS,EAAU,KAAK,GAAK,CAAClB,EAAM,UAAYsB,EAC1D,QAAS,IAAM,CACRE,EAAa,CACpB,EACD,kBAED,GACF,CAEJ,CAAC,EAEM6B,GAAQvD",
6
6
  "names": ["host", "read", "name", "makeStyles", "read", "tss_react_default", "host", "useStyles", "makeStyles", "theme", "CHOICES", "GraphAddTrackWidget", "observer", "model", "classes", "session", "getSession", "rootModel", "getRoot", "choice", "setChoice", "useState", "loc", "setLoc", "indexLoc", "setIndexLoc", "readsLoc", "setReadsLoc", "gbz", "sample", "setSample", "trackName", "setTrackName", "error", "setError", "submitting", "setSubmitting", "handleSubmit", "manifest", "readGraphManifest", "getEnv", "isAlive", "name", "addTrackFromWidget", "buildTrackConfig", "makeTrackId", "e", "jsxs", "Paper", "jsx", "ErrorMessage", "FormControl", "FormLabel", "RadioGroup", "event", "option", "FormControlLabel", "Radio", "GRAPH_FILE_LABELS", "FileSelector", "GRAPH_FILE_FIELDS", "GRAPH_INDEX_FIELDS", "TextField", "AssemblySelector", "arg", "Button", "AddTrackWorkflow_default"]
7
7
  }
@@ -1,2 +1,2 @@
1
- import{a as i}from"./chunk-AHFQHYJI.js";import"./chunk-6OPB7FOJ.js";import"./chunk-E3WYW7Z6.js";import{b as a,d as s}from"./chunk-JMBDTNUR.js";var t=s(function({model:o,onClose:p}){let{samples:l}=o.walkRowSampleChoices;return a(i,{title:"Choose samples",label:"Samples",testId:"graph-choose-samples",caption:`${l.length.toLocaleString()} samples in this cut. Their walks show in pairs, in the order picked.`,options:l,initial:o.walkRowSamples??[],confirmLabel:e=>e.length>0?"Show these samples":"Show every sample",onPick:e=>{o.setWalkRowSamples(e.length>0?e:void 0)},onClose:p})}),h=t;export{h as default};
2
- //# sourceMappingURL=ChooseSamplesDialog-BNVYBGJQ.js.map
1
+ import{a as i}from"./chunk-IBWORBGY.js";import"./chunk-YU7T5WOQ.js";import"./chunk-UGCNPRKC.js";import{b as a,d as s}from"./chunk-JMBDTNUR.js";var t=s(function({model:o,onClose:p}){let{samples:l}=o.walkRowSampleChoices;return a(i,{title:"Choose samples",label:"Samples",testId:"graph-choose-samples",caption:`${l.length.toLocaleString()} samples in this cut. Their walks show in pairs, in the order picked.`,options:l,initial:o.walkRowSamples??[],confirmLabel:e=>e.length>0?"Show these samples":"Show every sample",onPick:e=>{o.setWalkRowSamples(e.length>0?e:void 0)},onClose:p})}),h=t;export{h as default};
2
+ //# sourceMappingURL=ChooseSamplesDialog-JRN4YLBY.js.map
@@ -1,2 +1,2 @@
1
- import{a as t}from"./chunk-AHFQHYJI.js";import"./chunk-6OPB7FOJ.js";import"./chunk-E3WYW7Z6.js";import{b as l,d as i}from"./chunk-JMBDTNUR.js";var h=i(function({model:o,onClose:s}){let{walkChoices:a}=o;return l(t,{title:"Choose haplotypes",label:"Haplotypes",testId:"graph-choose-walks",caption:`${a.length.toLocaleString()} haplotypes in this cut. The ones picked are highlighted in this order, and the rest fades.`,options:a.map(e=>e.name),initial:o.walkLayers.map(e=>e.walk),optionLabel:e=>o.walkLabel(e),confirmLabel:e=>e.length>0?"Highlight these":"Clear highlights",onPick:e=>{o.liftWalks(e)},onClose:s})}),k=h;export{k as default};
2
- //# sourceMappingURL=ChooseWalksDialog-EN3C6AQA.js.map
1
+ import{a as t}from"./chunk-IBWORBGY.js";import"./chunk-YU7T5WOQ.js";import"./chunk-UGCNPRKC.js";import{b as l,d as i}from"./chunk-JMBDTNUR.js";var h=i(function({model:o,onClose:s}){let{walkChoices:a}=o;return l(t,{title:"Choose haplotypes",label:"Haplotypes",testId:"graph-choose-walks",caption:`${a.length.toLocaleString()} haplotypes in this cut. The ones picked are highlighted in this order, and the rest fades.`,options:a.map(e=>e.name),initial:o.walkLayers.map(e=>e.walk),optionLabel:e=>o.walkLabel(e),confirmLabel:e=>e.length>0?"Highlight these":"Clear highlights",onPick:e=>{o.liftWalks(e)},onClose:s})}),k=h;export{k as default};
2
+ //# sourceMappingURL=ChooseWalksDialog-UXLFZTZJ.js.map
@@ -0,0 +1,2 @@
1
+ import{a as ae,b as re,c as de,d as oe,e as se,f as q,g as ie}from"./chunk-EGDLEQ47.js";import{a as te}from"./chunk-KOB7YIMR.js";import{a as W,b as Z,c as U,d as u,e as M,f as $,g as k,h as E,i as z,k as X}from"./chunk-S75W27NT.js";import"./chunk-5O6YRJHL.js";import{a as O,b as H}from"./chunk-N6PQZSAX.js";import"./chunk-TUKJRVBS.js";import"./chunk-T4HOC6O3.js";import{a as ne}from"./chunk-FVK2THE4.js";import{a as V,c as L,m as ee}from"./chunk-FHZXW7IG.js";import"./chunk-YBLNN5FU.js";import"./chunk-MVVXRZXV.js";import{e as F,f as P,g as C}from"./chunk-ZMOMFU2B.js";import"./chunk-E7J6SB5P.js";import{b as G,e as Y,f as J,g as K}from"./chunk-5OI6LO4I.js";import"./chunk-7WIBGKBU.js";import{a as N}from"./chunk-OGQ6OQDN.js";import"./chunk-NDRLCGPC.js";import{e as j,o as I}from"./chunk-TENPGIPX.js";function A(n){let s=/\.gbz\.db$/i;return s.test(n)?n.replace(s,".haplotype-index.db"):void 0}function he(n){return J(n)&&A(Y(n))?K(n,s=>A(s)):void 0}function ue(n){return n instanceof Error&&(/\bHTTP 40[34]\b/.test(n.message)||n.code==="ENOENT")}async function ce(n,s){let t=he(n);if(!t)return{};try{return await N(t,s).stat(),{location:t}}catch(e){return ue(e)?{}:{unreadable:String(e)}}}var ge=5e3,ye=50*2**20,T=class extends Error{constructor(t){super(`this GAF is ${(t/2**20).toFixed(0)} MB without an index; bgzip it and index it with tabix -p gaf`);this.name="UnindexedGafTooLargeError"}};function be(n,s){if(n.length<=s)return n;let t=n.length/s;return Array.from({length:s},(e,a)=>n[Math.floor(a*t)])}var pe=(n,s)=>n.path.some(t=>s.has(t.name)),v=class{constructor(s,t){this.file=s;this.tabix=t;this.whole=L({setup:async({signal:s})=>{let{size:t}=await this.file.stat();if(t>ye)throw new T(t);let e=await this.file.readFile({signal:s}),a=new TextDecoder().decode(e[0]===31&&e[1]===139?await V(e):e);return H(a)}})}async readIndexed(s,t,e){let a=1/0,r=-1/0;for(let i of t){let c=Number(i);Number.isInteger(c)&&(a=Math.min(a,c),r=Math.max(r,c))}let o=[];return a<=r&&await s.getLines("{node}",a,r+1,{signal:e,lineCallback:i=>{let c=O(i);c&&pe(c,t)&&o.push(c)}}),o}async readsOver(s,{signal:t}={}){let e=this.tabix?await this.readIndexed(this.tabix,s,t):(await this.whole({signal:t})).filter(a=>pe(a,s));return{records:be(e,ge),total:e.length}}};var B=class extends Error{constructor(t){super(t===void 0?"no haplotype index is set, so the walks cannot be named; build one with gbz-haplotype-index (cargo install gbz-haplotype-index) and set it as haplotypeIndexLocation":`the haplotype index beside the graph database could not be checked (${t}), so the walks cannot be named; if it is there, reload to try again, or set haplotypeIndexLocation`);this.name="NoHaplotypeIndexError"}},Q=class extends Error{constructor(t,e){super(`the graph is cut on its reference, ${e}; a window on ${t} has no reference coordinates to cut at. Open the graph from a ${e} view, and find this haplotype's lane there`);this.name="HaplotypeWindowError"}};function we(n,s){return`${u(n)}#${n.contig}@${s.start.node}.${s.start.offset}`}function xe({alignment:n,assemblyName:s,refName:t,lane:e}){let{refStart:a,refEnd:r,strand:o,cigar:i}=n;if(!(r<=a)){let{name:c,hapStart:p,hapEnd:l}=n,d=we(c,n),g={uniqueId:d,assemblyName:s,refName:t,start:a,end:r,type:"match",strand:o==="-"?-1:1,CIGAR:i,syntenyId:d,mate:{refName:c.contig,start:p,end:l,assemblyName:e}};return new ae(g)}}function Se({pair:n,lane:s,mateLane:t}){return ie({target:{name:`${u(n.target)}#${n.target.contig}`,assemblyName:s,refName:n.target.contig,start:n.targetStart,end:n.targetEnd},query:{name:`${u(n.query)}#${n.query.contig}`,assemblyName:t,refName:n.query.contig,start:n.queryStart,end:n.queryEnd},strand:n.strand,cigar:n.cigar,matches:n.matches,columns:n.columns})}var Pe=12e4;function Ne(n,s){let t=new Map;for(let a of n)for(let{name:r,start:o,end:i}of a.walkSpans()){let c=u(r);if(s.has(c)){let p=`${c}#${r.contig}`,l=t.get(p);l?l.push([o,i]):t.set(p,[[o,i]])}}let e=0;for(let a of t.values()){a.sort((r,o)=>r[0]-o[0]);for(let r=1;r<a.length;r++){let o=a[r][0]-a[r-1][1];o>0&&o<=Pe&&(e=Math.max(e,o))}}return e}function _(n,s){return n instanceof U&&(n.message=`Zoom in to about ${j(n.fitsBp)} to see ${s}`),n}var le=n=>({...n,start:Math.floor(n.start),end:Math.ceil(n.end)}),D=class extends re{constructor(){super(...arguments);this.graph=L({label:"Opening pangenome database",setup:async()=>{let t=this.getConf("gbzDbLocation"),e=this.getConf("haplotypeIndexLocation"),a=G(e)?{location:e}:await ce(t,this.pluginManager),r=a.location,o=await Z.open({source:N(t,this.pluginManager),haplotypeIndex:r?N(r,this.pluginManager):void 0}),i=this.getConf("assemblyNames")[0];if(i===void 0)throw new Error("GbzBaseSyntenyAdapter needs assemblyNames: its first entry is the assembly the reference sample is loaded as");let c=await $(o),p=k({configured:this.getConf("referenceSample"),anchorPrefix:P(this,i),referenceSamples:c});return{db:o,anchor:i,referenceSample:p,referenceSamples:c,unreadableIndex:a.unreadable}}});this.gaf=L({label:"Opening reads",setup:async()=>{let t=this.getConf("readsLocation"),e=this.getConf(["readsIndex","location"]);return G(t)?new v(N(t,this.pluginManager),G(e)?ee(this,["readsLocation"],["readsIndex"]):void 0):void 0}})}async referenceQuery(t,e){let{db:a,referenceSample:r}=await this.graph(e);return E(a,r,t)}async getHaplotypes(t={}){let{db:e,referenceSamples:a}=await this.graph(t),r=new Map;for(let o of await e.paths()){let i=u(o.name),c=r.get(i);c?c.contigs.includes(o.name.contig)||c.contigs.push(o.name.contig):r.set(i,{prefix:i,sample:o.name.sample,haplotype:o.name.haplotype,contigs:[o.name.contig],isReference:a.includes(o.name.sample)})}return[...r.values()]}async getHeader(t={}){let{anchor:e,referenceSample:a}=await this.graph(t),r=C(this),o=[];for(let i of await this.getHaplotypes(t))i.sample!==a&&o.push({name:q(r,i.prefix),label:i.prefix,group:i.sample});return{hasCoarseTier:!1,anchorAssemblyName:e,lanes:o}}async getRefNames(t={}){let{db:e,anchor:a,referenceSample:r}=await this.graph(t),{assemblyName:o}=t,i=o===a?void 0:P(this,o),c=(await e.paths()).filter(p=>o===a?p.isIndexed&&p.name.sample===r:F(u(p.name),i)).map(p=>p.name.contig);return[...new Set(c)]}keepPredicate(t,e){let a=t===void 0||t.length===0?void 0:t.map(r=>P(this,r));return a===void 0&&e===void 0?void 0:r=>M(r,a)&&(e===void 0||F(u(r),e))}async getSubgraph(t,e={}){let{signal:a}=e,{db:r,anchor:o,referenceSample:i,unreadableIndex:c}=await this.graph({signal:a}),{assemblyName:p,refName:l,start:d,end:g}=le(t);if(p!==o)throw new Q(p,o);let y=await this.referenceQuery(l,{signal:a});if(!y)throw new Error(`${i} has no indexed path ${l} in this graph`);let w=this.keepPredicate(e.haplotypes);if(w!==void 0&&!r.hasHaplotypeIndex)throw new B(c);let f=await X(r,y,d,g,{context:this.getConf("context"),snarls:e.snarls??this.getConf("subgraphSnarls"),limit:this.getConf("nodeLimit"),signal:e.signal,keep:w}).catch(x=>{throw _(x,"the graph")});return typeof f=="string"?f:ne(f,{start:d,end:g},e.walkRows)}async getReads(t,e={}){let a=await this.gaf(e);return a?a.readsOver(new Set(t),e):{records:[],total:0}}async laneHaplotypes(t,e){let a=P(this,t);return(await this.getHaplotypes(e)).filter(r=>F(r.prefix,a))}async anchorFeatures(t,e){let{db:a}=await this.graph(e),{assemblyName:r,refName:o,start:i,end:c}=t,p=P(this,e.targetAssemblyName),l=C(this),d=this.keepPredicate(e.haplotypes,p),g=await this.referenceQuery(o,e),y=this.getConf("nodeLimit");return(g?await I(`Reading graph ${o}:${i.toLocaleString()}-${c.toLocaleString()}`,e.statusCallback,()=>a.getAlignments({path:g,start:i,end:c,context:this.getConf("context"),haplotypes:"all",limit:y,signal:e.signal,keep:d}).catch(f=>{throw _(z(f,y,c-i)??f,"lanes")})):[]).flatMap(f=>{let x=f.resolved?xe({alignment:f,assemblyName:r,refName:o,lane:q(l,u(f.name))}):void 0;return x===void 0?[]:[x]})}async pairFeatures(t,e,a){let{db:r}=await this.graph(a),{refName:o,start:i,end:c}=t,l=(await Promise.all(e.map(async m=>({pair:m,featureSide:await this.laneHaplotypes(m.queryAssemblyName,a),mateSide:await this.laneHaplotypes(m.targetAssemblyName,a)})))).filter(({featureSide:m,mateSide:h})=>m.length>0&&h.length>0),d=new Set(l.flatMap(({featureSide:m,mateSide:h})=>[...m,...h].map(u))),g=await this.referenceQuery(o,a),y=this.getConf("nodeLimit"),w=this.getConf("context"),f=(m,h)=>r.getSubgraphs({path:m,start:i,end:c,context:h,snarls:this.getConf("subgraphSnarls"),haplotypes:"all",limit:y,signal:a.signal,keep:S=>d.has(u(S))}),x=async m=>{let h=await f(m,w).catch(b=>{throw _(z(b,y,c-i)??b,"lanes")}),S=Ne(h,d),R=Math.ceil(S/2)+2e3;return S>0&&R>w?f(m,R).catch(b=>{if(b instanceof W)return h;throw b}):h};return(g&&l.length>0?await I(se(e),a.statusCallback,()=>x(g)):[]).flatMap(m=>l.flatMap(({pair:h,featureSide:S,mateSide:R})=>S.flatMap(b=>R.flatMap(me=>m.pairAlignments({target:b,query:me,bases:!1}).map(fe=>Se({pair:fe,lane:h.queryAssemblyName,mateLane:h.targetAssemblyName}))))))}getFeatures(t,e={}){let a=le(t);return te(async r=>{let{db:o,anchor:i,unreadableIndex:c}=await this.graph(e);if(!o.hasHaplotypeIndex)throw new B(c);if(a.assemblyName===i){let p=oe(e),l=p.length>0?await this.pairFeatures(a,p,e):await this.anchorFeatures(a,e);for(let d of l)r.next(d)}r.complete()},e.signal)}};export{B as NoHaplotypeIndexError,de as PairTargetError,D as default,Ne as widestWalkGap};
2
+ //# sourceMappingURL=GbzBaseSyntenyAdapter-LRECKL74.js.map
@@ -0,0 +1,7 @@
1
+ {
2
+ "version": 3,
3
+ "sources": ["../../node_modules/.pnpm/@jbrowse+bandage-core@10.2.6/node_modules/@jbrowse/bandage-core/dist/chunks/chunk-VQWHMWLF.js", "../../src/GbzBaseSyntenyAdapter/companion.ts", "../../src/gaf/gafFile.ts", "../../src/GbzBaseSyntenyAdapter/GbzBaseSyntenyAdapter.ts"],
4
+ "sourcesContent": ["import {\n GBZ_CUT_DEFAULTS,\n cutWindowGFA,\n haplotypeWanted,\n referencePathQuery,\n referenceSamplesOf,\n resolveReferenceSample\n} from \"./chunk-Q7LYHGGQ.js\";\n\n// src/gbzCut.ts\nimport { GBZBase } from \"@gmod/gbz-base\";\nvar HPRC_GBZ = {\n db: \"https://s3-us-west-2.amazonaws.com/human-pangenomics/pangenomes/freeze/release2/minigraph-cactus/v2.1/hprc-v2.1-mc-grch38/hprc-v2.1-mc-grch38.gbz.db\",\n index: \"https://jbrowse.org/demos/hprc/hprc-v2.1-mc-grch38.haplotype-index.f3.db\"\n};\nfunction haplotypeIndexBeside(db) {\n const GBZ_DB = /\\.gbz\\.db$/i;\n return GBZ_DB.test(db) ? db.replace(GBZ_DB, \".haplotype-index.db\") : void 0;\n}\nfunction parseRegion(text) {\n const m = /^\\s*([^:\\s]+):([\\d,]+)-([\\d,]+)\\s*$/.exec(text);\n if (!m) {\n throw new Error(\n `\"${text}\" is not a region like chr6:160,614,798-160,647,758`\n );\n }\n const n = (s) => Number(s.replaceAll(\",\", \"\"));\n const region = { refName: m[1], start: n(m[2]), end: n(m[3]) };\n if (region.end <= region.start) {\n throw new Error(`\"${text}\" ends at or before it starts`);\n }\n return region;\n}\nasync function openGbz(db, index) {\n const base = await GBZBase.open({ source: db, haplotypeIndex: index });\n return { base, referenceSamples: await referenceSamplesOf(base) };\n}\nfunction referenceSample(ref, samples) {\n try {\n return resolveReferenceSample({\n configured: \"\",\n anchorPrefix: ref ?? samples[0] ?? \"\",\n referenceSamples: samples\n });\n } catch (e) {\n if (ref) {\n return ref;\n }\n throw e;\n }\n}\nasync function cutGbzRegion({ base, referenceSamples }, src, signal) {\n const region = parseRegion(src.region);\n const sample = referenceSample(src.referenceSample, referenceSamples);\n const query = await referencePathQuery(base, sample, region.refName);\n if (!query) {\n throw new Error(`${sample} has no indexed path named ${region.refName}`);\n }\n signal?.throwIfAborted();\n const wanted = src.haplotypes?.length ? src.haplotypes : void 0;\n const text = await cutWindowGFA(base, query, region.start, region.end, {\n context: src.context ?? GBZ_CUT_DEFAULTS.context,\n snarls: src.snarls ?? GBZ_CUT_DEFAULTS.snarls,\n limit: GBZ_CUT_DEFAULTS.limit,\n signal,\n keep: wanted ? (name) => haplotypeWanted(name, wanted) : void 0\n });\n return { text, region, sample };\n}\n\nexport {\n HPRC_GBZ,\n haplotypeIndexBeside,\n parseRegion,\n openGbz,\n cutGbzRegion\n};\n", "import { haplotypeIndexBeside } from '@jbrowse/bandage-core/gbzCut'\nimport { openLocation } from '@jbrowse/core/util/io'\n\nimport { locationName, readsSiblings, renamed } from '../locationName'\n\nimport type PluginManager from '@jbrowse/core/PluginManager'\nimport type { FileLocation } from '@jbrowse/core/util/types'\n\nexport function siblingCompanion(gbzDb: FileLocation) {\n return readsSiblings(gbzDb) && haplotypeIndexBeside(locationName(gbzDb))\n ? renamed(gbzDb, name => haplotypeIndexBeside(name)!)\n : undefined\n}\n\n// 404 is a missing file, and so is 403 from a bucket that does not let a\n// reader list it; anything else leaves the file's presence unknown\nfunction isMissing(error: unknown) {\n return (\n error instanceof Error &&\n (/\\bHTTP 40[34]\\b/.test(error.message) ||\n (error as { code?: unknown }).code === 'ENOENT')\n )\n}\n\n// `unreadable` is why a companion that may be there could not be checked, such\n// as a network error, or a server that answers a missing file without CORS\nexport async function findCompanion(\n gbzDb: FileLocation,\n pluginManager?: PluginManager,\n): Promise<{ location?: FileLocation; unreadable?: string }> {\n const candidate = siblingCompanion(gbzDb)\n if (!candidate) {\n return {}\n }\n try {\n await openLocation(candidate, pluginManager).stat()\n return { location: candidate }\n } catch (error) {\n return isMissing(error) ? {} : { unreadable: String(error) }\n }\n}\n", "import { unzip } from '@gmod/bgzf-filehandle'\nimport { parseGaf, parseGafLine } from '@jbrowse/bandage-core/gaf/parseGaf'\nimport { cachedSetup } from '@jbrowse/core/data_adapters/BaseAdapter'\n\nimport type { TabixIndexedFile } from '@gmod/tabix'\nimport type { GafRecord } from '@jbrowse/bandage-core/gaf/parseGaf'\n\ninterface WholeFile {\n stat(): Promise<{ size: number }>\n readFile(opts?: { signal?: AbortSignal }): Promise<Uint8Array>\n}\n\nexport interface GafReads {\n records: GafRecord[]\n // reads over the nodes before sampling\n total: number\n}\n\n// The tube map lays reads out on the main thread at about 40 ms per thousand\nexport const MAX_READS = 5000\n\n// An unindexed GAF is read whole, once, into the worker\nconst MAX_UNINDEXED_BYTES = 50 * 2 ** 20\n\nexport class UnindexedGafTooLargeError extends Error {\n override name = 'UnindexedGafTooLargeError'\n\n constructor(size: number) {\n super(\n `this GAF is ${(size / 2 ** 20).toFixed(0)} MB without an index; bgzip it and index it with tabix -p gaf`,\n )\n }\n}\n\n// Evenly spaced, so a pileup thins out rather than losing its right half\nfunction sample(records: GafRecord[], max: number) {\n if (records.length <= max) {\n return records\n }\n const stride = records.length / max\n return Array.from({ length: max }, (_, k) => records[Math.floor(k * stride)]!)\n}\n\nconst touches = (record: GafRecord, names: ReadonlySet<string>) =>\n record.path.some(step => names.has(step.name))\n\n// A GAF over a graph's segments: bgzipped and indexed with `tabix -p gaf`,\n// which keys each line by its lowest and highest node id and so needs numeric\n// segment names, or plain or gzipped and read whole.\nexport class GafFile {\n // Shared by every window, so the signal is cachedSetup's own, which aborts\n // once no caller is left waiting\n private whole = cachedSetup({\n setup: async ({ signal }) => {\n const { size } = await this.file.stat()\n if (size > MAX_UNINDEXED_BYTES) {\n throw new UnindexedGafTooLargeError(size)\n }\n const bytes = await this.file.readFile({ signal })\n const text = new TextDecoder().decode(\n bytes[0] === 0x1f && bytes[1] === 0x8b ? await unzip(bytes) : bytes,\n )\n return parseGaf(text)\n },\n })\n\n constructor(\n private file: WholeFile,\n private tabix?: TabixIndexedFile,\n ) {}\n\n private async readIndexed(\n tabix: TabixIndexedFile,\n names: ReadonlySet<string>,\n signal?: AbortSignal,\n ) {\n let lo = Infinity\n let hi = -Infinity\n for (const name of names) {\n const id = Number(name)\n if (Number.isInteger(id)) {\n lo = Math.min(lo, id)\n hi = Math.max(hi, id)\n }\n }\n const records: GafRecord[] = []\n if (lo <= hi) {\n // a GAF index ignores the refName, and getLines keeps a read whose node\n // ids span min..max when min < end && max >= start, hence hi + 1\n await tabix.getLines('{node}', lo, hi + 1, {\n signal,\n lineCallback: line => {\n const record = parseGafLine(line)\n if (record && touches(record, names)) {\n records.push(record)\n }\n },\n })\n }\n return records\n }\n\n async readsOver(\n names: ReadonlySet<string>,\n { signal }: { signal?: AbortSignal } = {},\n ): Promise<GafReads> {\n const records = this.tabix\n ? await this.readIndexed(this.tabix, names, signal)\n : (await this.whole({ signal })).filter(r => touches(r, names))\n return { records: sample(records, MAX_READS), total: records.length }\n }\n}\n", "import { GBZBase, SubgraphLimitError } from '@gmod/gbz-base'\nimport {\n NodeLimitError,\n cutWindowTables,\n haplotypePrefix,\n haplotypeWanted,\n nodeLimitError,\n referencePathQuery,\n referenceSamplesOf,\n resolveReferenceSample,\n} from '@jbrowse/bandage-core/gbzWindow'\nimport {\n assemblyByPanSNPrefix,\n panSNMatchesPrefix,\n resolvePanSNPrefix,\n} from '@jbrowse/bandage-core/pansn'\nimport { cachedSetup } from '@jbrowse/core/data_adapters/BaseAdapter'\nimport { getBpDisplayStr, updateStatus } from '@jbrowse/core/util'\nimport { openLocation } from '@jbrowse/core/util/io'\nimport { ObservableCreate } from '@jbrowse/core/util/rxjs'\n\nimport { findCompanion } from './companion.ts'\nimport { walkCutFor } from '../RgfaTabixAdapter/walkRowRuns.ts'\nimport { GafFile } from '../gaf/gafFile.ts'\nimport { locationIsSet } from '../locationName.ts'\nimport { openTabixSlot } from '../panSNTabix.ts'\nimport { ComparativeAdapterBase } from '../synteny/ComparativeAdapterBase.ts'\nimport SyntenyFeature from '../synteny/SyntenyFeature.ts'\nimport {\n laneAssemblyName,\n lanePairsStatus,\n requestedPairs,\n} from '../synteny/lanePairs.ts'\nimport { matchFeature } from '../synteny/matchFeature.ts'\n\nimport type { GbzBaseSyntenyAdapterConfig } from './configSchema.ts'\nimport type { SubgraphAdapterOptions } from '../GetSubgraph.ts'\nimport type { GafReads } from '../gaf/gafFile.ts'\nimport type { LaneFeatureOptions, LanePair } from '../synteny/lanePairs.ts'\nimport type {\n HaplotypeAlignment,\n HaplotypeRef,\n PairAlignment,\n PathName,\n PathQuery,\n Subgraph,\n} from '@gmod/gbz-base'\nimport type { BaseOptions } from '@jbrowse/core/data_adapters/BaseAdapter'\nimport type { Feature, SimpleFeatureSerialized } from '@jbrowse/core/util'\nimport type { FileLocation, Region } from '@jbrowse/core/util/types'\n\nexport class NoHaplotypeIndexError extends Error {\n override name = 'NoHaplotypeIndexError'\n\n constructor(unreadable?: string) {\n super(\n unreadable === undefined\n ? 'no haplotype index is set, so the walks cannot be named; build one with gbz-haplotype-index (cargo install gbz-haplotype-index) and set it as haplotypeIndexLocation'\n : `the haplotype index beside the graph database could not be checked (${unreadable}), so the walks cannot be named; if it is there, reload to try again, or set haplotypeIndexLocation`,\n )\n }\n}\n\ninterface GbzHaplotype {\n prefix: string\n sample: string\n haplotype: number\n contigs: string[]\n isReference: boolean\n}\n\n/**\n * One lane the header declares, in the shape MultiWaySyntenyDisplay's lane\n * picker reads: `name` is the lane's assembly name, the same string the\n * fetched features' mates carry, `label` its PanSN prefix and `group` its\n * sample, so a diploid sample's two lanes sit together\n */\ninterface GbzHeaderLane {\n name: string\n label: string\n group: string\n}\n\nexport { PairTargetError } from '../synteny/lanePairs.ts'\n\nclass HaplotypeWindowError extends Error {\n override name = 'HaplotypeWindowError'\n\n constructor(assemblyName: string, anchor: string) {\n super(\n `the graph is cut on its reference, ${anchor}; a window on ${assemblyName} has no reference coordinates to cut at. Open the graph from a ${anchor} view, and find this haplotype's lane there`,\n )\n }\n}\n\ntype ResolvedAlignment = Extract<HaplotypeAlignment, { resolved: true }>\n\n/**\n * The id of one haplotype fragment: the walk's GBWT position at its first node\n * in the window, which is a property of the graph, so the same fragment gets\n * the same id on a refetch of the same window and the display's hover and\n * selection survive it. `clipToRegion` suffixes the window.\n */\nfunction fragmentId(name: PathName, alignment: ResolvedAlignment) {\n return `${haplotypePrefix(name)}#${name.contig}@${alignment.start.node}.${alignment.start.offset}`\n}\n\nfunction fragmentFeature({\n alignment,\n assemblyName,\n refName,\n lane,\n}: {\n alignment: ResolvedAlignment\n assemblyName: string\n refName: string\n lane: string\n}) {\n const { refStart, refEnd, strand, cigar } = alignment\n if (refEnd <= refStart) {\n return undefined\n } else {\n const { name, hapStart, hapEnd } = alignment\n const id = fragmentId(name, alignment)\n const data: SimpleFeatureSerialized = {\n uniqueId: id,\n assemblyName,\n refName,\n start: refStart,\n end: refEnd,\n type: 'match',\n strand: strand === '-' ? -1 : 1,\n CIGAR: cigar,\n syntenyId: id,\n mate: {\n refName: name.contig,\n start: hapStart,\n end: hapEnd,\n assemblyName: lane,\n },\n }\n return new SyntenyFeature(data)\n }\n}\n\n// gbz-base writes the CIGAR along the target, which is the side a JBrowse\n// feature is, so it passes through unchanged\nfunction pairFeature({\n pair,\n lane,\n mateLane,\n}: {\n pair: PairAlignment\n lane: string\n mateLane: string\n}) {\n return matchFeature({\n target: {\n name: `${haplotypePrefix(pair.target)}#${pair.target.contig}`,\n assemblyName: lane,\n refName: pair.target.contig,\n start: pair.targetStart,\n end: pair.targetEnd,\n },\n query: {\n name: `${haplotypePrefix(pair.query)}#${pair.query.contig}`,\n assemblyName: mateLane,\n refName: pair.query.contig,\n start: pair.queryStart,\n end: pair.queryEnd,\n },\n strand: pair.strand,\n cigar: pair.cigar,\n matches: pair.matches,\n columns: pair.columns,\n })\n}\n\n// The longest detour a pair cut widens its context to take in\nconst BRIDGE_MAX_BP = 120_000\n\n/**\n * The widest gap, up to `BRIDGE_MAX_BP`, between two pieces of one kept walk:\n * a haplotype that left the cut's nodes for sequence the reference lacks and\n * came back. Two kept haplotypes that share that sequence only align across it\n * once its nodes are in the cut.\n */\nexport function widestWalkGap(subgraphs: Subgraph[], kept: Set<string>) {\n const byPath = new Map<string, [number, number][]>()\n for (const subgraph of subgraphs) {\n for (const { name, start, end } of subgraph.walkSpans()) {\n const prefix = haplotypePrefix(name)\n if (kept.has(prefix)) {\n const key = `${prefix}#${name.contig}`\n const spans = byPath.get(key)\n if (spans) {\n spans.push([start, end])\n } else {\n byPath.set(key, [[start, end]])\n }\n }\n }\n }\n let widest = 0\n for (const spans of byPath.values()) {\n spans.sort((a, b) => a[0] - b[0])\n for (let i = 1; i < spans.length; i++) {\n const gap = spans[i]![0] - spans[i - 1]![1]\n if (gap > 0 && gap <= BRIDGE_MAX_BP) {\n widest = Math.max(widest, gap)\n }\n }\n }\n return widest\n}\n\n// A window past the node limit fails as a zoom-in notice naming a span that\n// fits\nfunction zoomInNotice(error: unknown, what: string) {\n if (error instanceof NodeLimitError) {\n error.message = `Zoom in to about ${getBpDisplayStr(error.fitsBp)} to see ${what}`\n }\n return error\n}\n\n// gbz-base reads whole bp; a view asks at fractional ones\nconst wholeBp = (region: Region) => ({\n ...region,\n start: Math.floor(region.start),\n end: Math.ceil(region.end),\n})\n\nexport default class GbzBaseSyntenyAdapter extends ComparativeAdapterBase<GbzBaseSyntenyAdapterConfig> {\n private graph = cachedSetup({\n label: 'Opening pangenome database',\n setup: async () => {\n const gbzDb: FileLocation = this.getConf('gbzDbLocation')\n const configured: FileLocation = this.getConf('haplotypeIndexLocation')\n const companion = locationIsSet(configured)\n ? { location: configured }\n : await findCompanion(gbzDb, this.pluginManager)\n const indexLocation = companion.location\n const db = await GBZBase.open({\n source: openLocation(gbzDb, this.pluginManager),\n haplotypeIndex: indexLocation\n ? openLocation(indexLocation, this.pluginManager)\n : undefined,\n })\n const anchor = this.getConf('assemblyNames')[0]\n if (anchor === undefined) {\n throw new Error(\n 'GbzBaseSyntenyAdapter needs assemblyNames: its first entry is the assembly the reference sample is loaded as',\n )\n }\n const referenceSamples = await referenceSamplesOf(db)\n const referenceSample = resolveReferenceSample({\n configured: this.getConf('referenceSample'),\n anchorPrefix: resolvePanSNPrefix(this, anchor),\n referenceSamples,\n })\n return {\n db,\n anchor,\n referenceSample,\n referenceSamples,\n unreadableIndex: companion.unreadable,\n }\n },\n })\n\n /**\n * The indexed reference path a window on `refName` resolves against, or\n * undefined when the reference sample has no indexed path by that contig.\n */\n private async referenceQuery(\n refName: string,\n opts: BaseOptions,\n ): Promise<PathQuery | undefined> {\n const { db, referenceSample } = await this.graph(opts)\n return referencePathQuery(db, referenceSample, refName)\n }\n\n /**\n * Every haplotype the graph carries a walk for, `sample#haplotype` with its\n * contigs, from the Paths scan the adapter already makes; the lane picker's\n * source list. Reference samples are included and flagged.\n */\n async getHaplotypes(opts: BaseOptions = {}): Promise<GbzHaplotype[]> {\n const { db, referenceSamples } = await this.graph(opts)\n const byPrefix = new Map<string, GbzHaplotype>()\n for (const path of await db.paths()) {\n const prefix = haplotypePrefix(path.name)\n const entry = byPrefix.get(prefix)\n if (entry) {\n if (!entry.contigs.includes(path.name.contig)) {\n entry.contigs.push(path.name.contig)\n }\n } else {\n byPrefix.set(prefix, {\n prefix,\n sample: path.name.sample,\n haplotype: path.name.haplotype,\n contigs: [path.name.contig],\n isReference: referenceSamples.includes(path.name.sample),\n })\n }\n }\n return [...byPrefix.values()]\n }\n\n /**\n * `lanes` is every haplotype but the reference sample's own, which is the\n * anchor rather than a lane, declared up front so the display can offer the\n * whole graph's haplotypes before a fetch has placed any of them.\n */\n async getHeader(opts: BaseOptions = {}) {\n const { anchor, referenceSample } = await this.graph(opts)\n const asmByPrefix = assemblyByPanSNPrefix(this)\n const lanes: GbzHeaderLane[] = []\n for (const haplotype of await this.getHaplotypes(opts)) {\n if (haplotype.sample !== referenceSample) {\n lanes.push({\n name: laneAssemblyName(asmByPrefix, haplotype.prefix),\n label: haplotype.prefix,\n group: haplotype.sample,\n })\n }\n }\n return {\n hasCoarseTier: false,\n anchorAssemblyName: anchor,\n lanes,\n }\n }\n\n async getRefNames(opts: BaseOptions = {}) {\n const { db, anchor, referenceSample } = await this.graph(opts)\n const { assemblyName } = opts\n const prefix =\n assemblyName === anchor\n ? undefined\n : resolvePanSNPrefix(this, assemblyName)\n const contigs = (await db.paths())\n .filter(path =>\n assemblyName === anchor\n ? path.isIndexed && path.name.sample === referenceSample\n : panSNMatchesPrefix(haplotypePrefix(path.name), prefix),\n )\n .map(path => path.name.contig)\n return [...new Set(contigs)]\n }\n\n /**\n * The predicate gbz-base builds a window for, from the lanes a fetch asks\n * for. With a companion whose stray rows cover the window gbz-base walks\n * only the wanted haplotypes; otherwise it names every walk and keeps the\n * ones the predicate accepts, and throws rather than guess when a walk has\n * no name. Either way a cut holds those walks, the reference, and the nodes\n * they visit. Undefined when every haplotype is wanted.\n */\n //\n // `targetPrefix` is a pairwise view's one lane, and it belongs in here: left\n // to filter what came back, it handed gbz-base no predicate, and every\n // haplotype in the graph was walked for one to be kept.\n private keepPredicate(\n haplotypes: string[] | undefined,\n targetPrefix?: string,\n ) {\n const wantedPrefixes =\n haplotypes === undefined || haplotypes.length === 0\n ? undefined\n : haplotypes.map(lane => resolvePanSNPrefix(this, lane))\n return wantedPrefixes === undefined && targetPrefix === undefined\n ? undefined\n : (name: PathName) =>\n haplotypeWanted(name, wantedPrefixes) &&\n (targetPrefix === undefined ||\n panSNMatchesPrefix(haplotypePrefix(name), targetPrefix))\n }\n\n /**\n * The graph view's cut of the window, as tables: the reference walk, every\n * top-level snarl contained in it, and one walk per haplotype, PanSN-named\n * when the database (or its companion) carries the haplotype index and\n * `unknown#N` otherwise. The reference walk is the first, which is the one\n * the view anchors on by default. With `haplotypes` the walks are that\n * set's and the nodes those walks visit, the reference walk kept. With\n * `walkRows` the rows' runs come along, as on a walk-file cut.\n *\n * Only a window on the anchor can be cut: a haplotype lane's coordinates\n * are its own contig's, and the graph is indexed for random access on the\n * reference sample's paths alone.\n */\n async getSubgraph(region: Region, opts: SubgraphAdapterOptions = {}) {\n const { signal } = opts\n const { db, anchor, referenceSample, unreadableIndex } = await this.graph({\n signal,\n })\n const { assemblyName, refName, start, end } = wholeBp(region)\n if (assemblyName !== anchor) {\n throw new HaplotypeWindowError(assemblyName, anchor)\n }\n const query = await this.referenceQuery(refName, { signal })\n if (!query) {\n throw new Error(\n `${referenceSample} has no indexed path ${refName} in this graph`,\n )\n }\n const keep = this.keepPredicate(opts.haplotypes)\n if (keep !== undefined && !db.hasHaplotypeIndex) {\n throw new NoHaplotypeIndexError(unreadableIndex)\n }\n const cut = await cutWindowTables(db, query, start, end, {\n context: this.getConf('context'),\n snarls: opts.snarls ?? this.getConf('subgraphSnarls'),\n limit: this.getConf('nodeLimit'),\n signal: opts.signal,\n keep,\n }).catch((error: unknown) => {\n throw zoomInNotice(error, 'the graph')\n })\n return typeof cut === 'string'\n ? cut\n : walkCutFor(cut, { start, end }, opts.walkRows)\n }\n\n private gaf = cachedSetup({\n label: 'Opening reads',\n setup: async () => {\n const location: FileLocation = this.getConf('readsLocation')\n const index: FileLocation = this.getConf(['readsIndex', 'location'])\n return locationIsSet(location)\n ? new GafFile(\n openLocation(location, this.pluginManager),\n locationIsSet(index)\n ? openTabixSlot(this, ['readsLocation'], ['readsIndex'])\n : undefined,\n )\n : undefined\n },\n })\n\n /**\n * The reads over a cut's segments, from `readsLocation`, sampled down to\n * what the tube map lays out at interactive speed\n */\n async getReads(\n nodeNames: string[],\n opts: BaseOptions = {},\n ): Promise<GafReads> {\n const gaf = await this.gaf(opts)\n return gaf\n ? gaf.readsOver(new Set(nodeNames), opts)\n : { records: [], total: 0 }\n }\n\n private async laneHaplotypes(\n assemblyName: string,\n opts: BaseOptions,\n ): Promise<HaplotypeRef[]> {\n const prefix = resolvePanSNPrefix(this, assemblyName)\n return (await this.getHaplotypes(opts)).filter(haplotype =>\n panSNMatchesPrefix(haplotype.prefix, prefix),\n )\n }\n\n private async anchorFeatures(region: Region, opts: LaneFeatureOptions) {\n const { db } = await this.graph(opts)\n const { assemblyName, refName, start, end } = region\n const targetPrefix = resolvePanSNPrefix(this, opts.targetAssemblyName)\n const asmByPrefix = assemblyByPanSNPrefix(this)\n const keep = this.keepPredicate(opts.haplotypes, targetPrefix)\n const query = await this.referenceQuery(refName, opts)\n const nodeLimit: number = this.getConf('nodeLimit')\n const alignments = query\n ? await updateStatus(\n `Reading graph ${refName}:${start.toLocaleString()}-${end.toLocaleString()}`,\n opts.statusCallback,\n () =>\n db\n .getAlignments({\n path: query,\n start,\n end,\n context: this.getConf('context'),\n haplotypes: 'all',\n limit: nodeLimit,\n signal: opts.signal,\n keep,\n })\n .catch((error: unknown) => {\n throw zoomInNotice(\n nodeLimitError(error, nodeLimit, end - start) ?? error,\n 'lanes',\n )\n }),\n )\n : []\n return alignments.flatMap(alignment => {\n const feature = alignment.resolved\n ? fragmentFeature({\n alignment,\n assemblyName,\n refName,\n lane: laneAssemblyName(\n asmByPrefix,\n haplotypePrefix(alignment.name),\n ),\n })\n : undefined\n return feature === undefined ? [] : [feature]\n })\n }\n\n /**\n * Each pair's two lanes cut out of the anchor window alone, all pairs in one\n * cut, and each walk of a pair's query lane aligned to each walk of its\n * target lane on the nodes both visit, with no base compared. A pair's\n * records sit on its query lane's contigs, the lane the display draws on top.\n *\n * Whether the index has anchor rows changes only how fast gbz-base finds\n * the walks: both routes cut the same pieces.\n */\n private async pairFeatures(\n region: Region,\n pairs: LanePair[],\n opts: LaneFeatureOptions,\n ) {\n const { db } = await this.graph(opts)\n const { refName, start, end } = region\n const sides = await Promise.all(\n pairs.map(async pair => ({\n pair,\n featureSide: await this.laneHaplotypes(pair.queryAssemblyName, opts),\n mateSide: await this.laneHaplotypes(pair.targetAssemblyName, opts),\n })),\n )\n const answerable = sides.filter(\n ({ featureSide, mateSide }) =>\n featureSide.length > 0 && mateSide.length > 0,\n )\n const kept = new Set(\n answerable.flatMap(({ featureSide, mateSide }) =>\n [...featureSide, ...mateSide].map(haplotypePrefix),\n ),\n )\n const query = await this.referenceQuery(refName, opts)\n const nodeLimit: number = this.getConf('nodeLimit')\n const baseContext: number = this.getConf('context')\n const cut = (path: PathQuery, context: number) =>\n db.getSubgraphs({\n path,\n start,\n end,\n context,\n snarls: this.getConf('subgraphSnarls'),\n haplotypes: 'all',\n limit: nodeLimit,\n signal: opts.signal,\n keep: name => kept.has(haplotypePrefix(name)),\n })\n // a cut past the node limit for the wider context keeps the first one\n const bridged = async (path: PathQuery) => {\n const subgraphs = await cut(path, baseContext).catch((error: unknown) => {\n throw zoomInNotice(\n nodeLimitError(error, nodeLimit, end - start) ?? error,\n 'lanes',\n )\n })\n const gap = widestWalkGap(subgraphs, kept)\n const context = Math.ceil(gap / 2) + 2000\n return gap > 0 && context > baseContext\n ? cut(path, context).catch((error: unknown) => {\n if (error instanceof SubgraphLimitError) {\n return subgraphs\n }\n throw error\n })\n : subgraphs\n }\n const subgraphs =\n query && answerable.length > 0\n ? await updateStatus(lanePairsStatus(pairs), opts.statusCallback, () =>\n bridged(query),\n )\n : []\n return subgraphs.flatMap(subgraph =>\n answerable.flatMap(({ pair, featureSide, mateSide }) =>\n featureSide.flatMap(target =>\n mateSide.flatMap(mate =>\n subgraph\n .pairAlignments({ target, query: mate, bases: false })\n .map(alignment =>\n pairFeature({\n pair: alignment,\n lane: pair.queryAssemblyName,\n mateLane: pair.targetAssemblyName,\n }),\n ),\n ),\n ),\n ),\n )\n }\n\n getFeatures(asked: Region, opts: LaneFeatureOptions = {}) {\n const region = wholeBp(asked)\n return ObservableCreate<Feature>(async observer => {\n const { db, anchor, unreadableIndex } = await this.graph(opts)\n if (!db.hasHaplotypeIndex) {\n throw new NoHaplotypeIndexError(unreadableIndex)\n }\n // the graph is indexed on its reference alone, so a window on a\n // haplotype lane has no answer: a lane pair is read inside the anchor's\n if (region.assemblyName === anchor) {\n const pairs = requestedPairs(opts)\n const features =\n pairs.length > 0\n ? await this.pairFeatures(region, pairs, opts)\n : await this.anchorFeatures(region, opts)\n for (const feature of features) {\n observer.next(feature)\n }\n }\n observer.complete()\n }, opts.signal)\n }\n}\n"],
5
+ "mappings": "iyBAeA,SAASA,EAAqBC,EAAI,CAChC,IAAMC,EAAS,cACf,OAAOA,EAAO,KAAKD,CAAE,EAAIA,EAAG,QAAQC,EAAQ,qBAAqB,EAAI,MACvE,CCVO,SAASC,GAAiBC,EAAqB,CACpD,OAAOC,EAAcD,CAAK,GAAKE,EAAqBC,EAAaH,CAAK,CAAC,EACnEI,EAAQJ,EAAOK,GAAQH,EAAqBG,CAAI,CAAE,EAClD,MACN,CAIA,SAASC,GAAUC,EAAgB,CACjC,OACEA,aAAiB,QAChB,kBAAkB,KAAKA,EAAM,OAAO,GAClCA,EAA6B,OAAS,SAE7C,CAIA,eAAsBC,GACpBR,EACAS,EAC2D,CAC3D,IAAMC,EAAYX,GAAiBC,CAAK,EACxC,GAAI,CAACU,EACH,MAAO,CAAC,EAEV,GAAI,CACF,aAAMC,EAAaD,EAAWD,CAAa,EAAE,KAAK,EAC3C,CAAE,SAAUC,CAAU,CAC/B,OAASH,EAAO,CACd,OAAOD,GAAUC,CAAK,EAAI,CAAC,EAAI,CAAE,WAAY,OAAOA,CAAK,CAAE,CAC7D,CACF,CCrBO,IAAMK,GAAY,IAGnBC,GAAsB,GAAK,GAAK,GAEzBC,EAAN,cAAwC,KAAM,CAGnD,YAAYC,EAAc,CACxB,MACE,gBAAgBA,EAAO,GAAK,IAAI,QAAQ,CAAC,CAAC,+DAC5C,EALF,KAAS,KAAO,2BAMhB,CACF,EAGA,SAASC,GAAOC,EAAsBC,EAAa,CACjD,GAAID,EAAQ,QAAUC,EACpB,OAAOD,EAET,IAAME,EAASF,EAAQ,OAASC,EAChC,OAAO,MAAM,KAAK,CAAE,OAAQA,CAAI,EAAG,CAACE,EAAGC,IAAMJ,EAAQ,KAAK,MAAMI,EAAIF,CAAM,CAAC,CAAE,CAC/E,CAEA,IAAMG,GAAU,CAACC,EAAmBC,IAClCD,EAAO,KAAK,KAAKE,GAAQD,EAAM,IAAIC,EAAK,IAAI,CAAC,EAKlCC,EAAN,KAAc,CAiBnB,YACUC,EACAC,EACR,CAFQ,UAAAD,EACA,WAAAC,EAhBV,KAAQ,MAAQC,EAAY,CAC1B,MAAO,MAAO,CAAE,OAAAC,CAAO,IAAM,CAC3B,GAAM,CAAE,KAAAf,CAAK,EAAI,MAAM,KAAK,KAAK,KAAK,EACtC,GAAIA,EAAOF,GACT,MAAM,IAAIC,EAA0BC,CAAI,EAE1C,IAAMgB,EAAQ,MAAM,KAAK,KAAK,SAAS,CAAE,OAAAD,CAAO,CAAC,EAC3CE,EAAO,IAAI,YAAY,EAAE,OAC7BD,EAAM,CAAC,IAAM,IAAQA,EAAM,CAAC,IAAM,IAAO,MAAME,EAAMF,CAAK,EAAIA,CAChE,EACA,OAAOG,EAASF,CAAI,CACtB,CACF,CAAC,CAKE,CAEH,MAAc,YACZJ,EACAJ,EACAM,EACA,CACA,IAAIK,EAAK,IACLC,EAAK,KACT,QAAWC,KAAQb,EAAO,CACxB,IAAMc,EAAK,OAAOD,CAAI,EAClB,OAAO,UAAUC,CAAE,IACrBH,EAAK,KAAK,IAAIA,EAAIG,CAAE,EACpBF,EAAK,KAAK,IAAIA,EAAIE,CAAE,EAExB,CACA,IAAMrB,EAAuB,CAAC,EAC9B,OAAIkB,GAAMC,GAGR,MAAMR,EAAM,SAAS,SAAUO,EAAIC,EAAK,EAAG,CACzC,OAAAN,EACA,aAAcS,GAAQ,CACpB,IAAMhB,EAASiB,EAAaD,CAAI,EAC5BhB,GAAUD,GAAQC,EAAQC,CAAK,GACjCP,EAAQ,KAAKM,CAAM,CAEvB,CACF,CAAC,EAEIN,CACT,CAEA,MAAM,UACJO,EACA,CAAE,OAAAM,CAAO,EAA8B,CAAC,EACrB,CACnB,IAAMb,EAAU,KAAK,MACjB,MAAM,KAAK,YAAY,KAAK,MAAOO,EAAOM,CAAM,GAC/C,MAAM,KAAK,MAAM,CAAE,OAAAA,CAAO,CAAC,GAAG,OAAOW,GAAKnB,GAAQmB,EAAGjB,CAAK,CAAC,EAChE,MAAO,CAAE,QAASR,GAAOC,EAASL,EAAS,EAAG,MAAOK,EAAQ,MAAO,CACtE,CACF,EC5DO,IAAMyB,EAAN,cAAoC,KAAM,CAG/C,YAAYC,EAAqB,CAC/B,MACEA,IAAe,OACX,uKACA,uEAAuEA,CAAU,qGACvF,EAPF,KAAS,KAAO,uBAQhB,CACF,EAwBMC,EAAN,cAAmC,KAAM,CAGvC,YAAYC,EAAsBC,EAAgB,CAChD,MACE,sCAAsCA,CAAM,iBAAiBD,CAAY,kEAAkEC,CAAM,6CACnJ,EALF,KAAS,KAAO,sBAMhB,CACF,EAUA,SAASC,GAAWC,EAAgBC,EAA8B,CAChE,MAAO,GAAGC,EAAgBF,CAAI,CAAC,IAAIA,EAAK,MAAM,IAAIC,EAAU,MAAM,IAAI,IAAIA,EAAU,MAAM,MAAM,EAClG,CAEA,SAASE,GAAgB,CACvB,UAAAF,EACA,aAAAJ,EACA,QAAAO,EACA,KAAAC,CACF,EAKG,CACD,GAAM,CAAE,SAAAC,EAAU,OAAAC,EAAQ,OAAAC,EAAQ,MAAAC,CAAM,EAAIR,EAC5C,GAAI,EAAAM,GAAUD,GAEP,CACL,GAAM,CAAE,KAAAN,EAAM,SAAAU,EAAU,OAAAC,CAAO,EAAIV,EAC7BW,EAAKb,GAAWC,EAAMC,CAAS,EAC/BY,EAAgC,CACpC,SAAUD,EACV,aAAAf,EACA,QAAAO,EACA,MAAOE,EACP,IAAKC,EACL,KAAM,QACN,OAAQC,IAAW,IAAM,GAAK,EAC9B,MAAOC,EACP,UAAWG,EACX,KAAM,CACJ,QAASZ,EAAK,OACd,MAAOU,EACP,IAAKC,EACL,aAAcN,CAChB,CACF,EACA,OAAO,IAAIS,GAAeD,CAAI,CAChC,CACF,CAIA,SAASE,GAAY,CACnB,KAAAC,EACA,KAAAX,EACA,SAAAY,CACF,EAIG,CACD,OAAOC,GAAa,CAClB,OAAQ,CACN,KAAM,GAAGhB,EAAgBc,EAAK,MAAM,CAAC,IAAIA,EAAK,OAAO,MAAM,GAC3D,aAAcX,EACd,QAASW,EAAK,OAAO,OACrB,MAAOA,EAAK,YACZ,IAAKA,EAAK,SACZ,EACA,MAAO,CACL,KAAM,GAAGd,EAAgBc,EAAK,KAAK,CAAC,IAAIA,EAAK,MAAM,MAAM,GACzD,aAAcC,EACd,QAASD,EAAK,MAAM,OACpB,MAAOA,EAAK,WACZ,IAAKA,EAAK,QACZ,EACA,OAAQA,EAAK,OACb,MAAOA,EAAK,MACZ,QAASA,EAAK,QACd,QAASA,EAAK,OAChB,CAAC,CACH,CAGA,IAAMG,GAAgB,KAQf,SAASC,GAAcC,EAAuBC,EAAmB,CACtE,IAAMC,EAAS,IAAI,IACnB,QAAWC,KAAYH,EACrB,OAAW,CAAE,KAAArB,EAAM,MAAAyB,EAAO,IAAAC,CAAI,IAAKF,EAAS,UAAU,EAAG,CACvD,IAAMG,EAASzB,EAAgBF,CAAI,EACnC,GAAIsB,EAAK,IAAIK,CAAM,EAAG,CACpB,IAAMC,EAAM,GAAGD,CAAM,IAAI3B,EAAK,MAAM,GAC9B6B,EAAQN,EAAO,IAAIK,CAAG,EACxBC,EACFA,EAAM,KAAK,CAACJ,EAAOC,CAAG,CAAC,EAEvBH,EAAO,IAAIK,EAAK,CAAC,CAACH,EAAOC,CAAG,CAAC,CAAC,CAElC,CACF,CAEF,IAAII,EAAS,EACb,QAAWD,KAASN,EAAO,OAAO,EAAG,CACnCM,EAAM,KAAK,CAACE,EAAGC,IAAMD,EAAE,CAAC,EAAIC,EAAE,CAAC,CAAC,EAChC,QAASC,EAAI,EAAGA,EAAIJ,EAAM,OAAQI,IAAK,CACrC,IAAMC,EAAML,EAAMI,CAAC,EAAG,CAAC,EAAIJ,EAAMI,EAAI,CAAC,EAAG,CAAC,EACtCC,EAAM,GAAKA,GAAOf,KACpBW,EAAS,KAAK,IAAIA,EAAQI,CAAG,EAEjC,CACF,CACA,OAAOJ,CACT,CAIA,SAASK,EAAaC,EAAgBC,EAAc,CAClD,OAAID,aAAiBE,IACnBF,EAAM,QAAU,oBAAoBG,EAAgBH,EAAM,MAAM,CAAC,WAAWC,CAAI,IAE3ED,CACT,CAGA,IAAMI,GAAWC,IAAoB,CACnC,GAAGA,EACH,MAAO,KAAK,MAAMA,EAAO,KAAK,EAC9B,IAAK,KAAK,KAAKA,EAAO,GAAG,CAC3B,GAEqBC,EAArB,cAAmDC,EAAoD,CAAvG,kCACE,KAAQ,MAAQC,EAAY,CAC1B,MAAO,6BACP,MAAO,SAAY,CACjB,IAAMC,EAAsB,KAAK,QAAQ,eAAe,EAClDC,EAA2B,KAAK,QAAQ,wBAAwB,EAChEC,EAAYC,EAAcF,CAAU,EACtC,CAAE,SAAUA,CAAW,EACvB,MAAMG,GAAcJ,EAAO,KAAK,aAAa,EAC3CK,EAAgBH,EAAU,SAC1BI,EAAK,MAAMC,EAAQ,KAAK,CAC5B,OAAQC,EAAaR,EAAO,KAAK,aAAa,EAC9C,eAAgBK,EACZG,EAAaH,EAAe,KAAK,aAAa,EAC9C,MACN,CAAC,EACKpD,EAAS,KAAK,QAAQ,eAAe,EAAE,CAAC,EAC9C,GAAIA,IAAW,OACb,MAAM,IAAI,MACR,8GACF,EAEF,IAAMwD,EAAmB,MAAMC,EAAmBJ,CAAE,EAC9CK,EAAkBC,EAAuB,CAC7C,WAAY,KAAK,QAAQ,iBAAiB,EAC1C,aAAcC,EAAmB,KAAM5D,CAAM,EAC7C,iBAAAwD,CACF,CAAC,EACD,MAAO,CACL,GAAAH,EACA,OAAArD,EACA,gBAAA0D,EACA,iBAAAF,EACA,gBAAiBP,EAAU,UAC7B,CACF,CACF,CAAC,EA8JD,KAAQ,IAAMH,EAAY,CACxB,MAAO,gBACP,MAAO,SAAY,CACjB,IAAMe,EAAyB,KAAK,QAAQ,eAAe,EACrDC,EAAsB,KAAK,QAAQ,CAAC,aAAc,UAAU,CAAC,EACnE,OAAOZ,EAAcW,CAAQ,EACzB,IAAIE,EACFR,EAAaM,EAAU,KAAK,aAAa,EACzCX,EAAcY,CAAK,EACfE,GAAc,KAAM,CAAC,eAAe,EAAG,CAAC,YAAY,CAAC,EACrD,MACN,EACA,MACN,CACF,CAAC,EAtKD,MAAc,eACZ1D,EACA2D,EACgC,CAChC,GAAM,CAAE,GAAAZ,EAAI,gBAAAK,CAAgB,EAAI,MAAM,KAAK,MAAMO,CAAI,EACrD,OAAOC,EAAmBb,EAAIK,EAAiBpD,CAAO,CACxD,CAOA,MAAM,cAAc2D,EAAoB,CAAC,EAA4B,CACnE,GAAM,CAAE,GAAAZ,EAAI,iBAAAG,CAAiB,EAAI,MAAM,KAAK,MAAMS,CAAI,EAChDE,EAAW,IAAI,IACrB,QAAWC,KAAQ,MAAMf,EAAG,MAAM,EAAG,CACnC,IAAMxB,EAASzB,EAAgBgE,EAAK,IAAI,EAClCC,EAAQF,EAAS,IAAItC,CAAM,EAC7BwC,EACGA,EAAM,QAAQ,SAASD,EAAK,KAAK,MAAM,GAC1CC,EAAM,QAAQ,KAAKD,EAAK,KAAK,MAAM,EAGrCD,EAAS,IAAItC,EAAQ,CACnB,OAAAA,EACA,OAAQuC,EAAK,KAAK,OAClB,UAAWA,EAAK,KAAK,UACrB,QAAS,CAACA,EAAK,KAAK,MAAM,EAC1B,YAAaZ,EAAiB,SAASY,EAAK,KAAK,MAAM,CACzD,CAAC,CAEL,CACA,MAAO,CAAC,GAAGD,EAAS,OAAO,CAAC,CAC9B,CAOA,MAAM,UAAUF,EAAoB,CAAC,EAAG,CACtC,GAAM,CAAE,OAAAjE,EAAQ,gBAAA0D,CAAgB,EAAI,MAAM,KAAK,MAAMO,CAAI,EACnDK,EAAcC,EAAsB,IAAI,EACxCC,EAAyB,CAAC,EAChC,QAAWC,KAAa,MAAM,KAAK,cAAcR,CAAI,EAC/CQ,EAAU,SAAWf,GACvBc,EAAM,KAAK,CACT,KAAME,EAAiBJ,EAAaG,EAAU,MAAM,EACpD,MAAOA,EAAU,OACjB,MAAOA,EAAU,MACnB,CAAC,EAGL,MAAO,CACL,cAAe,GACf,mBAAoBzE,EACpB,MAAAwE,CACF,CACF,CAEA,MAAM,YAAYP,EAAoB,CAAC,EAAG,CACxC,GAAM,CAAE,GAAAZ,EAAI,OAAArD,EAAQ,gBAAA0D,CAAgB,EAAI,MAAM,KAAK,MAAMO,CAAI,EACvD,CAAE,aAAAlE,CAAa,EAAIkE,EACnBpC,EACJ9B,IAAiBC,EACb,OACA4D,EAAmB,KAAM7D,CAAY,EACrC4E,GAAW,MAAMtB,EAAG,MAAM,GAC7B,OAAOe,GACNrE,IAAiBC,EACboE,EAAK,WAAaA,EAAK,KAAK,SAAWV,EACvCkB,EAAmBxE,EAAgBgE,EAAK,IAAI,EAAGvC,CAAM,CAC3D,EACC,IAAIuC,GAAQA,EAAK,KAAK,MAAM,EAC/B,MAAO,CAAC,GAAG,IAAI,IAAIO,CAAO,CAAC,CAC7B,CAcQ,cACNE,EACAC,EACA,CACA,IAAMC,EACJF,IAAe,QAAaA,EAAW,SAAW,EAC9C,OACAA,EAAW,IAAItE,GAAQqD,EAAmB,KAAMrD,CAAI,CAAC,EAC3D,OAAOwE,IAAmB,QAAaD,IAAiB,OACpD,OACC5E,GACC8E,EAAgB9E,EAAM6E,CAAc,IACnCD,IAAiB,QAChBF,EAAmBxE,EAAgBF,CAAI,EAAG4E,CAAY,EAChE,CAeA,MAAM,YAAYnC,EAAgBsB,EAA+B,CAAC,EAAG,CACnE,GAAM,CAAE,OAAAgB,CAAO,EAAIhB,EACb,CAAE,GAAAZ,EAAI,OAAArD,EAAQ,gBAAA0D,EAAiB,gBAAAwB,CAAgB,EAAI,MAAM,KAAK,MAAM,CACxE,OAAAD,CACF,CAAC,EACK,CAAE,aAAAlF,EAAc,QAAAO,EAAS,MAAAqB,EAAO,IAAAC,CAAI,EAAIc,GAAQC,CAAM,EAC5D,GAAI5C,IAAiBC,EACnB,MAAM,IAAIF,EAAqBC,EAAcC,CAAM,EAErD,IAAMmF,EAAQ,MAAM,KAAK,eAAe7E,EAAS,CAAE,OAAA2E,CAAO,CAAC,EAC3D,GAAI,CAACE,EACH,MAAM,IAAI,MACR,GAAGzB,CAAe,wBAAwBpD,CAAO,gBACnD,EAEF,IAAM8E,EAAO,KAAK,cAAcnB,EAAK,UAAU,EAC/C,GAAImB,IAAS,QAAa,CAAC/B,EAAG,kBAC5B,MAAM,IAAIzD,EAAsBsF,CAAe,EAEjD,IAAMG,EAAM,MAAMC,EAAgBjC,EAAI8B,EAAOxD,EAAOC,EAAK,CACvD,QAAS,KAAK,QAAQ,SAAS,EAC/B,OAAQqC,EAAK,QAAU,KAAK,QAAQ,gBAAgB,EACpD,MAAO,KAAK,QAAQ,WAAW,EAC/B,OAAQA,EAAK,OACb,KAAAmB,CACF,CAAC,EAAE,MAAO9C,GAAmB,CAC3B,MAAMD,EAAaC,EAAO,WAAW,CACvC,CAAC,EACD,OAAO,OAAO+C,GAAQ,SAClBA,EACAE,GAAWF,EAAK,CAAE,MAAA1D,EAAO,IAAAC,CAAI,EAAGqC,EAAK,QAAQ,CACnD,CAsBA,MAAM,SACJuB,EACAvB,EAAoB,CAAC,EACF,CACnB,IAAMwB,EAAM,MAAM,KAAK,IAAIxB,CAAI,EAC/B,OAAOwB,EACHA,EAAI,UAAU,IAAI,IAAID,CAAS,EAAGvB,CAAI,EACtC,CAAE,QAAS,CAAC,EAAG,MAAO,CAAE,CAC9B,CAEA,MAAc,eACZlE,EACAkE,EACyB,CACzB,IAAMpC,EAAS+B,EAAmB,KAAM7D,CAAY,EACpD,OAAQ,MAAM,KAAK,cAAckE,CAAI,GAAG,OAAOQ,GAC7CG,EAAmBH,EAAU,OAAQ5C,CAAM,CAC7C,CACF,CAEA,MAAc,eAAec,EAAgBsB,EAA0B,CACrE,GAAM,CAAE,GAAAZ,CAAG,EAAI,MAAM,KAAK,MAAMY,CAAI,EAC9B,CAAE,aAAAlE,EAAc,QAAAO,EAAS,MAAAqB,EAAO,IAAAC,CAAI,EAAIe,EACxCmC,EAAelB,EAAmB,KAAMK,EAAK,kBAAkB,EAC/DK,EAAcC,EAAsB,IAAI,EACxCa,EAAO,KAAK,cAAcnB,EAAK,WAAYa,CAAY,EACvDK,EAAQ,MAAM,KAAK,eAAe7E,EAAS2D,CAAI,EAC/CyB,EAAoB,KAAK,QAAQ,WAAW,EAyBlD,OAxBmBP,EACf,MAAMQ,EACJ,iBAAiBrF,CAAO,IAAIqB,EAAM,eAAe,CAAC,IAAIC,EAAI,eAAe,CAAC,GAC1EqC,EAAK,eACL,IACEZ,EACG,cAAc,CACb,KAAM8B,EACN,MAAAxD,EACA,IAAAC,EACA,QAAS,KAAK,QAAQ,SAAS,EAC/B,WAAY,MACZ,MAAO8D,EACP,OAAQzB,EAAK,OACb,KAAAmB,CACF,CAAC,EACA,MAAO9C,GAAmB,CACzB,MAAMD,EACJuD,EAAetD,EAAOoD,EAAW9D,EAAMD,CAAK,GAAKW,EACjD,OACF,CACF,CAAC,CACP,EACA,CAAC,GACa,QAAQnC,GAAa,CACrC,IAAM0F,EAAU1F,EAAU,SACtBE,GAAgB,CACd,UAAAF,EACA,aAAAJ,EACA,QAAAO,EACA,KAAMoE,EACJJ,EACAlE,EAAgBD,EAAU,IAAI,CAChC,CACF,CAAC,EACD,OACJ,OAAO0F,IAAY,OAAY,CAAC,EAAI,CAACA,CAAO,CAC9C,CAAC,CACH,CAWA,MAAc,aACZlD,EACAmD,EACA7B,EACA,CACA,GAAM,CAAE,GAAAZ,CAAG,EAAI,MAAM,KAAK,MAAMY,CAAI,EAC9B,CAAE,QAAA3D,EAAS,MAAAqB,EAAO,IAAAC,CAAI,EAAIe,EAQ1BoD,GAPQ,MAAM,QAAQ,IAC1BD,EAAM,IAAI,MAAM5E,IAAS,CACvB,KAAAA,EACA,YAAa,MAAM,KAAK,eAAeA,EAAK,kBAAmB+C,CAAI,EACnE,SAAU,MAAM,KAAK,eAAe/C,EAAK,mBAAoB+C,CAAI,CACnE,EAAE,CACJ,GACyB,OACvB,CAAC,CAAE,YAAA+B,EAAa,SAAAC,CAAS,IACvBD,EAAY,OAAS,GAAKC,EAAS,OAAS,CAChD,EACMzE,EAAO,IAAI,IACfuE,EAAW,QAAQ,CAAC,CAAE,YAAAC,EAAa,SAAAC,CAAS,IAC1C,CAAC,GAAGD,EAAa,GAAGC,CAAQ,EAAE,IAAI7F,CAAe,CACnD,CACF,EACM+E,EAAQ,MAAM,KAAK,eAAe7E,EAAS2D,CAAI,EAC/CyB,EAAoB,KAAK,QAAQ,WAAW,EAC5CQ,EAAsB,KAAK,QAAQ,SAAS,EAC5Cb,EAAM,CAACjB,EAAiB+B,IAC5B9C,EAAG,aAAa,CACd,KAAAe,EACA,MAAAzC,EACA,IAAAC,EACA,QAAAuE,EACA,OAAQ,KAAK,QAAQ,gBAAgB,EACrC,WAAY,MACZ,MAAOT,EACP,OAAQzB,EAAK,OACb,KAAM/D,GAAQsB,EAAK,IAAIpB,EAAgBF,CAAI,CAAC,CAC9C,CAAC,EAEGkG,EAAU,MAAOhC,GAAoB,CACzC,IAAM7C,EAAY,MAAM8D,EAAIjB,EAAM8B,CAAW,EAAE,MAAO5D,GAAmB,CACvE,MAAMD,EACJuD,EAAetD,EAAOoD,EAAW9D,EAAMD,CAAK,GAAKW,EACjD,OACF,CACF,CAAC,EACKF,EAAMd,GAAcC,EAAWC,CAAI,EACnC2E,EAAU,KAAK,KAAK/D,EAAM,CAAC,EAAI,IACrC,OAAOA,EAAM,GAAK+D,EAAUD,EACxBb,EAAIjB,EAAM+B,CAAO,EAAE,MAAO7D,GAAmB,CAC3C,GAAIA,aAAiB+D,EACnB,OAAO9E,EAET,MAAMe,CACR,CAAC,EACDf,CACN,EAOA,OALE4D,GAASY,EAAW,OAAS,EACzB,MAAMJ,EAAaW,GAAgBR,CAAK,EAAG7B,EAAK,eAAgB,IAC9DmC,EAAQjB,CAAK,CACf,EACA,CAAC,GACU,QAAQzD,GACvBqE,EAAW,QAAQ,CAAC,CAAE,KAAA7E,EAAM,YAAA8E,EAAa,SAAAC,CAAS,IAChDD,EAAY,QAAQO,GAClBN,EAAS,QAAQO,IACf9E,EACG,eAAe,CAAE,OAAA6E,EAAQ,MAAOC,GAAM,MAAO,EAAM,CAAC,EACpD,IAAIrG,IACHc,GAAY,CACV,KAAMd,GACN,KAAMe,EAAK,kBACX,SAAUA,EAAK,kBACjB,CAAC,CACH,CACJ,CACF,CACF,CACF,CACF,CAEA,YAAYuF,EAAexC,EAA2B,CAAC,EAAG,CACxD,IAAMtB,EAASD,GAAQ+D,CAAK,EAC5B,OAAOC,GAA0B,MAAMC,GAAY,CACjD,GAAM,CAAE,GAAAtD,EAAI,OAAArD,EAAQ,gBAAAkF,CAAgB,EAAI,MAAM,KAAK,MAAMjB,CAAI,EAC7D,GAAI,CAACZ,EAAG,kBACN,MAAM,IAAIzD,EAAsBsF,CAAe,EAIjD,GAAIvC,EAAO,eAAiB3C,EAAQ,CAClC,IAAM8F,EAAQc,GAAe3C,CAAI,EAC3B4C,EACJf,EAAM,OAAS,EACX,MAAM,KAAK,aAAanD,EAAQmD,EAAO7B,CAAI,EAC3C,MAAM,KAAK,eAAetB,EAAQsB,CAAI,EAC5C,QAAW4B,KAAWgB,EACpBF,EAAS,KAAKd,CAAO,CAEzB,CACAc,EAAS,SAAS,CACpB,EAAG1C,EAAK,MAAM,CAChB,CACF",
6
+ "names": ["haplotypeIndexBeside", "db", "GBZ_DB", "siblingCompanion", "gbzDb", "readsSiblings", "haplotypeIndexBeside", "locationName", "renamed", "name", "isMissing", "error", "findCompanion", "pluginManager", "candidate", "openLocation", "MAX_READS", "MAX_UNINDEXED_BYTES", "UnindexedGafTooLargeError", "size", "sample", "records", "max", "stride", "_", "k", "touches", "record", "names", "step", "GafFile", "file", "tabix", "cachedSetup", "signal", "bytes", "text", "unzip", "parseGaf", "lo", "hi", "name", "id", "line", "parseGafLine", "r", "NoHaplotypeIndexError", "unreadable", "HaplotypeWindowError", "assemblyName", "anchor", "fragmentId", "name", "alignment", "haplotypePrefix", "fragmentFeature", "refName", "lane", "refStart", "refEnd", "strand", "cigar", "hapStart", "hapEnd", "id", "data", "SyntenyFeature", "pairFeature", "pair", "mateLane", "matchFeature", "BRIDGE_MAX_BP", "widestWalkGap", "subgraphs", "kept", "byPath", "subgraph", "start", "end", "prefix", "key", "spans", "widest", "a", "b", "i", "gap", "zoomInNotice", "error", "what", "NodeLimitError", "getBpDisplayStr", "wholeBp", "region", "GbzBaseSyntenyAdapter", "ComparativeAdapterBase", "cachedSetup", "gbzDb", "configured", "companion", "locationIsSet", "findCompanion", "indexLocation", "db", "GBZBase", "openLocation", "referenceSamples", "referenceSamplesOf", "referenceSample", "resolveReferenceSample", "resolvePanSNPrefix", "location", "index", "GafFile", "openTabixSlot", "opts", "referencePathQuery", "byPrefix", "path", "entry", "asmByPrefix", "assemblyByPanSNPrefix", "lanes", "haplotype", "laneAssemblyName", "contigs", "panSNMatchesPrefix", "haplotypes", "targetPrefix", "wantedPrefixes", "haplotypeWanted", "signal", "unreadableIndex", "query", "keep", "cut", "cutWindowTables", "walkCutFor", "nodeNames", "gaf", "nodeLimit", "updateStatus", "nodeLimitError", "feature", "pairs", "answerable", "featureSide", "mateSide", "baseContext", "context", "bridged", "SubgraphLimitError", "lanePairsStatus", "target", "mate", "asked", "ObservableCreate", "observer", "requestedPairs", "features"]
7
+ }
@@ -1,4 +1,4 @@
1
- import{a as C}from"./chunk-FBGWMI6M.js";import"./chunk-YAADZ2AI.js";import"./chunk-G63ZHV6F.js";import"./chunk-6CAAGNOP.js";import{k as d}from"./chunk-GUNY323E.js";import"./chunk-NKBE5KYL.js";import"./chunk-SIGBYUXP.js";import"./chunk-CPCTMFXV.js";import"./chunk-ABNOMXT3.js";import"./chunk-S3QFRJR3.js";import"./chunk-JT3RXYIL.js";import"./chunk-77GS3NVS.js";import"./chunk-TXI43PEN.js";import"./chunk-YBLNN5FU.js";import"./chunk-MVVXRZXV.js";import"./chunk-IIGYYZI5.js";import"./chunk-5OI6LO4I.js";import"./chunk-NDRLCGPC.js";import"./chunk-TENPGIPX.js";import{g as x}from"./chunk-6OPB7FOJ.js";import{A as l,F as L,H as p,d as s,i as v,o as n,x as G}from"./chunk-E3WYW7Z6.js";import{a as y,b as a,c as r,d as i}from"./chunk-JMBDTNUR.js";var w=d()({paper:{padding:16,margin:8,maxWidth:560,marginInline:"auto",display:"flex",flexDirection:"column",gap:12},row:{display:"flex",alignItems:"center",justifyContent:"space-between",gap:8}}),E=i(function({model:e}){let{classes:t}=w();return e.isLoading?r(l,{className:t.paper,"data-testid":"graph-genome-loading",children:[r("div",{className:t.row,children:[a(v,{variant:"h6",message:e.statusMessage}),a(n,{size:"small","data-testid":"graph-genome-cancel",onClick:()=>{e.cancelLoad()},children:"Cancel"})]}),a(G,{variant:"indeterminate"})]}):e.error?a(l,{className:t.paper,children:a(s,{error:e.error,onReset:()=>{e.retryLoad()}})}):a(l,{className:t.paper,"data-testid":"graph-genome-load-canceled",children:r("div",{className:t.row,children:[a(p,{variant:"h6",children:"Loading canceled"}),a(n,{variant:"contained",size:"small","data-testid":"graph-genome-retry",onClick:()=>{e.retryLoad()},children:"Retry"})]})})}),g=E;var F=`H VN:Z:1.0
1
+ import{a as C}from"./chunk-MURN3JPK.js";import"./chunk-BUY2PV7D.js";import"./chunk-2MCOY3YG.js";import"./chunk-NWJYLP6Z.js";import{k as d}from"./chunk-VJ3SO4CH.js";import"./chunk-WEGWTC52.js";import"./chunk-FITPB5X6.js";import"./chunk-WAMBI7OX.js";import"./chunk-W6EQLQ6H.js";import"./chunk-5O6YRJHL.js";import"./chunk-N6PQZSAX.js";import"./chunk-TUKJRVBS.js";import"./chunk-FVK2THE4.js";import"./chunk-YBLNN5FU.js";import"./chunk-MVVXRZXV.js";import"./chunk-ZMOMFU2B.js";import"./chunk-5OI6LO4I.js";import"./chunk-NDRLCGPC.js";import"./chunk-TENPGIPX.js";import{h as x}from"./chunk-YU7T5WOQ.js";import{A as l,F as L,H as p,e as s,j as v,o as n,x as G}from"./chunk-UGCNPRKC.js";import{a as y,b as a,c as r,d as i}from"./chunk-JMBDTNUR.js";var w=d()({paper:{padding:16,margin:8,maxWidth:560,marginInline:"auto",display:"flex",flexDirection:"column",gap:12},row:{display:"flex",alignItems:"center",justifyContent:"space-between",gap:8}}),E=i(function({model:e}){let{classes:t}=w();return e.isLoading?r(l,{className:t.paper,"data-testid":"graph-genome-loading",children:[r("div",{className:t.row,children:[a(v,{variant:"h6",message:e.statusMessage}),a(n,{size:"small","data-testid":"graph-genome-cancel",onClick:()=>{e.cancelLoad()},children:"Cancel"})]}),a(G,{variant:"indeterminate"})]}):e.error?a(l,{className:t.paper,children:a(s,{error:e.error,onReset:()=>{e.retryLoad()}})}):a(l,{className:t.paper,"data-testid":"graph-genome-load-canceled",children:r("div",{className:t.row,children:[a(p,{variant:"h6",children:"Loading canceled"}),a(n,{variant:"contained",size:"small","data-testid":"graph-genome-retry",onClick:()=>{e.retryLoad()},children:"Retry"})]})})}),g=E;var F=`H VN:Z:1.0
2
2
  S 1 ACGT
3
3
  S 2 GGCC
4
4
  S 3 TTAA
@@ -7,4 +7,4 @@ L 1 + 2 + 0M
7
7
  L 1 + 3 + 0M
8
8
  L 2 + 4 + 0M
9
9
  L 3 + 4 + 0M`,S=d()({column:{display:"flex",flexDirection:"column",gap:8},rowEnd:{display:"flex",gap:8,alignItems:"flex-end"},rowCenter:{display:"flex",gap:8,alignItems:"center"},flex1:{flex:1},paper:{padding:16,margin:8,maxWidth:560,marginInline:"auto"},header:{marginBottom:12},footer:{marginTop:12,display:"flex",justifyContent:"flex-end"}}),T=i(function({model:e}){let{classes:t}=S(),[m,M]=x("");function u(){let o=m.trim();o&&e.loadGFAFromLocation({uri:o,locationType:"UriLocation"})}function N(o){let c=o.target.files?.[0];o.target.value="",c&&c.text().then(f=>e.loadGFA(f,c.name)).catch(f=>{e.setError(f)})}return r(l,{className:t.paper,children:[a("div",{className:t.header,children:a(p,{variant:"h6",children:"Load a GFA graph"})}),r("div",{className:t.column,children:[r("div",{className:t.rowCenter,children:[r(n,{variant:"outlined",component:"label",size:"small",children:["Choose file",a("input",{type:"file",accept:".gfa,.gfa1,.gfa2",hidden:!0,onChange:o=>{N(o)}})]}),a(p,{variant:"caption",color:"text.secondary",children:"Whole-file GFA; best for small/medium graphs."})]}),r("div",{className:t.rowEnd,children:[a(L,{size:"small",label:"URL",placeholder:"https://example.com/graph.gfa",value:m,onChange:o=>{M(o.target.value)},onKeyDown:o=>{o.key==="Enter"&&u()},className:t.flex1}),a(n,{variant:"contained",onClick:()=>{u()},disabled:!m.trim(),children:"Open"})]})]}),a("div",{className:t.footer,children:a(n,{size:"small",onClick:()=>{e.loadGFA(F,"Example graph")},children:"Load 4-node example"})}),e.error?a(s,{error:e.error}):null]})}),P=T;var A=i(function({model:e}){return e.hasGraph?a(C,{model:e}):e.canRetryLoad&&(e.error||e.loadCanceled)?a(g,{model:e}):r(y,{children:[e.isLoading?a(g,{model:e}):null,a("div",{hidden:e.isLoading,children:a(P,{model:e})})]})}),J=A;export{J as default};
10
- //# sourceMappingURL=GraphGenomeView-BG7UHKCT.js.map
10
+ //# sourceMappingURL=GraphGenomeView-NJ5XCREB.js.map
@@ -0,0 +1,6 @@
1
+ import{b as G}from"./chunk-BUY2PV7D.js";import{b as v,c as x}from"./chunk-E2QPSDQL.js";import{k as h}from"./chunk-VJ3SO4CH.js";import"./chunk-WAMBI7OX.js";import"./chunk-TUKJRVBS.js";import"./chunk-MVVXRZXV.js";import"./chunk-ZMOMFU2B.js";import"./chunk-5OI6LO4I.js";import{h as u}from"./chunk-YU7T5WOQ.js";import{D as H,F as y,H as c,n as g,s as b,w as S,z as C}from"./chunk-UGCNPRKC.js";import{b as o,c as l,d as i}from"./chunk-JMBDTNUR.js";var F=h()({section:{marginBottom:24},formControl:{minWidth:200}}),M=[{value:0,label:"None"},{value:1,label:"1 hop"},{value:2,label:"2 hops"}],I=i(function({model:t}){let{classes:n}=F(),a=t.sourceAdapter;return a?.type==="RgfaTabixAdapter"&&!v(a)?l("div",{className:n.section,children:[l(b,{className:n.formControl,children:[o(S,{children:"Graph context"}),o(H,{value:t.subgraphContext,label:"Graph context","data-testid":"graph-context-select",onChange:p=>{t.setSubgraphContext(p.target.value),t.recut()},children:M.map(({value:p,label:e})=>o(C,{value:p,children:e},p))})]}),o(c,{variant:"caption",component:"div",color:"text.secondary",children:"How far the cut follows links out of the region. At none, a detour that leaves the window draws as two stubs. No effect on gfa-to-tabix's default anchored index."})]}):null}),L=I;var N=h()({section:{marginBottom:24}});function P(r){return r.split(/[\s,]+/).filter(t=>t!=="")}var T=i(function({model:t}){let[n,a]=u(t.cutHaplotypes?.join(", ")??""),p=()=>{let e=P(n);e.join(`
2
+ `)===(t.cutHaplotypes??[]).join(`
3
+ `)||(t.setSubgraphHaplotypes(e),t.recut())};return o(y,{fullWidth:!0,label:"Haplotypes",placeholder:"every haplotype",value:n,slotProps:{htmlInput:{"data-testid":"graph-haplotypes-field"}},onChange:e=>{a(e.target.value)},onBlur:()=>{p()},onKeyDown:e=>{e.key==="Enter"&&p()}})}),w=i(function({model:t,names:n}){let[a,p]=u(t.cutHaplotypes??[]),[e,d]=u(!1),m=s=>{s.join(`
4
+ `)!==(t.cutHaplotypes??[]).join(`
5
+ `)&&(t.setSubgraphHaplotypes(s),t.recut())};return o(g,{multiple:!0,freeSolo:!0,disableCloseOnSelect:!0,options:n,value:a,open:e,onOpen:()=>{d(!0)},onClose:()=>{d(!1),m(a)},onChange:(s,f)=>{p(f),e||m(f)},renderInput:s=>o(y,{...s,label:"Haplotypes",placeholder:a.length===0?"every haplotype":void 0,slotProps:{...s.slotProps,htmlInput:{...s.slotProps.htmlInput,"data-testid":"graph-haplotypes-field"}}})})}),A=i(function({model:t}){let{classes:n}=N(),a=t.haplotypeNames,p=t.cutHaplotypes?.join(",")??"",e=t.sourceAdapter;return e&&x(e)?l("div",{className:n.section,children:[a?o(w,{model:t,names:a},p):o(T,{model:t},p),l(c,{variant:"caption",component:"div",color:"text.secondary",children:[a?`Which of the graph's ${a.length.toLocaleString()} haplotypes to cut, or a PanSN prefix (HG002 for both).`:"Which haplotypes to cut, as lane assembly names or PanSN prefixes (HG002#1, or HG002 for both) separated by commas."," ","Empty is every haplotype."]})]}):null}),k=A;function D({model:r,open:t,onClose:n}){return l(G,{model:r,open:t,onClose:n,children:[o(L,{model:r}),o(k,{model:r})]})}export{D as default};
6
+ //# sourceMappingURL=GraphTrackSettingsDialog-47PQ6LM7.js.map
@@ -0,0 +1,7 @@
1
+ {
2
+ "version": 3,
3
+ "sources": ["../../src/LinearGraphDisplay/components/SubgraphContextSelect.tsx", "../../src/LinearGraphDisplay/components/SubgraphHaplotypesField.tsx", "../../src/LinearGraphDisplay/components/GraphTrackSettingsDialog.tsx"],
4
+ "sourcesContent": ["import {\n FormControl,\n InputLabel,\n MenuItem,\n Select,\n Typography,\n} from '@mui/material'\nimport { observer } from 'mobx-react'\nimport { makeStyles } from 'tss-react/mui'\n\nimport { walkIndexed } from '../../graphTrackConfig'\n\nimport type { LinearGraphCutModel } from '../model'\n\nconst useStyles = makeStyles()({\n section: {\n marginBottom: 24,\n },\n formControl: {\n minWidth: 200,\n },\n})\n\n// Only an rGFA cut from segment rows reads hops: a GBZ cut widens by its\n// adapter's bp `context`, and a walk-indexed cut follows its walks.\n// Hops past the region's own segments, each costing a tabix query per\n// off-reference segment already reached. One is the default: at none a single\n// detour draws as two unrelated stubs, which is a wrong picture of the graph\n// rather than a cheaper one. Two exists for a graph whose alleles have alleles\n// (HPRC's amylase window keeps growing at 2, the E. coli paa locus does not), and\n// it stops there because a frontier is still not a bubble decomposition. Cut an\n// exact slice with gfatools when that is what is wanted.\nconst SUBGRAPH_CONTEXTS = [\n { value: 0, label: 'None' },\n { value: 1, label: '1 hop' },\n { value: 2, label: '2 hops' },\n]\n\nconst SubgraphContextSelect = observer(function SubgraphContextSelect({\n model,\n}: {\n model: LinearGraphCutModel\n}) {\n const { classes } = useStyles()\n const adapter = model.sourceAdapter\n return adapter?.type === 'RgfaTabixAdapter' && !walkIndexed(adapter) ? (\n <div className={classes.section}>\n <FormControl className={classes.formControl}>\n <InputLabel>Graph context</InputLabel>\n <Select\n value={model.subgraphContext}\n label=\"Graph context\"\n data-testid=\"graph-context-select\"\n onChange={e => {\n model.setSubgraphContext(e.target.value)\n model.recut()\n }}\n >\n {SUBGRAPH_CONTEXTS.map(({ value, label }) => (\n <MenuItem key={value} value={value}>\n {label}\n </MenuItem>\n ))}\n </Select>\n </FormControl>\n <Typography variant=\"caption\" component=\"div\" color=\"text.secondary\">\n How far the cut follows links out of the region. At none, a detour that\n leaves the window draws as two stubs. No effect on gfa-to-tabix&apos;s\n default anchored index.\n </Typography>\n </div>\n ) : null\n})\n\nexport default SubgraphContextSelect\n", "import { useState } from 'react'\n\nimport { Autocomplete, TextField, Typography } from '@mui/material'\nimport { observer } from 'mobx-react'\nimport { makeStyles } from 'tss-react/mui'\n\nimport { cutsByHaplotype } from '../../graphTrackConfig'\n\nimport type { LinearGraphCutModel } from '../model'\n\nconst useStyles = makeStyles()({\n section: {\n marginBottom: 24,\n },\n})\n\nfunction parseHaplotypeList(text: string) {\n return text.split(/[\\s,]+/).filter(name => name !== '')\n}\n\nconst HaplotypeListField = observer(function HaplotypeListField({\n model,\n}: {\n model: LinearGraphCutModel\n}) {\n const [draft, setDraft] = useState(model.cutHaplotypes?.join(', ') ?? '')\n const apply = () => {\n const parsed = parseHaplotypeList(draft)\n const unchanged =\n parsed.join('\\n') === (model.cutHaplotypes ?? []).join('\\n')\n if (!unchanged) {\n model.setSubgraphHaplotypes(parsed)\n model.recut()\n }\n }\n return (\n <TextField\n fullWidth\n label=\"Haplotypes\"\n placeholder=\"every haplotype\"\n value={draft}\n slotProps={{ htmlInput: { 'data-testid': 'graph-haplotypes-field' } }}\n onChange={e => {\n setDraft(e.target.value)\n }}\n onBlur={() => {\n apply()\n }}\n onKeyDown={e => {\n if (e.key === 'Enter') {\n apply()\n }\n }}\n />\n )\n})\n\n// The haplotypes a walk-indexed graph's header names, searched as typed. A\n// typed prefix the list lacks (HG002 for both haplotypes) is kept too. The\n// set applies when the list closes or, with the list closed, as it changes.\nconst HaplotypePickField = observer(function HaplotypePickField({\n model,\n names,\n}: {\n model: LinearGraphCutModel\n names: string[]\n}) {\n const [draft, setDraft] = useState(model.cutHaplotypes ?? [])\n const [open, setOpen] = useState(false)\n const apply = (haplotypes: string[]) => {\n if (haplotypes.join('\\n') !== (model.cutHaplotypes ?? []).join('\\n')) {\n model.setSubgraphHaplotypes(haplotypes)\n model.recut()\n }\n }\n return (\n <Autocomplete\n multiple\n freeSolo\n disableCloseOnSelect\n options={names}\n value={draft}\n open={open}\n onOpen={() => {\n setOpen(true)\n }}\n onClose={() => {\n setOpen(false)\n apply(draft)\n }}\n onChange={(_, value) => {\n setDraft(value)\n if (!open) {\n apply(value)\n }\n }}\n renderInput={params => (\n <TextField\n {...params}\n label=\"Haplotypes\"\n placeholder={draft.length === 0 ? 'every haplotype' : undefined}\n slotProps={{\n ...params.slotProps,\n htmlInput: {\n ...params.slotProps.htmlInput,\n 'data-testid': 'graph-haplotypes-field',\n },\n }}\n />\n )}\n />\n )\n})\n\n// Only a GBZ cut and a walk-indexed rGFA cut read the set. The inner field is\n// keyed on it so an outside change (a restored session, a launch) resets the\n// draft rather than fighting it.\nconst SubgraphHaplotypesField = observer(function SubgraphHaplotypesField({\n model,\n}: {\n model: LinearGraphCutModel\n}) {\n const { classes } = useStyles()\n const names = model.haplotypeNames\n const key = model.cutHaplotypes?.join(',') ?? ''\n const adapter = model.sourceAdapter\n return adapter && cutsByHaplotype(adapter) ? (\n <div className={classes.section}>\n {names ? (\n <HaplotypePickField key={key} model={model} names={names} />\n ) : (\n <HaplotypeListField key={key} model={model} />\n )}\n <Typography variant=\"caption\" component=\"div\" color=\"text.secondary\">\n {names\n ? `Which of the graph's ${names.length.toLocaleString()} haplotypes to cut, or a PanSN prefix (HG002 for both).`\n : 'Which haplotypes to cut, as lane assembly names or PanSN prefixes (HG002#1, or HG002 for both) separated by commas.'}{' '}\n Empty is every haplotype.\n </Typography>\n </div>\n ) : null\n})\n\nexport default SubgraphHaplotypesField\n", "import SubgraphContextSelect from './SubgraphContextSelect'\nimport SubgraphHaplotypesField from './SubgraphHaplotypesField'\nimport GraphSettingsDialog from '../../GraphGenomeView/components/GraphSettingsDialog'\n\nimport type { LinearGraphCutModel } from '../model'\n\nexport default function GraphTrackSettingsDialog({\n model,\n open,\n onClose,\n}: {\n model: LinearGraphCutModel\n open: boolean\n onClose: () => void\n}) {\n return (\n <GraphSettingsDialog model={model} open={open} onClose={onClose}>\n <SubgraphContextSelect model={model} />\n <SubgraphHaplotypesField model={model} />\n </GraphSettingsDialog>\n )\n}\n"],
5
+ "mappings": "2bAcA,IAAMA,EAAYC,EAAW,EAAE,CAC7B,QAAS,CACP,aAAc,EAChB,EACA,YAAa,CACX,SAAU,GACZ,CACF,CAAC,EAWKC,EAAoB,CACxB,CAAE,MAAO,EAAG,MAAO,MAAO,EAC1B,CAAE,MAAO,EAAG,MAAO,OAAQ,EAC3B,CAAE,MAAO,EAAG,MAAO,QAAS,CAC9B,EAEMC,EAAwBC,EAAS,SAA+B,CACpE,MAAAC,CACF,EAEG,CACD,GAAM,CAAE,QAAAC,CAAQ,EAAIN,EAAU,EACxBO,EAAUF,EAAM,cACtB,OAAOE,GAAS,OAAS,oBAAsB,CAACC,EAAYD,CAAO,EACjEE,EAAC,OAAI,UAAWH,EAAQ,QACtB,UAAAG,EAACC,EAAA,CAAY,UAAWJ,EAAQ,YAC9B,UAAAK,EAACC,EAAA,CAAW,yBAAa,EACzBD,EAACE,EAAA,CACC,MAAOR,EAAM,gBACb,MAAM,gBACN,cAAY,uBACZ,SAAUS,GAAK,CACbT,EAAM,mBAAmBS,EAAE,OAAO,KAAK,EACvCT,EAAM,MAAM,CACd,EAEC,SAAAH,EAAkB,IAAI,CAAC,CAAE,MAAAa,EAAO,MAAAC,CAAM,IACrCL,EAACM,EAAA,CAAqB,MAAOF,EAC1B,SAAAC,GADYD,CAEf,CACD,EACH,GACF,EACAJ,EAACO,EAAA,CAAW,QAAQ,UAAU,UAAU,MAAM,MAAM,iBAAiB,6KAIrE,GACF,EACE,IACN,CAAC,EAEMC,EAAQhB,EChEf,IAAMiB,EAAYC,EAAW,EAAE,CAC7B,QAAS,CACP,aAAc,EAChB,CACF,CAAC,EAED,SAASC,EAAmBC,EAAc,CACxC,OAAOA,EAAK,MAAM,QAAQ,EAAE,OAAOC,GAAQA,IAAS,EAAE,CACxD,CAEA,IAAMC,EAAqBC,EAAS,SAA4B,CAC9D,MAAAC,CACF,EAEG,CACD,GAAM,CAACC,EAAOC,CAAQ,EAAIC,EAASH,EAAM,eAAe,KAAK,IAAI,GAAK,EAAE,EAClEI,EAAQ,IAAM,CAClB,IAAMC,EAASV,EAAmBM,CAAK,EAErCI,EAAO,KAAK;AAAA,CAAI,KAAOL,EAAM,eAAiB,CAAC,GAAG,KAAK;AAAA,CAAI,IAE3DA,EAAM,sBAAsBK,CAAM,EAClCL,EAAM,MAAM,EAEhB,EACA,OACEM,EAACC,EAAA,CACC,UAAS,GACT,MAAM,aACN,YAAY,kBACZ,MAAON,EACP,UAAW,CAAE,UAAW,CAAE,cAAe,wBAAyB,CAAE,EACpE,SAAU,GAAK,CACbC,EAAS,EAAE,OAAO,KAAK,CACzB,EACA,OAAQ,IAAM,CACZE,EAAM,CACR,EACA,UAAW,GAAK,CACV,EAAE,MAAQ,SACZA,EAAM,CAEV,EACF,CAEJ,CAAC,EAKKI,EAAqBT,EAAS,SAA4B,CAC9D,MAAAC,EACA,MAAAS,CACF,EAGG,CACD,GAAM,CAACR,EAAOC,CAAQ,EAAIC,EAASH,EAAM,eAAiB,CAAC,CAAC,EACtD,CAACU,EAAMC,CAAO,EAAIR,EAAS,EAAK,EAChCC,EAASQ,GAAyB,CAClCA,EAAW,KAAK;AAAA,CAAI,KAAOZ,EAAM,eAAiB,CAAC,GAAG,KAAK;AAAA,CAAI,IACjEA,EAAM,sBAAsBY,CAAU,EACtCZ,EAAM,MAAM,EAEhB,EACA,OACEM,EAACO,EAAA,CACC,SAAQ,GACR,SAAQ,GACR,qBAAoB,GACpB,QAASJ,EACT,MAAOR,EACP,KAAMS,EACN,OAAQ,IAAM,CACZC,EAAQ,EAAI,CACd,EACA,QAAS,IAAM,CACbA,EAAQ,EAAK,EACbP,EAAMH,CAAK,CACb,EACA,SAAU,CAACa,EAAGC,IAAU,CACtBb,EAASa,CAAK,EACTL,GACHN,EAAMW,CAAK,CAEf,EACA,YAAaC,GACXV,EAACC,EAAA,CACE,GAAGS,EACJ,MAAM,aACN,YAAaf,EAAM,SAAW,EAAI,kBAAoB,OACtD,UAAW,CACT,GAAGe,EAAO,UACV,UAAW,CACT,GAAGA,EAAO,UAAU,UACpB,cAAe,wBACjB,CACF,EACF,EAEJ,CAEJ,CAAC,EAKKC,EAA0BlB,EAAS,SAAiC,CACxE,MAAAC,CACF,EAEG,CACD,GAAM,CAAE,QAAAkB,CAAQ,EAAIzB,EAAU,EACxBgB,EAAQT,EAAM,eACdmB,EAAMnB,EAAM,eAAe,KAAK,GAAG,GAAK,GACxCoB,EAAUpB,EAAM,cACtB,OAAOoB,GAAWC,EAAgBD,CAAO,EACvCE,EAAC,OAAI,UAAWJ,EAAQ,QACrB,UAAAT,EACCH,EAACE,EAAA,CAA6B,MAAOR,EAAO,MAAOS,GAA1BU,CAAiC,EAE1Db,EAACR,EAAA,CAA6B,MAAOE,GAAZmB,CAAmB,EAE9CG,EAACC,EAAA,CAAW,QAAQ,UAAU,UAAU,MAAM,MAAM,iBACjD,UAAAd,EACG,wBAAwBA,EAAM,OAAO,eAAe,CAAC,0DACrD,sHAAuH,IAAI,6BAEjI,GACF,EACE,IACN,CAAC,EAEMe,EAAQP,ECzIA,SAARQ,EAA0C,CAC/C,MAAAC,EACA,KAAAC,EACA,QAAAC,CACF,EAIG,CACD,OACEC,EAACC,EAAA,CAAoB,MAAOJ,EAAO,KAAMC,EAAM,QAASC,EACtD,UAAAG,EAACC,EAAA,CAAsB,MAAON,EAAO,EACrCK,EAACE,EAAA,CAAwB,MAAOP,EAAO,GACzC,CAEJ",
6
+ "names": ["useStyles", "makeStyles", "SUBGRAPH_CONTEXTS", "SubgraphContextSelect", "observer", "model", "classes", "adapter", "walkIndexed", "jsxs", "FormControl", "jsx", "InputLabel", "Select", "e", "value", "label", "MenuItem", "Typography", "SubgraphContextSelect_default", "useStyles", "makeStyles", "parseHaplotypeList", "text", "name", "HaplotypeListField", "observer", "model", "draft", "setDraft", "useState", "apply", "parsed", "jsx", "TextField", "HaplotypePickField", "names", "open", "setOpen", "haplotypes", "Autocomplete", "_", "value", "params", "SubgraphHaplotypesField", "classes", "key", "adapter", "cutsByHaplotype", "jsxs", "Typography", "SubgraphHaplotypesField_default", "GraphTrackSettingsDialog", "model", "open", "onClose", "jsxs", "GraphSettingsDialog_default", "jsx", "SubgraphContextSelect_default", "SubgraphHaplotypesField_default"]
7
+ }
@@ -1,2 +1,2 @@
1
- import{o as v,p as C}from"./chunk-SIGBYUXP.js";import"./chunk-CPCTMFXV.js";import"./chunk-77GS3NVS.js";import"./chunk-IIGYYZI5.js";import{F as g,H as t,c as u,o as y,q as c,r as d,z as s}from"./chunk-E3WYW7Z6.js";import{a as h,b as l,c as p,d as m}from"./chunk-JMBDTNUR.js";var b=m(function({model:o,onClose:r}){return p(u,{open:!0,onClose:r,title:"Color highlighted haplotypes",children:[p(d,{sx:{display:"grid",gridTemplateColumns:"auto 1fr 1fr",alignItems:"center",gap:1.5},children:[l(t,{variant:"subtitle2",children:"Haplotype"}),l(t,{variant:"subtitle2",children:"Color by"}),l(t,{variant:"subtitle2",children:"Palette"}),o.drawnWalks.map(({name:a,encoding:n})=>l(f,{label:o.walkLabel(a),field:n.field,scheme:n.scheme,onField:e=>{o.setWalkColor(a,{field:e})},onScheme:e=>{o.setWalkColor(a,{scheme:e})}},a))]}),l(c,{children:l(y,{variant:"contained",onClick:r,children:"Done"})})]})});function f({label:i,field:o,scheme:r,onField:a,onScheme:n}){return p(h,{children:[l(t,{variant:"body2",children:i}),l(g,{select:!0,size:"small",value:o,slotProps:{select:{SelectDisplayProps:{"aria-label":`Color ${i} by`}}},onChange:e=>{a(e.target.value)},children:v.map(e=>l(s,{value:e.value,children:e.label},e.value))}),l(g,{select:!0,size:"small",value:r,slotProps:{select:{SelectDisplayProps:{"aria-label":`Palette for ${i}`}}},onChange:e=>{n(e.target.value)},children:C.filter(e=>e.value!=="rainbow"||o==="reference").map(e=>l(s,{value:e.value,children:e.label},e.value))})]})}var P=b;export{P as default};
2
- //# sourceMappingURL=HighlightColorDialog-YT55H62Y.js.map
1
+ import{p as v,q as C}from"./chunk-FITPB5X6.js";import"./chunk-WAMBI7OX.js";import"./chunk-TUKJRVBS.js";import"./chunk-ZMOMFU2B.js";import{F as g,H as t,d as u,o as y,q as c,r as d,z as s}from"./chunk-UGCNPRKC.js";import{a as h,b as l,c as p,d as m}from"./chunk-JMBDTNUR.js";var b=m(function({model:o,onClose:r}){return p(u,{open:!0,onClose:r,title:"Color highlighted haplotypes",children:[p(d,{sx:{display:"grid",gridTemplateColumns:"auto 1fr 1fr",alignItems:"center",gap:1.5},children:[l(t,{variant:"subtitle2",children:"Haplotype"}),l(t,{variant:"subtitle2",children:"Color by"}),l(t,{variant:"subtitle2",children:"Palette"}),o.drawnWalks.map(({name:a,encoding:n})=>l(f,{label:o.walkLabel(a),field:n.field,scheme:n.scheme,onField:e=>{o.setWalkColor(a,{field:e})},onScheme:e=>{o.setWalkColor(a,{scheme:e})}},a))]}),l(c,{children:l(y,{variant:"contained",onClick:r,children:"Done"})})]})});function f({label:i,field:o,scheme:r,onField:a,onScheme:n}){return p(h,{children:[l(t,{variant:"body2",children:i}),l(g,{select:!0,size:"small",value:o,slotProps:{select:{SelectDisplayProps:{"aria-label":`Color ${i} by`}}},onChange:e=>{a(e.target.value)},children:v.map(e=>l(s,{value:e.value,children:e.label},e.value))}),l(g,{select:!0,size:"small",value:r,slotProps:{select:{SelectDisplayProps:{"aria-label":`Palette for ${i}`}}},onChange:e=>{n(e.target.value)},children:C.filter(e=>e.value!=="rainbow"||o==="reference").map(e=>l(s,{value:e.value,children:e.label},e.value))})]})}var P=b;export{P as default};
2
+ //# sourceMappingURL=HighlightColorDialog-UK5VDDV6.js.map
@@ -1,2 +1,2 @@
1
- import{a as g}from"./chunk-FBGWMI6M.js";import"./chunk-YAADZ2AI.js";import"./chunk-G63ZHV6F.js";import"./chunk-6CAAGNOP.js";import"./chunk-GUNY323E.js";import"./chunk-NKBE5KYL.js";import"./chunk-SIGBYUXP.js";import"./chunk-CPCTMFXV.js";import"./chunk-ABNOMXT3.js";import"./chunk-S3QFRJR3.js";import"./chunk-JT3RXYIL.js";import"./chunk-77GS3NVS.js";import"./chunk-TXI43PEN.js";import"./chunk-YBLNN5FU.js";import"./chunk-MVVXRZXV.js";import"./chunk-IIGYYZI5.js";import"./chunk-5OI6LO4I.js";import"./chunk-NDRLCGPC.js";import"./chunk-TENPGIPX.js";import{g as y}from"./chunk-6OPB7FOJ.js";import"./chunk-E3WYW7Z6.js";import{b as e,c as l,d as T}from"./chunk-JMBDTNUR.js";var s=JBrowseExports["@jbrowse/display-kit/DisplayBackgroundProgress"];var v=s&&s.__esModule?s.default:s;var i=JBrowseExports["@jbrowse/display-kit/DisplayErrorBar"];var h=i&&i.__esModule?i.default:i;var a=JBrowseExports["@jbrowse/display-kit/DisplayLoadingOverlay"];var O=a&&a.__esModule?a.default:a;var n=JBrowseExports["@jbrowse/display-kit/TooLargeMessage"];var f=n&&n.__esModule?n.default:n;var t=JBrowseExports["@jbrowse/display-ui"],p=r=>t!=null&&r in Object(t)?t[r]:Object.assign(function(){throw new Error("This JBrowse does not provide '"+r+"' from '@jbrowse/display-ui' to plugins: GraphGenomeView was built for a different JBrowse version, so update the plugin or JBrowse")},{jbrowseHostMissing:r});var w=p("BottomRightCornerContext");var c=p("TrackOverlayPortal");var L=p("useChromeOverlayOverride");var j=t&&t.__esModule?t.default:t;var D={TooLarge:f,ErrorBar:h,Loading:O,BackgroundProgress:v};function d({model:r,phase:o,drawn:u,testid:_,style:C,children:E,...B}){let{TooLarge:G,ErrorBar:S,Loading:b,BackgroundProgress:m}=L()??D,[P,A]=y(null),x={...B,"data-testid":_,"data-display-id":r.configuration.displayId,"data-display-drawn":u,"data-display-phase":o};return o==="tooLarge"?e("div",{style:{display:"contents"},...x,children:e(c,{children:e("div",{style:{pointerEvents:"auto"},children:e(G,{model:r})})})}):l("div",{style:{position:"relative",height:r.height,...C},...x,children:[e(w,{value:P,children:E}),l(c,{children:[e(S,{model:r,visible:o==="error"}),e(b,{model:r,visible:o==="loading"||o==="canceled",immediate:!u}),e("div",{ref:A,style:{position:"absolute",bottom:2,right:2,display:"flex",flexDirection:"column",alignItems:"flex-end",gap:4,pointerEvents:"none"},children:e(m,{model:r,visible:o==="ready"})})]})]})}var I=T(function({model:o}){return e(d,{model:o,phase:o.displayPhase,drawn:o.painted&&o.geometryPainted,testid:"linear-graph-display","data-layout":o.chosenLayoutMode,"data-cut-tier":o.cutTier,"data-recuts":o.recuts,"data-node-count":o.hasGraph?o.nodeCount:void 0,"data-loading":o.isLoading?"":void 0,style:{width:o.paneWidth,height:o.height,overflow:"hidden"},children:e(g,{model:o,ownChrome:!1})})}),Q=I;export{Q as default};
2
- //# sourceMappingURL=LinearGraphDisplay-SPILNH4J.js.map
1
+ import{a as g}from"./chunk-MURN3JPK.js";import"./chunk-BUY2PV7D.js";import"./chunk-2MCOY3YG.js";import"./chunk-NWJYLP6Z.js";import"./chunk-VJ3SO4CH.js";import"./chunk-WEGWTC52.js";import"./chunk-FITPB5X6.js";import"./chunk-WAMBI7OX.js";import"./chunk-W6EQLQ6H.js";import"./chunk-5O6YRJHL.js";import"./chunk-N6PQZSAX.js";import"./chunk-TUKJRVBS.js";import"./chunk-FVK2THE4.js";import"./chunk-YBLNN5FU.js";import"./chunk-MVVXRZXV.js";import"./chunk-ZMOMFU2B.js";import"./chunk-5OI6LO4I.js";import"./chunk-NDRLCGPC.js";import"./chunk-TENPGIPX.js";import{h as y}from"./chunk-YU7T5WOQ.js";import"./chunk-UGCNPRKC.js";import{b as e,c as l,d as T}from"./chunk-JMBDTNUR.js";var s=JBrowseExports["@jbrowse/display-kit/DisplayBackgroundProgress"];var v=s&&s.__esModule?s.default:s;var i=JBrowseExports["@jbrowse/display-kit/DisplayErrorBar"];var h=i&&i.__esModule?i.default:i;var a=JBrowseExports["@jbrowse/display-kit/DisplayLoadingOverlay"];var O=a&&a.__esModule?a.default:a;var n=JBrowseExports["@jbrowse/display-kit/TooLargeMessage"];var f=n&&n.__esModule?n.default:n;var t=JBrowseExports["@jbrowse/display-ui"],p=r=>t!=null&&r in Object(t)?t[r]:Object.assign(function(){throw new Error("This JBrowse does not provide '"+r+"' from '@jbrowse/display-ui' to plugins: GraphGenomeView was built for a different JBrowse version, so update the plugin or JBrowse")},{jbrowseHostMissing:r});var w=p("BottomRightCornerContext");var c=p("TrackOverlayPortal");var L=p("useChromeOverlayOverride");var j=t&&t.__esModule?t.default:t;var D={TooLarge:f,ErrorBar:h,Loading:O,BackgroundProgress:v};function d({model:r,phase:o,drawn:u,testid:_,style:C,children:E,...B}){let{TooLarge:G,ErrorBar:S,Loading:b,BackgroundProgress:m}=L()??D,[P,A]=y(null),x={...B,"data-testid":_,"data-display-id":r.configuration.displayId,"data-display-drawn":u,"data-display-phase":o};return o==="tooLarge"?e("div",{style:{display:"contents"},...x,children:e(c,{children:e("div",{style:{pointerEvents:"auto"},children:e(G,{model:r})})})}):l("div",{style:{position:"relative",height:r.height,...C},...x,children:[e(w,{value:P,children:E}),l(c,{children:[e(S,{model:r,visible:o==="error"}),e(b,{model:r,visible:o==="loading"||o==="canceled",immediate:!u}),e("div",{ref:A,style:{position:"absolute",bottom:2,right:2,display:"flex",flexDirection:"column",alignItems:"flex-end",gap:4,pointerEvents:"none"},children:e(m,{model:r,visible:o==="ready"})})]})]})}var I=T(function({model:o}){return e(d,{model:o,phase:o.displayPhase,drawn:o.painted&&o.geometryPainted,testid:"linear-graph-display","data-layout":o.chosenLayoutMode,"data-cut-tier":o.cutTier,"data-recuts":o.recuts,"data-node-count":o.hasGraph?o.nodeCount:void 0,"data-loading":o.isLoading?"":void 0,style:{width:o.paneWidth,height:o.height,overflow:"hidden"},children:e(g,{model:o,ownChrome:!1})})}),Q=I;export{Q as default};
2
+ //# sourceMappingURL=LinearGraphDisplay-TXTKAACH.js.map
@@ -0,0 +1,2 @@
1
+ import{a as f,d as c,f as h}from"./chunk-W6EQLQ6H.js";import{a as g}from"./chunk-KOB7YIMR.js";import{b,h as m,m as p,n as u}from"./chunk-FHZXW7IG.js";import"./chunk-YBLNN5FU.js";import"./chunk-MVVXRZXV.js";import"./chunk-ZMOMFU2B.js";import"./chunk-OGQ6OQDN.js";import"./chunk-NDRLCGPC.js";import{a as i,o as l}from"./chunk-TENPGIPX.js";var a=class extends b{constructor(){super(...arguments);this.bubbles=p(this,["bubblesLocation"],["index"]);this.refNames=new u(this.bubbles,this)}static{this.capabilities=["getFeatures","getRefNames"]}async getRefNames(t={}){return this.refNames.assemblyRefNames(t)}async getRegionByteSize(t,s={}){return this.refNames.regionBytes(t,s)}async hasDataForRefName(){return!0}getFeatures(t,s={}){let{signal:n,statusCallback:d}=s;return g(async r=>{let o=await this.refNames.resolve(t,s);o!==void 0&&await l("Downloading bubbles",d,()=>m(this.bubbles,o,t.start,t.end,{signal:n,lineCallback:(N,y)=>{let e=f(N);r.next(new i({uniqueId:`bubble-${y}`,refName:t.refName,start:e.start,end:e.end,name:c(e),description:h(e),type:"bubble",score:e.segmentCount,segmentCount:e.segmentCount,pathCount:e.pathCount,inversion:e.inversion,shortestAlleleLength:e.shortestAlleleLength,longestAlleleLength:e.longestAlleleLength,segments:e.segments,shortestAllele:e.shortestAllele,longestAllele:e.longestAllele}))}})),r.complete()},n)}};export{a as default};
2
+ //# sourceMappingURL=MinigraphBubbleAdapter-HW7XQMAI.js.map
@@ -0,0 +1,7 @@
1
+ {
2
+ "version": 3,
3
+ "sources": ["../../src/MinigraphBubbleAdapter/MinigraphBubbleAdapter.ts"],
4
+ "sourcesContent": ["import {\n bubbleDescription,\n bubbleLabel,\n parseBubbleLine,\n} from '@jbrowse/bandage-core/bubbles/bubbleLine'\nimport { BaseFeatureDataAdapter } from '@jbrowse/core/data_adapters/BaseAdapter'\nimport { SimpleFeature, updateStatus } from '@jbrowse/core/util'\nimport { ObservableCreate } from '@jbrowse/core/util/rxjs'\n\nimport { PanSNRefNames, openTabixSlot } from '../panSNTabix.ts'\nimport { getLines } from '../tabixRange.ts'\n\nimport type { MinigraphBubbleAdapterConfig } from './configSchema.ts'\nimport type { BaseOptions } from '@jbrowse/core/data_adapters/BaseAdapter'\nimport type { Feature } from '@jbrowse/core/util'\nimport type { Region } from '@jbrowse/core/util/types'\n\nexport default class MinigraphBubbleAdapter extends BaseFeatureDataAdapter<MinigraphBubbleAdapterConfig> {\n public static capabilities = ['getFeatures', 'getRefNames']\n\n private readonly bubbles = openTabixSlot(this, ['bubblesLocation'], ['index'])\n // `gfatools bubble` names each row after the graph's stable sequence, so a\n // Minigraph-Cactus graph produces PanSN rows (`GRCh38#0#chr6`) where a plain\n // minigraph graph produces bare ones (`chr6`). Same resolution the segment\n // adapter does, through the same `assemblyNameToPanSN` slot.\n private readonly refNames = new PanSNRefNames(this.bubbles, this)\n\n async getRefNames(opts: BaseOptions = {}) {\n return this.refNames.assemblyRefNames(opts)\n }\n\n async getRegionByteSize(regions: Region[], opts: BaseOptions = {}) {\n return this.refNames.regionBytes(regions, opts)\n }\n\n public async hasDataForRefName() {\n return true\n }\n\n getFeatures(query: Region, opts: BaseOptions = {}) {\n const { signal, statusCallback } = opts\n return ObservableCreate<Feature>(async observer => {\n const tabixRefName = await this.refNames.resolve(query, opts)\n if (tabixRefName !== undefined) {\n await updateStatus('Downloading bubbles', statusCallback, () =>\n getLines(this.bubbles, tabixRefName, query.start, query.end, {\n signal,\n lineCallback: (line, fileOffset) => {\n const bubble = parseBubbleLine(line)\n observer.next(\n new SimpleFeature({\n uniqueId: `bubble-${fileOffset}`,\n refName: query.refName,\n start: bubble.start,\n end: bubble.end,\n // How variable this spot is, in two words, since that is the\n // whole point of the track.\n name: bubbleLabel(bubble),\n description: bubbleDescription(bubble),\n type: 'bubble',\n score: bubble.segmentCount,\n segmentCount: bubble.segmentCount,\n pathCount: bubble.pathCount,\n inversion: bubble.inversion,\n shortestAlleleLength: bubble.shortestAlleleLength,\n longestAlleleLength: bubble.longestAlleleLength,\n segments: bubble.segments,\n shortestAllele: bubble.shortestAllele,\n longestAllele: bubble.longestAllele,\n }),\n )\n },\n }),\n )\n }\n observer.complete()\n }, signal)\n }\n}\n"],
5
+ "mappings": "iVAiBA,IAAqBA,EAArB,cAAoDC,CAAqD,CAAzG,kCAGE,KAAiB,QAAUC,EAAc,KAAM,CAAC,iBAAiB,EAAG,CAAC,OAAO,CAAC,EAK7E,KAAiB,SAAW,IAAIC,EAAc,KAAK,QAAS,IAAI,EAPhE,YAAc,aAAe,CAAC,cAAe,aAAa,EAS1D,MAAM,YAAYC,EAAoB,CAAC,EAAG,CACxC,OAAO,KAAK,SAAS,iBAAiBA,CAAI,CAC5C,CAEA,MAAM,kBAAkBC,EAAmBD,EAAoB,CAAC,EAAG,CACjE,OAAO,KAAK,SAAS,YAAYC,EAASD,CAAI,CAChD,CAEA,MAAa,mBAAoB,CAC/B,MAAO,EACT,CAEA,YAAYE,EAAeF,EAAoB,CAAC,EAAG,CACjD,GAAM,CAAE,OAAAG,EAAQ,eAAAC,CAAe,EAAIJ,EACnC,OAAOK,EAA0B,MAAMC,GAAY,CACjD,IAAMC,EAAe,MAAM,KAAK,SAAS,QAAQL,EAAOF,CAAI,EACxDO,IAAiB,QACnB,MAAMC,EAAa,sBAAuBJ,EAAgB,IACxDK,EAAS,KAAK,QAASF,EAAcL,EAAM,MAAOA,EAAM,IAAK,CAC3D,OAAAC,EACA,aAAc,CAACO,EAAMC,IAAe,CAClC,IAAMC,EAASC,EAAgBH,CAAI,EACnCJ,EAAS,KACP,IAAIQ,EAAc,CAChB,SAAU,UAAUH,CAAU,GAC9B,QAAST,EAAM,QACf,MAAOU,EAAO,MACd,IAAKA,EAAO,IAGZ,KAAMG,EAAYH,CAAM,EACxB,YAAaI,EAAkBJ,CAAM,EACrC,KAAM,SACN,MAAOA,EAAO,aACd,aAAcA,EAAO,aACrB,UAAWA,EAAO,UAClB,UAAWA,EAAO,UAClB,qBAAsBA,EAAO,qBAC7B,oBAAqBA,EAAO,oBAC5B,SAAUA,EAAO,SACjB,eAAgBA,EAAO,eACvB,cAAeA,EAAO,aACxB,CAAC,CACH,CACF,CACF,CAAC,CACH,EAEFN,EAAS,SAAS,CACpB,EAAGH,CAAM,CACX,CACF",
6
+ "names": ["MinigraphBubbleAdapter", "BaseFeatureDataAdapter", "openTabixSlot", "PanSNRefNames", "opts", "regions", "query", "signal", "statusCallback", "ObservableCreate", "observer", "tabixRefName", "updateStatus", "getLines", "line", "fileOffset", "bubble", "parseBubbleLine", "SimpleFeature", "bubbleLabel", "bubbleDescription"]
7
+ }
@@ -0,0 +1,3 @@
1
+ import{a as F}from"./chunk-MND7GAOG.js";import{a as w}from"./chunk-KOB7YIMR.js";import"./chunk-FVK2THE4.js";import{b as L,h as d,m as p,n as y}from"./chunk-FHZXW7IG.js";import"./chunk-YBLNN5FU.js";import"./chunk-MVVXRZXV.js";import"./chunk-ZMOMFU2B.js";import{b as x}from"./chunk-5OI6LO4I.js";import"./chunk-OGQ6OQDN.js";import"./chunk-NDRLCGPC.js";import{a as h,o as k}from"./chunk-TENPGIPX.js";function u(s){let n=s.split(" "),t=n.length>=8,a=t?5:0;return{anchored:t,segment:{refName:n[a],start:+n[a+1],end:+n[a+2],id:n[3],rank:+n[4],tags:n[t?8:5]??""}}}function C(s){let n=s.tags.split(" ").find(a=>a.startsWith("SM:Z:"));if(n===void 0)return;let t=n.slice(5).split(",").map(a=>a.trim()).filter(Boolean);return t.length>0?t:void 0}function I(s){return{id:s.slice(0,-1),strand:s.slice(-1)}}function v(s){let n=s.split(" "),t=I(n[3]),a=I(n[4]);return{source:t.id,sourceStrand:t.strand,target:a.id,targetStrand:a.strand,sourceSegment:{id:t.id,refName:n[5],start:+n[6],end:+n[7],rank:+n[8],tags:n[13]??""},targetSegment:{id:a.id,refName:n[9],start:+n[10],end:+n[11],rank:+n[12],tags:n[14]??""}}}function M(s){return`${s.source}${s.sourceStrand}${s.target}${s.targetStrand}`}function S(s,n,t){n.set(M(t),t);let a=[];for(let e of[t.sourceSegment,t.targetSegment])s.has(e.id)||(s.set(e.id,e),a.push(e));return a}function B(s,n,t){let a=i=>i.rank!==0||i.refName!==t.refName?0:i.end<=t.start?-1:i.start>=t.end?1:0,e=new Map(s),r=new Map;for(let i of n)(s.has(i.source)||s.has(i.target)||a(i.sourceSegment)*a(i.targetSegment)===-1)&&S(e,r,i);for(let i of n)e.has(i.source)&&e.has(i.target)&&r.set(M(i),i);return{segments:e,links:r}}function O(s,n,t){let a=e=>e.refName!==t.refName||e.end<=t.start||e.start>=t.end;return[...$(s.filter(e=>e.rank>0),e=>n.has(e.id),t.end-t.start),...$(s.filter(e=>e.rank===0&&a(e)),()=>!0,0)]}function $(s,n,t){let a=[...s].sort((r,i)=>l(r.refName,i.refName)||r.start-i.start),e=[];for(let r of a){let i=e.at(-1);i?.refName===r.refName&&r.start-i.end<=t?(i.needed||=r.start>i.end||n(r),i.end=Math.max(i.end,r.end)):e.push({refName:r.refName,start:r.start,end:r.end,needed:n(r)})}return e.filter(r=>r.needed).map(({refName:r,start:i,end:o})=>({refName:r,start:i,end:o}))}function T(s,n,t){let a=e=>n.has(e.id)||t.some(r=>r.refName===e.refName&&e.start>=r.start&&e.end<=r.end);return s.filter(e=>a(e.sourceSegment)&&a(e.targetSegment))}var P=/^[A-Za-z][A-Za-z0-9]:[AifZJHB]:/;function Z(s){let{id:n,refName:t,start:a,end:e,rank:r,tags:i}=s,o=i.split(" ").filter(f=>P.test(f)),g=o.length===0?"":` ${o.join(" ")}`;return`S ${n} * LN:i:${e-a} SN:Z:${t} SO:i:${a} SR:i:${r}${g}`}function W(s){let{source:n,sourceStrand:t,target:a,targetStrand:e}=s;return`L ${n} ${t} ${a} ${e} 0M`}function l(s,n){return s<n?-1:s>n?1:0}function N(s,n){let t=[...s.values()].sort((e,r)=>l(e.refName,r.refName)||e.start-r.start||l(e.id,r.id)),a=[...n.keys()].sort((e,r)=>l(e,r));return["H VN:Z:1.0",...t.map(e=>Z(e)),...a.map(e=>W(n.get(e)))].join(`
2
+ `)}function A(s){return s.filter(n=>n.rank>0)}async function R(s,n,t){let a=[];return await d(s,n.refName,n.start,n.end,{signal:t,lineCallback:e=>{a.push(v(e))}}),a}async function G(s,n,{hops:t=0,signal:a}={}){let e=new Map,r=new Map,i={anchored:!1};if(await d(s.segments,n.refName,n.start,n.end,{signal:a,lineCallback:f=>{let m=u(f);i.anchored||=m.anchored,e.set(m.segment.id,m.segment)}}),i.anchored){let f=B(e,await R(s.links,n,a),n);return N(f.segments,f.links)}let o=async f=>(await R(s.links,f,a)).flatMap(m=>S(e,r,m)),g=A(await o(n));for(let f=0;f<t;f++){let m=await Promise.all(g.map(o));g=A(m.flat())}if(t>0){let f=O([...e.values()],new Set(g.map(c=>c.id)),n),m=await Promise.all(f.map(c=>R(s.links,c,a)));for(let c of T(m.flat(),e,f))S(e,r,c)}return N(e,r)}var b=class extends L{constructor(){super(...arguments);this.fine=this.openIndex([])}static{this.capabilities=["getFeatures","getRefNames"]}openIndex(t){let a=p(this,[...t,"segmentsLocation"],[...t,"segmentsIndex"]),e=p(this,[...t,"linksLocation"],[...t,"linksIndex"]),r=this.getConf("walksLocation"),i=t.length===0&&x(r)?new F(p(this,["walksLocation"],["walksIndex"]),this):void 0;return{segments:a,links:e,walks:i,refNames:new y(i?.file??a,this)}}async getHaplotypeNames(t={}){return this.fine.walks?.haplotypeNames(t)}index(t="fine"){return t==="fine"?this.fine:(this.coarse??=this.openIndex(["coarse"]),this.coarse)}async getRefNames(t={}){return this.fine.refNames.assemblyRefNames(t)}async getRegionByteSize(t,a={}){let{walks:e,refNames:r}=this.fine;return r.regionBytes(t,a,e&&(i=>e.reader.queryStart(i,a)))}async hasDataForRefName(){return!0}getFeatures(t,a={}){return w(async e=>{let r=await this.fine.refNames.resolve(t,a);if(r!==void 0){let{walks:i}=this.fine,o=g=>{e.next(g)};await(i?this.walkFeatures(i,r,t,a,o):this.segmentFeatures(r,t,a,o))}e.complete()},a.signal)}async walkFeatures(t,a,e,r,i){let o=await k("Downloading walks",r.statusCallback,()=>t.referenceNodes(a,e,r));for(let g of o){let f=String(g.id);i(new h({uniqueId:f,refName:e.refName,start:g.start,end:g.end,name:f,type:"segment",rank:0,stableName:g.refName}))}}async segmentFeatures(t,a,e,r){await k("Downloading segments",e.statusCallback,()=>d(this.fine.segments,t,a.start,a.end,{signal:e.signal,lineCallback:i=>{let{segment:o,anchored:g}=u(i);if(g&&o.refName!==t)return;let f=C(o);r(new h({uniqueId:o.id,refName:a.refName,start:o.start,end:o.end,name:o.id,type:"segment",rank:o.rank,stableName:o.refName,...f&&{samples:f,sampleCount:f.length}}))}}))}async getSubgraph(t,a={}){let{hops:e,signal:r,tier:i}=a,o=this.index(i),g=await o.refNames.resolve(t,{signal:r});if(g===void 0)throw new Error(`${t.assemblyName} ${t.refName} is not in this graph's index; a graph with PanSN names (GRCh38#0#chr1) needs ${t.assemblyName} mapped to its prefix in assemblyNameToPanSN`);return o.walks?o.walks.subgraph(g,t,a):G(o,{refName:g,start:t.start,end:t.end},{hops:e,signal:r})}};export{b as default};
3
+ //# sourceMappingURL=RgfaTabixAdapter-PSKYWD2X.js.map
@@ -0,0 +1,7 @@
1
+ {
2
+ "version": 3,
3
+ "sources": ["../../src/RgfaTabixAdapter/rgfaBed.ts", "../../src/RgfaTabixAdapter/rgfaCut.ts", "../../src/RgfaTabixAdapter/RgfaTabixAdapter.ts"],
4
+ "sourcesContent": ["// A segment as it appears in segs.bed.gz, and as it is repeated inside every\n// links.bed.gz row: `stableName start end segmentId rank [tags]`. The first\n// five are the SN/SO/SR tags rGFA already carries, projected to BED by\n// `gfatools gfa2bed -m`.\n//\n// `tags` is a space-separated list of GFA tags, written verbatim onto the\n// S-line formatSegment synthesizes, where the GFA parser reads them into\n// `GraphNode.tags` with no work here. It is the extension point for everything\n// the five columns cannot say: carriage on a path-derived graph (`SM:Z:`), the\n// collapse summary on a level-of-detail tier (`ct:Z:`, `cn:i:`, ...), a\n// precomputed layout position later. rGFA files have no sixth column at all,\n// so this is empty for them and nothing changes.\nexport interface RgfaSegment {\n refName: string\n start: number\n end: number\n id: string\n rank: number\n tags: string\n}\n\nexport interface RgfaLink {\n source: string\n sourceStrand: string\n target: string\n targetStrand: string\n sourceSegment: RgfaSegment\n targetSegment: RgfaSegment\n}\n\nexport function parseSegmentLine(line: string): RgfaSegment {\n return parseSegmentRow(line).segment\n}\n\n// gfa-to-tabix's anchored layout files a segment under the reference interval\n// its bubble hangs from and states the segment's own coordinate after the rank:\n// `anchorName anchorStart anchorEnd segmentId rank stableName start end [tags]`.\n// One query over a region then returns every segment of every bubble under it,\n// and `anchored` tells a cut it has nothing left to follow.\nexport function parseSegmentRow(line: string) {\n const cols = line.split('\\t')\n const anchored = cols.length >= 8\n const own = anchored ? 5 : 0\n return {\n anchored,\n segment: {\n refName: cols[own]!,\n start: +cols[own + 1]!,\n end: +cols[own + 2]!,\n id: cols[3]!,\n rank: +cols[4]!,\n tags: cols[anchored ? 8 : 5] ?? '',\n } satisfies RgfaSegment,\n }\n}\n\n// Carriage as an ordinary feature attribute, for the linear track rather than\n// the graph. The tag column reaches the graph view through the synthesized\n// S-line (formatSegment below, then gfaConverter reads SM into\n// `GraphNode.samples`), and that route ends at the graph: `getFeatures` builds\n// its own features and had no way to say who carries a segment, so a lane\n// colored by carriage was not expressible.\n//\n// `SM:Z:` is its own grammar check here. build_pggb_tabix.sh wrote a bare\n// comma-separated list in this column before the tag column existed, and those\n// files are still hosted; a bare list has no prefix to match, so it is dropped\n// rather than read as one sample named `K12,Sakai,CFT073`.\n//\n// Undefined rather than [] when the tag is absent, matching gfaConverter: an\n// rGFA has no sixth column at all, so a segment that says nothing about\n// carriage stays distinguishable from one carried by nobody.\nexport function segmentSamples(segment: RgfaSegment) {\n const tag = segment.tags.split(' ').find(t => t.startsWith('SM:Z:'))\n if (tag === undefined) {\n return undefined\n }\n const samples = tag\n .slice('SM:Z:'.length)\n .split(',')\n .map(s => s.trim())\n .filter(Boolean)\n return samples.length > 0 ? samples : undefined\n}\n\n// `s322+` \u2014 the L-line's segment id with its orientation appended, mirroring\n// the GFA record it came from.\nfunction parseEndpoint(field: string) {\n return { id: field.slice(0, -1), strand: field.slice(-1) }\n}\n\nexport function parseLinkLine(line: string): RgfaLink {\n const cols = line.split('\\t')\n const source = parseEndpoint(cols[3]!)\n const target = parseEndpoint(cols[4]!)\n return {\n source: source.id,\n sourceStrand: source.strand,\n target: target.id,\n targetStrand: target.strand,\n sourceSegment: {\n id: source.id,\n refName: cols[5]!,\n start: +cols[6]!,\n end: +cols[7]!,\n rank: +cols[8]!,\n tags: cols[13] ?? '',\n },\n targetSegment: {\n id: target.id,\n refName: cols[9]!,\n start: +cols[10]!,\n end: +cols[11]!,\n rank: +cols[12]!,\n tags: cols[14] ?? '',\n },\n }\n}\n\nexport function linkKey(link: RgfaLink) {\n return `${link.source}${link.sourceStrand}${link.target}${link.targetStrand}`\n}\n\n// gfa-to-tabix writes both endpoints of a link into every row and files the\n// link under each endpoint's stable sequence, so a read covering both returns\n// it twice; linkKey dedupes it. Returns the endpoints the cut lacked.\nexport function addLink(\n segments: Map<string, RgfaSegment>,\n links: Map<string, RgfaLink>,\n link: RgfaLink,\n) {\n links.set(linkKey(link), link)\n const added: RgfaSegment[] = []\n for (const segment of [link.sourceSegment, link.targetSegment]) {\n if (!segments.has(segment.id)) {\n segments.set(segment.id, segment)\n added.push(segment)\n }\n }\n return added\n}\n\n// The cut an anchored index gives: `under` is what the segment file returned\n// for the window, every segment of every bubble under it, and `found` is what\n// the link file returned. gfa-to-tabix files a link widely enough that `found`\n// holds every link touching a segment in `under` and every link between the\n// segments those lead to, plus links of neither kind, which are dropped here.\n// A backbone link that jumps the whole window, as a deletion or inversion\n// spanning it does, touches nothing in `under` and is kept by its own rule.\nexport function anchoredCut(\n under: Map<string, RgfaSegment>,\n found: RgfaLink[],\n window: StableSpan,\n) {\n const side = (segment: RgfaSegment) =>\n segment.rank !== 0 || segment.refName !== window.refName\n ? 0\n : segment.end <= window.start\n ? -1\n : segment.start >= window.end\n ? 1\n : 0\n const segments = new Map(under)\n const links = new Map<string, RgfaLink>()\n for (const link of found) {\n if (\n under.has(link.source) ||\n under.has(link.target) ||\n side(link.sourceSegment) * side(link.targetSegment) === -1\n ) {\n addLink(segments, links, link)\n }\n }\n for (const link of found) {\n if (segments.has(link.source) && segments.has(link.target)) {\n links.set(linkKey(link), link)\n }\n }\n return { segments, links }\n}\n\nexport interface StableSpan {\n refName: string\n start: number\n end: number\n}\n\n// Where a cut reads links once more to close itself. A hop reads the links of\n// the segments the round before it reached, so a link between two segments of\n// the last round is never read and an allele reached from both ends draws as\n// two stubs. rGFA lays an allele out as one stretch of the stable sequence that\n// introduced it, so pieces of one sequence within the window's length of each\n// other merge into a span, and reading it brings the pieces between too.\n// Backbone reached past the window is unread for the same reason; its spans\n// join only abutting pieces, since the reference past the window is not the\n// cut's to fill. A span whose segments were all read and abut is skipped.\nexport function closingSpans(\n held: RgfaSegment[],\n unread: Set<string>,\n window: StableSpan,\n) {\n const outside = (segment: RgfaSegment) =>\n segment.refName !== window.refName ||\n segment.end <= window.start ||\n segment.start >= window.end\n return [\n ...mergedSpans(\n held.filter(segment => segment.rank > 0),\n segment => unread.has(segment.id),\n window.end - window.start,\n ),\n ...mergedSpans(\n held.filter(segment => segment.rank === 0 && outside(segment)),\n () => true,\n 0,\n ),\n ]\n}\n\nfunction mergedSpans(\n segments: RgfaSegment[],\n isUnread: (segment: RgfaSegment) => boolean,\n maxGap: number,\n) {\n const sorted = [...segments].sort(\n (a, b) => compare(a.refName, b.refName) || a.start - b.start,\n )\n const spans: (StableSpan & { needed: boolean })[] = []\n for (const segment of sorted) {\n const span = spans.at(-1)\n if (\n span?.refName === segment.refName &&\n segment.start - span.end <= maxGap\n ) {\n span.needed ||= segment.start > span.end || isUnread(segment)\n span.end = Math.max(span.end, segment.end)\n } else {\n spans.push({\n refName: segment.refName,\n start: segment.start,\n end: segment.end,\n needed: isUnread(segment),\n })\n }\n }\n return spans\n .filter(span => span.needed)\n .map(({ refName, start, end }): StableSpan => ({ refName, start, end }))\n}\n\n// The links a closing read keeps: each end is a segment the cut holds or one\n// lying inside a span. Following anything else is what another hop is for.\nexport function closingLinks(\n found: RgfaLink[],\n held: Map<string, RgfaSegment>,\n spans: StableSpan[],\n) {\n const inside = (segment: RgfaSegment) =>\n held.has(segment.id) ||\n spans.some(\n span =>\n span.refName === segment.refName &&\n segment.start >= span.start &&\n segment.end <= span.end,\n )\n return found.filter(\n link => inside(link.sourceSegment) && inside(link.targetSegment),\n )\n}\n\n// Segments carry no sequence here \u2014 the BED records only their span \u2014 so every\n// S-line is written with `*` and an LN tag, which is what the GFA spec asks for\n// and what packages/graph-core's parser reads back as the node length.\n// `NM:i:0`, `SM:Z:HG002.1,HG002.2` \u2014 the GFA tag grammar, checked rather than\n// trusted. Files built before the tag column existed put a bare comma-separated\n// sample list here, and passing that through would put a field that is not a\n// tag on an S-line, which is a malformed GFA rather than a missing annotation.\n// Dropping it degrades to the pre-tag behaviour instead.\nconst GFA_TAG = /^[A-Za-z][A-Za-z0-9]:[AifZJHB]:/\n\nfunction formatSegment(segment: RgfaSegment) {\n const { id, refName, start, end, rank, tags } = segment\n const valid = tags.split(' ').filter(t => GFA_TAG.test(t))\n const extra = valid.length === 0 ? '' : `\\t${valid.join('\\t')}`\n return `S\\t${id}\\t*\\tLN:i:${end - start}\\tSN:Z:${refName}\\tSO:i:${start}\\tSR:i:${rank}${extra}`\n}\n\nfunction formatLink(link: RgfaLink) {\n const { source, sourceStrand, target, targetStrand } = link\n return `L\\t${source}\\t${sourceStrand}\\t${target}\\t${targetStrand}\\t0M`\n}\n\n// Codepoint order, not localeCompare: the byte-identical-output guarantee below\n// has to hold across locales.\nfunction compare(a: string, b: string) {\n return a < b ? -1 : a > b ? 1 : 0\n}\n\n// Deterministic ordering: the same region always produces byte-identical GFA,\n// so a layout computed from it is reproducible.\nexport function formatSubgraph(\n segments: Map<string, RgfaSegment>,\n links: Map<string, RgfaLink>,\n) {\n const sortedSegments = [...segments.values()].sort(\n (a, b) =>\n compare(a.refName, b.refName) || a.start - b.start || compare(a.id, b.id),\n )\n const sortedLinks = [...links.keys()].sort((a, b) => compare(a, b))\n return [\n 'H\\tVN:Z:1.0',\n ...sortedSegments.map(s => formatSegment(s)),\n ...sortedLinks.map(k => formatLink(links.get(k)!)),\n ].join('\\n')\n}\n", "import {\n addLink,\n anchoredCut,\n closingLinks,\n closingSpans,\n formatSubgraph,\n parseLinkLine,\n parseSegmentRow,\n} from './rgfaBed.ts'\nimport { getLines } from '../tabixRange.ts'\n\nimport type { RgfaLink, RgfaSegment, StableSpan } from './rgfaBed.ts'\nimport type { TabixIndexedFile } from '@gmod/tabix'\n\n// What a hop follows: alleles, never the backbone. A rank-0 segment reached\n// through a link is flanking backbone, usually outside the window, so querying\n// its links walks the reference away from the region, brings back more backbone,\n// and closes nothing. Every segment that closes a bubble is an allele interior,\n// and an allele interior is rank > 0 in both index formats (rGFA's SR tag, and\n// pggb's rank 0-or-1 against the anchor path). Reached backbone still becomes a\n// node and its edge still draws; it is only not walked through.\n//\n// This is gfabase's rule, arrived at from the other end. `gfabase sub` expands\n// from a seed either to the whole connected component (`--connected`, which its\n// own README calls overkill) or with `--cutpoints 1`, which takes the biconnected\n// component by *stopping at cutpoints*, the segments every end-to-end walk of\n// the chromosome must traverse. In a reference-anchored pangenome the cutpoints\n// are the backbone, and rank is what the index already states about it, so\n// `rank > 0` is that stopping rule without a connectivity index to build. The\n// third tool is the same shape again: `gfatools view -R` maps the region onto its\n// precomputed bubble decomposition and returns whole overlapping bubbles\n// (gfa_query_by_reg), which is complete for the same reason. Bandage's own\n// `--distance` is a node-step radius, but Bandage holds the entire graph in\n// memory, so its scope is a drawing filter rather than a fetch strategy.\n//\n// Measured with this rule: the E. coli paa window converges at one hop (17\n// segments at 1 and at 2), so the cut closes and stays closed. HPRC's amylase\n// window keeps growing (63 at 0, 78 at 1, 92 at 2) because its alleles have\n// alleles, and wall-clock is flat across all three, since a hop's queries go out\n// together and the remote index dominates.\n//\n// The exact version, when a graph ships a bubble index beside its segments:\n// gfatools' bubble rows carry the member segment ids of every bubble, so the\n// bubbles overlapping the region state which segments a complete cut must\n// contain. That turns the hop count into a termination condition rather than a\n// setting, and MinigraphBubbleAdapter already parses those rows.\nfunction offReference(segments: RgfaSegment[]) {\n return segments.filter(segment => segment.rank > 0)\n}\n\nexport interface RgfaFiles {\n segments: TabixIndexedFile\n links: TabixIndexedFile\n}\n\nasync function linksOver(\n file: TabixIndexedFile,\n span: StableSpan,\n signal?: AbortSignal,\n) {\n const found: RgfaLink[] = []\n await getLines(file, span.refName, span.start, span.end, {\n signal,\n lineCallback: line => {\n found.push(parseLinkLine(line))\n },\n })\n return found\n}\n\n/**\n * The graph around a window of one stable sequence, as GFA text. The window's\n * own segments come from the segment file; the link file supplies the edges\n * incident to them *and* the coordinates of the segments on the other end,\n * which typically sit on a different stable sequence (a rank>0 bubble) and so\n * are not reachable by any coordinate query on this window.\n *\n * `hops` is extra rounds of link-following past the window's own segments and\n * their immediate neighbours, each costing one tabix query per off-reference\n * segment newly reached. After the last round, one more read per stretch of\n * stable sequence the cut holds closes it (closingSpans), so an allele reached\n * from both ends draws whole. The view asks for one hop; a hopless cut stays\n * as cheap as it was. `signal` goes to every one of those queries, so a cut\n * the view has replaced stops.\n */\nexport async function rgfaCut(\n files: RgfaFiles,\n window: StableSpan,\n { hops = 0, signal }: { hops?: number; signal?: AbortSignal } = {},\n) {\n const segments = new Map<string, RgfaSegment>()\n const links = new Map<string, RgfaLink>()\n const layout = { anchored: false }\n await getLines(files.segments, window.refName, window.start, window.end, {\n signal,\n lineCallback: line => {\n const row = parseSegmentRow(line)\n layout.anchored ||= row.anchored\n segments.set(row.segment.id, row.segment)\n },\n })\n // An anchored index files every allele and every link under the reference\n // interval its bubble hangs from, so one read of each file holds the whole\n // graph under the window and there is nothing to hop to or close.\n if (layout.anchored) {\n const cut = anchoredCut(\n segments,\n await linksOver(files.links, window, signal),\n window,\n )\n return formatSubgraph(cut.segments, cut.links)\n }\n const addLinksOver = async (span: StableSpan) =>\n (await linksOver(files.links, span, signal)).flatMap(link =>\n addLink(segments, links, link),\n )\n let frontier = offReference(await addLinksOver(window))\n for (let hop = 0; hop < hops; hop++) {\n // One hop's queries are independent of each other, so they go out\n // together rather than one round-trip at a time, and the output is sorted\n // at the end regardless of the order they arrive in.\n const reached = await Promise.all(frontier.map(addLinksOver))\n frontier = offReference(reached.flat())\n }\n if (hops > 0) {\n const spans = closingSpans(\n [...segments.values()],\n new Set(frontier.map(segment => segment.id)),\n window,\n )\n const found = await Promise.all(\n spans.map(span => linksOver(files.links, span, signal)),\n )\n for (const link of closingLinks(found.flat(), segments, spans)) {\n addLink(segments, links, link)\n }\n }\n return formatSubgraph(segments, links)\n}\n", "import { BaseFeatureDataAdapter } from '@jbrowse/core/data_adapters/BaseAdapter'\nimport { SimpleFeature, updateStatus } from '@jbrowse/core/util'\nimport { ObservableCreate } from '@jbrowse/core/util/rxjs'\n\nimport { locationIsSet } from '../locationName.ts'\nimport { PanSNRefNames, openTabixSlot } from '../panSNTabix.ts'\nimport { getLines } from '../tabixRange.ts'\nimport { parseSegmentRow, segmentSamples } from './rgfaBed.ts'\nimport { rgfaCut } from './rgfaCut.ts'\nimport { WalkIndex } from './walkReader.ts'\n\nimport type { RgfaTabixAdapterConfig } from './configSchema.ts'\nimport type { SubgraphAdapterOptions, SubgraphTier } from '../GetSubgraph.ts'\nimport type { WalkCut } from './walkRowRuns.ts'\nimport type { TabixIndexedFile } from '@gmod/tabix'\nimport type { BaseOptions } from '@jbrowse/core/data_adapters/BaseAdapter'\nimport type { Feature } from '@jbrowse/core/util'\nimport type { FileLocation, Region } from '@jbrowse/core/util/types'\n\ninterface GraphIndex {\n segments: TabixIndexedFile\n links: TabixIndexedFile\n walks?: WalkIndex\n refNames: PanSNRefNames\n}\n\nexport default class RgfaTabixAdapter extends BaseFeatureDataAdapter<RgfaTabixAdapterConfig> {\n public static capabilities = ['getFeatures', 'getRefNames']\n\n private readonly fine = this.openIndex([])\n private coarse: GraphIndex | undefined\n\n // refNames off the SEGMENT index: it is the one that names every stable\n // sequence the graph is anchored to, where the link index only names those a\n // link touches. A walk-indexed graph's walk file names the same sequences,\n // and it has no segment file.\n private openIndex(under: string[]): GraphIndex {\n const segments = openTabixSlot(\n this,\n [...under, 'segmentsLocation'],\n [...under, 'segmentsIndex'],\n )\n const links = openTabixSlot(\n this,\n [...under, 'linksLocation'],\n [...under, 'linksIndex'],\n )\n const walksLocation: FileLocation = this.getConf('walksLocation')\n const walks =\n under.length === 0 && locationIsSet(walksLocation)\n ? new WalkIndex(\n openTabixSlot(this, ['walksLocation'], ['walksIndex']),\n this,\n )\n : undefined\n return {\n segments,\n links,\n walks,\n refNames: new PanSNRefNames(walks?.file ?? segments, this),\n }\n }\n\n /**\n * The haplotypes a walk-indexed graph's header names, `sample#haplotype`\n * in byte order, or undefined when it names none\n */\n async getHaplotypeNames(opts: BaseOptions = {}) {\n return this.fine.walks?.haplotypeNames(opts)\n }\n\n private index(tier: SubgraphTier = 'fine') {\n if (tier === 'fine') {\n return this.fine\n }\n this.coarse ??= this.openIndex(['coarse'])\n return this.coarse\n }\n\n async getRefNames(opts: BaseOptions = {}) {\n return this.fine.refNames.assemblyRefNames(opts)\n }\n\n // getFeatures reads a walk file from the chunk before the window\n async getRegionByteSize(regions: Region[], opts: BaseOptions = {}) {\n const { walks, refNames } = this.fine\n return refNames.regionBytes(\n regions,\n opts,\n walks && (start => walks.reader.queryStart(start, opts)),\n )\n }\n\n public async hasDataForRefName() {\n return true\n }\n\n getFeatures(query: Region, opts: BaseOptions = {}) {\n return ObservableCreate<Feature>(async observer => {\n const tabixRefName = await this.fine.refNames.resolve(query, opts)\n if (tabixRefName !== undefined) {\n const { walks } = this.fine\n const emit = (feature: Feature) => {\n observer.next(feature)\n }\n await (walks\n ? this.walkFeatures(walks, tabixRefName, query, opts, emit)\n : this.segmentFeatures(tabixRefName, query, opts, emit))\n }\n observer.complete()\n }, opts.signal)\n }\n\n // A walk-indexed graph's reference nodes are its ref rows' steps\n private async walkFeatures(\n walks: WalkIndex,\n tabixRefName: string,\n query: Region,\n opts: BaseOptions,\n emit: (feature: Feature) => void,\n ) {\n const nodes = await updateStatus(\n 'Downloading walks',\n opts.statusCallback,\n () => walks.referenceNodes(tabixRefName, query, opts),\n )\n for (const node of nodes) {\n const id = String(node.id)\n emit(\n new SimpleFeature({\n uniqueId: id,\n refName: query.refName,\n start: node.start,\n end: node.end,\n name: id,\n type: 'segment',\n rank: 0,\n stableName: node.refName,\n }),\n )\n }\n }\n\n private async segmentFeatures(\n tabixRefName: string,\n query: Region,\n opts: BaseOptions,\n emit: (feature: Feature) => void,\n ) {\n await updateStatus('Downloading segments', opts.statusCallback, () =>\n getLines(this.fine.segments, tabixRefName, query.start, query.end, {\n signal: opts.signal,\n lineCallback: line => {\n const { segment, anchored } = parseSegmentRow(line)\n // An anchored index also returns the alleles under the region,\n // which lie on other stable sequences and have no place on this one\n if (anchored && segment.refName !== tabixRefName) {\n return\n }\n // `samples` lists the haplotypes whose paths visit the segment and\n // `sampleCount` counts them. A lane's color reads the count, and\n // a count in jexl would otherwise rely on `.length` resolving\n // through a member access on an array.\n const samples = segmentSamples(segment)\n emit(\n new SimpleFeature({\n uniqueId: segment.id,\n refName: query.refName,\n start: segment.start,\n end: segment.end,\n name: segment.id,\n type: 'segment',\n rank: segment.rank,\n stableName: segment.refName,\n ...(samples && { samples, sampleCount: samples.length }),\n }),\n )\n },\n }),\n )\n }\n\n // The graph around a region for GraphGenomeView: GFA text (rgfaCut), or a\n // walk-indexed graph's tables (WalkIndex.subgraph). `tier: 'coarse'` reads\n // the `coarse` pair instead, the same way.\n async getSubgraph(\n region: Region,\n opts: SubgraphAdapterOptions = {},\n ): Promise<string | WalkCut> {\n const { hops, signal, tier } = opts\n const index = this.index(tier)\n const tabixRefName = await index.refNames.resolve(region, { signal })\n if (tabixRefName === undefined) {\n throw new Error(\n `${region.assemblyName} ${region.refName} is not in this graph's index; a graph with PanSN names (GRCh38#0#chr1) needs ${region.assemblyName} mapped to its prefix in assemblyNameToPanSN`,\n )\n }\n if (index.walks) {\n return index.walks.subgraph(tabixRefName, region, opts)\n }\n return rgfaCut(\n index,\n { refName: tabixRefName, start: region.start, end: region.end },\n { hops, signal },\n )\n }\n}\n"],
5
+ "mappings": "4YAuCO,SAASA,EAAgBC,EAAc,CAC5C,IAAMC,EAAOD,EAAK,MAAM,GAAI,EACtBE,EAAWD,EAAK,QAAU,EAC1BE,EAAMD,EAAW,EAAI,EAC3B,MAAO,CACL,SAAAA,EACA,QAAS,CACP,QAASD,EAAKE,CAAG,EACjB,MAAO,CAACF,EAAKE,EAAM,CAAC,EACpB,IAAK,CAACF,EAAKE,EAAM,CAAC,EAClB,GAAIF,EAAK,CAAC,EACV,KAAM,CAACA,EAAK,CAAC,EACb,KAAMA,EAAKC,EAAW,EAAI,CAAC,GAAK,EAClC,CACF,CACF,CAiBO,SAASE,EAAeC,EAAsB,CACnD,IAAMC,EAAMD,EAAQ,KAAK,MAAM,GAAG,EAAE,KAAKE,GAAKA,EAAE,WAAW,OAAO,CAAC,EACnE,GAAID,IAAQ,OACV,OAEF,IAAME,EAAUF,EACb,MAAM,CAAc,EACpB,MAAM,GAAG,EACT,IAAIG,GAAKA,EAAE,KAAK,CAAC,EACjB,OAAO,OAAO,EACjB,OAAOD,EAAQ,OAAS,EAAIA,EAAU,MACxC,CAIA,SAASE,EAAcC,EAAe,CACpC,MAAO,CAAE,GAAIA,EAAM,MAAM,EAAG,EAAE,EAAG,OAAQA,EAAM,MAAM,EAAE,CAAE,CAC3D,CAEO,SAASC,EAAcZ,EAAwB,CACpD,IAAMC,EAAOD,EAAK,MAAM,GAAI,EACtBa,EAASH,EAAcT,EAAK,CAAC,CAAE,EAC/Ba,EAASJ,EAAcT,EAAK,CAAC,CAAE,EACrC,MAAO,CACL,OAAQY,EAAO,GACf,aAAcA,EAAO,OACrB,OAAQC,EAAO,GACf,aAAcA,EAAO,OACrB,cAAe,CACb,GAAID,EAAO,GACX,QAASZ,EAAK,CAAC,EACf,MAAO,CAACA,EAAK,CAAC,EACd,IAAK,CAACA,EAAK,CAAC,EACZ,KAAM,CAACA,EAAK,CAAC,EACb,KAAMA,EAAK,EAAE,GAAK,EACpB,EACA,cAAe,CACb,GAAIa,EAAO,GACX,QAASb,EAAK,CAAC,EACf,MAAO,CAACA,EAAK,EAAE,EACf,IAAK,CAACA,EAAK,EAAE,EACb,KAAM,CAACA,EAAK,EAAE,EACd,KAAMA,EAAK,EAAE,GAAK,EACpB,CACF,CACF,CAEO,SAASc,EAAQC,EAAgB,CACtC,MAAO,GAAGA,EAAK,MAAM,GAAGA,EAAK,YAAY,GAAGA,EAAK,MAAM,GAAGA,EAAK,YAAY,EAC7E,CAKO,SAASC,EACdC,EACAC,EACAH,EACA,CACAG,EAAM,IAAIJ,EAAQC,CAAI,EAAGA,CAAI,EAC7B,IAAMI,EAAuB,CAAC,EAC9B,QAAWf,IAAW,CAACW,EAAK,cAAeA,EAAK,aAAa,EACtDE,EAAS,IAAIb,EAAQ,EAAE,IAC1Ba,EAAS,IAAIb,EAAQ,GAAIA,CAAO,EAChCe,EAAM,KAAKf,CAAO,GAGtB,OAAOe,CACT,CASO,SAASC,EACdC,EACAC,EACAC,EACA,CACA,IAAMC,EAAQpB,GACZA,EAAQ,OAAS,GAAKA,EAAQ,UAAYmB,EAAO,QAC7C,EACAnB,EAAQ,KAAOmB,EAAO,MACpB,GACAnB,EAAQ,OAASmB,EAAO,IACtB,EACA,EACJN,EAAW,IAAI,IAAII,CAAK,EACxBH,EAAQ,IAAI,IAClB,QAAWH,KAAQO,GAEfD,EAAM,IAAIN,EAAK,MAAM,GACrBM,EAAM,IAAIN,EAAK,MAAM,GACrBS,EAAKT,EAAK,aAAa,EAAIS,EAAKT,EAAK,aAAa,IAAM,KAExDC,EAAQC,EAAUC,EAAOH,CAAI,EAGjC,QAAWA,KAAQO,EACbL,EAAS,IAAIF,EAAK,MAAM,GAAKE,EAAS,IAAIF,EAAK,MAAM,GACvDG,EAAM,IAAIJ,EAAQC,CAAI,EAAGA,CAAI,EAGjC,MAAO,CAAE,SAAAE,EAAU,MAAAC,CAAM,CAC3B,CAiBO,SAASO,EACdC,EACAC,EACAJ,EACA,CACA,IAAMK,EAAWxB,GACfA,EAAQ,UAAYmB,EAAO,SAC3BnB,EAAQ,KAAOmB,EAAO,OACtBnB,EAAQ,OAASmB,EAAO,IAC1B,MAAO,CACL,GAAGM,EACDH,EAAK,OAAOtB,GAAWA,EAAQ,KAAO,CAAC,EACvCA,GAAWuB,EAAO,IAAIvB,EAAQ,EAAE,EAChCmB,EAAO,IAAMA,EAAO,KACtB,EACA,GAAGM,EACDH,EAAK,OAAOtB,GAAWA,EAAQ,OAAS,GAAKwB,EAAQxB,CAAO,CAAC,EAC7D,IAAM,GACN,CACF,CACF,CACF,CAEA,SAASyB,EACPZ,EACAa,EACAC,EACA,CACA,IAAMC,EAAS,CAAC,GAAGf,CAAQ,EAAE,KAC3B,CAACgB,EAAGC,IAAMC,EAAQF,EAAE,QAASC,EAAE,OAAO,GAAKD,EAAE,MAAQC,EAAE,KACzD,EACME,EAA8C,CAAC,EACrD,QAAWhC,KAAW4B,EAAQ,CAC5B,IAAMK,EAAOD,EAAM,GAAG,EAAE,EAEtBC,GAAM,UAAYjC,EAAQ,SAC1BA,EAAQ,MAAQiC,EAAK,KAAON,GAE5BM,EAAK,SAAWjC,EAAQ,MAAQiC,EAAK,KAAOP,EAAS1B,CAAO,EAC5DiC,EAAK,IAAM,KAAK,IAAIA,EAAK,IAAKjC,EAAQ,GAAG,GAEzCgC,EAAM,KAAK,CACT,QAAShC,EAAQ,QACjB,MAAOA,EAAQ,MACf,IAAKA,EAAQ,IACb,OAAQ0B,EAAS1B,CAAO,CAC1B,CAAC,CAEL,CACA,OAAOgC,EACJ,OAAOC,GAAQA,EAAK,MAAM,EAC1B,IAAI,CAAC,CAAE,QAAAC,EAAS,MAAAC,EAAO,IAAAC,CAAI,KAAmB,CAAE,QAAAF,EAAS,MAAAC,EAAO,IAAAC,CAAI,EAAE,CAC3E,CAIO,SAASC,EACdnB,EACAI,EACAU,EACA,CACA,IAAMM,EAAUtC,GACdsB,EAAK,IAAItB,EAAQ,EAAE,GACnBgC,EAAM,KACJC,GACEA,EAAK,UAAYjC,EAAQ,SACzBA,EAAQ,OAASiC,EAAK,OACtBjC,EAAQ,KAAOiC,EAAK,GACxB,EACF,OAAOf,EAAM,OACXP,GAAQ2B,EAAO3B,EAAK,aAAa,GAAK2B,EAAO3B,EAAK,aAAa,CACjE,CACF,CAUA,IAAM4B,EAAU,kCAEhB,SAASC,EAAcxC,EAAsB,CAC3C,GAAM,CAAE,GAAAyC,EAAI,QAAAP,EAAS,MAAAC,EAAO,IAAAC,EAAK,KAAAM,EAAM,KAAAC,CAAK,EAAI3C,EAC1C4C,EAAQD,EAAK,MAAM,GAAG,EAAE,OAAOzC,GAAKqC,EAAQ,KAAKrC,CAAC,CAAC,EACnD2C,EAAQD,EAAM,SAAW,EAAI,GAAK,IAAKA,EAAM,KAAK,GAAI,CAAC,GAC7D,MAAO,KAAMH,CAAE,WAAaL,EAAMD,CAAK,SAAUD,CAAO,SAAUC,CAAK,SAAUO,CAAI,GAAGG,CAAK,EAC/F,CAEA,SAASC,EAAWnC,EAAgB,CAClC,GAAM,CAAE,OAAAH,EAAQ,aAAAuC,EAAc,OAAAtC,EAAQ,aAAAuC,CAAa,EAAIrC,EACvD,MAAO,KAAMH,CAAM,IAAKuC,CAAY,IAAKtC,CAAM,IAAKuC,CAAY,KAClE,CAIA,SAASjB,EAAQF,EAAWC,EAAW,CACrC,OAAOD,EAAIC,EAAI,GAAKD,EAAIC,EAAI,EAAI,CAClC,CAIO,SAASmB,EACdpC,EACAC,EACA,CACA,IAAMoC,EAAiB,CAAC,GAAGrC,EAAS,OAAO,CAAC,EAAE,KAC5C,CAACgB,EAAGC,IACFC,EAAQF,EAAE,QAASC,EAAE,OAAO,GAAKD,EAAE,MAAQC,EAAE,OAASC,EAAQF,EAAE,GAAIC,EAAE,EAAE,CAC5E,EACMqB,EAAc,CAAC,GAAGrC,EAAM,KAAK,CAAC,EAAE,KAAK,CAACe,EAAGC,IAAMC,EAAQF,EAAGC,CAAC,CAAC,EAClE,MAAO,CACL,aACA,GAAGoB,EAAe,IAAI9C,GAAKoC,EAAcpC,CAAC,CAAC,EAC3C,GAAG+C,EAAY,IAAIC,GAAKN,EAAWhC,EAAM,IAAIsC,CAAC,CAAE,CAAC,CACnD,EAAE,KAAK;AAAA,CAAI,CACb,CC3QA,SAASC,EAAaC,EAAyB,CAC7C,OAAOA,EAAS,OAAOC,GAAWA,EAAQ,KAAO,CAAC,CACpD,CAOA,eAAeC,EACbC,EACAC,EACAC,EACA,CACA,IAAMC,EAAoB,CAAC,EAC3B,aAAMC,EAASJ,EAAMC,EAAK,QAASA,EAAK,MAAOA,EAAK,IAAK,CACvD,OAAAC,EACA,aAAcG,GAAQ,CACpBF,EAAM,KAAKG,EAAcD,CAAI,CAAC,CAChC,CACF,CAAC,EACMF,CACT,CAiBA,eAAsBI,EACpBC,EACAC,EACA,CAAE,KAAAC,EAAO,EAAG,OAAAR,CAAO,EAA6C,CAAC,EACjE,CACA,IAAML,EAAW,IAAI,IACfc,EAAQ,IAAI,IACZC,EAAS,CAAE,SAAU,EAAM,EAYjC,GAXA,MAAMR,EAASI,EAAM,SAAUC,EAAO,QAASA,EAAO,MAAOA,EAAO,IAAK,CACvE,OAAAP,EACA,aAAcG,GAAQ,CACpB,IAAMQ,EAAMC,EAAgBT,CAAI,EAChCO,EAAO,WAAaC,EAAI,SACxBhB,EAAS,IAAIgB,EAAI,QAAQ,GAAIA,EAAI,OAAO,CAC1C,CACF,CAAC,EAIGD,EAAO,SAAU,CACnB,IAAMG,EAAMC,EACVnB,EACA,MAAME,EAAUS,EAAM,MAAOC,EAAQP,CAAM,EAC3CO,CACF,EACA,OAAOQ,EAAeF,EAAI,SAAUA,EAAI,KAAK,CAC/C,CACA,IAAMG,EAAe,MAAOjB,IACzB,MAAMF,EAAUS,EAAM,MAAOP,EAAMC,CAAM,GAAG,QAAQiB,GACnDC,EAAQvB,EAAUc,EAAOQ,CAAI,CAC/B,EACEE,EAAWzB,EAAa,MAAMsB,EAAaT,CAAM,CAAC,EACtD,QAASa,EAAM,EAAGA,EAAMZ,EAAMY,IAAO,CAInC,IAAMC,EAAU,MAAM,QAAQ,IAAIF,EAAS,IAAIH,CAAY,CAAC,EAC5DG,EAAWzB,EAAa2B,EAAQ,KAAK,CAAC,CACxC,CACA,GAAIb,EAAO,EAAG,CACZ,IAAMc,EAAQC,EACZ,CAAC,GAAG5B,EAAS,OAAO,CAAC,EACrB,IAAI,IAAIwB,EAAS,IAAIvB,GAAWA,EAAQ,EAAE,CAAC,EAC3CW,CACF,EACMN,EAAQ,MAAM,QAAQ,IAC1BqB,EAAM,IAAIvB,GAAQF,EAAUS,EAAM,MAAOP,EAAMC,CAAM,CAAC,CACxD,EACA,QAAWiB,KAAQO,EAAavB,EAAM,KAAK,EAAGN,EAAU2B,CAAK,EAC3DJ,EAAQvB,EAAUc,EAAOQ,CAAI,CAEjC,CACA,OAAOF,EAAepB,EAAUc,CAAK,CACvC,CChHA,IAAqBgB,EAArB,cAA8CC,CAA+C,CAA7F,kCAGE,KAAiB,KAAO,KAAK,UAAU,CAAC,CAAC,EAFzC,YAAc,aAAe,CAAC,cAAe,aAAa,EASlD,UAAUC,EAA6B,CAC7C,IAAMC,EAAWC,EACf,KACA,CAAC,GAAGF,EAAO,kBAAkB,EAC7B,CAAC,GAAGA,EAAO,eAAe,CAC5B,EACMG,EAAQD,EACZ,KACA,CAAC,GAAGF,EAAO,eAAe,EAC1B,CAAC,GAAGA,EAAO,YAAY,CACzB,EACMI,EAA8B,KAAK,QAAQ,eAAe,EAC1DC,EACJL,EAAM,SAAW,GAAKM,EAAcF,CAAa,EAC7C,IAAIG,EACFL,EAAc,KAAM,CAAC,eAAe,EAAG,CAAC,YAAY,CAAC,EACrD,IACF,EACA,OACN,MAAO,CACL,SAAAD,EACA,MAAAE,EACA,MAAAE,EACA,SAAU,IAAIG,EAAcH,GAAO,MAAQJ,EAAU,IAAI,CAC3D,CACF,CAMA,MAAM,kBAAkBQ,EAAoB,CAAC,EAAG,CAC9C,OAAO,KAAK,KAAK,OAAO,eAAeA,CAAI,CAC7C,CAEQ,MAAMC,EAAqB,OAAQ,CACzC,OAAIA,IAAS,OACJ,KAAK,MAEd,KAAK,SAAW,KAAK,UAAU,CAAC,QAAQ,CAAC,EAClC,KAAK,OACd,CAEA,MAAM,YAAYD,EAAoB,CAAC,EAAG,CACxC,OAAO,KAAK,KAAK,SAAS,iBAAiBA,CAAI,CACjD,CAGA,MAAM,kBAAkBE,EAAmBF,EAAoB,CAAC,EAAG,CACjE,GAAM,CAAE,MAAAJ,EAAO,SAAAO,CAAS,EAAI,KAAK,KACjC,OAAOA,EAAS,YACdD,EACAF,EACAJ,IAAUQ,GAASR,EAAM,OAAO,WAAWQ,EAAOJ,CAAI,EACxD,CACF,CAEA,MAAa,mBAAoB,CAC/B,MAAO,EACT,CAEA,YAAYK,EAAeL,EAAoB,CAAC,EAAG,CACjD,OAAOM,EAA0B,MAAMC,GAAY,CACjD,IAAMC,EAAe,MAAM,KAAK,KAAK,SAAS,QAAQH,EAAOL,CAAI,EACjE,GAAIQ,IAAiB,OAAW,CAC9B,GAAM,CAAE,MAAAZ,CAAM,EAAI,KAAK,KACjBa,EAAQC,GAAqB,CACjCH,EAAS,KAAKG,CAAO,CACvB,EACA,MAAOd,EACH,KAAK,aAAaA,EAAOY,EAAcH,EAAOL,EAAMS,CAAI,EACxD,KAAK,gBAAgBD,EAAcH,EAAOL,EAAMS,CAAI,EAC1D,CACAF,EAAS,SAAS,CACpB,EAAGP,EAAK,MAAM,CAChB,CAGA,MAAc,aACZJ,EACAY,EACAH,EACAL,EACAS,EACA,CACA,IAAME,EAAQ,MAAMC,EAClB,oBACAZ,EAAK,eACL,IAAMJ,EAAM,eAAeY,EAAcH,EAAOL,CAAI,CACtD,EACA,QAAWa,KAAQF,EAAO,CACxB,IAAMG,EAAK,OAAOD,EAAK,EAAE,EACzBJ,EACE,IAAIM,EAAc,CAChB,SAAUD,EACV,QAAST,EAAM,QACf,MAAOQ,EAAK,MACZ,IAAKA,EAAK,IACV,KAAMC,EACN,KAAM,UACN,KAAM,EACN,WAAYD,EAAK,OACnB,CAAC,CACH,CACF,CACF,CAEA,MAAc,gBACZL,EACAH,EACAL,EACAS,EACA,CACA,MAAMG,EAAa,uBAAwBZ,EAAK,eAAgB,IAC9DgB,EAAS,KAAK,KAAK,SAAUR,EAAcH,EAAM,MAAOA,EAAM,IAAK,CACjE,OAAQL,EAAK,OACb,aAAciB,GAAQ,CACpB,GAAM,CAAE,QAAAC,EAAS,SAAAC,CAAS,EAAIC,EAAgBH,CAAI,EAGlD,GAAIE,GAAYD,EAAQ,UAAYV,EAClC,OAMF,IAAMa,EAAUC,EAAeJ,CAAO,EACtCT,EACE,IAAIM,EAAc,CAChB,SAAUG,EAAQ,GAClB,QAASb,EAAM,QACf,MAAOa,EAAQ,MACf,IAAKA,EAAQ,IACb,KAAMA,EAAQ,GACd,KAAM,UACN,KAAMA,EAAQ,KACd,WAAYA,EAAQ,QACpB,GAAIG,GAAW,CAAE,QAAAA,EAAS,YAAaA,EAAQ,MAAO,CACxD,CAAC,CACH,CACF,CACF,CAAC,CACH,CACF,CAKA,MAAM,YACJE,EACAvB,EAA+B,CAAC,EACL,CAC3B,GAAM,CAAE,KAAAwB,EAAM,OAAAC,EAAQ,KAAAxB,CAAK,EAAID,EACzB0B,EAAQ,KAAK,MAAMzB,CAAI,EACvBO,EAAe,MAAMkB,EAAM,SAAS,QAAQH,EAAQ,CAAE,OAAAE,CAAO,CAAC,EACpE,GAAIjB,IAAiB,OACnB,MAAM,IAAI,MACR,GAAGe,EAAO,YAAY,IAAIA,EAAO,OAAO,iFAAiFA,EAAO,YAAY,8CAC9I,EAEF,OAAIG,EAAM,MACDA,EAAM,MAAM,SAASlB,EAAce,EAAQvB,CAAI,EAEjD2B,EACLD,EACA,CAAE,QAASlB,EAAc,MAAOe,EAAO,MAAO,IAAKA,EAAO,GAAI,EAC9D,CAAE,KAAAC,EAAM,OAAAC,CAAO,CACjB,CACF,CACF",
6
+ "names": ["parseSegmentRow", "line", "cols", "anchored", "own", "segmentSamples", "segment", "tag", "t", "samples", "s", "parseEndpoint", "field", "parseLinkLine", "source", "target", "linkKey", "link", "addLink", "segments", "links", "added", "anchoredCut", "under", "found", "window", "side", "closingSpans", "held", "unread", "outside", "mergedSpans", "isUnread", "maxGap", "sorted", "a", "b", "compare", "spans", "span", "refName", "start", "end", "closingLinks", "inside", "GFA_TAG", "formatSegment", "id", "rank", "tags", "valid", "extra", "formatLink", "sourceStrand", "targetStrand", "formatSubgraph", "sortedSegments", "sortedLinks", "k", "offReference", "segments", "segment", "linksOver", "file", "span", "signal", "found", "getLines", "line", "parseLinkLine", "rgfaCut", "files", "window", "hops", "links", "layout", "row", "parseSegmentRow", "cut", "anchoredCut", "formatSubgraph", "addLinksOver", "link", "addLink", "frontier", "hop", "reached", "spans", "closingSpans", "closingLinks", "RgfaTabixAdapter", "BaseFeatureDataAdapter", "under", "segments", "openTabixSlot", "links", "walksLocation", "walks", "locationIsSet", "WalkIndex", "PanSNRefNames", "opts", "tier", "regions", "refNames", "start", "query", "ObservableCreate", "observer", "tabixRefName", "emit", "feature", "nodes", "updateStatus", "node", "id", "SimpleFeature", "getLines", "line", "segment", "anchored", "parseSegmentRow", "samples", "segmentSamples", "region", "hops", "signal", "index", "rgfaCut"]
7
+ }
@@ -1,2 +1,2 @@
1
- import{H as i,g as o,m as a}from"./chunk-E3WYW7Z6.js";import{c as r,d as e}from"./chunk-JMBDTNUR.js";var u=e(function({model:n,title:s,children:l,handleClose:m}){return r(o,{open:!0,title:s,onClose:()=>{m()},children:[l,n.trackWarnings.flatMap(({name:f,warnings:g},p)=>g.map((t,c)=>r(a,{severity:"warning",style:{marginBottom:8},children:[r(i,{variant:"subtitle2",children:[f,": ",t.message]}),t.effect]},`${p}_${c}`)))]})}),x=u;export{x as default};
2
- //# sourceMappingURL=TrackWarningsDialog-QT4YEXGD.js.map
1
+ import{H as i,h as o,m as a}from"./chunk-UGCNPRKC.js";import{c as r,d as e}from"./chunk-JMBDTNUR.js";var u=e(function({model:n,title:s,children:l,handleClose:m}){return r(o,{open:!0,title:s,onClose:()=>{m()},children:[l,n.trackWarnings.flatMap(({name:f,warnings:g},p)=>g.map((t,c)=>r(a,{severity:"warning",style:{marginBottom:8},children:[r(i,{variant:"subtitle2",children:[f,": ",t.message]}),t.effect]},`${p}_${c}`)))]})}),x=u;export{x as default};
2
+ //# sourceMappingURL=TrackWarningsDialog-TRYBOIHX.js.map
@@ -0,0 +1,2 @@
1
+ import{a as G}from"./chunk-MND7GAOG.js";import{a as z}from"./chunk-WCGVUZWP.js";import{b as j,d as H,e as D,f as E,g as X}from"./chunk-EGDLEQ47.js";import{a as $}from"./chunk-KOB7YIMR.js";import{p as W,q as R,r as I}from"./chunk-T4HOC6O3.js";import"./chunk-FVK2THE4.js";import{e as q,m as T,n as _}from"./chunk-FHZXW7IG.js";import"./chunk-YBLNN5FU.js";import"./chunk-MVVXRZXV.js";import{a as F,c as v,d as S,e as P,f as k,g as x}from"./chunk-ZMOMFU2B.js";import"./chunk-E7J6SB5P.js";import"./chunk-7WIBGKBU.js";import"./chunk-OGQ6OQDN.js";import"./chunk-NDRLCGPC.js";import{o as C}from"./chunk-TENPGIPX.js";var ee=50;function te(i){let o=[];for(let e=0;e<i.length;e++){let[t,s]=i[e],[n,a]=i[e+1]??[];s===a&&s<=ee&&(t==="I"&&n==="D"||t==="D"&&n==="I")?(o.push(["X",s]),e++):o.push([t,s])}return o}function M(i){return Array.from(i.ids,(o,e)=>2*o+i.rev[e])}function K(i,o){let e=i.get(o);return e===void 0?0:e.end-e.start}function U(i){let o=new Map,e=n=>{let a=K(i,n),r=o.get(a);return r===void 0&&(r="N".repeat(a),o.set(a,r)),r},t=new Map,s=n=>{let a=t.get(n);return a===void 0&&(a=M(n),t.set(n,a)),a};return(n,a)=>R({query:s(n),target:s(a),sequenceOf:e,bases:!1})}var ae=100;function O(i,o,e){let t=i.at(-1);e>0&&(t?.[0]===o?t[1]+=e:i.push([o,e]))}function J(i){let o=t=>K(i,t>>1),e=t=>{let s=new Float64Array(t.length+1);return t.forEach((n,a)=>{s[a+1]=s[a]+o(n)}),s};return(t,s)=>{let n=M(s),a=M(t),[r,f]=W(a,n,o),[g,y]=W(a.map(u=>u^1).reverse(),n,o),l=y>f,p=l?g:r,N=Math.max(f,y);if(N<ae)return[];let h=e(n),c=e(a),b=u=>l?a.length-1-p[u][0]:p[u][0],m=[];p.forEach(([,u],w)=>{if(w>0){let[Q,V]=[b(w-1),b(w)].sort((Y,Z)=>Y-Z);O(m,"I",c[V]-c[Q+1]),O(m,"D",h[u]-h[p[w-1][1]+1])}O(m,"=",o(n[u]))});let[d,A]=[b(0),b(p.length-1)].sort((u,w)=>u-w);return[{queryStart:c[d],queryEnd:c[A+1],targetStart:h[p[0][1]],targetEnd:h[p.at(-1)[1]+1],strand:l?"-":"+",edits:m,sharedBases:N}]}}function B({chain:i,target:o,query:e,assemblyName:t,refName:s,mateAssemblyName:n}){let a=te(i.edits);return X({target:{name:o.name,assemblyName:t,refName:s,start:o.hapStart+i.targetStart,end:o.hapStart+i.targetEnd},query:{name:e.name,assemblyName:n,refName:S(e.name),start:e.hapStart+i.queryStart,end:e.hapStart+i.queryEnd},strand:i.strand,cigar:I(a),matches:a.reduce((r,[f,g])=>r+(f==="="?g:0),0),columns:a.reduce((r,[,f])=>r+f,0)})}var ne=i=>v(i.name)??i.name,L=class extends j{constructor(){super(...arguments);this.walks=new G(T(this,["walksLocation"],["walksIndex"]),this);this.refNames=new _(this.walks.file,this)}get anchor(){let e=this.getConf("assemblyNames")[0];if(e===void 0)throw new Error("WalkTabixSyntenyAdapter needs assemblyNames: its first entry is the assembly the reference sample is loaded as");return e}async getHaplotypeNames(e={}){return this.walks.haplotypeNames(e)}async getHeader(e={}){let{haplotypes:t}=await this.walks.reader.header(e),{anchor:s}=this,n=k(this,s),a=x(this);return{hasCoarseTier:!1,anchorAssemblyName:s,lanes:t.filter(r=>r!==n&&F(r)!==n).map(r=>({name:E(a,r),label:r,group:F(r)}))}}async getRefNames(e={}){let{assemblyName:t}=e;return t===void 0||t===this.anchor?this.refNames.assemblyRefNames({...e,assemblyName:this.anchor}):[]}async getSubgraph(e,t={}){let s=await this.refNames.resolve(e,{signal:t.signal});if(s===void 0)throw new Error(`${e.assemblyName} ${e.refName} is not in this graph's index; the graph is cut on ${this.anchor}, mapped to its PanSN prefix in assemblyNameToPanSN`);return this.walks.subgraph(s,e,t)}async laneWalks(e,t,s,n,a){return C(a,n.statusCallback,()=>this.walks.reader.cut(e,t,{...this.walks.budgets(),signal:n.signal,keep:q(s,e),links:!1}).catch(r=>{throw z(r)&&r instanceof Error&&(r.message=r.message.replace("to see the graph","to see lanes")),r}))}async anchorFeatures(e,t,s){let n=s.haplotypes?.length?s.haplotypes.map(m=>k(this,m)):void 0,a=k(this,s.targetAssemblyName),r=a===void 0?n:n===void 0||n.some(m=>P(a,m)||P(m,a))?[a]:[];if(r?.length===0)return[];let{start:f,end:g}=t,{graph:y,fragments:l}=await this.laneWalks(e,t,r,s,`Reading walks ${t.refName}:${f.toLocaleString()}-${g.toLocaleString()}`),p=J(y.nodes),N=x(this),h=l.filter(m=>m.name===e),c=m=>r===void 0?0:r.findIndex(d=>P(m.name,d));return l.filter(m=>m.name!==e&&c(m)>=0).sort((m,d)=>c(m)-c(d)).flatMap(m=>h.flatMap(d=>p(m,d).map(A=>B({chain:A,target:d,query:m,assemblyName:t.assemblyName,refName:t.refName,mateAssemblyName:E(N,ne(m))}))))}async pairFeatures(e,t,s,n){let a=s.map(l=>({pair:l,featurePrefix:k(this,l.queryAssemblyName),matePrefix:k(this,l.targetAssemblyName)})),{graph:r,fragments:f}=await this.laneWalks(e,t,a.flatMap(l=>[l.featurePrefix,l.matePrefix]),n,D(s)),g=U(r.nodes),y=l=>f.filter(p=>P(p.name,l));return a.flatMap(({pair:l,featurePrefix:p,matePrefix:N})=>y(p).flatMap(h=>y(N).flatMap(c=>g(c,h).map(b=>B({chain:b,target:h,query:c,assemblyName:l.queryAssemblyName,refName:S(h.name),mateAssemblyName:l.targetAssemblyName})))))}getFeatures(e,t={}){return $(async s=>{let n=e.assemblyName===this.anchor?await this.refNames.resolve(e,t):void 0;if(n!==void 0){let a=H(t),r=a.length>0?await this.pairFeatures(n,e,a,t):await this.anchorFeatures(n,e,t);for(let f of r)s.next(f)}s.complete()},t.signal)}};export{L as default};
2
+ //# sourceMappingURL=WalkTabixSyntenyAdapter-U53JLYVU.js.map
@@ -0,0 +1,7 @@
1
+ {
2
+ "version": 3,
3
+ "sources": ["../../src/WalkTabixSyntenyAdapter/walkLanes.ts", "../../src/WalkTabixSyntenyAdapter/WalkTabixSyntenyAdapter.ts"],
4
+ "sourcesContent": ["import { weightedLcs } from '@gmod/gbz-base/lcs'\nimport { pairAlignments, pairCigar } from '@gmod/gbz-base/pairAlignment'\nimport { panSNContig } from '@jbrowse/bandage-core/pansn'\n\nimport { matchFeature } from '../synteny/matchFeature.ts'\n\nimport type { PairChain, PairEdit } from '@gmod/gbz-base/pairAlignment'\nimport type { WalkFragment, WalkNode } from '@jbrowse/graph-walks'\n\n// The longest bubble with equal bases on both sides that is written as\n// mismatches; a longer one stays an insertion and a deletion\nconst BALANCED_BUBBLE_MAX_BP = 50\n\n/**\n * A chain's edits with each short balanced bubble as `X`. The walks carry no\n * bases, so two walks differ wherever they leave a shared node, and a SNP is\n * one base inserted and one deleted until it is read as the mismatch it is.\n */\nexport function bubblesAsMismatches(edits: PairEdit[]) {\n const out: PairEdit[] = []\n for (let i = 0; i < edits.length; i++) {\n const [op, len] = edits[i]!\n const [nextOp, nextLen] = edits[i + 1] ?? []\n if (\n len === nextLen &&\n len <= BALANCED_BUBBLE_MAX_BP &&\n ((op === 'I' && nextOp === 'D') || (op === 'D' && nextOp === 'I'))\n ) {\n out.push(['X', len])\n i++\n } else {\n out.push([op, len])\n }\n }\n return out\n}\n\n// A walk's steps as the handles gbz-base reads, 2 * id + reversed\nfunction handlesOf(fragment: WalkFragment) {\n return Array.from(fragment.ids, (id, i) => 2 * id + fragment.rev[i]!)\n}\n\nfunction nodeLength(nodes: Map<number, WalkNode>, id: number) {\n const node = nodes.get(id)\n return node === undefined ? 0 : node.end - node.start\n}\n\n/**\n * Aligns one walk fragment to another on the nodes both visit, from node\n * lengths alone: the fragments' steps as the handles gbz-base's chainer reads,\n * and a stand-in sequence of each node's length.\n */\nexport function walkAligner(nodes: Map<number, WalkNode>) {\n const standIn = new Map<number, string>()\n const sequenceOf = (id: number) => {\n const length = nodeLength(nodes, id)\n let sequence = standIn.get(length)\n if (sequence === undefined) {\n sequence = 'N'.repeat(length)\n standIn.set(length, sequence)\n }\n return sequence\n }\n const handles = new Map<WalkFragment, number[]>()\n const cachedHandles = (fragment: WalkFragment) => {\n let steps = handles.get(fragment)\n if (steps === undefined) {\n steps = handlesOf(fragment)\n handles.set(fragment, steps)\n }\n return steps\n }\n return (query: WalkFragment, target: WalkFragment) =>\n pairAlignments({\n query: cachedHandles(query),\n target: cachedHandles(target),\n sequenceOf,\n bases: false,\n })\n}\n\n// The fewest bases on shared nodes a record against the reference keeps,\n// gbz-base's floor for a pair record\nconst MIN_MATCH_BP = 100\n\nfunction appendEdit(edits: PairEdit[], op: PairEdit[0], len: number) {\n const last = edits.at(-1)\n if (len > 0) {\n if (last?.[0] === op) {\n last[1] += len\n } else {\n edits.push([op, len])\n }\n }\n}\n\n/**\n * Aligns a haplotype's walk fragment to the reference's as gbz-base aligns a\n * walk to its reference path: the heaviest common subsequence of their steps,\n * each node weighing its length, read in whichever direction shares more. A\n * pair chain leaves out every node either walk passes twice, which is the\n * whole of a collapsed copy-number array (amylase, KIV-2); the subsequence\n * matches one pass over it, and the other copies are insertions.\n */\nexport function referenceAligner(nodes: Map<number, WalkNode>) {\n const lengthOf = (handle: number) => nodeLength(nodes, handle >> 1)\n const offsets = (handles: number[]) => {\n const out = new Float64Array(handles.length + 1)\n handles.forEach((handle, i) => {\n out[i + 1] = out[i]! + lengthOf(handle)\n })\n return out\n }\n\n return (query: WalkFragment, target: WalkFragment): PairChain[] => {\n const ref = handlesOf(target)\n const hap = handlesOf(query)\n const [same, sameWeight] = weightedLcs(hap, ref, lengthOf)\n const [turned, turnedWeight] = weightedLcs(\n hap.map(handle => handle ^ 1).reverse(),\n ref,\n lengthOf,\n )\n const reverse = turnedWeight > sameWeight\n const pairs = reverse ? turned : same\n const sharedBases = Math.max(sameWeight, turnedWeight)\n if (sharedBases < MIN_MATCH_BP) {\n return []\n }\n const refAt = offsets(ref)\n const hapAt = offsets(hap)\n // a pair's haplotype step, counted along the haplotype\n const stepOf = (k: number) =>\n reverse ? hap.length - 1 - pairs[k]![0] : pairs[k]![0]\n const edits: PairEdit[] = []\n pairs.forEach(([, r], k) => {\n if (k > 0) {\n const [lo, hi] = [stepOf(k - 1), stepOf(k)].sort((x, y) => x - y)\n appendEdit(edits, 'I', hapAt[hi!]! - hapAt[lo! + 1]!)\n appendEdit(edits, 'D', refAt[r]! - refAt[pairs[k - 1]![1] + 1]!)\n }\n appendEdit(edits, '=', lengthOf(ref[r]!))\n })\n const [lo, hi] = [stepOf(0), stepOf(pairs.length - 1)].sort((x, y) => x - y)\n return [\n {\n queryStart: hapAt[lo!]!,\n queryEnd: hapAt[hi! + 1]!,\n targetStart: refAt[pairs[0]![1]]!,\n targetEnd: refAt[pairs.at(-1)![1] + 1]!,\n strand: reverse ? '-' : '+',\n edits,\n sharedBases,\n },\n ]\n }\n}\n\n/**\n * One record of a walk aligned to another, with its short balanced bubbles as\n * mismatches\n */\nexport function chainFeature({\n chain,\n target,\n query,\n assemblyName,\n refName,\n mateAssemblyName,\n}: {\n chain: PairChain\n target: WalkFragment\n query: WalkFragment\n assemblyName: string\n refName: string\n mateAssemblyName: string\n}) {\n const edits = bubblesAsMismatches(chain.edits)\n return matchFeature({\n target: {\n name: target.name,\n assemblyName,\n refName,\n start: target.hapStart + chain.targetStart,\n end: target.hapStart + chain.targetEnd,\n },\n query: {\n name: query.name,\n assemblyName: mateAssemblyName,\n refName: panSNContig(query.name),\n start: query.hapStart + chain.queryStart,\n end: query.hapStart + chain.queryEnd,\n },\n strand: chain.strand,\n cigar: pairCigar(edits),\n matches: edits.reduce((sum, [op, n]) => sum + (op === '=' ? n : 0), 0),\n columns: edits.reduce((sum, [, n]) => sum + n, 0),\n })\n}\n", "import {\n assemblyByPanSNPrefix,\n panSNContig,\n panSNHaplotype,\n panSNMatchesPrefix,\n panSNSample,\n resolvePanSNPrefix,\n} from '@jbrowse/bandage-core/pansn'\nimport { updateStatus } from '@jbrowse/core/util'\nimport { ObservableCreate } from '@jbrowse/core/util/rxjs'\nimport { walksKept } from '@jbrowse/graph-walks'\n\nimport { chainFeature, referenceAligner, walkAligner } from './walkLanes.ts'\nimport { isNodeLimitError } from '../LinearGraphDisplay/denseWindow.ts'\nimport { WalkIndex } from '../RgfaTabixAdapter/walkReader.ts'\nimport { PanSNRefNames, openTabixSlot } from '../panSNTabix.ts'\nimport { ComparativeAdapterBase } from '../synteny/ComparativeAdapterBase.ts'\nimport {\n laneAssemblyName,\n lanePairsStatus,\n requestedPairs,\n} from '../synteny/lanePairs.ts'\n\nimport type { WalkTabixSyntenyAdapterConfig } from './configSchema.ts'\nimport type { SubgraphAdapterOptions } from '../GetSubgraph.ts'\nimport type { LaneFeatureOptions, LanePair } from '../synteny/lanePairs.ts'\nimport type { BaseOptions } from '@jbrowse/core/data_adapters/BaseAdapter'\nimport type { Feature } from '@jbrowse/core/util'\nimport type { Region } from '@jbrowse/core/util/types'\nimport type { WalkFragment } from '@jbrowse/graph-walks'\n\nconst haplotypeOf = (fragment: WalkFragment) =>\n panSNHaplotype(fragment.name) ?? fragment.name\n\nexport default class WalkTabixSyntenyAdapter extends ComparativeAdapterBase<WalkTabixSyntenyAdapterConfig> {\n private readonly walks: WalkIndex = new WalkIndex(\n openTabixSlot(this, ['walksLocation'], ['walksIndex']),\n this,\n )\n private readonly refNames: PanSNRefNames = new PanSNRefNames(\n this.walks.file,\n this,\n )\n\n private get anchor() {\n const anchor = this.getConf('assemblyNames')[0]\n if (anchor === undefined) {\n throw new Error(\n 'WalkTabixSyntenyAdapter needs assemblyNames: its first entry is the assembly the reference sample is loaded as',\n )\n }\n return anchor\n }\n\n /**\n * The haplotypes the walk file's header names, `sample#haplotype` in byte\n * order, or undefined when it names none\n */\n async getHaplotypeNames(opts: BaseOptions = {}) {\n return this.walks.haplotypeNames(opts)\n }\n\n /**\n * `lanes` is every haplotype the header names but the anchor's own,\n * declared up front so the display can offer the whole graph's haplotypes\n * before a fetch has placed any of them\n */\n async getHeader(opts: BaseOptions = {}) {\n const { haplotypes } = await this.walks.reader.header(opts)\n const { anchor } = this\n const anchorPrefix = resolvePanSNPrefix(this, anchor)\n const asmByPrefix = assemblyByPanSNPrefix(this)\n return {\n hasCoarseTier: false,\n anchorAssemblyName: anchor,\n lanes: haplotypes\n .filter(\n prefix =>\n prefix !== anchorPrefix && panSNSample(prefix) !== anchorPrefix,\n )\n .map(prefix => ({\n name: laneAssemblyName(asmByPrefix, prefix),\n label: prefix,\n group: panSNSample(prefix),\n })),\n }\n }\n\n // A lane's contigs are named only in its walk rows, so only the anchor's\n // are listed\n async getRefNames(opts: BaseOptions = {}) {\n const { assemblyName } = opts\n return assemblyName === undefined || assemblyName === this.anchor\n ? this.refNames.assemblyRefNames({ ...opts, assemblyName: this.anchor })\n : []\n }\n\n /**\n * The graph display's cut of the window, as RgfaTabixAdapter cuts a\n * walk-indexed graph: the reference walk first, then the walks `haplotypes`\n * asks for, and the nodes and links they visit\n */\n async getSubgraph(region: Region, opts: SubgraphAdapterOptions = {}) {\n const refName = await this.refNames.resolve(region, { signal: opts.signal })\n if (refName === undefined) {\n throw new Error(\n `${region.assemblyName} ${region.refName} is not in this graph's index; the graph is cut on ${this.anchor}, mapped to its PanSN prefix in assemblyNameToPanSN`,\n )\n }\n return this.walks.subgraph(refName, region, opts)\n }\n\n // The window's walks for the haplotypes named and the reference, with no\n // link read; a window past a budget fails as a zoom-in notice for lanes\n private async laneWalks(\n refName: string,\n region: Region,\n prefixes: string[] | undefined,\n opts: LaneFeatureOptions,\n status: string,\n ) {\n return updateStatus(status, opts.statusCallback, () =>\n this.walks.reader\n .cut(refName, region, {\n ...this.walks.budgets(),\n signal: opts.signal,\n keep: walksKept(prefixes, refName),\n links: false,\n })\n .catch((error: unknown) => {\n if (isNodeLimitError(error) && error instanceof Error) {\n error.message = error.message.replace(\n 'to see the graph',\n 'to see lanes',\n )\n }\n throw error\n }),\n )\n }\n\n /**\n * One record per stretch a haplotype's walk shares with the reference's,\n * on the reference with the haplotype as its mate\n */\n private async anchorFeatures(\n refName: string,\n region: Region,\n opts: LaneFeatureOptions,\n ) {\n const wanted = opts.haplotypes?.length\n ? opts.haplotypes.map(lane => resolvePanSNPrefix(this, lane))\n : undefined\n const target = resolvePanSNPrefix(this, opts.targetAssemblyName)\n const prefixes =\n target === undefined\n ? wanted\n : wanted === undefined ||\n wanted.some(\n prefix =>\n panSNMatchesPrefix(target, prefix) ||\n panSNMatchesPrefix(prefix, target),\n )\n ? [target]\n : []\n if (prefixes?.length === 0) {\n return []\n }\n const { start, end } = region\n const { graph, fragments } = await this.laneWalks(\n refName,\n region,\n prefixes,\n opts,\n `Reading walks ${region.refName}:${start.toLocaleString()}-${end.toLocaleString()}`,\n )\n const align = referenceAligner(graph.nodes)\n const asmByPrefix = assemblyByPanSNPrefix(this)\n const references = fragments.filter(f => f.name === refName)\n // lanes come back in the order the fetch names them, which is the order\n // the display stacks lanes of equal weight in\n const rank = (fragment: WalkFragment) =>\n prefixes === undefined\n ? 0\n : prefixes.findIndex(prefix =>\n panSNMatchesPrefix(fragment.name, prefix),\n )\n const shown = fragments\n .filter(fragment => fragment.name !== refName && rank(fragment) >= 0)\n .sort((a, b) => rank(a) - rank(b))\n return shown.flatMap(query =>\n references.flatMap(reference =>\n align(query, reference).map(chain =>\n chainFeature({\n chain,\n target: reference,\n query,\n assemblyName: region.assemblyName,\n refName: region.refName,\n mateAssemblyName: laneAssemblyName(asmByPrefix, haplotypeOf(query)),\n }),\n ),\n ),\n )\n }\n\n /**\n * Each pair's two lanes read out of the anchor window, all pairs in one\n * read, and each walk of a pair's query lane aligned to each walk of its\n * target lane. A pair's records sit on its query lane's contigs, the lane\n * the display draws on top.\n */\n private async pairFeatures(\n refName: string,\n region: Region,\n pairs: LanePair[],\n opts: LaneFeatureOptions,\n ) {\n const sides = pairs.map(pair => ({\n pair,\n featurePrefix: resolvePanSNPrefix(this, pair.queryAssemblyName),\n matePrefix: resolvePanSNPrefix(this, pair.targetAssemblyName),\n }))\n const { graph, fragments } = await this.laneWalks(\n refName,\n region,\n sides.flatMap(side => [side.featurePrefix, side.matePrefix]),\n opts,\n lanePairsStatus(pairs),\n )\n const align = walkAligner(graph.nodes)\n const walksOf = (prefix: string) =>\n fragments.filter(f => panSNMatchesPrefix(f.name, prefix))\n return sides.flatMap(({ pair, featurePrefix, matePrefix }) =>\n walksOf(featurePrefix).flatMap(target =>\n walksOf(matePrefix).flatMap(query =>\n align(query, target).map(chain =>\n chainFeature({\n chain,\n target,\n query,\n assemblyName: pair.queryAssemblyName,\n refName: panSNContig(target.name),\n mateAssemblyName: pair.targetAssemblyName,\n }),\n ),\n ),\n ),\n )\n }\n\n // The files are filed under the reference alone, so a window on a\n // haplotype lane has no answer: a lane pair is read inside the anchor's\n getFeatures(region: Region, opts: LaneFeatureOptions = {}) {\n return ObservableCreate<Feature>(async observer => {\n const refName =\n region.assemblyName === this.anchor\n ? await this.refNames.resolve(region, opts)\n : undefined\n if (refName !== undefined) {\n const pairs = requestedPairs(opts)\n const features =\n pairs.length > 0\n ? await this.pairFeatures(refName, region, pairs, opts)\n : await this.anchorFeatures(refName, region, opts)\n for (const feature of features) {\n observer.next(feature)\n }\n }\n observer.complete()\n }, opts.signal)\n }\n}\n"],
5
+ "mappings": "+lBAWA,IAAMA,GAAyB,GAOxB,SAASC,GAAoBC,EAAmB,CACrD,IAAMC,EAAkB,CAAC,EACzB,QAASC,EAAI,EAAGA,EAAIF,EAAM,OAAQE,IAAK,CACrC,GAAM,CAACC,EAAIC,CAAG,EAAIJ,EAAME,CAAC,EACnB,CAACG,EAAQC,CAAO,EAAIN,EAAME,EAAI,CAAC,GAAK,CAAC,EAEzCE,IAAQE,GACRF,GAAON,KACLK,IAAO,KAAOE,IAAW,KAASF,IAAO,KAAOE,IAAW,MAE7DJ,EAAI,KAAK,CAAC,IAAKG,CAAG,CAAC,EACnBF,KAEAD,EAAI,KAAK,CAACE,EAAIC,CAAG,CAAC,CAEtB,CACA,OAAOH,CACT,CAGA,SAASM,EAAUC,EAAwB,CACzC,OAAO,MAAM,KAAKA,EAAS,IAAK,CAACC,EAAIP,IAAM,EAAIO,EAAKD,EAAS,IAAIN,CAAC,CAAE,CACtE,CAEA,SAASQ,EAAWC,EAA8BF,EAAY,CAC5D,IAAMG,EAAOD,EAAM,IAAIF,CAAE,EACzB,OAAOG,IAAS,OAAY,EAAIA,EAAK,IAAMA,EAAK,KAClD,CAOO,SAASC,EAAYF,EAA8B,CACxD,IAAMG,EAAU,IAAI,IACdC,EAAcN,GAAe,CACjC,IAAMO,EAASN,EAAWC,EAAOF,CAAE,EAC/BQ,EAAWH,EAAQ,IAAIE,CAAM,EACjC,OAAIC,IAAa,SACfA,EAAW,IAAI,OAAOD,CAAM,EAC5BF,EAAQ,IAAIE,EAAQC,CAAQ,GAEvBA,CACT,EACMC,EAAU,IAAI,IACdC,EAAiBX,GAA2B,CAChD,IAAIY,EAAQF,EAAQ,IAAIV,CAAQ,EAChC,OAAIY,IAAU,SACZA,EAAQb,EAAUC,CAAQ,EAC1BU,EAAQ,IAAIV,EAAUY,CAAK,GAEtBA,CACT,EACA,MAAO,CAACC,EAAqBC,IAC3BC,EAAe,CACb,MAAOJ,EAAcE,CAAK,EAC1B,OAAQF,EAAcG,CAAM,EAC5B,WAAAP,EACA,MAAO,EACT,CAAC,CACL,CAIA,IAAMS,GAAe,IAErB,SAASC,EAAWzB,EAAmBG,EAAiBC,EAAa,CACnE,IAAMsB,EAAO1B,EAAM,GAAG,EAAE,EACpBI,EAAM,IACJsB,IAAO,CAAC,IAAMvB,EAChBuB,EAAK,CAAC,GAAKtB,EAEXJ,EAAM,KAAK,CAACG,EAAIC,CAAG,CAAC,EAG1B,CAUO,SAASuB,EAAiBhB,EAA8B,CAC7D,IAAMiB,EAAYC,GAAmBnB,EAAWC,EAAOkB,GAAU,CAAC,EAC5DC,EAAWZ,GAAsB,CACrC,IAAMjB,EAAM,IAAI,aAAaiB,EAAQ,OAAS,CAAC,EAC/C,OAAAA,EAAQ,QAAQ,CAACW,EAAQ3B,IAAM,CAC7BD,EAAIC,EAAI,CAAC,EAAID,EAAIC,CAAC,EAAK0B,EAASC,CAAM,CACxC,CAAC,EACM5B,CACT,EAEA,MAAO,CAACoB,EAAqBC,IAAsC,CACjE,IAAMS,EAAMxB,EAAUe,CAAM,EACtBU,EAAMzB,EAAUc,CAAK,EACrB,CAACY,EAAMC,CAAU,EAAIC,EAAYH,EAAKD,EAAKH,CAAQ,EACnD,CAACQ,EAAQC,CAAY,EAAIF,EAC7BH,EAAI,IAAIH,GAAUA,EAAS,CAAC,EAAE,QAAQ,EACtCE,EACAH,CACF,EACMU,EAAUD,EAAeH,EACzBK,EAAQD,EAAUF,EAASH,EAC3BO,EAAc,KAAK,IAAIN,EAAYG,CAAY,EACrD,GAAIG,EAAchB,GAChB,MAAO,CAAC,EAEV,IAAMiB,EAAQX,EAAQC,CAAG,EACnBW,EAAQZ,EAAQE,CAAG,EAEnBW,EAAUC,GACdN,EAAUN,EAAI,OAAS,EAAIO,EAAMK,CAAC,EAAG,CAAC,EAAIL,EAAMK,CAAC,EAAG,CAAC,EACjD5C,EAAoB,CAAC,EAC3BuC,EAAM,QAAQ,CAAC,CAAC,CAAEM,CAAC,EAAGD,IAAM,CAC1B,GAAIA,EAAI,EAAG,CACT,GAAM,CAACE,EAAIC,CAAE,EAAI,CAACJ,EAAOC,EAAI,CAAC,EAAGD,EAAOC,CAAC,CAAC,EAAE,KAAK,CAACI,EAAGC,IAAMD,EAAIC,CAAC,EAChExB,EAAWzB,EAAO,IAAK0C,EAAMK,CAAG,EAAKL,EAAMI,EAAM,CAAC,CAAE,EACpDrB,EAAWzB,EAAO,IAAKyC,EAAMI,CAAC,EAAKJ,EAAMF,EAAMK,EAAI,CAAC,EAAG,CAAC,EAAI,CAAC,CAAE,CACjE,CACAnB,EAAWzB,EAAO,IAAK4B,EAASG,EAAIc,CAAC,CAAE,CAAC,CAC1C,CAAC,EACD,GAAM,CAACC,EAAIC,CAAE,EAAI,CAACJ,EAAO,CAAC,EAAGA,EAAOJ,EAAM,OAAS,CAAC,CAAC,EAAE,KAAK,CAACS,EAAGC,IAAMD,EAAIC,CAAC,EAC3E,MAAO,CACL,CACE,WAAYP,EAAMI,CAAG,EACrB,SAAUJ,EAAMK,EAAM,CAAC,EACvB,YAAaN,EAAMF,EAAM,CAAC,EAAG,CAAC,CAAC,EAC/B,UAAWE,EAAMF,EAAM,GAAG,EAAE,EAAG,CAAC,EAAI,CAAC,EACrC,OAAQD,EAAU,IAAM,IACxB,MAAAtC,EACA,YAAAwC,CACF,CACF,CACF,CACF,CAMO,SAASU,EAAa,CAC3B,MAAAC,EACA,OAAA7B,EACA,MAAAD,EACA,aAAA+B,EACA,QAAAC,EACA,iBAAAC,CACF,EAOG,CACD,IAAMtD,EAAQD,GAAoBoD,EAAM,KAAK,EAC7C,OAAOI,EAAa,CAClB,OAAQ,CACN,KAAMjC,EAAO,KACb,aAAA8B,EACA,QAAAC,EACA,MAAO/B,EAAO,SAAW6B,EAAM,YAC/B,IAAK7B,EAAO,SAAW6B,EAAM,SAC/B,EACA,MAAO,CACL,KAAM9B,EAAM,KACZ,aAAciC,EACd,QAASE,EAAYnC,EAAM,IAAI,EAC/B,MAAOA,EAAM,SAAW8B,EAAM,WAC9B,IAAK9B,EAAM,SAAW8B,EAAM,QAC9B,EACA,OAAQA,EAAM,OACd,MAAOM,EAAUzD,CAAK,EACtB,QAASA,EAAM,OAAO,CAAC0D,EAAK,CAACvD,EAAIwD,CAAC,IAAMD,GAAOvD,IAAO,IAAMwD,EAAI,GAAI,CAAC,EACrE,QAAS3D,EAAM,OAAO,CAAC0D,EAAK,CAAC,CAAEC,CAAC,IAAMD,EAAMC,EAAG,CAAC,CAClD,CAAC,CACH,CCvKA,IAAMC,GAAeC,GACnBC,EAAeD,EAAS,IAAI,GAAKA,EAAS,KAEvBE,EAArB,cAAqDC,CAAsD,CAA3G,kCACE,KAAiB,MAAmB,IAAIC,EACtCC,EAAc,KAAM,CAAC,eAAe,EAAG,CAAC,YAAY,CAAC,EACrD,IACF,EACA,KAAiB,SAA0B,IAAIC,EAC7C,KAAK,MAAM,KACX,IACF,EAEA,IAAY,QAAS,CACnB,IAAMC,EAAS,KAAK,QAAQ,eAAe,EAAE,CAAC,EAC9C,GAAIA,IAAW,OACb,MAAM,IAAI,MACR,gHACF,EAEF,OAAOA,CACT,CAMA,MAAM,kBAAkBC,EAAoB,CAAC,EAAG,CAC9C,OAAO,KAAK,MAAM,eAAeA,CAAI,CACvC,CAOA,MAAM,UAAUA,EAAoB,CAAC,EAAG,CACtC,GAAM,CAAE,WAAAC,CAAW,EAAI,MAAM,KAAK,MAAM,OAAO,OAAOD,CAAI,EACpD,CAAE,OAAAD,CAAO,EAAI,KACbG,EAAeC,EAAmB,KAAMJ,CAAM,EAC9CK,EAAcC,EAAsB,IAAI,EAC9C,MAAO,CACL,cAAe,GACf,mBAAoBN,EACpB,MAAOE,EACJ,OACCK,GACEA,IAAWJ,GAAgBK,EAAYD,CAAM,IAAMJ,CACvD,EACC,IAAII,IAAW,CACd,KAAME,EAAiBJ,EAAaE,CAAM,EAC1C,MAAOA,EACP,MAAOC,EAAYD,CAAM,CAC3B,EAAE,CACN,CACF,CAIA,MAAM,YAAYN,EAAoB,CAAC,EAAG,CACxC,GAAM,CAAE,aAAAS,CAAa,EAAIT,EACzB,OAAOS,IAAiB,QAAaA,IAAiB,KAAK,OACvD,KAAK,SAAS,iBAAiB,CAAE,GAAGT,EAAM,aAAc,KAAK,MAAO,CAAC,EACrE,CAAC,CACP,CAOA,MAAM,YAAYU,EAAgBV,EAA+B,CAAC,EAAG,CACnE,IAAMW,EAAU,MAAM,KAAK,SAAS,QAAQD,EAAQ,CAAE,OAAQV,EAAK,MAAO,CAAC,EAC3E,GAAIW,IAAY,OACd,MAAM,IAAI,MACR,GAAGD,EAAO,YAAY,IAAIA,EAAO,OAAO,sDAAsD,KAAK,MAAM,qDAC3G,EAEF,OAAO,KAAK,MAAM,SAASC,EAASD,EAAQV,CAAI,CAClD,CAIA,MAAc,UACZW,EACAD,EACAE,EACAZ,EACAa,EACA,CACA,OAAOC,EAAaD,EAAQb,EAAK,eAAgB,IAC/C,KAAK,MAAM,OACR,IAAIW,EAASD,EAAQ,CACpB,GAAG,KAAK,MAAM,QAAQ,EACtB,OAAQV,EAAK,OACb,KAAMe,EAAUH,EAAUD,CAAO,EACjC,MAAO,EACT,CAAC,EACA,MAAOK,GAAmB,CACzB,MAAIC,EAAiBD,CAAK,GAAKA,aAAiB,QAC9CA,EAAM,QAAUA,EAAM,QAAQ,QAC5B,mBACA,cACF,GAEIA,CACR,CAAC,CACL,CACF,CAMA,MAAc,eACZL,EACAD,EACAV,EACA,CACA,IAAMkB,EAASlB,EAAK,YAAY,OAC5BA,EAAK,WAAW,IAAImB,GAAQhB,EAAmB,KAAMgB,CAAI,CAAC,EAC1D,OACEC,EAASjB,EAAmB,KAAMH,EAAK,kBAAkB,EACzDY,EACJQ,IAAW,OACPF,EACAA,IAAW,QACTA,EAAO,KACLZ,GACEe,EAAmBD,EAAQd,CAAM,GACjCe,EAAmBf,EAAQc,CAAM,CACrC,EACA,CAACA,CAAM,EACP,CAAC,EACT,GAAIR,GAAU,SAAW,EACvB,MAAO,CAAC,EAEV,GAAM,CAAE,MAAAU,EAAO,IAAAC,CAAI,EAAIb,EACjB,CAAE,MAAAc,EAAO,UAAAC,CAAU,EAAI,MAAM,KAAK,UACtCd,EACAD,EACAE,EACAZ,EACA,iBAAiBU,EAAO,OAAO,IAAIY,EAAM,eAAe,CAAC,IAAIC,EAAI,eAAe,CAAC,EACnF,EACMG,EAAQC,EAAiBH,EAAM,KAAK,EACpCpB,EAAcC,EAAsB,IAAI,EACxCuB,EAAaH,EAAU,OAAOI,GAAKA,EAAE,OAASlB,CAAO,EAGrDmB,EAAQtC,GACZoB,IAAa,OACT,EACAA,EAAS,UAAUN,GACjBe,EAAmB7B,EAAS,KAAMc,CAAM,CAC1C,EAIN,OAHcmB,EACX,OAAOjC,GAAYA,EAAS,OAASmB,GAAWmB,EAAKtC,CAAQ,GAAK,CAAC,EACnE,KAAK,CAACuC,EAAGC,IAAMF,EAAKC,CAAC,EAAID,EAAKE,CAAC,CAAC,EACtB,QAAQC,GACnBL,EAAW,QAAQM,GACjBR,EAAMO,EAAOC,CAAS,EAAE,IAAIC,GAC1BC,EAAa,CACX,MAAAD,EACA,OAAQD,EACR,MAAAD,EACA,aAAcvB,EAAO,aACrB,QAASA,EAAO,QAChB,iBAAkBF,EAAiBJ,EAAab,GAAY0C,CAAK,CAAC,CACpE,CAAC,CACH,CACF,CACF,CACF,CAQA,MAAc,aACZtB,EACAD,EACA2B,EACArC,EACA,CACA,IAAMsC,EAAQD,EAAM,IAAIE,IAAS,CAC/B,KAAAA,EACA,cAAepC,EAAmB,KAAMoC,EAAK,iBAAiB,EAC9D,WAAYpC,EAAmB,KAAMoC,EAAK,kBAAkB,CAC9D,EAAE,EACI,CAAE,MAAAf,EAAO,UAAAC,CAAU,EAAI,MAAM,KAAK,UACtCd,EACAD,EACA4B,EAAM,QAAQE,GAAQ,CAACA,EAAK,cAAeA,EAAK,UAAU,CAAC,EAC3DxC,EACAyC,EAAgBJ,CAAK,CACvB,EACMX,EAAQgB,EAAYlB,EAAM,KAAK,EAC/BmB,EAAWrC,GACfmB,EAAU,OAAOI,GAAKR,EAAmBQ,EAAE,KAAMvB,CAAM,CAAC,EAC1D,OAAOgC,EAAM,QAAQ,CAAC,CAAE,KAAAC,EAAM,cAAAK,EAAe,WAAAC,CAAW,IACtDF,EAAQC,CAAa,EAAE,QAAQxB,GAC7BuB,EAAQE,CAAU,EAAE,QAAQZ,GAC1BP,EAAMO,EAAOb,CAAM,EAAE,IAAIe,GACvBC,EAAa,CACX,MAAAD,EACA,OAAAf,EACA,MAAAa,EACA,aAAcM,EAAK,kBACnB,QAASO,EAAY1B,EAAO,IAAI,EAChC,iBAAkBmB,EAAK,kBACzB,CAAC,CACH,CACF,CACF,CACF,CACF,CAIA,YAAY7B,EAAgBV,EAA2B,CAAC,EAAG,CACzD,OAAO+C,EAA0B,MAAMC,GAAY,CACjD,IAAMrC,EACJD,EAAO,eAAiB,KAAK,OACzB,MAAM,KAAK,SAAS,QAAQA,EAAQV,CAAI,EACxC,OACN,GAAIW,IAAY,OAAW,CACzB,IAAM0B,EAAQY,EAAejD,CAAI,EAC3BkD,EACJb,EAAM,OAAS,EACX,MAAM,KAAK,aAAa1B,EAASD,EAAQ2B,EAAOrC,CAAI,EACpD,MAAM,KAAK,eAAeW,EAASD,EAAQV,CAAI,EACrD,QAAWmD,KAAWD,EACpBF,EAAS,KAAKG,CAAO,CAEzB,CACAH,EAAS,SAAS,CACpB,EAAGhD,EAAK,MAAM,CAChB,CACF",
6
+ "names": ["BALANCED_BUBBLE_MAX_BP", "bubblesAsMismatches", "edits", "out", "i", "op", "len", "nextOp", "nextLen", "handlesOf", "fragment", "id", "nodeLength", "nodes", "node", "walkAligner", "standIn", "sequenceOf", "length", "sequence", "handles", "cachedHandles", "steps", "query", "target", "pairAlignments", "MIN_MATCH_BP", "appendEdit", "last", "referenceAligner", "lengthOf", "handle", "offsets", "ref", "hap", "same", "sameWeight", "weightedLcs", "turned", "turnedWeight", "reverse", "pairs", "sharedBases", "refAt", "hapAt", "stepOf", "k", "r", "lo", "hi", "x", "y", "chainFeature", "chain", "assemblyName", "refName", "mateAssemblyName", "matchFeature", "panSNContig", "pairCigar", "sum", "n", "haplotypeOf", "fragment", "panSNHaplotype", "WalkTabixSyntenyAdapter", "ComparativeAdapterBase", "WalkIndex", "openTabixSlot", "PanSNRefNames", "anchor", "opts", "haplotypes", "anchorPrefix", "resolvePanSNPrefix", "asmByPrefix", "assemblyByPanSNPrefix", "prefix", "panSNSample", "laneAssemblyName", "assemblyName", "region", "refName", "prefixes", "status", "updateStatus", "walksKept", "error", "isNodeLimitError", "wanted", "lane", "target", "panSNMatchesPrefix", "start", "end", "graph", "fragments", "align", "referenceAligner", "references", "f", "rank", "a", "b", "query", "reference", "chain", "chainFeature", "pairs", "sides", "pair", "side", "lanePairsStatus", "walkAligner", "walksOf", "featurePrefix", "matePrefix", "panSNContig", "ObservableCreate", "observer", "requestedPairs", "features", "feature"]
7
+ }