jbrowse-plugin-graphgenomeviewer 7.3.1 → 7.3.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.
- package/dist/chunks/{GbzBaseSyntenyAdapter-NYI2YVNG.js → GbzBaseSyntenyAdapter-3MELVQRM.js} +2 -2
- package/dist/chunks/GbzBaseSyntenyAdapter-3MELVQRM.js.map +7 -0
- package/dist/chunks/{GraphGenomeView-WJ3FIDY2.js → GraphGenomeView-6EDKLBOJ.js} +2 -2
- package/dist/chunks/{GraphTrackSettingsDialog-ZW7ALMDV.js → GraphTrackSettingsDialog-KYEZWI4Z.js} +2 -2
- package/dist/chunks/{HighlightColorDialog-VBZRUKD5.js → HighlightColorDialog-IMLYI64R.js} +2 -2
- package/dist/chunks/{LinearGraphDisplay-SPOPTJEN.js → LinearGraphDisplay-CRIZLJDE.js} +2 -2
- package/dist/chunks/{MinigraphBubbleAdapter-ABZPIZYV.js → MinigraphBubbleAdapter-2HPTGZ5N.js} +2 -2
- package/dist/chunks/RgfaTabixAdapter-GNH3QCPZ.js +2 -0
- package/dist/chunks/RgfaTabixAdapter-GNH3QCPZ.js.map +7 -0
- package/dist/chunks/{WalkTabixSyntenyAdapter-ZJJ4SQHH.js → WalkTabixSyntenyAdapter-CUPHY36S.js} +2 -2
- package/dist/chunks/{chunk-44CJ2N3N.js → chunk-4EKOQHUX.js} +2 -2
- package/dist/chunks/{chunk-44CJ2N3N.js.map → chunk-4EKOQHUX.js.map} +4 -4
- package/dist/chunks/{chunk-RFDSEB2X.js → chunk-6DRBCLWY.js} +4 -4
- package/dist/chunks/chunk-B5SVBPXN.js +2 -0
- package/dist/chunks/chunk-B5SVBPXN.js.map +7 -0
- package/dist/chunks/chunk-CJMBJJLO.js +3 -0
- package/dist/chunks/chunk-CJMBJJLO.js.map +7 -0
- package/dist/chunks/{chunk-BXRTG7FU.js → chunk-F4T7JI7S.js} +2 -2
- package/dist/chunks/{chunk-BXRTG7FU.js.map → chunk-F4T7JI7S.js.map} +1 -1
- package/dist/chunks/chunk-GN4HSVLB.js +4 -0
- package/dist/chunks/{chunk-VQ77X5QU.js.map → chunk-GN4HSVLB.js.map} +4 -4
- package/dist/chunks/{chunk-WUSKVEA6.js → chunk-HOSCIP4U.js} +2 -2
- package/dist/chunks/{chunk-WUSKVEA6.js.map → chunk-HOSCIP4U.js.map} +2 -2
- package/dist/chunks/chunk-HTSG4T32.js +3 -0
- package/dist/chunks/chunk-HTSG4T32.js.map +7 -0
- package/dist/chunks/{chunk-YXVEGNFQ.js → chunk-IFBKAC2H.js} +2 -2
- package/dist/chunks/{chunk-RW775LQN.js → chunk-MNL7IPTH.js} +2 -2
- package/dist/chunks/{chunk-RW775LQN.js.map → chunk-MNL7IPTH.js.map} +2 -2
- package/dist/chunks/chunk-N5VWQHTK.js +2 -0
- package/dist/chunks/{chunk-QTOA4RB7.js.map → chunk-N5VWQHTK.js.map} +2 -2
- package/dist/chunks/{chunk-DBPEUVUC.js → chunk-NGR473EL.js} +2 -2
- package/dist/chunks/{chunk-DBPEUVUC.js.map → chunk-NGR473EL.js.map} +2 -2
- package/dist/chunks/chunk-ORHUPW3G.js +3 -0
- package/dist/chunks/chunk-ORHUPW3G.js.map +7 -0
- package/dist/chunks/{chunk-DUR437TH.js → chunk-PUFH6DMJ.js} +2 -2
- package/dist/chunks/{chunk-DUR437TH.js.map → chunk-PUFH6DMJ.js.map} +1 -1
- package/dist/chunks/{chunk-FIRQAKMX.js → chunk-Q55PHRHV.js} +1 -1
- package/dist/chunks/{chunk-FIRQAKMX.js.map → chunk-Q55PHRHV.js.map} +1 -1
- package/dist/chunks/chunk-QQCTNB2Z.js +2 -0
- package/dist/chunks/chunk-QQCTNB2Z.js.map +7 -0
- package/dist/chunks/chunk-R25GOTED.js +2 -0
- package/dist/chunks/chunk-R25GOTED.js.map +7 -0
- package/dist/chunks/chunk-TQKSF3NI.js +2 -0
- package/dist/chunks/{chunk-7HWHWPRW.js.map → chunk-TQKSF3NI.js.map} +3 -3
- package/dist/chunks/{chunk-6MUBJVBG.js → chunk-ULLEUORS.js} +2 -2
- package/dist/chunks/{chunk-6MUBJVBG.js.map → chunk-ULLEUORS.js.map} +4 -4
- package/dist/chunks/chunk-VW4PPXH6.js +2 -0
- package/dist/chunks/{renderSvg-WD6YDV2I.js → renderSvg-HEEQJD5B.js} +2 -2
- package/dist/jbrowse-plugin-graphgenomeviewer.esm.js +2 -2
- package/dist/jbrowse-plugin-graphgenomeviewer.esm.js.map +2 -2
- package/package.json +3 -3
- package/src/RgfaTabixAdapter/RgfaTabixAdapter.ts +3 -1
- package/src/RgfaTabixAdapter/walkReader.ts +5 -1
- package/src/panSNTabix.ts +3 -3
- package/src/version.ts +1 -1
- package/dist/chunks/GbzBaseSyntenyAdapter-NYI2YVNG.js.map +0 -7
- package/dist/chunks/RgfaTabixAdapter-3NQC5QPT.js +0 -2
- package/dist/chunks/RgfaTabixAdapter-3NQC5QPT.js.map +0 -7
- package/dist/chunks/chunk-2M4XYSLG.js +0 -2
- package/dist/chunks/chunk-2M4XYSLG.js.map +0 -7
- package/dist/chunks/chunk-6KUHGRHU.js +0 -2
- package/dist/chunks/chunk-6KUHGRHU.js.map +0 -7
- package/dist/chunks/chunk-7HWHWPRW.js +0 -2
- package/dist/chunks/chunk-BKBH3WGW.js +0 -3
- package/dist/chunks/chunk-BKBH3WGW.js.map +0 -7
- package/dist/chunks/chunk-JCAZTNYG.js +0 -2
- package/dist/chunks/chunk-N74HGPAW.js +0 -3
- package/dist/chunks/chunk-N74HGPAW.js.map +0 -7
- package/dist/chunks/chunk-PEDWEPNK.js +0 -2
- package/dist/chunks/chunk-PEDWEPNK.js.map +0 -7
- package/dist/chunks/chunk-PYQG3JHP.js +0 -3
- package/dist/chunks/chunk-PYQG3JHP.js.map +0 -7
- package/dist/chunks/chunk-QTOA4RB7.js +0 -2
- package/dist/chunks/chunk-VQ77X5QU.js +0 -4
- /package/dist/chunks/{GraphGenomeView-WJ3FIDY2.js.map → GraphGenomeView-6EDKLBOJ.js.map} +0 -0
- /package/dist/chunks/{GraphTrackSettingsDialog-ZW7ALMDV.js.map → GraphTrackSettingsDialog-KYEZWI4Z.js.map} +0 -0
- /package/dist/chunks/{HighlightColorDialog-VBZRUKD5.js.map → HighlightColorDialog-IMLYI64R.js.map} +0 -0
- /package/dist/chunks/{LinearGraphDisplay-SPOPTJEN.js.map → LinearGraphDisplay-CRIZLJDE.js.map} +0 -0
- /package/dist/chunks/{MinigraphBubbleAdapter-ABZPIZYV.js.map → MinigraphBubbleAdapter-2HPTGZ5N.js.map} +0 -0
- /package/dist/chunks/{WalkTabixSyntenyAdapter-ZJJ4SQHH.js.map → WalkTabixSyntenyAdapter-CUPHY36S.js.map} +0 -0
- /package/dist/chunks/{chunk-RFDSEB2X.js.map → chunk-6DRBCLWY.js.map} +0 -0
- /package/dist/chunks/{chunk-YXVEGNFQ.js.map → chunk-IFBKAC2H.js.map} +0 -0
- /package/dist/chunks/{chunk-JCAZTNYG.js.map → chunk-VW4PPXH6.js.map} +0 -0
- /package/dist/chunks/{renderSvg-WD6YDV2I.js.map → renderSvg-HEEQJD5B.js.map} +0 -0
package/package.json
CHANGED
|
@@ -1,6 +1,6 @@
|
|
|
1
1
|
{
|
|
2
2
|
"name": "jbrowse-plugin-graphgenomeviewer",
|
|
3
|
-
"version": "7.3.
|
|
3
|
+
"version": "7.3.2",
|
|
4
4
|
"description": "JBrowse 2 graph genome (pangenome) view: rGFA reference-anchored layouts, the Bandage force-directed engine, and haplotype lanes read from a gbz-base database",
|
|
5
5
|
"license": "GPL-3.0-or-later",
|
|
6
6
|
"repository": {
|
|
@@ -32,7 +32,7 @@
|
|
|
32
32
|
"@jbrowse/core": "5.0.0-beta.11",
|
|
33
33
|
"@jbrowse/display-kit": "5.0.0-beta.11",
|
|
34
34
|
"@jbrowse/display-ui": "5.0.0-beta.11",
|
|
35
|
-
"@jbrowse/graph-walks": "^3.
|
|
35
|
+
"@jbrowse/graph-walks": "^3.1.0",
|
|
36
36
|
"@jbrowse/render-core": "5.0.0-beta.11",
|
|
37
37
|
"@jbrowse/synteny-core": "5.0.0-beta.11"
|
|
38
38
|
},
|
|
@@ -40,7 +40,7 @@
|
|
|
40
40
|
"@emotion/react": "^11.14.0",
|
|
41
41
|
"@emotion/styled": "^11.14.1",
|
|
42
42
|
"@eslint/js": "^10.0.1",
|
|
43
|
-
"@jbrowse/bandage-core": "^11.
|
|
43
|
+
"@jbrowse/bandage-core": "^11.6.1",
|
|
44
44
|
"@jbrowse/cli": "^4.3.0",
|
|
45
45
|
"@jbrowse/mobx-state-tree": "^6.9.0",
|
|
46
46
|
"@jbrowse/plugin-linear-genome-view": "5.0.0-beta.11",
|
|
@@ -90,7 +90,9 @@ export default class RgfaTabixAdapter extends BaseFeatureDataAdapter<RgfaTabixAd
|
|
|
90
90
|
return refNames.regionBytes(
|
|
91
91
|
regions,
|
|
92
92
|
opts,
|
|
93
|
-
walks &&
|
|
93
|
+
walks &&
|
|
94
|
+
((refName, { start, end }) =>
|
|
95
|
+
walks.reader.queryStart(start, { ...opts, refName, end })),
|
|
94
96
|
)
|
|
95
97
|
}
|
|
96
98
|
|
|
@@ -87,7 +87,11 @@ export class WalkIndex {
|
|
|
87
87
|
await getLines(
|
|
88
88
|
this.file,
|
|
89
89
|
refName,
|
|
90
|
-
await this.reader.queryStart(region.start,
|
|
90
|
+
await this.reader.queryStart(region.start, {
|
|
91
|
+
...opts,
|
|
92
|
+
refName,
|
|
93
|
+
end: region.end,
|
|
94
|
+
}),
|
|
91
95
|
region.end,
|
|
92
96
|
{
|
|
93
97
|
signal: opts.signal,
|
package/src/panSNTabix.ts
CHANGED
|
@@ -182,12 +182,12 @@ export class PanSNRefNames {
|
|
|
182
182
|
/**
|
|
183
183
|
* The compressed bytes a read of `regions` fetches from the file, as its
|
|
184
184
|
* index estimates them, for the display's byte gate. `queryStart` moves a
|
|
185
|
-
*
|
|
185
|
+
* window's start to where the read begins.
|
|
186
186
|
*/
|
|
187
187
|
async regionBytes(
|
|
188
188
|
regions: Region[],
|
|
189
189
|
opts: BaseOptions = {},
|
|
190
|
-
queryStart = async (
|
|
190
|
+
queryStart = async (_refName: string, region: Region) => region.start,
|
|
191
191
|
) {
|
|
192
192
|
const spans = []
|
|
193
193
|
for (const region of regions) {
|
|
@@ -195,7 +195,7 @@ export class PanSNRefNames {
|
|
|
195
195
|
if (refName !== undefined) {
|
|
196
196
|
spans.push({
|
|
197
197
|
refName,
|
|
198
|
-
start: await queryStart(region
|
|
198
|
+
start: await queryStart(refName, region),
|
|
199
199
|
end: region.end,
|
|
200
200
|
})
|
|
201
201
|
}
|
package/src/version.ts
CHANGED
|
@@ -1 +1 @@
|
|
|
1
|
-
export const version = '7.3.
|
|
1
|
+
export const version = '7.3.2'
|
|
@@ -1,7 +0,0 @@
|
|
|
1
|
-
{
|
|
2
|
-
"version": 3,
|
|
3
|
-
"sources": ["../../node_modules/.pnpm/@jbrowse+bandage-core@11.4.0/node_modules/@jbrowse/bandage-core/dist/chunks/chunk-TUQKGAK3.js", "../../src/GbzBaseSyntenyAdapter/companion.ts", "../../src/gaf/gafFile.ts", "../../src/GbzBaseSyntenyAdapter/GbzBaseSyntenyAdapter.ts"],
|
|
4
|
-
"sourcesContent": ["import {\n GBZ_CUT_DEFAULTS,\n NoReferenceSampleError,\n aliasedSample,\n cutWindowGFA,\n haplotypeWanted,\n referencePathQuery,\n referenceSamplesOf\n} from \"./chunk-NMWQUFPF.js\";\nimport {\n panSNSample\n} from \"./chunk-4XWZNPQ2.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 gbzLocations({ db, index }) {\n return db === \"hprc\" ? { db: HPRC_GBZ.db, index: index ?? HPRC_GBZ.index } : { db, index };\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 if (ref) {\n return aliasedSample(panSNSample(ref), samples) ?? ref;\n }\n const first = samples[0];\n if (first === void 0) {\n throw new NoReferenceSampleError(\"\", samples);\n }\n return first;\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 gbzLocations,\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/walkCut.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 { laneZoomNotice } from '../synteny/laneZoomNotice.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\nconst laneLimitNotice = (error: unknown, nodeLimit: number, windowBp: number) =>\n laneZoomNotice(\n zoomInNotice(nodeLimitError(error, nodeLimit, windowBp) ?? error, 'lanes'),\n windowBp,\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 laneLimitNotice(error, nodeLimit, end - start)\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 laneLimitNotice(error, nodeLimit, end - start)\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": "g2BAsBA,SAASA,EAAqBC,EAAI,CAChC,IAAMC,EAAS,cACf,OAAOA,EAAO,KAAKD,CAAE,EAAIA,EAAG,QAAQC,EAAQ,qBAAqB,EAAI,MACvE,CCjBO,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,EC3DO,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,GAAaC,EAAgBC,EAAc,CAClD,OAAID,aAAiBE,IACnBF,EAAM,QAAU,oBAAoBG,EAAgBH,EAAM,MAAM,CAAC,WAAWC,CAAI,IAE3ED,CACT,CAEA,IAAMI,GAAkB,CAACJ,EAAgBK,EAAmBC,IAC1DC,GACER,GAAaS,EAAeR,EAAOK,EAAWC,CAAQ,GAAKN,EAAO,OAAO,EACzEM,CACF,EAGIG,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,EACKzD,EAAS,KAAK,QAAQ,eAAe,EAAE,CAAC,EAC9C,GAAIA,IAAW,OACb,MAAM,IAAI,MACR,8GACF,EAEF,IAAM6D,EAAmB,MAAMC,EAAmBJ,CAAE,EAC9CK,EAAkBC,EAAuB,CAC7C,WAAY,KAAK,QAAQ,iBAAiB,EAC1C,aAAcC,EAAmB,KAAMjE,CAAM,EAC7C,iBAAA6D,CACF,CAAC,EACD,MAAO,CACL,GAAAH,EACA,OAAA1D,EACA,gBAAA+D,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,EAAc,KAAM,CAAC,eAAe,EAAG,CAAC,YAAY,CAAC,EACrD,MACN,EACA,MACN,CACF,CAAC,EAtKD,MAAc,eACZ/D,EACAgE,EACgC,CAChC,GAAM,CAAE,GAAAZ,EAAI,gBAAAK,CAAgB,EAAI,MAAM,KAAK,MAAMO,CAAI,EACrD,OAAOC,EAAmBb,EAAIK,EAAiBzD,CAAO,CACxD,CAOA,MAAM,cAAcgE,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,IAAM7B,EAASzB,EAAgBqE,EAAK,IAAI,EAClCC,EAAQF,EAAS,IAAI3C,CAAM,EAC7B6C,EACGA,EAAM,QAAQ,SAASD,EAAK,KAAK,MAAM,GAC1CC,EAAM,QAAQ,KAAKD,EAAK,KAAK,MAAM,EAGrCD,EAAS,IAAI3C,EAAQ,CACnB,OAAAA,EACA,OAAQ4C,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,OAAAtE,EAAQ,gBAAA+D,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,mBAAoB9E,EACpB,MAAA6E,CACF,CACF,CAEA,MAAM,YAAYP,EAAoB,CAAC,EAAG,CACxC,GAAM,CAAE,GAAAZ,EAAI,OAAA1D,EAAQ,gBAAA+D,CAAgB,EAAI,MAAM,KAAK,MAAMO,CAAI,EACvD,CAAE,aAAAvE,CAAa,EAAIuE,EACnBzC,EACJ9B,IAAiBC,EACb,OACAiE,EAAmB,KAAMlE,CAAY,EACrCiF,GAAW,MAAMtB,EAAG,MAAM,GAC7B,OAAOe,GACN1E,IAAiBC,EACbyE,EAAK,WAAaA,EAAK,KAAK,SAAWV,EACvCkB,EAAmB7E,EAAgBqE,EAAK,IAAI,EAAG5C,CAAM,CAC3D,EACC,IAAI4C,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,IAAI3E,GAAQ0D,EAAmB,KAAM1D,CAAI,CAAC,EAC3D,OAAO6E,IAAmB,QAAaD,IAAiB,OACpD,OACCjF,GACCmF,EAAgBnF,EAAMkF,CAAc,IACnCD,IAAiB,QAChBF,EAAmB7E,EAAgBF,CAAI,EAAGiF,CAAY,EAChE,CAeA,MAAM,YAAYnC,EAAgBsB,EAA+B,CAAC,EAAG,CACnE,GAAM,CAAE,OAAAgB,CAAO,EAAIhB,EACb,CAAE,GAAAZ,EAAI,OAAA1D,EAAQ,gBAAA+D,EAAiB,gBAAAwB,CAAgB,EAAI,MAAM,KAAK,MAAM,CACxE,OAAAD,CACF,CAAC,EACK,CAAE,aAAAvF,EAAc,QAAAO,EAAS,MAAAqB,EAAO,IAAAC,CAAI,EAAImB,GAAQC,CAAM,EAC5D,GAAIjD,IAAiBC,EACnB,MAAM,IAAIF,EAAqBC,EAAcC,CAAM,EAErD,IAAMwF,EAAQ,MAAM,KAAK,eAAelF,EAAS,CAAE,OAAAgF,CAAO,CAAC,EAC3D,GAAI,CAACE,EACH,MAAM,IAAI,MACR,GAAGzB,CAAe,wBAAwBzD,CAAO,gBACnD,EAEF,IAAMmF,EAAO,KAAK,cAAcnB,EAAK,UAAU,EAC/C,GAAImB,IAAS,QAAa,CAAC/B,EAAG,kBAC5B,MAAM,IAAI9D,EAAsB2F,CAAe,EAEjD,IAAMG,EAAM,MAAMC,EAAgBjC,EAAI8B,EAAO7D,EAAOC,EAAK,CACvD,QAAS,KAAK,QAAQ,SAAS,EAC/B,OAAQ0C,EAAK,QAAU,KAAK,QAAQ,gBAAgB,EACpD,MAAO,KAAK,QAAQ,WAAW,EAC/B,OAAQA,EAAK,OACb,KAAAmB,CACF,CAAC,EAAE,MAAOnD,GAAmB,CAC3B,MAAMD,GAAaC,EAAO,WAAW,CACvC,CAAC,EACD,OAAO,OAAOoD,GAAQ,SAClBA,EACAE,GAAWF,EAAK,CAAE,MAAA/D,EAAO,IAAAC,CAAI,EAAG0C,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,eACZvE,EACAuE,EACyB,CACzB,IAAMzC,EAASoC,EAAmB,KAAMlE,CAAY,EACpD,OAAQ,MAAM,KAAK,cAAcuE,CAAI,GAAG,OAAOQ,GAC7CG,EAAmBH,EAAU,OAAQjD,CAAM,CAC7C,CACF,CAEA,MAAc,eAAemB,EAAgBsB,EAA0B,CACrE,GAAM,CAAE,GAAAZ,CAAG,EAAI,MAAM,KAAK,MAAMY,CAAI,EAC9B,CAAE,aAAAvE,EAAc,QAAAO,EAAS,MAAAqB,EAAO,IAAAC,CAAI,EAAIoB,EACxCmC,EAAelB,EAAmB,KAAMK,EAAK,kBAAkB,EAC/DK,EAAcC,EAAsB,IAAI,EACxCa,EAAO,KAAK,cAAcnB,EAAK,WAAYa,CAAY,EACvDK,EAAQ,MAAM,KAAK,eAAelF,EAASgE,CAAI,EAC/C3B,EAAoB,KAAK,QAAQ,WAAW,EAsBlD,OArBmB6C,EACf,MAAMO,EACJ,iBAAiBzF,CAAO,IAAIqB,EAAM,eAAe,CAAC,IAAIC,EAAI,eAAe,CAAC,GAC1E0C,EAAK,eACL,IACEZ,EACG,cAAc,CACb,KAAM8B,EACN,MAAA7D,EACA,IAAAC,EACA,QAAS,KAAK,QAAQ,SAAS,EAC/B,WAAY,MACZ,MAAOe,EACP,OAAQ2B,EAAK,OACb,KAAAmB,CACF,CAAC,EACA,MAAOnD,GAAmB,CACzB,MAAMI,GAAgBJ,EAAOK,EAAWf,EAAMD,CAAK,CACrD,CAAC,CACP,EACA,CAAC,GACa,QAAQxB,GAAa,CACrC,IAAM6F,EAAU7F,EAAU,SACtBE,GAAgB,CACd,UAAAF,EACA,aAAAJ,EACA,QAAAO,EACA,KAAMyE,EACJJ,EACAvE,EAAgBD,EAAU,IAAI,CAChC,CACF,CAAC,EACD,OACJ,OAAO6F,IAAY,OAAY,CAAC,EAAI,CAACA,CAAO,CAC9C,CAAC,CACH,CAWA,MAAc,aACZhD,EACAiD,EACA3B,EACA,CACA,GAAM,CAAE,GAAAZ,CAAG,EAAI,MAAM,KAAK,MAAMY,CAAI,EAC9B,CAAE,QAAAhE,EAAS,MAAAqB,EAAO,IAAAC,CAAI,EAAIoB,EAQ1BkD,GAPQ,MAAM,QAAQ,IAC1BD,EAAM,IAAI,MAAM/E,IAAS,CACvB,KAAAA,EACA,YAAa,MAAM,KAAK,eAAeA,EAAK,kBAAmBoD,CAAI,EACnE,SAAU,MAAM,KAAK,eAAepD,EAAK,mBAAoBoD,CAAI,CACnE,EAAE,CACJ,GACyB,OACvB,CAAC,CAAE,YAAA6B,EAAa,SAAAC,CAAS,IACvBD,EAAY,OAAS,GAAKC,EAAS,OAAS,CAChD,EACM5E,EAAO,IAAI,IACf0E,EAAW,QAAQ,CAAC,CAAE,YAAAC,EAAa,SAAAC,CAAS,IAC1C,CAAC,GAAGD,EAAa,GAAGC,CAAQ,EAAE,IAAIhG,CAAe,CACnD,CACF,EACMoF,EAAQ,MAAM,KAAK,eAAelF,EAASgE,CAAI,EAC/C3B,EAAoB,KAAK,QAAQ,WAAW,EAC5C0D,EAAsB,KAAK,QAAQ,SAAS,EAC5CX,EAAM,CAACjB,EAAiB6B,IAC5B5C,EAAG,aAAa,CACd,KAAAe,EACA,MAAA9C,EACA,IAAAC,EACA,QAAA0E,EACA,OAAQ,KAAK,QAAQ,gBAAgB,EACrC,WAAY,MACZ,MAAO3D,EACP,OAAQ2B,EAAK,OACb,KAAMpE,GAAQsB,EAAK,IAAIpB,EAAgBF,CAAI,CAAC,CAC9C,CAAC,EAEGqG,EAAU,MAAO9B,GAAoB,CACzC,IAAMlD,EAAY,MAAMmE,EAAIjB,EAAM4B,CAAW,EAAE,MAAO/D,GAAmB,CACvE,MAAMI,GAAgBJ,EAAOK,EAAWf,EAAMD,CAAK,CACrD,CAAC,EACKS,EAAMd,GAAcC,EAAWC,CAAI,EACnC8E,EAAU,KAAK,KAAKlE,EAAM,CAAC,EAAI,IACrC,OAAOA,EAAM,GAAKkE,EAAUD,EACxBX,EAAIjB,EAAM6B,CAAO,EAAE,MAAOhE,GAAmB,CAC3C,GAAIA,aAAiBkE,EACnB,OAAOjF,EAET,MAAMe,CACR,CAAC,EACDf,CACN,EAOA,OALEiE,GAASU,EAAW,OAAS,EACzB,MAAMH,EAAaU,GAAgBR,CAAK,EAAG3B,EAAK,eAAgB,IAC9DiC,EAAQf,CAAK,CACf,EACA,CAAC,GACU,QAAQ9D,GACvBwE,EAAW,QAAQ,CAAC,CAAE,KAAAhF,EAAM,YAAAiF,EAAa,SAAAC,CAAS,IAChDD,EAAY,QAAQO,GAClBN,EAAS,QAAQO,IACfjF,EACG,eAAe,CAAE,OAAAgF,EAAQ,MAAOC,GAAM,MAAO,EAAM,CAAC,EACpD,IAAIxG,IACHc,GAAY,CACV,KAAMd,GACN,KAAMe,EAAK,kBACX,SAAUA,EAAK,kBACjB,CAAC,CACH,CACJ,CACF,CACF,CACF,CACF,CAEA,YAAY0F,EAAetC,EAA2B,CAAC,EAAG,CACxD,IAAMtB,EAASD,GAAQ6D,CAAK,EAC5B,OAAOC,GAA0B,MAAMC,GAAY,CACjD,GAAM,CAAE,GAAApD,EAAI,OAAA1D,EAAQ,gBAAAuF,CAAgB,EAAI,MAAM,KAAK,MAAMjB,CAAI,EAC7D,GAAI,CAACZ,EAAG,kBACN,MAAM,IAAI9D,EAAsB2F,CAAe,EAIjD,GAAIvC,EAAO,eAAiBhD,EAAQ,CAClC,IAAMiG,EAAQc,GAAezC,CAAI,EAC3B0C,EACJf,EAAM,OAAS,EACX,MAAM,KAAK,aAAajD,EAAQiD,EAAO3B,CAAI,EAC3C,MAAM,KAAK,eAAetB,EAAQsB,CAAI,EAC5C,QAAW0B,KAAWgB,EACpBF,EAAS,KAAKd,CAAO,CAEzB,CACAc,EAAS,SAAS,CACpB,EAAGxC,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", "laneLimitNotice", "nodeLimit", "windowBp", "laneZoomNotice", "nodeLimitError", "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", "updateStatus", "feature", "pairs", "answerable", "featureSide", "mateSide", "baseContext", "context", "bridged", "SubgraphLimitError", "lanePairsStatus", "target", "mate", "asked", "ObservableCreate", "observer", "requestedPairs", "features"]
|
|
7
|
-
}
|
|
@@ -1,2 +0,0 @@
|
|
|
1
|
-
import{a as b}from"./chunk-QTOA4RB7.js";import"./chunk-JCAZTNYG.js";import{a as N}from"./chunk-KOB7YIMR.js";import"./chunk-YBLNN5FU.js";import"./chunk-6KUHGRHU.js";import"./chunk-6MUBJVBG.js";import{b as g,h as c,i as u,j as h,k as w,p as m,q as k,r as x}from"./chunk-BKBH3WGW.js";import"./chunk-YBLNN5FU.js";import"./chunk-MVVXRZXV.js";import"./chunk-7HWHWPRW.js";import{d}from"./chunk-IKJZPPOA.js";import"./chunk-OGQ6OQDN.js";import"./chunk-NDRLCGPC.js";import{a as p,o as f}from"./chunk-TENPGIPX.js";var l=class extends g{constructor(){super(...arguments);this.fine=this.openIndex([])}static{this.capabilities=["getFeatures","getRefNames"]}openIndex(e){let t=m(this,[...e,"segmentsLocation"],[...e,"segmentsIndex"]),s=m(this,[...e,"linksLocation"],[...e,"linksIndex"]),n=this.getConf("walksLocation"),i=e.length===0&&d(n)?new b(m(this,["walksLocation"],["walksIndex"]),this):void 0;return{segments:t,links:s,walks:i,refNames:new x(i?.file??t,this)}}async getHaplotypeNames(e={}){return this.fine.walks?.haplotypeNames(e)}index(e="fine"){return e==="fine"?this.fine:(this.coarse??=this.openIndex(["coarse"]),this.coarse)}async getRefNames(e={}){return this.fine.refNames.assemblyRefNames(e)}async getRegionByteSize(e,t={}){let{walks:s,refNames:n}=this.fine;return n.regionBytes(e,t,s&&(i=>s.reader.queryStart(i,t)))}async hasDataForRefName(){return!0}getFeatures(e,t={}){return N(async s=>{let n=await this.fine.refNames.resolveMapped(e,t);if(n!==void 0){let{walks:i}=this.fine,a=r=>{s.next(r)};await(i?this.walkFeatures(i,n,e,t,a):this.segmentFeatures(n,e,t,a))}s.complete()},t.signal)}async walkFeatures(e,t,s,n,i){let a=await f("Downloading walks",n.statusCallback,()=>e.referenceNodes(t,s,n));for(let r of a){let o=String(r.id);i(new p({uniqueId:o,refName:s.refName,start:r.start,end:r.end,name:o,type:"segment",rank:0,stableName:r.refName}))}}async segmentFeatures(e,t,s,n){await f("Downloading segments",s.statusCallback,()=>w(this.fine.segments,e,t.start,t.end,{signal:s.signal,lineCallback:i=>{let{segment:a,anchored:r}=c(i);if(r&&a.refName!==e)return;let o=u(a);n(new p({uniqueId:a.id,refName:t.refName,start:a.start,end:a.end,name:a.id,type:"segment",rank:a.rank,stableName:a.refName,...o&&{samples:o,sampleCount:o.length}}))}}))}async getSubgraph(e,t={}){let{hops:s,signal:n,tier:i}=t,a=this.index(i),r=await a.refNames.resolve(e,{signal:n});if(r===void 0)throw new Error(k(e));return a.walks?a.walks.subgraph(r,e,t):h(a,{refName:r,start:e.start,end:e.end},{hops:s,signal:n})}};export{l as default};
|
|
2
|
-
//# sourceMappingURL=RgfaTabixAdapter-3NQC5QPT.js.map
|
|
@@ -1,7 +0,0 @@
|
|
|
1
|
-
{
|
|
2
|
-
"version": 3,
|
|
3
|
-
"sources": ["../../src/RgfaTabixAdapter/RgfaTabixAdapter.ts"],
|
|
4
|
-
"sourcesContent": ["import { BaseFeatureDataAdapter } from '@jbrowse/core/data_adapters/BaseAdapter'\nimport { SimpleFeature, updateStatus } from '@jbrowse/core/util'\nimport { ObservableCreate } from '@jbrowse/core/util/rxjs'\nimport { parseSegmentRow, rgfaCut, segmentSamples } from '@jbrowse/graph-walks'\n\nimport { locationIsSet } from '../locationName.ts'\nimport {\n PanSNRefNames,\n openTabixSlot,\n unmappedRegionMessage,\n} from '../panSNTabix.ts'\nimport { getLines } from '../tabixRange.ts'\nimport { WalkIndex } from './walkReader.ts'\n\nimport type { RgfaTabixAdapterConfig } from './configSchema.ts'\nimport type { SubgraphAdapterOptions, SubgraphTier } from '../GetSubgraph.ts'\nimport type { WalkCut } from './walkCut.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.resolveMapped(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(unmappedRegionMessage(region))\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": "ufA6BA,IAAqBA,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,cAAcH,EAAOL,CAAI,EACvE,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,MAAMmB,EAAsBJ,CAAM,CAAC,EAE/C,OAAIG,EAAM,MACDA,EAAM,MAAM,SAASlB,EAAce,EAAQvB,CAAI,EAEjD4B,EACLF,EACA,CAAE,QAASlB,EAAc,MAAOe,EAAO,MAAO,IAAKA,EAAO,GAAI,EAC9D,CAAE,KAAAC,EAAM,OAAAC,CAAO,CACjB,CACF,CACF",
|
|
6
|
-
"names": ["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", "unmappedRegionMessage", "rgfaCut"]
|
|
7
|
-
}
|
|
@@ -1,2 +0,0 @@
|
|
|
1
|
-
import{B as Yt,C as jt,D as Et,E as xt,M as Vt,N as lt,a as Gt,k as Wt,l as gt,n as Bt,o as $t,p as Xt,q as Kt,u as Ut}from"./chunk-PEDWEPNK.js";import{a as J,g as Q}from"./chunk-6MUBJVBG.js";function p(t,e,n,o){return(o<<24|n<<16|e<<8|t)>>>0}function Mt(t){return t&255}function _t(t){return t>>>8&255}function Rt(t){return t>>>16&255}function yt(t){return t>>>24&255}function zt(t){return`rgba(${Mt(t)},${_t(t)},${Rt(t)},${yt(t)/255})`}function me(t,e){return`rgba(${Math.round(t[0]*255)},${Math.round(t[1]*255)},${Math.round(t[2]*255)},${e})`}function j(t,e){return p(Mt(t),_t(t),Rt(t),Math.round(yt(t)*e))}function pe(t,e){let n=Math.max(0,1-1/e),o=r=>Math.round(r+(255-r)*n);return p(o(Mt(t)),o(_t(t)),o(Rt(t)),yt(t))}var qt=300,Z=.7,tt=.5;function An(t){return`hsl(${t}, ${Z*100}%, ${tt*100}%)`}function nt(t){return Math.max(0,Math.min(1,t))*qt}function At(t,e){return nt((t-e.start)/(e.end-e.start))}var bn=30;function Ee({start:t,bp:e,reversed:n},o){let r=At(t,o),s=At(t+e,o);if(s-r<1)return[(r+s)/2];let a=Math.ceil((s-r)/bn);return Array.from({length:a+1},(c,i)=>At(t+e*(n?a-i:i)/a,o))}var vn=Array.from({length:7},(t,e)=>({offset:e/6,color:An(e/6*qt)})),xe=`linear-gradient(to right, ${vn.map(t=>`${t.color} ${t.offset*100}%`).join(", ")})`;var Cn=Gt();function Jt(t,e,n,o,r,s,a=!1){return Cn(t,[e,n.scaleX,n.scaleY,o,s,a,r],()=>{let c=new Map;for(let i=0;i<e.edges.length;i++){let l=o?.get(i),d=Ut(t,e.edges[i],n,l,l&&s?.[i],a);d&&c.set(i,d)}return c})}var Sn=4;function In(t){let e=new Map;for(let n of t.edges)J(e,n.from,n.to),J(e,n.to,n.from);return e}function wn({nodeId:t,nodeById:e,neighbors:n}){let o=new Set([t]),r=[t],s=1/0,a=-1/0;for(let c=0;c<Sn&&r.length>0;c++){let i=[];for(let l of r){let d=n.get(l);if(d){for(let m of d)if(!o.has(m)){o.add(m);let h=e.get(m);h&&Q(h)?(s=Math.min(s,h.stable.start+h.length),a=Math.max(a,h.stable.start)):i.push(m)}}}r=i}return s===1/0?void 0:{start:Math.min(s,a),end:Math.max(s,a)}}function Tn({nodeId:t,nodeById:e,neighbors:n}){let o=e.get(t);return o===void 0?void 0:Q(o)?{start:o.stable.start,end:o.stable.start+o.length}:wn({nodeId:t,nodeById:e,neighbors:n})}function Qt(t){let e=new Map(t.nodes.map(r=>[r.id,r])),n=In(t),o=new Map;for(let r of t.nodes){let s=Tn({nodeId:r.id,nodeById:e,neighbors:n});s&&o.set(r.id,(s.start+s.end)/2)}return o}var Nn=p(119,119,119,217),Zt=p(136,136,136,217),Ln=p(24,24,28,240),On=2.2,cn=p(160,160,160,255),et=.18,Pe=zt(j(cn,et)),kn=11,Pn=1.2,tn=4,Dn=1e-6,Fn=7,Hn=3,Gn=3,Wn=1.5,nn=3,Bn=.5,$n=.45,bt=12,Xn=1e-6,w=new Float64Array(64),D=new Float64Array(64),F=new Float64Array(64);function Kn(t,e,n,o){let r=t.length*bt+1;w.length<r&&(w=new Float64Array(r),D=new Float64Array(r),F=new Float64Array(r));let s=t[t.length-1];w[0]=s.x1,D[0]=s.y1,F[0]=0;let a=1,c=0;for(let i=t.length-1;i>=0&&c<o;i--){let l=t[i];for(let d=1;d<=bt&&c<o;d++){let m=1-(d/bt)**2,h=1-m,_=h*h*h,M=3*h*h*m,T=3*h*m*m,A=m*m*m,N=_*l.x0+M*l.cx0+T*l.cx1+A*l.x1,L=_*l.y0+M*l.cy0+T*l.cy1+A*l.y1;c+=Math.hypot(N-w[a-1],(L-D[a-1])*n)*e,w[a]=N,D[a]=L,F[a]=c,a++}}return a}function en(t,e){let n=1;for(;n<t-1&&F[n]<e;)n++;let o=F[n-1],r=F[n],s=r>o?Math.min(1,Math.max(0,(e-o)/(r-o))):0;return{x:w[n-1]+(w[n]-w[n-1])*s,y:D[n-1]+(D[n]-D[n-1])*s}}function Un({curves:t,edgeHalfWidth:e,nodeHalfWidth:n,scale:o,yToX:r,color:s,fromSegments:a,toSegments:c}){let i=Fn+Hn*e,l=Kn(t,o,r,n+nn*i),d=F[l-1],m=d<Xn,h=Math.min(n,d),_=m?Math.min(1,Kt(a,o,o*r)/i):1;if(m?_<Bn:d<nn*i)return;let M=en(l,h),T=en(l,h+i*_),A=M.x-T.x,N=(M.y-T.y)*r,L=t[t.length-1];return{x:M.x,y:M.y,angle:Math.hypot(A,N)>0?Math.atan2(N,A):te(L,r,jn(L,c)),length:i*_,halfWidth:(Gn+Wn*e)*_,color:s}}function Yn(t,e,n,o,r){if(o)return Bt(t).from==="start";let s=a=>Math.hypot(a.x-n.x0,(a.y-n.y0)*r);return s(e[0])<s(e.at(-1))}function jn(t,e){if(e.length<2)return;let n=e[0],o=e[e.length-1];return Math.hypot(t.x1-n.x,t.y1-n.y)<=Math.hypot(t.x1-o.x,t.y1-o.y)?{from:n,to:e[1]}:{from:o,to:e[e.length-2]}}function on(t,e){if(!t||!e)return;let n=e.x-t.x,o=e.y-t.y,r=Math.hypot(n,o);return r>0?{nx:-o/r,ny:n/r}:void 0}function vt(t,e=1){let n=e===1?t:t.map(o=>({x:o.x,y:o.y*e}));return n.map((o,r)=>{let s=on(n[r-1],o),a=on(o,n[r+1]);if(s&&a){let c=(s.nx+a.nx)/2,i=(s.ny+a.ny)/2,l=s.nx*c+s.ny*i;return l>.1?{nx:c/l,ny:i/l}:{nx:c,ny:i}}return s??a??{nx:0,ny:0}})}function rn(t,e,n,o=1){return t.map((r,s)=>({x:r.x+e[s].nx*n,y:r.y+e[s].ny*n/o}))}function St(t,e,n){let o=(1-Math.abs(2*n-1))*e,r=o*(1-Math.abs(t/60%2-1)),s=n-o/2,a=0,c=0,i=0;return t<60?(a=o,c=r):t<120?(a=r,c=o):t<180?(c=o,i=r):t<240?(c=r,i=o):t<300?(a=r,i=o):(a=o,i=r),[a+s,c+s,i+s]}function ft(t,e,n,o=1){let[r,s,a]=St(t,e,n);return p(Math.round(r*255),Math.round(s*255),Math.round(a*255),Math.round(o*255))}var Vn=[[68,1,84],[59,82,139],[33,145,140],[94,201,98],[253,231,37]],zn=[[237,137,44],[158,42,122]],qn=[[220,50,50],[50,120,220]];function Ct(t,e){let n=Math.max(0,Math.min(1,e))*(t.length-1),o=Math.min(t.length-2,Math.floor(n)),r=n-o,s=t[o],a=t[o+1];return p(Math.round(s[0]+(a[0]-s[0])*r),Math.round(s[1]+(a[1]-s[1])*r),Math.round(s[2]+(a[2]-s[2])*r),255)}function Jn(t,e){return ft(Yt(t),Et,xt,e)}var ln=[60,65,72],It=p(...ln,255),De=`rgb(${ln.join(", ")})`;function Fe(t,e){let n=Qt(t);if(e&&e.end>e.start)return{start:e.start,span:e.end-e.start,midpoints:n};let o=1/0,r=-1/0;for(let s of n.values())o=Math.min(o,s),r=Math.max(r,s);return{start:o===1/0?0:o,span:r>o?r-o:1,midpoints:n}}function Qn(t,e){let n=1/0,o=-1/0,r=1/0,s=-1/0,a=0;for(let c of t.nodes)c.stable!==void 0&&c.stable.rank>a&&(a=c.stable.rank),c.depth<n&&(n=c.depth),c.depth>o&&(o=c.depth),c.length<r&&(r=c.length),c.length>s&&(s=c.length);return{minDepth:n,maxDepth:o,minLength:r,maxLength:s,maxRank:a,nodeCount:t.nodes.length,referenceRamp:e}}function Zn(t,e,n,o){switch(n){case"random":return Jn(t.id,1);case"depth":return Ct(Vn,o.maxDepth>o.minDepth?(t.depth-o.minDepth)/(o.maxDepth-o.minDepth):.5);case"node-length":return Ct(qn,o.maxLength>o.minLength?(t.length-o.minLength)/(o.maxLength-o.minLength):.5);case"stable-rank":return t.stable===void 0?p(160,160,160,255):t.stable.rank===0?p(52,152,219,255):Ct(zn,o.maxRank>1?(t.stable.rank-1)/(o.maxRank-1):0);case"reference-position":{let r=o.referenceRamp;if(r===void 0)return p(160,160,160,255);let s=r.midpoints.get(t.id);return!Q(t)||s===void 0?It:ft(nt((s-r.start)/r.span),Z,tt)}case"grey":return p(160,160,160,255);case"rainbow":{let r=o.nodeCount>1?e/o.nodeCount*360:0;return ft(r,.75,.5)}default:return p(52,152,219,255)}}function te(t,e=1,n){let o=[[t.cx1,t.cy1],[t.cx0,t.cy0],[t.x0,t.y0]];for(let[r,s]of o){let a=t.x1-r,c=(t.y1-s)*e;if(Math.hypot(a,c)>0)return Math.atan2(c,a)}return n?Math.atan2((n.to.y-n.from.y)*e,n.to.x-n.from.x):0}function ne(t,e){let n=!1;for(let o=0;o<t.length-1&&!n;o++){let r=t[o],s=t[o+1];n=Math.max(r.x,s.x)>=e.minX&&Math.min(r.x,s.x)<=e.maxX&&Math.max(r.y,s.y)>=e.minY&&Math.min(r.y,s.y)<=e.maxY}return n}function ee(t,e){for(let n of t){let o=Math.min(n.x0,n.cx0,n.cx1,n.x1),r=Math.max(n.x0,n.cx0,n.cx1,n.x1),s=Math.min(n.y0,n.cy0,n.cy1,n.y1),a=Math.max(n.y0,n.cy0,n.cy1,n.y1);if(r>=e.minX&&o<=e.maxX&&a>=e.minY&&s<=e.maxY)return!0}return!1}var sn=new WeakMap;function oe(t){let e=sn.get(t);if(!e){let n=[],o=new Map,r=t.paths??[];for(let s=0;s<r.length;s++)n.push(ft(jt(s,r.length),Et,xt,.85)),o.has(r[s].name)||o.set(r[s].name,s);e={pathColorByIndex:n,pathIndexByName:o},sn.set(t,e)}return e}var an=new WeakMap,re=new Map;function se(t){let e=an.get(t);if(!e){e=new Map;for(let[n,o]of(t.paths??[]).entries())for(let r of new Set(o.nodeIds))J(e,r,n);an.set(t,e)}return e}function He(t){let{nodePositions:e,graph:n,nodeById:o,colorScheme:r,contigThickness:s,connectorThickness:a,drawPaths:c,nodeWidth:i="uniform",highlight:l,axis:d,linearLayout:m,viewportBounds:h,referenceRamp:_,deletions:M,deletionRoutes:T,stranded:A,hiddenEdges:N,version:L=0}=t,ut=i==="depth"?Vt(n):0,ot=d.scaleX,dt=1/Math.max(ot,Dn),b=$t(d),gn=Jt(e,n,d,M,L,T,A),H=[],Nt=new Map,rt=[],Lt=new Map,G=[],Ot=new Map,En=Qn(n,_),kt=new Map;if(r==="rainbow")for(let f=0;f<n.nodes.length;f++)kt.set(n.nodes[f].id,f);let mt=n.paths?.length??0,{pathColorByIndex:Pt,pathIndexByName:xn}=oe(n),Mn=f=>Pt[xn.get(f)??-1]??Zt,_n=c&&mt>1?se(n):re,Rn=ot>(m?1:$n);for(let f=0;f<n.edges.length;f++){let u=n.edges[f],R=e[u.from],v=e[u.to],g=gn.get(f);if(!R?.length||!v?.length||!g||N?.has(f))continue;let W=u.pathIds?.length??0,V=M?.get(f)!==void 0,B=l?.edgeIndexes.has(f)??!1,$=a/2*(V?On:1),X=c&&W>0,O=V?Ln:Nn,C=!l||B?O:j(O,et),z=o.get(u.to),K=(z?lt(z,s,i,ut):s)/2,S=x=>!h||ee(x,h),k=(x,y,E)=>{if(V)for(let P of Wt(x,kn/ot))G.push({curves:P,thickness:y,color:E});else G.push({curves:x,thickness:y,color:E})},I=x=>{let y=Rn?Un({curves:x,edgeHalfWidth:a/2,nodeHalfWidth:K,scale:ot,yToX:b,color:C,fromSegments:R,toSegments:v}):void 0;y&&rt.push(y)},U=G.length,at=rt.length,ct=l?.walks.length&&B&&!X?l.walks:void 0;if(ct){let x=o.get(u.from),y=x?lt(x,s,i,ut):s,E=Math.max(y/ct.length,tn),P=vt(R,b),{nx:it,ny:q}=(Yn(u,R,g[0],A,b)?P[0]:P.at(-1))??{nx:0,ny:0};ct.forEach((Y,yn)=>{if(!Y.edgeIndexes.has(f))return;let pt=(yn-(ct.length-1)/2)*E*dt,Ht=pt===0?g:gt(g,it*pt,q*pt/b);S(Ht)&&k(Ht,E/2,Y.colors.get(u.from)??Y.colors.get(u.to)??O)}),G.length>U&&I(g)}else if(!X)S(g)&&(k(g,$,C),I(g));else{let x=Xt(R,v,W,d),y=Math.floor((W-1)/2);for(let E=0;E<W;E++){let{x:P,y:it}=x[E],q=P===0&&it===0?g:gt(g,P,it);if(!S(q))continue;let Y=Mn(u.pathIds[E]);k(q,$,l&&(l.walks.length?!l.names.has(u.pathIds[E]):!B)?j(Y,et):Y),E===y&&I(q)}}let Dt=G.length-U;Dt>0&&Ot.set(f,{start:U,count:Dt});let Ft=rt.length-at;Ft>0&&Lt.set(f,{start:at,count:Ft})}for(let[f,u]of Object.entries(e)){let R=o.get(f);if(!R||h&&!ne(u,h))continue;let v=l!==void 0&&!l.nodeIds.has(f),g=v?cn:Zn(R,kt.get(f)??0,r,En),W=v?j(g,et):g,st=lt(R,s,i,ut),V=st/2,B=H.length,$=_n.get(f),X=st/mt,O=u.length>=2,C=l?.walks.length&&!v&&!$?.length?l.walks:void 0;if(O&&C){let K=Math.max(st/C.length,tn),S=vt(u,b);C.forEach((k,I)=>{let U=k.colors.get(f);if(U===void 0)return;let at=(I-(C.length-1)/2)*K*dt;H.push({points:C.length===1?u:rn(u,S,at,b),thickness:K/2,color:U})})}else if(O&&$?.length&&X>=Pn){let K=vt(u,b);for(let S of $){let k=(S-(mt-1)/2)*X*dt,I=Pt[S]??Zt;H.push({points:rn(u,K,k,b),thickness:X/2,color:v?j(I,et):I})}}else O&&H.push({points:u,thickness:V,color:W});let z=H.length-B;z>0&&Nt.set(f,{start:B,count:z})}return{nodeStrokes:H,nodeStrokeRuns:Nt,arrows:rt,arrowRuns:Lt,edgeCurves:G,edgeCurveRuns:Ot}}var un=[{value:"progress",label:"Progress along the walk"},{value:"walk",label:"One colour for the walk"},{value:"reference",label:"Reference position"}],wt=[{value:"grey",label:"Grey to black",ends:["#a3a3a3","#151515"]},{value:"blue",label:"Cyan to navy",ends:["#6ee0f0","#0b2c7a"]},{value:"red",label:"Yellow to red",ends:["#ffd24a","#b3121b"]},{value:"purple",label:"Pink to plum",ends:["#f5a3d8","#5e1a6e"]},{value:"green",label:"Lime to forest",ends:["#b6e36c","#0b5d3b"]},{value:"rainbow",label:"Rainbow"}],fn=["blue","red","purple","green"],ae=t=>un.some(e=>e.value===t),ce=t=>wt.some(e=>e.value===t);function ie(t,e,n,o=!0){let{field:r,scheme:s}=t.color??{},a=ae(r)?r:o?"progress":"walk",c=e?"grey":fn[n%fn.length],i=ce(s)?s==="rainbow"&&a!=="reference"?c:s:c;return{field:a,scheme:i}}var dn=.5,le=t=>t<=.04045?t/12.92:((t+.055)/1.055)**2.4,fe=t=>t<=.0031308?t*12.92:1.055*t**(1/2.4)-.055;function hn(t){let[e,n,o]=[1,3,5].map(c=>le(Number.parseInt(t.slice(c,c+2),16)/255)),r=Math.cbrt(.4122214708*e+.5363325363*n+.0514459929*o),s=Math.cbrt(.2119034982*e+.6806995451*n+.1073969566*o),a=Math.cbrt(.0883024619*e+.2817188376*n+.6299787005*o);return[.2104542553*r+.793617785*s-.0040720468*a,1.9779984951*r-2.428592205*s+.4505937099*a,.0259040371*r+.7827717662*s-.808675766*a]}function he([t,e,n]){let o=(t+.3963377774*e+.2158037573*n)**3,r=(t-.1055613458*e-.0638541728*n)**3,s=(t-.0894841775*e-1.291485548*n)**3;return[4.0767416621*o-3.3077115913*r+.2309699292*s,-1.2684380046*o+2.6097574011*r-.3413193965*s,-.0041960863*o-.7034186147*r+1.707614701*s].map(a=>fe(Math.max(0,Math.min(1,a)))*255)}function Tt(t,e){let n=Math.max(0,Math.min(1,e)),o=wt.find(a=>a.value===t);if(!o||!("ends"in o))return St(nt(n),Z,tt).map(a=>a*255);let r=hn(o.ends[0]),s=hn(o.ends[1]);return he(r.map((a,c)=>a+(s[c]-a)*n))}function mn(t,e){let[n,o,r]=Tt(t,e).map(Math.round);return`rgb(${n}, ${o}, ${r})`}function ht(t,e){let[n,o,r]=Tt(t,e).map(Math.round);return p(n,o,r,255)}function ue({field:t,scheme:e},n){return t==="walk"?ht(e,dn):n===void 0?pn:ht(e,n)}var pn=It;function de({field:t,scheme:e}){let n=(o,r)=>({offset:o,color:mn(e,r),abgr:ht(e,r)});return t==="walk"?[n(0,dn)]:Array.from({length:7},(o,r)=>n(r/6,r/6))}export{In as a,Tn as b,p as c,zt as d,me as e,j as f,pe as g,An as h,Ee as i,vn as j,xe as k,Jt as l,Pe as m,St as n,De as o,Fe as p,Qn as q,Zn as r,He as s,un as t,wt as u,ie as v,ue as w,pn as x,de as y};
|
|
2
|
-
//# sourceMappingURL=chunk-2M4XYSLG.js.map
|
|
@@ -1,7 +0,0 @@
|
|
|
1
|
-
{
|
|
2
|
-
"version": 3,
|
|
3
|
-
"sources": ["../../node_modules/.pnpm/@jbrowse+bandage-core@11.4.0/node_modules/@jbrowse/bandage-core/dist/chunks/chunk-G2FD7S3C.js", "../../node_modules/.pnpm/@jbrowse+bandage-core@11.4.0/node_modules/@jbrowse/bandage-core/dist/chunks/chunk-GQZVR5DW.js", "../../node_modules/.pnpm/@jbrowse+bandage-core@11.4.0/node_modules/@jbrowse/bandage-core/dist/chunks/chunk-WA5WNP6V.js", "../../node_modules/.pnpm/@jbrowse+bandage-core@11.4.0/node_modules/@jbrowse/bandage-core/dist/chunks/chunk-VX656OAZ.js", "../../node_modules/.pnpm/@jbrowse+bandage-core@11.4.0/node_modules/@jbrowse/bandage-core/dist/chunks/chunk-XSF4RQRP.js", "../../node_modules/.pnpm/@jbrowse+bandage-core@11.4.0/node_modules/@jbrowse/bandage-core/dist/chunks/chunk-GM7M25NW.js"],
|
|
4
|
-
"sourcesContent": ["// src/renderer/colorBits.ts\nfunction packAbgr(r, g, b, a) {\n return (a << 24 | b << 16 | g << 8 | r) >>> 0;\n}\nfunction abgrRed(c) {\n return c & 255;\n}\nfunction abgrGreen(c) {\n return c >>> 8 & 255;\n}\nfunction abgrBlue(c) {\n return c >>> 16 & 255;\n}\nfunction abgrAlpha(c) {\n return c >>> 24 & 255;\n}\nfunction abgrToCssRgba(c) {\n return `rgba(${abgrRed(c)},${abgrGreen(c)},${abgrBlue(c)},${abgrAlpha(c) / 255})`;\n}\nfunction normalizedRgbToCssRgba(c, alpha) {\n return `rgba(${Math.round(c[0] * 255)},${Math.round(c[1] * 255)},${Math.round(c[2] * 255)},${alpha})`;\n}\nfunction fadeAbgr(c, alpha) {\n return packAbgr(\n abgrRed(c),\n abgrGreen(c),\n abgrBlue(c),\n Math.round(abgrAlpha(c) * alpha)\n );\n}\nfunction brightenAbgr(c, factor) {\n const t = Math.max(0, 1 - 1 / factor);\n const lift = (v) => Math.round(v + (255 - v) * t);\n return packAbgr(\n lift(abgrRed(c)),\n lift(abgrGreen(c)),\n lift(abgrBlue(c)),\n abgrAlpha(c)\n );\n}\n\nexport {\n packAbgr,\n abgrRed,\n abgrGreen,\n abgrBlue,\n abgrAlpha,\n abgrToCssRgba,\n normalizedRgbToCssRgba,\n fadeAbgr,\n brightenAbgr\n};\n", "// src/referenceRampCss.ts\nvar REFERENCE_RAMP_MAX_HUE = 300;\nvar REFERENCE_RAMP_SATURATION = 0.7;\nvar REFERENCE_RAMP_LIGHTNESS = 0.5;\nfunction rampHueCss(hue) {\n return `hsl(${hue}, ${REFERENCE_RAMP_SATURATION * 100}%, ${REFERENCE_RAMP_LIGHTNESS * 100}%)`;\n}\nfunction rampHueAt(frac) {\n return Math.max(0, Math.min(1, frac)) * REFERENCE_RAMP_MAX_HUE;\n}\nfunction rampHue(bp, domain) {\n return rampHueAt((bp - domain.start) / (domain.end - domain.start));\n}\nvar HUE_STOP_DEG = 30;\nfunction rampStops({ start, bp, reversed }, domain) {\n const low = rampHue(start, domain);\n const high = rampHue(start + bp, domain);\n if (high - low < 1) {\n return [(low + high) / 2];\n }\n const n = Math.ceil((high - low) / HUE_STOP_DEG);\n return Array.from(\n { length: n + 1 },\n (_, i) => rampHue(start + bp * (reversed ? n - i : i) / n, domain)\n );\n}\nvar RAMP_GRADIENT_STOPS = Array.from({ length: 7 }, (_, i) => ({\n offset: i / 6,\n color: rampHueCss(i / 6 * REFERENCE_RAMP_MAX_HUE)\n}));\nvar RAMP_GRADIENT_CSS = `linear-gradient(to right, ${RAMP_GRADIENT_STOPS.map(\n (s) => `${s.color} ${s.offset * 100}%`\n).join(\", \")})`;\n\nexport {\n REFERENCE_RAMP_MAX_HUE,\n REFERENCE_RAMP_SATURATION,\n REFERENCE_RAMP_LIGHTNESS,\n rampHueCss,\n rampHueAt,\n rampHue,\n rampStops,\n RAMP_GRADIENT_STOPS,\n RAMP_GRADIENT_CSS\n};\n", "import {\n weakMemo\n} from \"./chunk-HMXCN6PG.js\";\nimport {\n edgeCurvesOf\n} from \"./chunk-CGMSKIEO.js\";\n\n// src/util/edgeCurves.ts\nvar cache = weakMemo();\nfunction baseEdgeCurves(nodePositions, graph, axis, deletions, version, routes, stranded = false) {\n return cache(\n nodePositions,\n [graph, axis.scaleX, axis.scaleY, deletions, routes, stranded, version],\n () => {\n const curves = /* @__PURE__ */ new Map();\n for (let ei = 0; ei < graph.edges.length; ei++) {\n const bypassed = deletions?.get(ei);\n const edgeCurves = edgeCurvesOf(\n nodePositions,\n graph.edges[ei],\n axis,\n bypassed,\n bypassed && routes?.[ei],\n stranded\n );\n if (edgeCurves) {\n curves.set(ei, edgeCurves);\n }\n }\n return curves;\n }\n );\n}\n\nexport {\n baseEdgeCurves\n};\n", "import {\n pushTo\n} from \"./chunk-BSOZW3FC.js\";\nimport {\n isBackbone\n} from \"./chunk-RJWPJCIB.js\";\n\n// src/referenceSpan.ts\nvar MAX_BACKBONE_HOPS = 4;\nfunction buildNeighbors(graph) {\n const neighbors = /* @__PURE__ */ new Map();\n for (const edge of graph.edges) {\n pushTo(neighbors, edge.from, edge.to);\n pushTo(neighbors, edge.to, edge.from);\n }\n return neighbors;\n}\nfunction flankingSpan({\n nodeId,\n nodeById,\n neighbors\n}) {\n const seen = /* @__PURE__ */ new Set([nodeId]);\n let frontier = [nodeId];\n let minEnd = Infinity;\n let maxStart = -Infinity;\n for (let hop = 0; hop < MAX_BACKBONE_HOPS && frontier.length > 0; hop++) {\n const next = [];\n for (const id of frontier) {\n const adjacent = neighbors.get(id);\n if (adjacent) {\n for (const neighborId of adjacent) {\n if (!seen.has(neighborId)) {\n seen.add(neighborId);\n const neighbor = nodeById.get(neighborId);\n if (neighbor && isBackbone(neighbor)) {\n minEnd = Math.min(minEnd, neighbor.stable.start + neighbor.length);\n maxStart = Math.max(maxStart, neighbor.stable.start);\n } else {\n next.push(neighborId);\n }\n }\n }\n }\n }\n frontier = next;\n }\n return minEnd === Infinity ? void 0 : {\n start: Math.min(minEnd, maxStart),\n end: Math.max(minEnd, maxStart)\n };\n}\nfunction nodeReferenceSpan({\n nodeId,\n nodeById,\n neighbors\n}) {\n const node = nodeById.get(nodeId);\n return node === void 0 ? void 0 : isBackbone(node) ? { start: node.stable.start, end: node.stable.start + node.length } : flankingSpan({ nodeId, nodeById, neighbors });\n}\nfunction referenceMidpoints(graph) {\n const nodeById = new Map(graph.nodes.map((node) => [node.id, node]));\n const neighbors = buildNeighbors(graph);\n const midpoints = /* @__PURE__ */ new Map();\n for (const node of graph.nodes) {\n const span = nodeReferenceSpan({ nodeId: node.id, nodeById, neighbors });\n if (span) {\n midpoints.set(node.id, (span.start + span.end) / 2);\n }\n }\n return midpoints;\n}\n\nexport {\n buildNeighbors,\n nodeReferenceSpan,\n referenceMidpoints\n};\n", "import {\n abgrToCssRgba,\n fadeAbgr,\n packAbgr\n} from \"./chunk-G2FD7S3C.js\";\nimport {\n REFERENCE_RAMP_LIGHTNESS,\n REFERENCE_RAMP_SATURATION,\n rampHueAt\n} from \"./chunk-GQZVR5DW.js\";\nimport {\n baseEdgeCurves\n} from \"./chunk-WA5WNP6V.js\";\nimport {\n PATH_LIGHTNESS,\n PATH_SATURATION,\n nameHue,\n pathHueAt\n} from \"./chunk-3GP3KT5L.js\";\nimport {\n meanDepth,\n nodeWidthPx\n} from \"./chunk-CRCD3ODN.js\";\nimport {\n referenceMidpoints\n} from \"./chunk-VX656OAZ.js\";\nimport {\n dashCurves,\n pathRibbonOffsets,\n polylinePx,\n strandSides,\n translateCurves,\n yToXOf\n} from \"./chunk-UBQSRCWR.js\";\nimport {\n pushTo\n} from \"./chunk-BSOZW3FC.js\";\nimport {\n isBackbone\n} from \"./chunk-RJWPJCIB.js\";\n\n// src/renderer/GeometryBuilder.ts\nvar EDGE_DEFAULT_COLOR = packAbgr(119, 119, 119, 217);\nvar EDGE_PATH_FALLBACK_COLOR = packAbgr(136, 136, 136, 217);\nvar EDGE_DELETION_COLOR = packAbgr(24, 24, 28, 240);\nvar DELETION_THICKNESS_FACTOR = 2.2;\nvar LIFT_BACKDROP_COLOR = packAbgr(160, 160, 160, 255);\nvar FADED_ALPHA = 0.18;\nvar LIFT_BACKDROP_CSS = abgrToCssRgba(\n fadeAbgr(LIFT_BACKDROP_COLOR, FADED_ALPHA)\n);\nvar DELETION_DASH_PX = 11;\nvar MIN_PATH_STRIPE_PX = 1.2;\nvar MIN_WALK_LANE_PX = 4;\nvar MIN_SCALE_FOR_OFFSET = 1e-6;\nvar ARROW_LENGTH_BASE_PX = 7;\nvar ARROW_LENGTH_PER_EDGE_PX = 3;\nvar ARROW_HALF_WIDTH_BASE_PX = 3;\nvar ARROW_HALF_WIDTH_PER_EDGE_PX = 1.5;\nvar MIN_LINK_HEADS = 3;\nvar MIN_ARROW_FIT = 0.5;\nvar MIN_ARROW_SCALE = 0.45;\nvar ARROW_CURVE_STEPS = 12;\nvar ABUTTING_PX = 1e-6;\nvar walkX = new Float64Array(64);\nvar walkY = new Float64Array(64);\nvar walkPx = new Float64Array(64);\nfunction backFromEnd(curves, scale, yToX, untilPx) {\n const most = curves.length * ARROW_CURVE_STEPS + 1;\n if (walkX.length < most) {\n walkX = new Float64Array(most);\n walkY = new Float64Array(most);\n walkPx = new Float64Array(most);\n }\n const last = curves[curves.length - 1];\n walkX[0] = last.x1;\n walkY[0] = last.y1;\n walkPx[0] = 0;\n let n = 1;\n let px = 0;\n for (let ci = curves.length - 1; ci >= 0 && px < untilPx; ci--) {\n const c = curves[ci];\n for (let k = 1; k <= ARROW_CURVE_STEPS && px < untilPx; k++) {\n const t = 1 - (k / ARROW_CURVE_STEPS) ** 2;\n const u = 1 - t;\n const a = u * u * u;\n const b = 3 * u * u * t;\n const d = 3 * u * t * t;\n const e = t * t * t;\n const x = a * c.x0 + b * c.cx0 + d * c.cx1 + e * c.x1;\n const y = a * c.y0 + b * c.cy0 + d * c.cy1 + e * c.y1;\n px += Math.hypot(x - walkX[n - 1], (y - walkY[n - 1]) * yToX) * scale;\n walkX[n] = x;\n walkY[n] = y;\n walkPx[n] = px;\n n++;\n }\n }\n return n;\n}\nfunction walkPointAt(n, px) {\n let i = 1;\n while (i < n - 1 && walkPx[i] < px) {\n i++;\n }\n const from = walkPx[i - 1];\n const to = walkPx[i];\n const f = to > from ? Math.min(1, Math.max(0, (px - from) / (to - from))) : 0;\n return {\n x: walkX[i - 1] + (walkX[i] - walkX[i - 1]) * f,\n y: walkY[i - 1] + (walkY[i] - walkY[i - 1]) * f\n };\n}\nfunction arrowheadFor({\n curves,\n edgeHalfWidth,\n nodeHalfWidth,\n scale,\n yToX,\n color,\n fromSegments,\n toSegments\n}) {\n const length = ARROW_LENGTH_BASE_PX + ARROW_LENGTH_PER_EDGE_PX * edgeHalfWidth;\n const n = backFromEnd(\n curves,\n scale,\n yToX,\n nodeHalfWidth + MIN_LINK_HEADS * length\n );\n const lengthPx = walkPx[n - 1];\n const abutting = lengthPx < ABUTTING_PX;\n const inset = Math.min(nodeHalfWidth, lengthPx);\n const fit = abutting ? Math.min(1, polylinePx(fromSegments, scale, scale * yToX) / length) : 1;\n if (abutting ? fit < MIN_ARROW_FIT : lengthPx < MIN_LINK_HEADS * length) {\n return void 0;\n }\n const tip = walkPointAt(n, inset);\n const base = walkPointAt(n, inset + length * fit);\n const dx = tip.x - base.x;\n const dy = (tip.y - base.y) * yToX;\n const last = curves[curves.length - 1];\n return {\n x: tip.x,\n y: tip.y,\n angle: Math.hypot(dx, dy) > 0 ? Math.atan2(dy, dx) : endTangent(last, yToX, intoDirection(last, toSegments)),\n length: length * fit,\n halfWidth: (ARROW_HALF_WIDTH_BASE_PX + ARROW_HALF_WIDTH_PER_EDGE_PX * edgeHalfWidth) * fit,\n color\n };\n}\nfunction leavesStart(edge, segments, first, stranded, yToX) {\n if (stranded) {\n return strandSides(edge).from === \"start\";\n }\n const gap = (p) => Math.hypot(p.x - first.x0, (p.y - first.y0) * yToX);\n return gap(segments[0]) < gap(segments.at(-1));\n}\nfunction intoDirection(last, toSegments) {\n if (toSegments.length < 2) {\n return void 0;\n }\n const first = toSegments[0];\n const final = toSegments[toSegments.length - 1];\n const atStart = Math.hypot(last.x1 - first.x, last.y1 - first.y) <= Math.hypot(last.x1 - final.x, last.y1 - final.y);\n return atStart ? { from: first, to: toSegments[1] } : { from: final, to: toSegments[toSegments.length - 2] };\n}\nfunction normal(a, b) {\n if (!a || !b) {\n return void 0;\n }\n const dx = b.x - a.x;\n const dy = b.y - a.y;\n const len = Math.hypot(dx, dy);\n return len > 0 ? { nx: -dy / len, ny: dx / len } : void 0;\n}\nfunction pointNormalsOf(rawPoints, yToX = 1) {\n const points = yToX === 1 ? rawPoints : rawPoints.map((p) => ({ x: p.x, y: p.y * yToX }));\n return points.map((p, i) => {\n const before = normal(points[i - 1], p);\n const after = normal(p, points[i + 1]);\n if (before && after) {\n const nx = (before.nx + after.nx) / 2;\n const ny = (before.ny + after.ny) / 2;\n const dot = before.nx * nx + before.ny * ny;\n return dot > 0.1 ? { nx: nx / dot, ny: ny / dot } : { nx, ny };\n }\n return before ?? after ?? { nx: 0, ny: 0 };\n });\n}\nfunction offsetPolyline(points, normals, distance, yToX = 1) {\n return points.map((p, i) => ({\n x: p.x + normals[i].nx * distance,\n y: p.y + normals[i].ny * distance / yToX\n }));\n}\nfunction hslToRgb(h, s, l) {\n const c = (1 - Math.abs(2 * l - 1)) * s;\n const x = c * (1 - Math.abs(h / 60 % 2 - 1));\n const m = l - c / 2;\n let r = 0;\n let g = 0;\n let b = 0;\n if (h < 60) {\n r = c;\n g = x;\n } else if (h < 120) {\n r = x;\n g = c;\n } else if (h < 180) {\n g = c;\n b = x;\n } else if (h < 240) {\n g = x;\n b = c;\n } else if (h < 300) {\n r = x;\n b = c;\n } else {\n r = c;\n b = x;\n }\n return [r + m, g + m, b + m];\n}\nfunction packHsl(h, s, l, a = 1) {\n const [r, g, b] = hslToRgb(h, s, l);\n return packAbgr(\n Math.round(r * 255),\n Math.round(g * 255),\n Math.round(b * 255),\n Math.round(a * 255)\n );\n}\nvar DEPTH_GRADIENT = [\n [68, 1, 84],\n [59, 82, 139],\n [33, 145, 140],\n [94, 201, 98],\n [253, 231, 37]\n];\nvar STABLE_RANK_GRADIENT = [\n [237, 137, 44],\n [158, 42, 122]\n];\nvar NODE_LENGTH_GRADIENT = [\n [220, 50, 50],\n [50, 120, 220]\n];\nfunction sampleGradient(stops, t) {\n const pos = Math.max(0, Math.min(1, t)) * (stops.length - 1);\n const i = Math.min(stops.length - 2, Math.floor(pos));\n const f = pos - i;\n const a = stops[i];\n const b = stops[i + 1];\n return packAbgr(\n Math.round(a[0] + (b[0] - a[0]) * f),\n Math.round(a[1] + (b[1] - a[1]) * f),\n Math.round(a[2] + (b[2] - a[2]) * f),\n 255\n );\n}\nfunction hashColor(str, alpha) {\n return packHsl(nameHue(str), PATH_SATURATION, PATH_LIGHTNESS, alpha);\n}\nvar REFERENCE_RAMP_ALT_RGB = [60, 65, 72];\nvar REFERENCE_RAMP_ALT_COLOR = packAbgr(...REFERENCE_RAMP_ALT_RGB, 255);\nvar REFERENCE_RAMP_ALT_CSS = `rgb(${REFERENCE_RAMP_ALT_RGB.join(\", \")})`;\nfunction computeReferenceRamp(graph, domain) {\n const midpoints = referenceMidpoints(graph);\n if (domain && domain.end > domain.start) {\n return {\n start: domain.start,\n span: domain.end - domain.start,\n midpoints\n };\n }\n let min = Infinity;\n let max = -Infinity;\n for (const mid of midpoints.values()) {\n min = Math.min(min, mid);\n max = Math.max(max, mid);\n }\n return {\n start: min === Infinity ? 0 : min,\n span: max > min ? max - min : 1,\n midpoints\n };\n}\nfunction computeColorSchemeRange(graph, referenceRamp) {\n let minDepth = Infinity;\n let maxDepth = -Infinity;\n let minLength = Infinity;\n let maxLength = -Infinity;\n let maxRank = 0;\n for (const n of graph.nodes) {\n if (n.stable !== void 0 && n.stable.rank > maxRank) {\n maxRank = n.stable.rank;\n }\n if (n.depth < minDepth) {\n minDepth = n.depth;\n }\n if (n.depth > maxDepth) {\n maxDepth = n.depth;\n }\n if (n.length < minLength) {\n minLength = n.length;\n }\n if (n.length > maxLength) {\n maxLength = n.length;\n }\n }\n return {\n minDepth,\n maxDepth,\n minLength,\n maxLength,\n maxRank,\n nodeCount: graph.nodes.length,\n referenceRamp\n };\n}\nfunction getNodeColor(node, nodeIndex, colorScheme, range) {\n switch (colorScheme) {\n case \"random\":\n return hashColor(node.id, 1);\n case \"depth\":\n return sampleGradient(\n DEPTH_GRADIENT,\n range.maxDepth > range.minDepth ? (node.depth - range.minDepth) / (range.maxDepth - range.minDepth) : 0.5\n );\n case \"node-length\":\n return sampleGradient(\n NODE_LENGTH_GRADIENT,\n range.maxLength > range.minLength ? (node.length - range.minLength) / (range.maxLength - range.minLength) : 0.5\n );\n // rGFA states the backbone rather than leaving it to be inferred: rank 0 is\n // the reference, higher ranks are the sequence that diverges from it. A\n // graph with no SR tags draws every node as unranked.\n case \"stable-rank\":\n return node.stable === void 0 ? packAbgr(160, 160, 160, 255) : node.stable.rank === 0 ? packAbgr(52, 152, 219, 255) : sampleGradient(\n STABLE_RANK_GRADIENT,\n range.maxRank > 1 ? (node.stable.rank - 1) / (range.maxRank - 1) : 0\n );\n // Position on the reference, the one quantity both panels of a\n // graph-over-linear figure can state, so only the backbone takes a hue\n case \"reference-position\": {\n const ramp = range.referenceRamp;\n if (ramp === void 0) {\n return packAbgr(160, 160, 160, 255);\n }\n const mid = ramp.midpoints.get(node.id);\n if (!isBackbone(node) || mid === void 0) {\n return REFERENCE_RAMP_ALT_COLOR;\n }\n return packHsl(\n rampHueAt((mid - ramp.start) / ramp.span),\n REFERENCE_RAMP_SATURATION,\n REFERENCE_RAMP_LIGHTNESS\n );\n }\n case \"grey\":\n return packAbgr(160, 160, 160, 255);\n case \"rainbow\": {\n const hue = range.nodeCount > 1 ? nodeIndex / range.nodeCount * 360 : 0;\n return packHsl(hue, 0.75, 0.5);\n }\n default:\n return packAbgr(52, 152, 219, 255);\n }\n}\nfunction endTangent(c, yToX = 1, into) {\n const candidates = [\n [c.cx1, c.cy1],\n [c.cx0, c.cy0],\n [c.x0, c.y0]\n ];\n for (const [x, y] of candidates) {\n const dx = c.x1 - x;\n const dy = (c.y1 - y) * yToX;\n if (Math.hypot(dx, dy) > 0) {\n return Math.atan2(dy, dx);\n }\n }\n return into ? Math.atan2((into.to.y - into.from.y) * yToX, into.to.x - into.from.x) : 0;\n}\nfunction isPolylineInBounds(segments, bounds) {\n let hit = false;\n for (let i = 0; i < segments.length - 1 && !hit; i++) {\n const a = segments[i];\n const b = segments[i + 1];\n hit = Math.max(a.x, b.x) >= bounds.minX && Math.min(a.x, b.x) <= bounds.maxX && Math.max(a.y, b.y) >= bounds.minY && Math.min(a.y, b.y) <= bounds.maxY;\n }\n return hit;\n}\nfunction isBezierInBounds(curves, bounds) {\n for (const c of curves) {\n const cMinX = Math.min(c.x0, c.cx0, c.cx1, c.x1);\n const cMaxX = Math.max(c.x0, c.cx0, c.cx1, c.x1);\n const cMinY = Math.min(c.y0, c.cy0, c.cy1, c.y1);\n const cMaxY = Math.max(c.y0, c.cy0, c.cy1, c.y1);\n if (cMaxX >= bounds.minX && cMinX <= bounds.maxX && cMaxY >= bounds.minY && cMinY <= bounds.maxY) {\n return true;\n }\n }\n return false;\n}\nvar paletteCache = /* @__PURE__ */ new WeakMap();\nfunction pathPaletteOf(graph) {\n let palette = paletteCache.get(graph);\n if (!palette) {\n const pathColorByIndex = [];\n const pathIndexByName = /* @__PURE__ */ new Map();\n const paths = graph.paths ?? [];\n for (let i = 0; i < paths.length; i++) {\n pathColorByIndex.push(\n packHsl(\n pathHueAt(i, paths.length),\n PATH_SATURATION,\n PATH_LIGHTNESS,\n 0.85\n )\n );\n if (!pathIndexByName.has(paths[i].name)) {\n pathIndexByName.set(paths[i].name, i);\n }\n }\n palette = { pathColorByIndex, pathIndexByName };\n paletteCache.set(graph, palette);\n }\n return palette;\n}\nvar slotsCache = /* @__PURE__ */ new WeakMap();\nvar NO_SLOTS = /* @__PURE__ */ new Map();\nfunction pathSlotsOf(graph) {\n let slots = slotsCache.get(graph);\n if (!slots) {\n slots = /* @__PURE__ */ new Map();\n for (const [pathIdx, path] of (graph.paths ?? []).entries()) {\n for (const nodeId of new Set(path.nodeIds)) {\n pushTo(slots, nodeId, pathIdx);\n }\n }\n slotsCache.set(graph, slots);\n }\n return slots;\n}\nfunction buildGeometry(options) {\n const {\n nodePositions,\n graph,\n nodeById,\n colorScheme,\n contigThickness,\n connectorThickness,\n drawPaths,\n nodeWidth = \"uniform\",\n highlight,\n axis,\n linearLayout,\n viewportBounds,\n referenceRamp,\n deletions,\n deletionRoutes,\n stranded,\n hiddenEdges,\n version = 0\n } = options;\n const depthNorm = nodeWidth === \"depth\" ? meanDepth(graph) : 0;\n const scale = axis.scaleX;\n const worldPerScreenPx = 1 / Math.max(scale, MIN_SCALE_FOR_OFFSET);\n const yToX = yToXOf(axis);\n const sharedCurves = baseEdgeCurves(\n nodePositions,\n graph,\n axis,\n deletions,\n version,\n deletionRoutes,\n stranded\n );\n const nodeStrokes = [];\n const nodeStrokeRuns = /* @__PURE__ */ new Map();\n const arrows = [];\n const arrowRuns = /* @__PURE__ */ new Map();\n const edgeCurves = [];\n const edgeCurveRuns = /* @__PURE__ */ new Map();\n const colorRange = computeColorSchemeRange(graph, referenceRamp);\n const nodeIndexMap = /* @__PURE__ */ new Map();\n if (colorScheme === \"rainbow\") {\n for (let i = 0; i < graph.nodes.length; i++) {\n nodeIndexMap.set(graph.nodes[i].id, i);\n }\n }\n const pathCount = graph.paths?.length ?? 0;\n const { pathColorByIndex, pathIndexByName } = pathPaletteOf(graph);\n const pathColorByName = (name) => pathColorByIndex[pathIndexByName.get(name) ?? -1] ?? EDGE_PATH_FALLBACK_COLOR;\n const nodePathSlots = drawPaths && pathCount > 1 ? pathSlotsOf(graph) : NO_SLOTS;\n const showArrows = scale > (linearLayout ? 1 : MIN_ARROW_SCALE);\n for (let ei = 0; ei < graph.edges.length; ei++) {\n const edge = graph.edges[ei];\n const fromSegments = nodePositions[edge.from];\n const toSegments = nodePositions[edge.to];\n const baseCurves = sharedCurves.get(ei);\n if (!fromSegments?.length || !toSegments?.length || !baseCurves || hiddenEdges?.has(ei)) {\n continue;\n }\n const numPaths = edge.pathIds?.length ?? 0;\n const bypassed = deletions?.get(ei);\n const isDeletion = bypassed !== void 0;\n const onWalk = highlight?.edgeIndexes.has(ei) ?? false;\n const edgeThickness = connectorThickness / 2 * (isDeletion ? DELETION_THICKNESS_FACTOR : 1);\n const ribbons = drawPaths && numPaths > 0;\n const plainEdgeColor = isDeletion ? EDGE_DELETION_COLOR : EDGE_DEFAULT_COLOR;\n const edgeColor = !highlight || onWalk ? plainEdgeColor : fadeAbgr(plainEdgeColor, FADED_ALPHA);\n const toNode = nodeById.get(edge.to);\n const intoHalfWidth = (toNode ? nodeWidthPx(toNode, contigThickness, nodeWidth, depthNorm) : contigThickness) / 2;\n const inView = (curves) => !viewportBounds || isBezierInBounds(curves, viewportBounds);\n const stroke = (curves, thickness, color) => {\n if (isDeletion) {\n for (const dash of dashCurves(curves, DELETION_DASH_PX / scale)) {\n edgeCurves.push({ curves: dash, thickness, color });\n }\n } else {\n edgeCurves.push({ curves, thickness, color });\n }\n };\n const arrowOn = (curves) => {\n const head = showArrows ? arrowheadFor({\n curves,\n edgeHalfWidth: connectorThickness / 2,\n nodeHalfWidth: intoHalfWidth,\n scale,\n yToX,\n color: edgeColor,\n fromSegments,\n toSegments\n }) : void 0;\n if (head) {\n arrows.push(head);\n }\n };\n const edgeCurveStart = edgeCurves.length;\n const arrowStart = arrows.length;\n const laneWalks = highlight?.walks.length && onWalk && !ribbons ? highlight.walks : void 0;\n if (laneWalks) {\n const fromNode = nodeById.get(edge.from);\n const width = fromNode ? nodeWidthPx(fromNode, contigThickness, nodeWidth, depthNorm) : contigThickness;\n const laneWidth = Math.max(width / laneWalks.length, MIN_WALK_LANE_PX);\n const normals = pointNormalsOf(fromSegments, yToX);\n const { nx, ny } = (leavesStart(\n edge,\n fromSegments,\n baseCurves[0],\n stranded,\n yToX\n ) ? normals[0] : normals.at(-1)) ?? {\n nx: 0,\n ny: 0\n };\n laneWalks.forEach((walk, lane) => {\n if (!walk.edgeIndexes.has(ei)) {\n return;\n }\n const offset = (lane - (laneWalks.length - 1) / 2) * laneWidth * worldPerScreenPx;\n const curves = offset === 0 ? baseCurves : translateCurves(baseCurves, nx * offset, ny * offset / yToX);\n if (inView(curves)) {\n stroke(\n curves,\n laneWidth / 2,\n walk.colors.get(edge.from) ?? walk.colors.get(edge.to) ?? plainEdgeColor\n );\n }\n });\n if (edgeCurves.length > edgeCurveStart) {\n arrowOn(baseCurves);\n }\n } else if (!ribbons) {\n if (inView(baseCurves)) {\n stroke(baseCurves, edgeThickness, edgeColor);\n arrowOn(baseCurves);\n }\n } else {\n const offsets = pathRibbonOffsets(\n fromSegments,\n toSegments,\n numPaths,\n axis\n );\n const arrowRibbon = Math.floor((numPaths - 1) / 2);\n for (let pathIdx = 0; pathIdx < numPaths; pathIdx++) {\n const { x, y } = offsets[pathIdx];\n const curves = x === 0 && y === 0 ? baseCurves : translateCurves(baseCurves, x, y);\n if (!inView(curves)) {\n continue;\n }\n const ribbonColor = pathColorByName(edge.pathIds[pathIdx]);\n stroke(\n curves,\n edgeThickness,\n highlight && (highlight.walks.length ? !highlight.names.has(edge.pathIds[pathIdx]) : !onWalk) ? fadeAbgr(ribbonColor, FADED_ALPHA) : ribbonColor\n );\n if (pathIdx === arrowRibbon) {\n arrowOn(curves);\n }\n }\n }\n const edgeCurveCount = edgeCurves.length - edgeCurveStart;\n if (edgeCurveCount > 0) {\n edgeCurveRuns.set(ei, { start: edgeCurveStart, count: edgeCurveCount });\n }\n const arrowCount = arrows.length - arrowStart;\n if (arrowCount > 0) {\n arrowRuns.set(ei, { start: arrowStart, count: arrowCount });\n }\n }\n for (const [nodeId, segments] of Object.entries(nodePositions)) {\n const node = nodeById.get(nodeId);\n if (!node) {\n continue;\n }\n if (viewportBounds && !isPolylineInBounds(segments, viewportBounds)) {\n continue;\n }\n const faded = highlight !== void 0 && !highlight.nodeIds.has(nodeId);\n const own = faded ? LIFT_BACKDROP_COLOR : getNodeColor(\n node,\n nodeIndexMap.get(nodeId) ?? 0,\n colorScheme,\n colorRange\n );\n const color = faded ? fadeAbgr(own, FADED_ALPHA) : own;\n const width = nodeWidthPx(node, contigThickness, nodeWidth, depthNorm);\n const nodeThickness = width / 2;\n const start = nodeStrokes.length;\n const slots = nodePathSlots.get(nodeId);\n const slotWidth = width / pathCount;\n const drawable = segments.length >= 2;\n const lanes = highlight?.walks.length && !faded && !slots?.length ? highlight.walks : void 0;\n if (drawable && lanes) {\n const laneWidth = Math.max(width / lanes.length, MIN_WALK_LANE_PX);\n const normals = pointNormalsOf(segments, yToX);\n lanes.forEach((walk, lane) => {\n const laneColor = walk.colors.get(nodeId);\n if (laneColor === void 0) {\n return;\n }\n const offset = (lane - (lanes.length - 1) / 2) * laneWidth * worldPerScreenPx;\n nodeStrokes.push({\n points: lanes.length === 1 ? segments : offsetPolyline(segments, normals, offset, yToX),\n thickness: laneWidth / 2,\n color: laneColor\n });\n });\n } else if (drawable && slots?.length && slotWidth >= MIN_PATH_STRIPE_PX) {\n const normals = pointNormalsOf(segments, yToX);\n for (const slot of slots) {\n const offset = (slot - (pathCount - 1) / 2) * slotWidth * worldPerScreenPx;\n const stripe = pathColorByIndex[slot] ?? EDGE_PATH_FALLBACK_COLOR;\n nodeStrokes.push({\n points: offsetPolyline(segments, normals, offset, yToX),\n thickness: slotWidth / 2,\n // a node off the lifted walk fades whichever way it is painted\n color: faded ? fadeAbgr(stripe, FADED_ALPHA) : stripe\n });\n }\n } else if (drawable) {\n nodeStrokes.push({ points: segments, thickness: nodeThickness, color });\n }\n const count = nodeStrokes.length - start;\n if (count > 0) {\n nodeStrokeRuns.set(nodeId, { start, count });\n }\n }\n return {\n nodeStrokes,\n nodeStrokeRuns,\n arrows,\n arrowRuns,\n edgeCurves,\n edgeCurveRuns\n };\n}\n\nexport {\n LIFT_BACKDROP_COLOR,\n FADED_ALPHA,\n LIFT_BACKDROP_CSS,\n hslToRgb,\n REFERENCE_RAMP_ALT_COLOR,\n REFERENCE_RAMP_ALT_CSS,\n computeReferenceRamp,\n computeColorSchemeRange,\n getNodeColor,\n endTangent,\n buildGeometry\n};\n", "import {\n REFERENCE_RAMP_ALT_COLOR,\n hslToRgb\n} from \"./chunk-XSF4RQRP.js\";\nimport {\n packAbgr\n} from \"./chunk-G2FD7S3C.js\";\nimport {\n REFERENCE_RAMP_LIGHTNESS,\n REFERENCE_RAMP_SATURATION,\n rampHueAt\n} from \"./chunk-GQZVR5DW.js\";\n\n// src/walkEncoding.ts\nvar WALK_FIELDS = [\n { value: \"progress\", label: \"Progress along the walk\" },\n { value: \"walk\", label: \"One colour for the walk\" },\n { value: \"reference\", label: \"Reference position\" }\n];\nvar WALK_SCHEMES = [\n { value: \"grey\", label: \"Grey to black\", ends: [\"#a3a3a3\", \"#151515\"] },\n { value: \"blue\", label: \"Cyan to navy\", ends: [\"#6ee0f0\", \"#0b2c7a\"] },\n { value: \"red\", label: \"Yellow to red\", ends: [\"#ffd24a\", \"#b3121b\"] },\n { value: \"purple\", label: \"Pink to plum\", ends: [\"#f5a3d8\", \"#5e1a6e\"] },\n { value: \"green\", label: \"Lime to forest\", ends: [\"#b6e36c\", \"#0b5d3b\"] },\n { value: \"rainbow\", label: \"Rainbow\" }\n];\nvar FAMILIES = [\"blue\", \"red\", \"purple\", \"green\"];\nvar isField = (v) => WALK_FIELDS.some((f) => f.value === v);\nvar isScheme = (v) => WALK_SCHEMES.some((s) => s.value === v);\nfunction resolveEncoding(layer, reference, picked, alone = true) {\n const { field: askedField, scheme: askedScheme } = layer.color ?? {};\n const field = isField(askedField) ? askedField : alone ? \"progress\" : \"walk\";\n const family = reference ? \"grey\" : FAMILIES[picked % FAMILIES.length];\n const scheme = isScheme(askedScheme) ? askedScheme === \"rainbow\" && field !== \"reference\" ? family : askedScheme : family;\n return { field, scheme };\n}\nvar BASE_T = 0.5;\nvar toLinear = (c) => c <= 0.04045 ? c / 12.92 : ((c + 0.055) / 1.055) ** 2.4;\nvar toGamma = (c) => c <= 31308e-7 ? c * 12.92 : 1.055 * c ** (1 / 2.4) - 0.055;\nfunction oklab(hex) {\n const [r, g, b] = [1, 3, 5].map(\n (i) => toLinear(Number.parseInt(hex.slice(i, i + 2), 16) / 255)\n );\n const l = Math.cbrt(0.4122214708 * r + 0.5363325363 * g + 0.0514459929 * b);\n const m = Math.cbrt(0.2119034982 * r + 0.6806995451 * g + 0.1073969566 * b);\n const s = Math.cbrt(0.0883024619 * r + 0.2817188376 * g + 0.6299787005 * b);\n return [\n 0.2104542553 * l + 0.793617785 * m - 0.0040720468 * s,\n 1.9779984951 * l - 2.428592205 * m + 0.4505937099 * s,\n 0.0259040371 * l + 0.7827717662 * m - 0.808675766 * s\n ];\n}\nfunction oklabToRgb([L, a, b]) {\n const l = (L + 0.3963377774 * a + 0.2158037573 * b) ** 3;\n const m = (L - 0.1055613458 * a - 0.0638541728 * b) ** 3;\n const s = (L - 0.0894841775 * a - 1.291485548 * b) ** 3;\n return [\n 4.0767416621 * l - 3.3077115913 * m + 0.2309699292 * s,\n -1.2684380046 * l + 2.6097574011 * m - 0.3413193965 * s,\n -0.0041960863 * l - 0.7034186147 * m + 1.707614701 * s\n ].map((c) => toGamma(Math.max(0, Math.min(1, c))) * 255);\n}\nfunction schemeRgb(scheme, t) {\n const x = Math.max(0, Math.min(1, t));\n const family = WALK_SCHEMES.find((s) => s.value === scheme);\n if (!family || !(\"ends\" in family)) {\n return hslToRgb(\n rampHueAt(x),\n REFERENCE_RAMP_SATURATION,\n REFERENCE_RAMP_LIGHTNESS\n ).map((c) => c * 255);\n }\n const from = oklab(family.ends[0]);\n const to = oklab(family.ends[1]);\n return oklabToRgb(from.map((v, i) => v + (to[i] - v) * x));\n}\nfunction schemeCss(scheme, t) {\n const [r, g, b] = schemeRgb(scheme, t).map(Math.round);\n return `rgb(${r}, ${g}, ${b})`;\n}\nfunction schemeColor(scheme, t) {\n const [r, g, b] = schemeRgb(scheme, t).map(Math.round);\n return packAbgr(r, g, b, 255);\n}\nfunction encodedColor({ field, scheme }, t) {\n return field === \"walk\" ? schemeColor(scheme, BASE_T) : t === void 0 ? NO_VALUE_COLOR : schemeColor(scheme, t);\n}\nvar NO_VALUE_COLOR = REFERENCE_RAMP_ALT_COLOR;\nfunction encodingStops({ field, scheme }) {\n const stop = (offset, t) => ({\n offset,\n color: schemeCss(scheme, t),\n abgr: schemeColor(scheme, t)\n });\n return field === \"walk\" ? [stop(0, BASE_T)] : Array.from({ length: 7 }, (_, i) => stop(i / 6, i / 6));\n}\nfunction encodingSwatchCss(encoding) {\n const stops = encodingStops(encoding);\n return stops.length === 1 ? stops[0].color : `linear-gradient(to right, ${stops.map((s) => `${s.color} ${s.offset * 100}%`).join(\", \")})`;\n}\n\nexport {\n WALK_FIELDS,\n WALK_SCHEMES,\n resolveEncoding,\n schemeRgb,\n schemeCss,\n schemeColor,\n encodedColor,\n NO_VALUE_COLOR,\n encodingStops,\n encodingSwatchCss\n};\n"],
|
|
5
|
-
"mappings": "gMACA,SAASA,EAASC,EAAGC,EAAGC,EAAGC,EAAG,CAC5B,OAAQA,GAAK,GAAKD,GAAK,GAAKD,GAAK,EAAID,KAAO,CAC9C,CACA,SAASI,GAAQC,EAAG,CAClB,OAAOA,EAAI,GACb,CACA,SAASC,GAAUD,EAAG,CACpB,OAAOA,IAAM,EAAI,GACnB,CACA,SAASE,GAASF,EAAG,CACnB,OAAOA,IAAM,GAAK,GACpB,CACA,SAASG,GAAUH,EAAG,CACpB,OAAOA,IAAM,GAAK,GACpB,CACA,SAASI,GAAcJ,EAAG,CACxB,MAAO,QAAQD,GAAQC,CAAC,CAAC,IAAIC,GAAUD,CAAC,CAAC,IAAIE,GAASF,CAAC,CAAC,IAAIG,GAAUH,CAAC,EAAI,GAAG,GAChF,CACA,SAASK,GAAuBL,EAAGM,EAAO,CACxC,MAAO,QAAQ,KAAK,MAAMN,EAAE,CAAC,EAAI,GAAG,CAAC,IAAI,KAAK,MAAMA,EAAE,CAAC,EAAI,GAAG,CAAC,IAAI,KAAK,MAAMA,EAAE,CAAC,EAAI,GAAG,CAAC,IAAIM,CAAK,GACpG,CACA,SAASC,EAASP,EAAGM,EAAO,CAC1B,OAAOZ,EACLK,GAAQC,CAAC,EACTC,GAAUD,CAAC,EACXE,GAASF,CAAC,EACV,KAAK,MAAMG,GAAUH,CAAC,EAAIM,CAAK,CACjC,CACF,CACA,SAASE,GAAaR,EAAGS,EAAQ,CAC/B,IAAMC,EAAI,KAAK,IAAI,EAAG,EAAI,EAAID,CAAM,EAC9BE,EAAQC,GAAM,KAAK,MAAMA,GAAK,IAAMA,GAAKF,CAAC,EAChD,OAAOhB,EACLiB,EAAKZ,GAAQC,CAAC,CAAC,EACfW,EAAKV,GAAUD,CAAC,CAAC,EACjBW,EAAKT,GAASF,CAAC,CAAC,EAChBG,GAAUH,CAAC,CACb,CACF,CCtCA,IAAIa,GAAyB,IACzBC,EAA4B,GAC5BC,GAA2B,GAC/B,SAASC,GAAWC,EAAK,CACvB,MAAO,OAAOA,CAAG,KAAKH,EAA4B,GAAG,MAAMC,GAA2B,GAAG,IAC3F,CACA,SAASG,GAAUC,EAAM,CACvB,OAAO,KAAK,IAAI,EAAG,KAAK,IAAI,EAAGA,CAAI,CAAC,EAAIN,EAC1C,CACA,SAASO,GAAQC,EAAIC,EAAQ,CAC3B,OAAOJ,IAAWG,EAAKC,EAAO,QAAUA,EAAO,IAAMA,EAAO,MAAM,CACpE,CACA,IAAIC,GAAe,GACnB,SAASC,GAAU,CAAE,MAAAC,EAAO,GAAAJ,EAAI,SAAAK,CAAS,EAAGJ,EAAQ,CAClD,IAAMK,EAAMP,GAAQK,EAAOH,CAAM,EAC3BM,EAAOR,GAAQK,EAAQJ,EAAIC,CAAM,EACvC,GAAIM,EAAOD,EAAM,EACf,MAAO,EAAEA,EAAMC,GAAQ,CAAC,EAE1B,IAAMC,EAAI,KAAK,MAAMD,EAAOD,GAAOJ,EAAY,EAC/C,OAAO,MAAM,KACX,CAAE,OAAQM,EAAI,CAAE,EAChB,CAACC,EAAG,IAAMV,GAAQK,EAAQJ,GAAMK,EAAWG,EAAI,EAAI,GAAKA,EAAGP,CAAM,CACnE,CACF,CACA,IAAIS,GAAsB,MAAM,KAAK,CAAE,OAAQ,CAAE,EAAG,CAACD,EAAGE,KAAO,CAC7D,OAAQA,EAAI,EACZ,MAAOhB,GAAWgB,EAAI,EAAInB,EAAsB,CAClD,EAAE,EACEoB,GAAoB,6BAA6BF,GAAoB,IACtEG,GAAM,GAAGA,EAAE,KAAK,IAAIA,EAAE,OAAS,GAAG,GACrC,EAAE,KAAK,IAAI,CAAC,ICxBZ,IAAIC,GAAQC,GAAS,EACrB,SAASC,GAAeC,EAAeC,EAAOC,EAAMC,EAAWC,EAASC,EAAQC,EAAW,GAAO,CAChG,OAAOT,GACLG,EACA,CAACC,EAAOC,EAAK,OAAQA,EAAK,OAAQC,EAAWE,EAAQC,EAAUF,CAAO,EACtE,IAAM,CACJ,IAAMG,EAAyB,IAAI,IACnC,QAASC,EAAK,EAAGA,EAAKP,EAAM,MAAM,OAAQO,IAAM,CAC9C,IAAMC,EAAWN,GAAW,IAAIK,CAAE,EAC5BE,EAAaC,GACjBX,EACAC,EAAM,MAAMO,CAAE,EACdN,EACAO,EACAA,GAAYJ,IAASG,CAAE,EACvBF,CACF,EACII,GACFH,EAAO,IAAIC,EAAIE,CAAU,CAE7B,CACA,OAAOH,CACT,CACF,CACF,CCxBA,IAAIK,GAAoB,EACxB,SAASC,GAAeC,EAAO,CAC7B,IAAMC,EAA4B,IAAI,IACtC,QAAWC,KAAQF,EAAM,MACvBG,EAAOF,EAAWC,EAAK,KAAMA,EAAK,EAAE,EACpCC,EAAOF,EAAWC,EAAK,GAAIA,EAAK,IAAI,EAEtC,OAAOD,CACT,CACA,SAASG,GAAa,CACpB,OAAAC,EACA,SAAAC,EACA,UAAAL,CACF,EAAG,CACD,IAAMM,EAAuB,IAAI,IAAI,CAACF,CAAM,CAAC,EACzCG,EAAW,CAACH,CAAM,EAClBI,EAAS,IACTC,EAAW,KACf,QAASC,EAAM,EAAGA,EAAMb,IAAqBU,EAAS,OAAS,EAAGG,IAAO,CACvE,IAAMC,EAAO,CAAC,EACd,QAAWC,KAAML,EAAU,CACzB,IAAMM,EAAWb,EAAU,IAAIY,CAAE,EACjC,GAAIC,GACF,QAAWC,KAAcD,EACvB,GAAI,CAACP,EAAK,IAAIQ,CAAU,EAAG,CACzBR,EAAK,IAAIQ,CAAU,EACnB,IAAMC,EAAWV,EAAS,IAAIS,CAAU,EACpCC,GAAYC,EAAWD,CAAQ,GACjCP,EAAS,KAAK,IAAIA,EAAQO,EAAS,OAAO,MAAQA,EAAS,MAAM,EACjEN,EAAW,KAAK,IAAIA,EAAUM,EAAS,OAAO,KAAK,GAEnDJ,EAAK,KAAKG,CAAU,CAExB,EAGN,CACAP,EAAWI,CACb,CACA,OAAOH,IAAW,IAAW,OAAS,CACpC,MAAO,KAAK,IAAIA,EAAQC,CAAQ,EAChC,IAAK,KAAK,IAAID,EAAQC,CAAQ,CAChC,CACF,CACA,SAASQ,GAAkB,CACzB,OAAAb,EACA,SAAAC,EACA,UAAAL,CACF,EAAG,CACD,IAAMkB,EAAOb,EAAS,IAAID,CAAM,EAChC,OAAOc,IAAS,OAAS,OAASF,EAAWE,CAAI,EAAI,CAAE,MAAOA,EAAK,OAAO,MAAO,IAAKA,EAAK,OAAO,MAAQA,EAAK,MAAO,EAAIf,GAAa,CAAE,OAAAC,EAAQ,SAAAC,EAAU,UAAAL,CAAU,CAAC,CACxK,CACA,SAASmB,GAAmBpB,EAAO,CACjC,IAAMM,EAAW,IAAI,IAAIN,EAAM,MAAM,IAAKmB,GAAS,CAACA,EAAK,GAAIA,CAAI,CAAC,CAAC,EAC7DlB,EAAYF,GAAeC,CAAK,EAChCqB,EAA4B,IAAI,IACtC,QAAWF,KAAQnB,EAAM,MAAO,CAC9B,IAAMsB,EAAOJ,GAAkB,CAAE,OAAQC,EAAK,GAAI,SAAAb,EAAU,UAAAL,CAAU,CAAC,EACnEqB,GACFD,EAAU,IAAIF,EAAK,IAAKG,EAAK,MAAQA,EAAK,KAAO,CAAC,CAEtD,CACA,OAAOD,CACT,CC7BA,IAAIE,GAAqBC,EAAS,IAAK,IAAK,IAAK,GAAG,EAChDC,GAA2BD,EAAS,IAAK,IAAK,IAAK,GAAG,EACtDE,GAAsBF,EAAS,GAAI,GAAI,GAAI,GAAG,EAC9CG,GAA4B,IAC5BC,GAAsBJ,EAAS,IAAK,IAAK,IAAK,GAAG,EACjDK,GAAc,IACdC,GAAoBC,GACtBC,EAASJ,GAAqBC,EAAW,CAC3C,EACII,GAAmB,GACnBC,GAAqB,IACrBC,GAAmB,EACnBC,GAAuB,KACvBC,GAAuB,EACvBC,GAA2B,EAC3BC,GAA2B,EAC3BC,GAA+B,IAC/BC,GAAiB,EACjBC,GAAgB,GAChBC,GAAkB,IAClBC,GAAoB,GACpBC,GAAc,KACdC,EAAQ,IAAI,aAAa,EAAE,EAC3BC,EAAQ,IAAI,aAAa,EAAE,EAC3BC,EAAS,IAAI,aAAa,EAAE,EAChC,SAASC,GAAYC,EAAQC,EAAOC,EAAMC,EAAS,CACjD,IAAMC,EAAOJ,EAAO,OAASN,GAAoB,EAC7CE,EAAM,OAASQ,IACjBR,EAAQ,IAAI,aAAaQ,CAAI,EAC7BP,EAAQ,IAAI,aAAaO,CAAI,EAC7BN,EAAS,IAAI,aAAaM,CAAI,GAEhC,IAAMC,EAAOL,EAAOA,EAAO,OAAS,CAAC,EACrCJ,EAAM,CAAC,EAAIS,EAAK,GAChBR,EAAM,CAAC,EAAIQ,EAAK,GAChBP,EAAO,CAAC,EAAI,EACZ,IAAIQ,EAAI,EACJC,EAAK,EACT,QAASC,EAAKR,EAAO,OAAS,EAAGQ,GAAM,GAAKD,EAAKJ,EAASK,IAAM,CAC9D,IAAMC,EAAIT,EAAOQ,CAAE,EACnB,QAASE,EAAI,EAAGA,GAAKhB,IAAqBa,EAAKJ,EAASO,IAAK,CAC3D,IAAMC,EAAI,GAAKD,EAAIhB,KAAsB,EACnCkB,EAAI,EAAID,EACRE,EAAID,EAAIA,EAAIA,EACZE,EAAI,EAAIF,EAAIA,EAAID,EAChBI,EAAI,EAAIH,EAAID,EAAIA,EAChBK,EAAIL,EAAIA,EAAIA,EACZM,EAAIJ,EAAIJ,EAAE,GAAKK,EAAIL,EAAE,IAAMM,EAAIN,EAAE,IAAMO,EAAIP,EAAE,GAC7CS,EAAIL,EAAIJ,EAAE,GAAKK,EAAIL,EAAE,IAAMM,EAAIN,EAAE,IAAMO,EAAIP,EAAE,GACnDF,GAAM,KAAK,MAAMU,EAAIrB,EAAMU,EAAI,CAAC,GAAIY,EAAIrB,EAAMS,EAAI,CAAC,GAAKJ,CAAI,EAAID,EAChEL,EAAMU,CAAC,EAAIW,EACXpB,EAAMS,CAAC,EAAIY,EACXpB,EAAOQ,CAAC,EAAIC,EACZD,GACF,CACF,CACA,OAAOA,CACT,CACA,SAASa,GAAYb,EAAGC,EAAI,CAC1B,IAAIa,EAAI,EACR,KAAOA,EAAId,EAAI,GAAKR,EAAOsB,CAAC,EAAIb,GAC9Ba,IAEF,IAAMC,EAAOvB,EAAOsB,EAAI,CAAC,EACnBE,EAAKxB,EAAOsB,CAAC,EACbG,EAAID,EAAKD,EAAO,KAAK,IAAI,EAAG,KAAK,IAAI,GAAId,EAAKc,IAASC,EAAKD,EAAK,CAAC,EAAI,EAC5E,MAAO,CACL,EAAGzB,EAAMwB,EAAI,CAAC,GAAKxB,EAAMwB,CAAC,EAAIxB,EAAMwB,EAAI,CAAC,GAAKG,EAC9C,EAAG1B,EAAMuB,EAAI,CAAC,GAAKvB,EAAMuB,CAAC,EAAIvB,EAAMuB,EAAI,CAAC,GAAKG,CAChD,CACF,CACA,SAASC,GAAa,CACpB,OAAAxB,EACA,cAAAyB,EACA,cAAAC,EACA,MAAAzB,EACA,KAAAC,EACA,MAAAyB,EACA,aAAAC,EACA,WAAAC,CACF,EAAG,CACD,IAAMC,EAAS3C,GAAuBC,GAA2BqC,EAC3DnB,EAAIP,GACRC,EACAC,EACAC,EACAwB,EAAgBnC,GAAiBuC,CACnC,EACMC,EAAWjC,EAAOQ,EAAI,CAAC,EACvB0B,EAAWD,EAAWpC,GACtBsC,EAAQ,KAAK,IAAIP,EAAeK,CAAQ,EACxCG,EAAMF,EAAW,KAAK,IAAI,EAAGG,GAAWP,EAAc3B,EAAOA,EAAQC,CAAI,EAAI4B,CAAM,EAAI,EAC7F,GAAIE,EAAWE,EAAM1C,GAAgBuC,EAAWxC,GAAiBuC,EAC/D,OAEF,IAAMM,EAAMjB,GAAYb,EAAG2B,CAAK,EAC1BI,EAAOlB,GAAYb,EAAG2B,EAAQH,EAASI,CAAG,EAC1CI,EAAKF,EAAI,EAAIC,EAAK,EAClBE,GAAMH,EAAI,EAAIC,EAAK,GAAKnC,EACxBG,EAAOL,EAAOA,EAAO,OAAS,CAAC,EACrC,MAAO,CACL,EAAGoC,EAAI,EACP,EAAGA,EAAI,EACP,MAAO,KAAK,MAAME,EAAIC,CAAE,EAAI,EAAI,KAAK,MAAMA,EAAID,CAAE,EAAIE,GAAWnC,EAAMH,EAAMuC,GAAcpC,EAAMwB,CAAU,CAAC,EAC3G,OAAQC,EAASI,EACjB,WAAY7C,GAA2BC,GAA+BmC,GAAiBS,EACvF,MAAAP,CACF,CACF,CACA,SAASe,GAAYC,EAAMC,EAAUC,EAAOC,EAAU5C,EAAM,CAC1D,GAAI4C,EACF,OAAOC,GAAYJ,CAAI,EAAE,OAAS,QAEpC,IAAMK,EAAOC,GAAM,KAAK,MAAMA,EAAE,EAAIJ,EAAM,IAAKI,EAAE,EAAIJ,EAAM,IAAM3C,CAAI,EACrE,OAAO8C,EAAIJ,EAAS,CAAC,CAAC,EAAII,EAAIJ,EAAS,GAAG,EAAE,CAAC,CAC/C,CACA,SAASH,GAAcpC,EAAMwB,EAAY,CACvC,GAAIA,EAAW,OAAS,EACtB,OAEF,IAAMgB,EAAQhB,EAAW,CAAC,EACpBqB,EAAQrB,EAAWA,EAAW,OAAS,CAAC,EAE9C,OADgB,KAAK,MAAMxB,EAAK,GAAKwC,EAAM,EAAGxC,EAAK,GAAKwC,EAAM,CAAC,GAAK,KAAK,MAAMxC,EAAK,GAAK6C,EAAM,EAAG7C,EAAK,GAAK6C,EAAM,CAAC,EAClG,CAAE,KAAML,EAAO,GAAIhB,EAAW,CAAC,CAAE,EAAI,CAAE,KAAMqB,EAAO,GAAIrB,EAAWA,EAAW,OAAS,CAAC,CAAE,CAC7G,CACA,SAASsB,GAAOtC,EAAGC,EAAG,CACpB,GAAI,CAACD,GAAK,CAACC,EACT,OAEF,IAAMwB,EAAKxB,EAAE,EAAID,EAAE,EACb0B,EAAKzB,EAAE,EAAID,EAAE,EACbuC,EAAM,KAAK,MAAMd,EAAIC,CAAE,EAC7B,OAAOa,EAAM,EAAI,CAAE,GAAI,CAACb,EAAKa,EAAK,GAAId,EAAKc,CAAI,EAAI,MACrD,CACA,SAASC,GAAeC,EAAWpD,EAAO,EAAG,CAC3C,IAAMqD,EAASrD,IAAS,EAAIoD,EAAYA,EAAU,IAAKL,IAAO,CAAE,EAAGA,EAAE,EAAG,EAAGA,EAAE,EAAI/C,CAAK,EAAE,EACxF,OAAOqD,EAAO,IAAI,CAACN,EAAG7B,IAAM,CAC1B,IAAMoC,EAASL,GAAOI,EAAOnC,EAAI,CAAC,EAAG6B,CAAC,EAChCQ,EAAQN,GAAOF,EAAGM,EAAOnC,EAAI,CAAC,CAAC,EACrC,GAAIoC,GAAUC,EAAO,CACnB,IAAMC,GAAMF,EAAO,GAAKC,EAAM,IAAM,EAC9BE,GAAMH,EAAO,GAAKC,EAAM,IAAM,EAC9BG,EAAMJ,EAAO,GAAKE,EAAKF,EAAO,GAAKG,EACzC,OAAOC,EAAM,GAAM,CAAE,GAAIF,EAAKE,EAAK,GAAID,EAAKC,CAAI,EAAI,CAAE,GAAAF,EAAI,GAAAC,CAAG,CAC/D,CACA,OAAOH,GAAUC,GAAS,CAAE,GAAI,EAAG,GAAI,CAAE,CAC3C,CAAC,CACH,CACA,SAASI,GAAeN,EAAQO,EAASC,EAAU7D,EAAO,EAAG,CAC3D,OAAOqD,EAAO,IAAI,CAACN,EAAG7B,KAAO,CAC3B,EAAG6B,EAAE,EAAIa,EAAQ1C,CAAC,EAAE,GAAK2C,EACzB,EAAGd,EAAE,EAAIa,EAAQ1C,CAAC,EAAE,GAAK2C,EAAW7D,CACtC,EAAE,CACJ,CACA,SAAS8D,GAASC,EAAGC,EAAGC,EAAG,CACzB,IAAM1D,GAAK,EAAI,KAAK,IAAI,EAAI0D,EAAI,CAAC,GAAKD,EAChCjD,EAAIR,GAAK,EAAI,KAAK,IAAIwD,EAAI,GAAK,EAAI,CAAC,GACpCG,EAAID,EAAI1D,EAAI,EACd4D,EAAI,EACJC,EAAI,EACJxD,EAAI,EACR,OAAImD,EAAI,IACNI,EAAI5D,EACJ6D,EAAIrD,GACKgD,EAAI,KACbI,EAAIpD,EACJqD,EAAI7D,GACKwD,EAAI,KACbK,EAAI7D,EACJK,EAAIG,GACKgD,EAAI,KACbK,EAAIrD,EACJH,EAAIL,GACKwD,EAAI,KACbI,EAAIpD,EACJH,EAAIL,IAEJ4D,EAAI5D,EACJK,EAAIG,GAEC,CAACoD,EAAID,EAAGE,EAAIF,EAAGtD,EAAIsD,CAAC,CAC7B,CACA,SAASG,GAAQN,EAAGC,EAAGC,EAAGtD,EAAI,EAAG,CAC/B,GAAM,CAAC,EAAGyD,EAAGxD,CAAC,EAAIkD,GAASC,EAAGC,EAAGC,CAAC,EAClC,OAAO7F,EACL,KAAK,MAAM,EAAI,GAAG,EAClB,KAAK,MAAMgG,EAAI,GAAG,EAClB,KAAK,MAAMxD,EAAI,GAAG,EAClB,KAAK,MAAMD,EAAI,GAAG,CACpB,CACF,CACA,IAAI2D,GAAiB,CACnB,CAAC,GAAI,EAAG,EAAE,EACV,CAAC,GAAI,GAAI,GAAG,EACZ,CAAC,GAAI,IAAK,GAAG,EACb,CAAC,GAAI,IAAK,EAAE,EACZ,CAAC,IAAK,IAAK,EAAE,CACf,EACIC,GAAuB,CACzB,CAAC,IAAK,IAAK,EAAE,EACb,CAAC,IAAK,GAAI,GAAG,CACf,EACIC,GAAuB,CACzB,CAAC,IAAK,GAAI,EAAE,EACZ,CAAC,GAAI,IAAK,GAAG,CACf,EACA,SAASC,GAAeC,EAAOjE,EAAG,CAChC,IAAMkE,EAAM,KAAK,IAAI,EAAG,KAAK,IAAI,EAAGlE,CAAC,CAAC,GAAKiE,EAAM,OAAS,GACpDxD,EAAI,KAAK,IAAIwD,EAAM,OAAS,EAAG,KAAK,MAAMC,CAAG,CAAC,EAC9CtD,EAAIsD,EAAMzD,EACVP,EAAI+D,EAAMxD,CAAC,EACXN,EAAI8D,EAAMxD,EAAI,CAAC,EACrB,OAAO9C,EACL,KAAK,MAAMuC,EAAE,CAAC,GAAKC,EAAE,CAAC,EAAID,EAAE,CAAC,GAAKU,CAAC,EACnC,KAAK,MAAMV,EAAE,CAAC,GAAKC,EAAE,CAAC,EAAID,EAAE,CAAC,GAAKU,CAAC,EACnC,KAAK,MAAMV,EAAE,CAAC,GAAKC,EAAE,CAAC,EAAID,EAAE,CAAC,GAAKU,CAAC,EACnC,GACF,CACF,CACA,SAASuD,GAAUC,EAAKC,EAAO,CAC7B,OAAOT,GAAQU,GAAQF,CAAG,EAAGG,GAAiBC,GAAgBH,CAAK,CACrE,CACA,IAAII,GAAyB,CAAC,GAAI,GAAI,EAAE,EACpCC,GAA2B/G,EAAS,GAAG8G,GAAwB,GAAG,EAClEE,GAAyB,OAAOF,GAAuB,KAAK,IAAI,CAAC,IACrE,SAASG,GAAqBC,EAAOC,EAAQ,CAC3C,IAAMC,EAAYC,GAAmBH,CAAK,EAC1C,GAAIC,GAAUA,EAAO,IAAMA,EAAO,MAChC,MAAO,CACL,MAAOA,EAAO,MACd,KAAMA,EAAO,IAAMA,EAAO,MAC1B,UAAAC,CACF,EAEF,IAAIE,EAAM,IACNC,EAAM,KACV,QAAWC,KAAOJ,EAAU,OAAO,EACjCE,EAAM,KAAK,IAAIA,EAAKE,CAAG,EACvBD,EAAM,KAAK,IAAIA,EAAKC,CAAG,EAEzB,MAAO,CACL,MAAOF,IAAQ,IAAW,EAAIA,EAC9B,KAAMC,EAAMD,EAAMC,EAAMD,EAAM,EAC9B,UAAAF,CACF,CACF,CACA,SAASK,GAAwBP,EAAOQ,EAAe,CACrD,IAAIC,EAAW,IACXC,EAAW,KACXC,EAAY,IACZC,EAAY,KACZC,EAAU,EACd,QAAW/F,KAAKkF,EAAM,MAChBlF,EAAE,SAAW,QAAUA,EAAE,OAAO,KAAO+F,IACzCA,EAAU/F,EAAE,OAAO,MAEjBA,EAAE,MAAQ2F,IACZA,EAAW3F,EAAE,OAEXA,EAAE,MAAQ4F,IACZA,EAAW5F,EAAE,OAEXA,EAAE,OAAS6F,IACbA,EAAY7F,EAAE,QAEZA,EAAE,OAAS8F,IACbA,EAAY9F,EAAE,QAGlB,MAAO,CACL,SAAA2F,EACA,SAAAC,EACA,UAAAC,EACA,UAAAC,EACA,QAAAC,EACA,UAAWb,EAAM,MAAM,OACvB,cAAAQ,CACF,CACF,CACA,SAASM,GAAaC,EAAMC,EAAWC,EAAaC,EAAO,CACzD,OAAQD,EAAa,CACnB,IAAK,SACH,OAAO3B,GAAUyB,EAAK,GAAI,CAAC,EAC7B,IAAK,QACH,OAAO5B,GACLH,GACAkC,EAAM,SAAWA,EAAM,UAAYH,EAAK,MAAQG,EAAM,WAAaA,EAAM,SAAWA,EAAM,UAAY,EACxG,EACF,IAAK,cACH,OAAO/B,GACLD,GACAgC,EAAM,UAAYA,EAAM,WAAaH,EAAK,OAASG,EAAM,YAAcA,EAAM,UAAYA,EAAM,WAAa,EAC9G,EAIF,IAAK,cACH,OAAOH,EAAK,SAAW,OAASjI,EAAS,IAAK,IAAK,IAAK,GAAG,EAAIiI,EAAK,OAAO,OAAS,EAAIjI,EAAS,GAAI,IAAK,IAAK,GAAG,EAAIqG,GACpHF,GACAiC,EAAM,QAAU,GAAKH,EAAK,OAAO,KAAO,IAAMG,EAAM,QAAU,GAAK,CACrE,EAGF,IAAK,qBAAsB,CACzB,IAAMC,EAAOD,EAAM,cACnB,GAAIC,IAAS,OACX,OAAOrI,EAAS,IAAK,IAAK,IAAK,GAAG,EAEpC,IAAMwH,EAAMa,EAAK,UAAU,IAAIJ,EAAK,EAAE,EACtC,MAAI,CAACK,EAAWL,CAAI,GAAKT,IAAQ,OACxBT,GAEFd,GACLsC,IAAWf,EAAMa,EAAK,OAASA,EAAK,IAAI,EACxCG,EACAC,EACF,CACF,CACA,IAAK,OACH,OAAOzI,EAAS,IAAK,IAAK,IAAK,GAAG,EACpC,IAAK,UAAW,CACd,IAAM0I,EAAMN,EAAM,UAAY,EAAIF,EAAYE,EAAM,UAAY,IAAM,EACtE,OAAOnC,GAAQyC,EAAK,IAAM,EAAG,CAC/B,CACA,QACE,OAAO1I,EAAS,GAAI,IAAK,IAAK,GAAG,CACrC,CACF,CACA,SAASkE,GAAW/B,EAAGP,EAAO,EAAG+G,EAAM,CACrC,IAAMC,EAAa,CACjB,CAACzG,EAAE,IAAKA,EAAE,GAAG,EACb,CAACA,EAAE,IAAKA,EAAE,GAAG,EACb,CAACA,EAAE,GAAIA,EAAE,EAAE,CACb,EACA,OAAW,CAACQ,EAAGC,CAAC,IAAKgG,EAAY,CAC/B,IAAM5E,EAAK7B,EAAE,GAAKQ,EACZsB,GAAM9B,EAAE,GAAKS,GAAKhB,EACxB,GAAI,KAAK,MAAMoC,EAAIC,CAAE,EAAI,EACvB,OAAO,KAAK,MAAMA,EAAID,CAAE,CAE5B,CACA,OAAO2E,EAAO,KAAK,OAAOA,EAAK,GAAG,EAAIA,EAAK,KAAK,GAAK/G,EAAM+G,EAAK,GAAG,EAAIA,EAAK,KAAK,CAAC,EAAI,CACxF,CACA,SAASE,GAAmBvE,EAAUwE,EAAQ,CAC5C,IAAIC,EAAM,GACV,QAASjG,EAAI,EAAGA,EAAIwB,EAAS,OAAS,GAAK,CAACyE,EAAKjG,IAAK,CACpD,IAAMP,EAAI+B,EAASxB,CAAC,EACdN,EAAI8B,EAASxB,EAAI,CAAC,EACxBiG,EAAM,KAAK,IAAIxG,EAAE,EAAGC,EAAE,CAAC,GAAKsG,EAAO,MAAQ,KAAK,IAAIvG,EAAE,EAAGC,EAAE,CAAC,GAAKsG,EAAO,MAAQ,KAAK,IAAIvG,EAAE,EAAGC,EAAE,CAAC,GAAKsG,EAAO,MAAQ,KAAK,IAAIvG,EAAE,EAAGC,EAAE,CAAC,GAAKsG,EAAO,IACpJ,CACA,OAAOC,CACT,CACA,SAASC,GAAiBtH,EAAQoH,EAAQ,CACxC,QAAW3G,KAAKT,EAAQ,CACtB,IAAMuH,EAAQ,KAAK,IAAI9G,EAAE,GAAIA,EAAE,IAAKA,EAAE,IAAKA,EAAE,EAAE,EACzC+G,EAAQ,KAAK,IAAI/G,EAAE,GAAIA,EAAE,IAAKA,EAAE,IAAKA,EAAE,EAAE,EACzCgH,EAAQ,KAAK,IAAIhH,EAAE,GAAIA,EAAE,IAAKA,EAAE,IAAKA,EAAE,EAAE,EACzCiH,EAAQ,KAAK,IAAIjH,EAAE,GAAIA,EAAE,IAAKA,EAAE,IAAKA,EAAE,EAAE,EAC/C,GAAI+G,GAASJ,EAAO,MAAQG,GAASH,EAAO,MAAQM,GAASN,EAAO,MAAQK,GAASL,EAAO,KAC1F,MAAO,EAEX,CACA,MAAO,EACT,CACA,IAAIO,GAA+B,IAAI,QACvC,SAASC,GAAcpC,EAAO,CAC5B,IAAIqC,EAAUF,GAAa,IAAInC,CAAK,EACpC,GAAI,CAACqC,EAAS,CACZ,IAAMC,EAAmB,CAAC,EACpBC,EAAkC,IAAI,IACtCC,EAAQxC,EAAM,OAAS,CAAC,EAC9B,QAASpE,EAAI,EAAGA,EAAI4G,EAAM,OAAQ5G,IAChC0G,EAAiB,KACfvD,GACE0D,GAAU7G,EAAG4G,EAAM,MAAM,EACzB9C,GACAC,GACA,GACF,CACF,EACK4C,EAAgB,IAAIC,EAAM5G,CAAC,EAAE,IAAI,GACpC2G,EAAgB,IAAIC,EAAM5G,CAAC,EAAE,KAAMA,CAAC,EAGxCyG,EAAU,CAAE,iBAAAC,EAAkB,gBAAAC,CAAgB,EAC9CJ,GAAa,IAAInC,EAAOqC,CAAO,CACjC,CACA,OAAOA,CACT,CACA,IAAIK,GAA6B,IAAI,QACjCC,GAA2B,IAAI,IACnC,SAASC,GAAY5C,EAAO,CAC1B,IAAI6C,EAAQH,GAAW,IAAI1C,CAAK,EAChC,GAAI,CAAC6C,EAAO,CACVA,EAAwB,IAAI,IAC5B,OAAW,CAACC,EAASC,CAAI,IAAM/C,EAAM,OAAS,CAAC,GAAG,QAAQ,EACxD,QAAWgD,KAAU,IAAI,IAAID,EAAK,OAAO,EACvCE,EAAOJ,EAAOG,EAAQF,CAAO,EAGjCJ,GAAW,IAAI1C,EAAO6C,CAAK,CAC7B,CACA,OAAOA,CACT,CACA,SAASK,GAAcC,EAAS,CAC9B,GAAM,CACJ,cAAAC,EACA,MAAApD,EACA,SAAAqD,EACA,YAAApC,EACA,gBAAAqC,EACA,mBAAAC,EACA,UAAAC,EACA,UAAAC,EAAY,UACZ,UAAAC,EACA,KAAAC,EACA,aAAAC,EACA,eAAAC,EACA,cAAArD,EACA,UAAAsD,EACA,eAAAC,EACA,SAAAzG,EACA,YAAA0G,EACA,QAAAC,EAAU,CACZ,EAAId,EACEe,GAAYT,IAAc,QAAUU,GAAUnE,CAAK,EAAI,EACvDvF,GAAQkJ,EAAK,OACbS,GAAmB,EAAI,KAAK,IAAI3J,GAAOf,EAAoB,EAC3DgB,EAAO2J,GAAOV,CAAI,EAClBW,GAAeC,GACnBnB,EACApD,EACA2D,EACAG,EACAG,EACAF,EACAzG,CACF,EACMkH,EAAc,CAAC,EACfC,GAAiC,IAAI,IACrCC,GAAS,CAAC,EACVC,GAA4B,IAAI,IAChCC,EAAa,CAAC,EACdC,GAAgC,IAAI,IACpCC,GAAavE,GAAwBP,EAAOQ,CAAa,EACzDuE,GAA+B,IAAI,IACzC,GAAI9D,IAAgB,UAClB,QAASrF,EAAI,EAAGA,EAAIoE,EAAM,MAAM,OAAQpE,IACtCmJ,GAAa,IAAI/E,EAAM,MAAMpE,CAAC,EAAE,GAAIA,CAAC,EAGzC,IAAMoJ,GAAYhF,EAAM,OAAO,QAAU,EACnC,CAAE,iBAAAsC,GAAkB,gBAAAC,EAAgB,EAAIH,GAAcpC,CAAK,EAC3DiF,GAAmBC,GAAS5C,GAAiBC,GAAgB,IAAI2C,CAAI,GAAK,EAAE,GAAKnM,GACjFoM,GAAgB3B,GAAawB,GAAY,EAAIpC,GAAY5C,CAAK,EAAI2C,GAClEyC,GAAa3K,IAASmJ,EAAe,EAAI3J,IAC/C,QAASoL,EAAK,EAAGA,EAAKrF,EAAM,MAAM,OAAQqF,IAAM,CAC9C,IAAMlI,EAAO6C,EAAM,MAAMqF,CAAE,EACrBjJ,EAAegH,EAAcjG,EAAK,IAAI,EACtCd,EAAa+G,EAAcjG,EAAK,EAAE,EAClCmI,EAAahB,GAAa,IAAIe,CAAE,EACtC,GAAI,CAACjJ,GAAc,QAAU,CAACC,GAAY,QAAU,CAACiJ,GAActB,GAAa,IAAIqB,CAAE,EACpF,SAEF,IAAME,EAAWpI,EAAK,SAAS,QAAU,EAEnCqI,EADW1B,GAAW,IAAIuB,CAAE,IACF,OAC1BI,EAAS/B,GAAW,YAAY,IAAI2B,CAAE,GAAK,GAC3CK,EAAgBnC,EAAqB,GAAKiC,EAAavM,GAA4B,GACnF0M,EAAUnC,GAAa+B,EAAW,EAClCK,EAAiBJ,EAAaxM,GAAsBH,GACpDgN,EAAY,CAACnC,GAAa+B,EAASG,EAAiBtM,EAASsM,EAAgBzM,EAAW,EACxF2M,EAASzC,EAAS,IAAIlG,EAAK,EAAE,EAC7B4I,GAAiBD,EAASE,GAAYF,EAAQxC,EAAiBG,EAAWS,EAAS,EAAIZ,GAAmB,EAC1G2C,EAAUzL,GAAW,CAACqJ,GAAkB/B,GAAiBtH,EAAQqJ,CAAc,EAC/EqC,EAAS,CAAC1L,EAAQ2L,EAAWhK,IAAU,CAC3C,GAAIqJ,EACF,QAAWY,KAAQC,GAAW7L,EAAQjB,GAAmBkB,EAAK,EAC5DmK,EAAW,KAAK,CAAE,OAAQwB,EAAM,UAAAD,EAAW,MAAAhK,CAAM,CAAC,OAGpDyI,EAAW,KAAK,CAAE,OAAApK,EAAQ,UAAA2L,EAAW,MAAAhK,CAAM,CAAC,CAEhD,EACMmK,EAAW9L,GAAW,CAC1B,IAAM+L,EAAOnB,GAAapJ,GAAa,CACrC,OAAAxB,EACA,cAAe+I,EAAqB,EACpC,cAAewC,EACf,MAAAtL,GACA,KAAAC,EACA,MAAOmL,EACP,aAAAzJ,EACA,WAAAC,CACF,CAAC,EAAI,OACDkK,GACF7B,GAAO,KAAK6B,CAAI,CAEpB,EACMC,EAAiB5B,EAAW,OAC5B6B,GAAa/B,GAAO,OACpBgC,GAAYhD,GAAW,MAAM,QAAU+B,GAAU,CAACE,EAAUjC,EAAU,MAAQ,OACpF,GAAIgD,GAAW,CACb,IAAMC,EAAWtD,EAAS,IAAIlG,EAAK,IAAI,EACjCyJ,EAAQD,EAAWX,GAAYW,EAAUrD,EAAiBG,EAAWS,EAAS,EAAIZ,EAClFuD,EAAY,KAAK,IAAID,EAAQF,GAAU,OAAQjN,EAAgB,EAC/D6E,EAAUT,GAAezB,EAAc1B,CAAI,EAC3C,CAAE,GAAAwD,GAAI,GAAAC,CAAG,GAAKjB,GAClBC,EACAf,EACAkJ,EAAW,CAAC,EACZhI,EACA5C,CACF,EAAI4D,EAAQ,CAAC,EAAIA,EAAQ,GAAG,EAAE,IAAM,CAClC,GAAI,EACJ,GAAI,CACN,EACAoI,GAAU,QAAQ,CAACI,EAAMC,KAAS,CAChC,GAAI,CAACD,EAAK,YAAY,IAAIzB,CAAE,EAC1B,OAEF,IAAM2B,IAAUD,IAAQL,GAAU,OAAS,GAAK,GAAKG,EAAYzC,GAC3D5J,GAASwM,KAAW,EAAI1B,EAAa2B,GAAgB3B,EAAYpH,GAAK8I,GAAQ7I,EAAK6I,GAAStM,CAAI,EAClGuL,EAAOzL,EAAM,GACf0L,EACE1L,GACAqM,EAAY,EACZC,EAAK,OAAO,IAAI3J,EAAK,IAAI,GAAK2J,EAAK,OAAO,IAAI3J,EAAK,EAAE,GAAKyI,CAC5D,CAEJ,CAAC,EACGhB,EAAW,OAAS4B,GACtBF,EAAQhB,CAAU,CAEtB,SAAW,CAACK,EACNM,EAAOX,CAAU,IACnBY,EAAOZ,EAAYI,EAAeG,CAAS,EAC3CS,EAAQhB,CAAU,OAEf,CACL,IAAM4B,EAAUC,GACd/K,EACAC,EACAkJ,EACA5B,CACF,EACMyD,EAAc,KAAK,OAAO7B,EAAW,GAAK,CAAC,EACjD,QAASzC,EAAU,EAAGA,EAAUyC,EAAUzC,IAAW,CACnD,GAAM,CAAE,EAAArH,EAAG,EAAAC,EAAE,EAAIwL,EAAQpE,CAAO,EAC1BtI,EAASiB,IAAM,GAAKC,KAAM,EAAI4J,EAAa2B,GAAgB3B,EAAY7J,EAAGC,EAAC,EACjF,GAAI,CAACuK,EAAOzL,CAAM,EAChB,SAEF,IAAM6M,EAAcpC,GAAgB9H,EAAK,QAAQ2F,CAAO,CAAC,EACzDoD,EACE1L,EACAkL,EACAhC,IAAcA,EAAU,MAAM,OAAS,CAACA,EAAU,MAAM,IAAIvG,EAAK,QAAQ2F,CAAO,CAAC,EAAI,CAAC2C,GAAUnM,EAAS+N,EAAalO,EAAW,EAAIkO,CACvI,EACIvE,IAAYsE,GACdd,EAAQ9L,CAAM,CAElB,CACF,CACA,IAAM8M,GAAiB1C,EAAW,OAAS4B,EACvCc,GAAiB,GACnBzC,GAAc,IAAIQ,EAAI,CAAE,MAAOmB,EAAgB,MAAOc,EAAe,CAAC,EAExE,IAAMC,GAAa7C,GAAO,OAAS+B,GAC/Bc,GAAa,GACf5C,GAAU,IAAIU,EAAI,CAAE,MAAOoB,GAAY,MAAOc,EAAW,CAAC,CAE9D,CACA,OAAW,CAACvE,EAAQ5F,CAAQ,IAAK,OAAO,QAAQgG,CAAa,EAAG,CAC9D,IAAMrC,EAAOsC,EAAS,IAAIL,CAAM,EAIhC,GAHI,CAACjC,GAGD8C,GAAkB,CAAClC,GAAmBvE,EAAUyG,CAAc,EAChE,SAEF,IAAM2D,EAAQ9D,IAAc,QAAU,CAACA,EAAU,QAAQ,IAAIV,CAAM,EAC7DyE,EAAMD,EAAQtO,GAAsB4H,GACxCC,EACAgE,GAAa,IAAI/B,CAAM,GAAK,EAC5B/B,EACA6D,EACF,EACM3I,EAAQqL,EAAQlO,EAASmO,EAAKtO,EAAW,EAAIsO,EAC7Cb,GAAQZ,GAAYjF,EAAMuC,EAAiBG,EAAWS,EAAS,EAC/DwD,EAAgBd,GAAQ,EACxBe,EAAQnD,EAAY,OACpB3B,EAAQsC,GAAc,IAAInC,CAAM,EAChC4E,EAAYhB,GAAQ5B,GACpB6C,EAAWzK,EAAS,QAAU,EAC9B0K,EAAQpE,GAAW,MAAM,QAAU,CAAC8D,GAAS,CAAC3E,GAAO,OAASa,EAAU,MAAQ,OACtF,GAAImE,GAAYC,EAAO,CACrB,IAAMjB,EAAY,KAAK,IAAID,GAAQkB,EAAM,OAAQrO,EAAgB,EAC3D6E,EAAUT,GAAeT,EAAU1C,CAAI,EAC7CoN,EAAM,QAAQ,CAAChB,EAAMC,IAAS,CAC5B,IAAMgB,EAAYjB,EAAK,OAAO,IAAI9D,CAAM,EACxC,GAAI+E,IAAc,OAChB,OAEF,IAAMf,IAAUD,GAAQe,EAAM,OAAS,GAAK,GAAKjB,EAAYzC,GAC7DI,EAAY,KAAK,CACf,OAAQsD,EAAM,SAAW,EAAI1K,EAAWiB,GAAejB,EAAUkB,EAAS0I,GAAQtM,CAAI,EACtF,UAAWmM,EAAY,EACvB,MAAOkB,CACT,CAAC,CACH,CAAC,CACH,SAAWF,GAAYhF,GAAO,QAAU+E,GAAapO,GAAoB,CACvE,IAAM8E,EAAUT,GAAeT,EAAU1C,CAAI,EAC7C,QAAWsN,KAAQnF,EAAO,CACxB,IAAMmE,GAAUgB,GAAQhD,GAAY,GAAK,GAAK4C,EAAYxD,GACpD6D,EAAS3F,GAAiB0F,CAAI,GAAKjP,GACzCyL,EAAY,KAAK,CACf,OAAQnG,GAAejB,EAAUkB,EAAS0I,EAAQtM,CAAI,EACtD,UAAWkN,EAAY,EAEvB,MAAOJ,EAAQlO,EAAS2O,EAAQ9O,EAAW,EAAI8O,CACjD,CAAC,CACH,CACF,MAAWJ,GACTrD,EAAY,KAAK,CAAE,OAAQpH,EAAU,UAAWsK,EAAe,MAAAvL,CAAM,CAAC,EAExE,IAAM+L,EAAQ1D,EAAY,OAASmD,EAC/BO,EAAQ,GACVzD,GAAe,IAAIzB,EAAQ,CAAE,MAAA2E,EAAO,MAAAO,CAAM,CAAC,CAE/C,CACA,MAAO,CACL,YAAA1D,EACA,eAAAC,GACA,OAAAC,GACA,UAAAC,GACA,WAAAC,EACA,cAAAC,EACF,CACF,CC3pBA,IAAIsD,GAAc,CAChB,CAAE,MAAO,WAAY,MAAO,yBAA0B,EACtD,CAAE,MAAO,OAAQ,MAAO,yBAA0B,EAClD,CAAE,MAAO,YAAa,MAAO,oBAAqB,CACpD,EACIC,GAAe,CACjB,CAAE,MAAO,OAAQ,MAAO,gBAAiB,KAAM,CAAC,UAAW,SAAS,CAAE,EACtE,CAAE,MAAO,OAAQ,MAAO,eAAgB,KAAM,CAAC,UAAW,SAAS,CAAE,EACrE,CAAE,MAAO,MAAO,MAAO,gBAAiB,KAAM,CAAC,UAAW,SAAS,CAAE,EACrE,CAAE,MAAO,SAAU,MAAO,eAAgB,KAAM,CAAC,UAAW,SAAS,CAAE,EACvE,CAAE,MAAO,QAAS,MAAO,iBAAkB,KAAM,CAAC,UAAW,SAAS,CAAE,EACxE,CAAE,MAAO,UAAW,MAAO,SAAU,CACvC,EACIC,GAAW,CAAC,OAAQ,MAAO,SAAU,OAAO,EAC5CC,GAAWC,GAAMJ,GAAY,KAAMK,GAAMA,EAAE,QAAUD,CAAC,EACtDE,GAAYF,GAAMH,GAAa,KAAMM,GAAMA,EAAE,QAAUH,CAAC,EAC5D,SAASI,GAAgBC,EAAOC,EAAWC,EAAQC,EAAQ,GAAM,CAC/D,GAAM,CAAE,MAAOC,EAAY,OAAQC,CAAY,EAAIL,EAAM,OAAS,CAAC,EAC7DM,EAAQZ,GAAQU,CAAU,EAAIA,EAAaD,EAAQ,WAAa,OAChEI,EAASN,EAAY,OAASR,GAASS,EAAST,GAAS,MAAM,EAC/De,EAASX,GAASQ,CAAW,EAAIA,IAAgB,WAAaC,IAAU,YAAcC,EAASF,EAAcE,EACnH,MAAO,CAAE,MAAAD,EAAO,OAAAE,CAAO,CACzB,CACA,IAAIC,GAAS,GACTC,GAAYC,GAAMA,GAAK,OAAUA,EAAI,QAAUA,EAAI,MAAS,QAAU,IACtEC,GAAWD,GAAMA,GAAK,SAAWA,EAAI,MAAQ,MAAQA,IAAM,EAAI,KAAO,KAC1E,SAASE,GAAMC,EAAK,CAClB,GAAM,CAACC,EAAGC,EAAGC,CAAC,EAAI,CAAC,EAAG,EAAG,CAAC,EAAE,IACzBC,GAAMR,GAAS,OAAO,SAASI,EAAI,MAAMI,EAAGA,EAAI,CAAC,EAAG,EAAE,EAAI,GAAG,CAChE,EACMC,EAAI,KAAK,KAAK,YAAeJ,EAAI,YAAeC,EAAI,YAAeC,CAAC,EACpEG,EAAI,KAAK,KAAK,YAAeL,EAAI,YAAeC,EAAI,YAAeC,CAAC,EACpEnB,EAAI,KAAK,KAAK,YAAeiB,EAAI,YAAeC,EAAI,YAAeC,CAAC,EAC1E,MAAO,CACL,YAAeE,EAAI,WAAcC,EAAI,YAAetB,EACpD,aAAeqB,EAAI,YAAcC,EAAI,YAAetB,EACpD,YAAeqB,EAAI,YAAeC,EAAI,WAActB,CACtD,CACF,CACA,SAASuB,GAAW,CAACC,EAAGC,EAAGN,CAAC,EAAG,CAC7B,IAAME,GAAKG,EAAI,YAAeC,EAAI,YAAeN,IAAM,EACjDG,GAAKE,EAAI,YAAeC,EAAI,YAAeN,IAAM,EACjD,GAAKK,EAAI,YAAeC,EAAI,YAAcN,IAAM,EACtD,MAAO,CACL,aAAeE,EAAI,aAAeC,EAAI,YAAe,EACrD,cAAgBD,EAAI,aAAeC,EAAI,YAAe,EACtD,aAAgBD,EAAI,YAAeC,EAAI,YAAc,CACvD,EAAE,IAAKT,GAAMC,GAAQ,KAAK,IAAI,EAAG,KAAK,IAAI,EAAGD,CAAC,CAAC,CAAC,EAAI,GAAG,CACzD,CACA,SAASa,GAAUhB,EAAQiB,EAAG,CAC5B,IAAMC,EAAI,KAAK,IAAI,EAAG,KAAK,IAAI,EAAGD,CAAC,CAAC,EAC9BlB,EAASf,GAAa,KAAMM,GAAMA,EAAE,QAAUU,CAAM,EAC1D,GAAI,CAACD,GAAU,EAAE,SAAUA,GACzB,OAAOoB,GACLC,GAAUF,CAAC,EACXG,EACAC,EACF,EAAE,IAAKnB,GAAMA,EAAI,GAAG,EAEtB,IAAMoB,EAAOlB,GAAMN,EAAO,KAAK,CAAC,CAAC,EAC3ByB,EAAKnB,GAAMN,EAAO,KAAK,CAAC,CAAC,EAC/B,OAAOc,GAAWU,EAAK,IAAI,CAACpC,EAAGuB,IAAMvB,GAAKqC,EAAGd,CAAC,EAAIvB,GAAK+B,CAAC,CAAC,CAC3D,CACA,SAASO,GAAUzB,EAAQiB,EAAG,CAC5B,GAAM,CAACV,EAAGC,EAAGC,CAAC,EAAIO,GAAUhB,EAAQiB,CAAC,EAAE,IAAI,KAAK,KAAK,EACrD,MAAO,OAAOV,CAAC,KAAKC,CAAC,KAAKC,CAAC,GAC7B,CACA,SAASiB,GAAY1B,EAAQiB,EAAG,CAC9B,GAAM,CAACV,EAAGC,EAAGC,CAAC,EAAIO,GAAUhB,EAAQiB,CAAC,EAAE,IAAI,KAAK,KAAK,EACrD,OAAOU,EAASpB,EAAGC,EAAGC,EAAG,GAAG,CAC9B,CACA,SAASmB,GAAa,CAAE,MAAA9B,EAAO,OAAAE,CAAO,EAAGiB,EAAG,CAC1C,OAAOnB,IAAU,OAAS4B,GAAY1B,EAAQC,EAAM,EAAIgB,IAAM,OAASY,GAAiBH,GAAY1B,EAAQiB,CAAC,CAC/G,CACA,IAAIY,GAAiBC,GACrB,SAASC,GAAc,CAAE,MAAAjC,EAAO,OAAAE,CAAO,EAAG,CACxC,IAAMgC,EAAO,CAACC,EAAQhB,KAAO,CAC3B,OAAAgB,EACA,MAAOR,GAAUzB,EAAQiB,CAAC,EAC1B,KAAMS,GAAY1B,EAAQiB,CAAC,CAC7B,GACA,OAAOnB,IAAU,OAAS,CAACkC,EAAK,EAAG/B,EAAM,CAAC,EAAI,MAAM,KAAK,CAAE,OAAQ,CAAE,EAAG,CAACiC,EAAGxB,IAAMsB,EAAKtB,EAAI,EAAGA,EAAI,CAAC,CAAC,CACtG",
|
|
6
|
-
"names": ["packAbgr", "r", "g", "b", "a", "abgrRed", "c", "abgrGreen", "abgrBlue", "abgrAlpha", "abgrToCssRgba", "normalizedRgbToCssRgba", "alpha", "fadeAbgr", "brightenAbgr", "factor", "t", "lift", "v", "REFERENCE_RAMP_MAX_HUE", "REFERENCE_RAMP_SATURATION", "REFERENCE_RAMP_LIGHTNESS", "rampHueCss", "hue", "rampHueAt", "frac", "rampHue", "bp", "domain", "HUE_STOP_DEG", "rampStops", "start", "reversed", "low", "high", "n", "_", "RAMP_GRADIENT_STOPS", "i", "RAMP_GRADIENT_CSS", "s", "cache", "weakMemo", "baseEdgeCurves", "nodePositions", "graph", "axis", "deletions", "version", "routes", "stranded", "curves", "ei", "bypassed", "edgeCurves", "edgeCurvesOf", "MAX_BACKBONE_HOPS", "buildNeighbors", "graph", "neighbors", "edge", "pushTo", "flankingSpan", "nodeId", "nodeById", "seen", "frontier", "minEnd", "maxStart", "hop", "next", "id", "adjacent", "neighborId", "neighbor", "isBackbone", "nodeReferenceSpan", "node", "referenceMidpoints", "midpoints", "span", "EDGE_DEFAULT_COLOR", "packAbgr", "EDGE_PATH_FALLBACK_COLOR", "EDGE_DELETION_COLOR", "DELETION_THICKNESS_FACTOR", "LIFT_BACKDROP_COLOR", "FADED_ALPHA", "LIFT_BACKDROP_CSS", "abgrToCssRgba", "fadeAbgr", "DELETION_DASH_PX", "MIN_PATH_STRIPE_PX", "MIN_WALK_LANE_PX", "MIN_SCALE_FOR_OFFSET", "ARROW_LENGTH_BASE_PX", "ARROW_LENGTH_PER_EDGE_PX", "ARROW_HALF_WIDTH_BASE_PX", "ARROW_HALF_WIDTH_PER_EDGE_PX", "MIN_LINK_HEADS", "MIN_ARROW_FIT", "MIN_ARROW_SCALE", "ARROW_CURVE_STEPS", "ABUTTING_PX", "walkX", "walkY", "walkPx", "backFromEnd", "curves", "scale", "yToX", "untilPx", "most", "last", "n", "px", "ci", "c", "k", "t", "u", "a", "b", "d", "e", "x", "y", "walkPointAt", "i", "from", "to", "f", "arrowheadFor", "edgeHalfWidth", "nodeHalfWidth", "color", "fromSegments", "toSegments", "length", "lengthPx", "abutting", "inset", "fit", "polylinePx", "tip", "base", "dx", "dy", "endTangent", "intoDirection", "leavesStart", "edge", "segments", "first", "stranded", "strandSides", "gap", "p", "final", "normal", "len", "pointNormalsOf", "rawPoints", "points", "before", "after", "nx", "ny", "dot", "offsetPolyline", "normals", "distance", "hslToRgb", "h", "s", "l", "m", "r", "g", "packHsl", "DEPTH_GRADIENT", "STABLE_RANK_GRADIENT", "NODE_LENGTH_GRADIENT", "sampleGradient", "stops", "pos", "hashColor", "str", "alpha", "nameHue", "PATH_SATURATION", "PATH_LIGHTNESS", "REFERENCE_RAMP_ALT_RGB", "REFERENCE_RAMP_ALT_COLOR", "REFERENCE_RAMP_ALT_CSS", "computeReferenceRamp", "graph", "domain", "midpoints", "referenceMidpoints", "min", "max", "mid", "computeColorSchemeRange", "referenceRamp", "minDepth", "maxDepth", "minLength", "maxLength", "maxRank", "getNodeColor", "node", "nodeIndex", "colorScheme", "range", "ramp", "isBackbone", "rampHueAt", "REFERENCE_RAMP_SATURATION", "REFERENCE_RAMP_LIGHTNESS", "hue", "into", "candidates", "isPolylineInBounds", "bounds", "hit", "isBezierInBounds", "cMinX", "cMaxX", "cMinY", "cMaxY", "paletteCache", "pathPaletteOf", "palette", "pathColorByIndex", "pathIndexByName", "paths", "pathHueAt", "slotsCache", "NO_SLOTS", "pathSlotsOf", "slots", "pathIdx", "path", "nodeId", "pushTo", "buildGeometry", "options", "nodePositions", "nodeById", "contigThickness", "connectorThickness", "drawPaths", "nodeWidth", "highlight", "axis", "linearLayout", "viewportBounds", "deletions", "deletionRoutes", "hiddenEdges", "version", "depthNorm", "meanDepth", "worldPerScreenPx", "yToXOf", "sharedCurves", "baseEdgeCurves", "nodeStrokes", "nodeStrokeRuns", "arrows", "arrowRuns", "edgeCurves", "edgeCurveRuns", "colorRange", "nodeIndexMap", "pathCount", "pathColorByName", "name", "nodePathSlots", "showArrows", "ei", "baseCurves", "numPaths", "isDeletion", "onWalk", "edgeThickness", "ribbons", "plainEdgeColor", "edgeColor", "toNode", "intoHalfWidth", "nodeWidthPx", "inView", "stroke", "thickness", "dash", "dashCurves", "arrowOn", "head", "edgeCurveStart", "arrowStart", "laneWalks", "fromNode", "width", "laneWidth", "walk", "lane", "offset", "translateCurves", "offsets", "pathRibbonOffsets", "arrowRibbon", "ribbonColor", "edgeCurveCount", "arrowCount", "faded", "own", "nodeThickness", "start", "slotWidth", "drawable", "lanes", "laneColor", "slot", "stripe", "count", "WALK_FIELDS", "WALK_SCHEMES", "FAMILIES", "isField", "v", "f", "isScheme", "s", "resolveEncoding", "layer", "reference", "picked", "alone", "askedField", "askedScheme", "field", "family", "scheme", "BASE_T", "toLinear", "c", "toGamma", "oklab", "hex", "r", "g", "b", "i", "l", "m", "oklabToRgb", "L", "a", "schemeRgb", "t", "x", "hslToRgb", "rampHueAt", "REFERENCE_RAMP_SATURATION", "REFERENCE_RAMP_LIGHTNESS", "from", "to", "schemeCss", "schemeColor", "packAbgr", "encodedColor", "NO_VALUE_COLOR", "REFERENCE_RAMP_ALT_COLOR", "encodingStops", "stop", "offset", "_"]
|
|
7
|
-
}
|
|
@@ -1,2 +0,0 @@
|
|
|
1
|
-
import{c as Lt,o as Wt,r as dt}from"./chunk-6MUBJVBG.js";import{a as At,d as Ct,e as Gt}from"./chunk-7HWHWPRW.js";var se="ws",ne="we",Dt=["elsewhere","break","end"],de={missing:"not in the graph",elsewhere:"on nodes filed elsewhere",unknown:"gap, cause unknown"},pe="outside the extracted subgraph";function Zt(t,e,s){if(!e&&!s)return[{bp:t}];let n=Math.min(e?.trimmed??0,t),r=Math.min(s?.trimmed??0,t-n),c=t-n-r,a=e?.kind==="break"?Math.min(e.bp,c):0;c-=a;let u=s?.kind==="break"?Math.min(s.bp,c):0;c-=u;let i=[e,s].find(b=>b?.kind==="elsewhere");return[{bp:n},{bp:a,cause:"missing"},i?.kind==="elsewhere"?{bp:c,cause:"elsewhere",chunk:i.chunk}:{bp:c,cause:"unknown"},{bp:u,cause:"missing"},{bp:r}].filter(b=>b.bp>0)}function Et(t,e){let s=t.bp[e],n=t.chunkRefs[e];return Ht(t.kindNames[t.kinds[e]],Number.isNaN(s)?void 0:s,t.trimmed[e],n<0?void 0:{refName:t.chunkRefNames[n],start:t.chunkStarts[e]})}function Ht(t,e,s,n){let r=s>0?{trimmed:s}:{};return t==="elsewhere"&&n?{kind:t,chunk:n,...e===void 0?{}:{bp:e},...r}:t==="break"&&e!==void 0?{kind:t,bp:e,...n?{chunk:n}:{},...r}:t==="end"?{kind:t,...r}:void 0}function Vt(t){let[e,...s]=t.split(","),n,r=0,c;for(let a of s)if(a.startsWith("bp="))n=+a.slice(3);else if(a.startsWith("trimmed="))r=+a.slice(8);else if(a.startsWith("chunk=")){let u=a.lastIndexOf(":");c={refName:a.slice(6,u),start:+a.slice(u+1)}}return Ht(e,n,r,c)}function Ot(t,e){return{...t?{before:t}:{},...e?{after:e}:{}}}function ge(t){let e=t[se],s=t[ne];return Ot(typeof e=="string"?Vt(e):void 0,typeof s=="string"?Vt(s):void 0)}var ve=class{constructor(){this.kinds=[],this.bp=[],this.trimmed=[],this.chunkRefs=[],this.chunkStarts=[],this.chunkRefNames=[],this.refIndex=new Map}push(t){this.kinds.push(Dt.indexOf(t.kind)),this.bp.push("bp"in t&&t.bp!==void 0?t.bp:NaN),this.trimmed.push(t.trimmed??0);let e="chunk"in t?t.chunk:void 0;if(e){let s=this.refIndex.get(e.refName);s===void 0&&(s=this.chunkRefNames.length,this.refIndex.set(e.refName,s),this.chunkRefNames.push(e.refName)),this.chunkRefs.push(s),this.chunkStarts.push(e.start)}else this.chunkRefs.push(-1),this.chunkStarts.push(NaN)}table(){return{kinds:Uint8Array.from(this.kinds),kindNames:[...Dt],bp:Float64Array.from(this.bp),trimmed:Float64Array.from(this.trimmed),chunkRefs:Int32Array.from(this.chunkRefs),chunkStarts:Float64Array.from(this.chunkStarts),chunkRefNames:this.chunkRefNames}}};function re(t,e,s){return e===0?t:`${t}:${e}-${s}`}function xt(t,e){let s=t.names[e],n=Gt(s),r=t.starts[e];return n?{name:`${n.sample}#${n.haplotype}#${n.contig}`,start:r,...n,...Ot(t.before&&Et(t.before,e),t.after&&Et(t.after,e))}:{name:re(s,r,t.ends[e]),start:r}}function oe(t,e,s,n){let r=new Int32Array(s+1);for(let i of t)r[i+1]++;for(let i=0;i<s;i++)r[i+1]+=r[i];let c=r.slice(0,s),a=new Int32Array(t.length);t.forEach((i,d)=>{a[c[i]++]=e[d]});let u=new Int32Array(n.length).fill(-1);return Array.from({length:s},(i,d)=>{let b=[];for(let I=r[d];I<r[d+1];I++){let U=a[I];u[U]!==d&&(u[U]=d,b.push(n[U]))}return b})}function Ft({nodes:t,links:e,walks:s}){let n=new Uint8Array(t.names.length),r=(a,u)=>{n[a]===0&&(n[a]=u?2:1)};for(let a=0;a<e.from.length;a++)r(e.from[a],e.strands[a]&1),r(e.to[a],e.strands[a]&2);let c=s.offsets[s.names.length]??0;for(let a=0;a<c;a++)r(s.steps[a],s.reversed[a]);return t.names.map((a,u)=>`${a}${n[u]===2?"-":"+"}`)}function Se({nodes:t,links:e,walks:s},n="Imported GFA"){let r=t.names.length,c=t.lengths.length,a=[];for(let o=0;o<s.names.length;o++)s.offsets[o+1]>s.offsets[o]&&a.push(o);let u=new Int32Array(r);for(let o of s.steps.subarray(0,s.offsets[s.names.length]))u[o]++;let i=Ft({nodes:t,links:e,walks:s}),d=[];for(let o=0;o<c;o++){let g=t.refs[o];d.push({id:i[o],name:t.names[o],length:t.lengths[o],depth:Math.max(u[o],1),stable:g<0?void 0:{refName:t.refNames[g],start:t.starts[o],rank:t.ranks[o]},samples:void 0})}let b=[],I=new Array(r),U=(o,g)=>{let h=I[o];if(h){for(let y=0;y<h.length;y+=2)if(h[y]===g)return h[y+1]}},B=(o,g,h)=>{let y=I[o];y?y.push(g,h):I[o]=[g,h]},st=[],X=0;for(let o=0;o<e.from.length;o++){let g=e.from[o],h=e.to[o];b.push({from:i[g],to:i[h],fromStrand:e.strands[o]&1?"-":"+",toStrand:e.strands[o]&2?"-":"+"});let y=U(g,h);y===void 0&&(y=X++,B(g,h,y),B(h,g,y)),st.push(y)}let F=new Int32Array(r+1);for(let o=0;o<r;o++)F[o+1]=F[o]+(I[o]?.length??0)/2;let nt=new Int32Array(F[r]),rt=new Int32Array(F[r]);I.forEach((o,g)=>{for(let h=0;h<o.length;h+=2)nt[F[g]+h/2]=o[h],rt[F[g]+h/2]=o[h+1]});let St=(o,g)=>{for(let h=F[o];h<F[o+1];h++)if(nt[h]===g)return rt[h];return-1},O=new Map,M=[],K=new Int32Array(s.offsets[s.names.length]??0),ot=new Int32Array(K.length),$=0,gt=new Int32Array(X).fill(-1),_=[],q=[],z=s.offsets[s.names.length]??0,vt=new Int32Array(z),G=new Float64Array(z),R=new Uint8Array(z),x=new Int32Array(r).fill(-1),at=[],L=[0],W=[],it=new Map,ft=[];for(let o of a){let g=xt(s,o),{start:h,name:y}=g,f=g.sample===void 0?s.names[o]:y,m=O.get(f);m===void 0&&(m=M.length,O.set(f,m),M.push(f),W.push(At(f)));let S=it.get(y);S===void 0&&(S=ft.length,it.set(y,S),ft.push(y));let P=W[m],k=s.offsets[o],lt=new Array(s.offsets[o+1]-k),T=h,D=-1;for(let C=k;C<s.offsets[o+1];C++){let E=s.steps[C];lt[C-k]=i[E],x[E]===-1&&(x[E]=L.at(-1),L.push(x[E]+u[E]),at.push(E));let V=x[E]++;if(vt[V]=m,G[V]=T,R[V]=s.reversed[C],T+=E<c?t.lengths[E]:0,D>=0){let J=St(D,E);J>=0&>[J]!==o&&(gt[J]=o,K[$]=J,ot[$++]=S)}D=E}_.push({...g,nodeIds:lt}),q.push({name:f,sample:P,start:h,length:T-h})}let It=oe(K.subarray(0,$),ot,X,ft);b.forEach((o,g)=>{let h=It[st[g]];h.length>0&&(o.pathIds=h)});let ct=a.length>0;return{name:n,nodes:d,edges:b,paths:ct?_:void 0,anchorPaths:ct?q:void 0,pathVisits:ct?new Wt({segments:at.map(o=>t.names[o]),offsets:Int32Array.from(L),path:vt.subarray(0,L.at(-1)),start:G.subarray(0,L.at(-1)),reversed:R.subarray(0,L.at(-1)),paths:M,samples:W}):void 0,anchoredBy:d.some(o=>o.stable)?"tags":void 0}}var Tt=1,Kt=2,wt=4,qt=[void 0,"missing","elsewhere","unknown"],zt=3,ae=t=>qt.indexOf(t)<<zt;function Jt(t,e,s){return e!==void 0&&t.includes(e)?e:t.find(n=>s.has(n))??t[0]}function ie(t){return Jt((t.paths??[]).map(e=>dt(e.name).name),t.referencePath,new Set(t.nodes.flatMap(e=>e.stable?.rank===0?[e.stable.refName]:[])))}function fe(t){let e=new Map,s=[],n=[],r=d=>{let b=e.get(d);return b===void 0&&(b=s.length,e.set(d,b),s.push(d),n.push(0)),b};for(let d of t.nodes)n[r(d.id)]=d.length;let c=t.paths??[],a=new Int32Array(c.length),u=new Int32Array(c.length),i=[];return c.forEach((d,b)=>{a[b]=i.length;for(let I of d.nodeIds)i.push(r(I));u[b]=i.length}),{heads:c,from:a,to:u,steps:i,lengths:n,names:s,reference:ie(t)}}function ce(t,e,s){let{nodes:n,walks:r}=t,c=[];for(let i=0;i<r.names.length;i++)r.offsets[i+1]>r.offsets[i]&&c.push(i);let a=c.map(i=>xt(r,i)),u=new Set;for(let i=0;i<n.lengths.length;i++)n.refs[i]>=0&&n.ranks[i]===0&&u.add(n.refNames[n.refs[i]]);return{heads:a,from:c.map(i=>r.offsets[i]),to:c.map(i=>r.offsets[i+1]),steps:r.steps,lengths:n.lengths,names:s?Ft(t):n.names,reference:Jt(a.map(i=>dt(i.name).name),e,u)}}function Qt(t){return t.sample??At(t.name)}function le(t){return t.haplotype!==void 0&&t.sample!==void 0?`${t.sample}#${t.haplotype}`:Qt(t)}var Xt=class{constructor(t){this.offsets=[0],this.values=Array.from({length:t},()=>[])}get length(){return this.values[0].length}push(...t){t.forEach((e,s)=>this.values[s].push(e))}endRow(){this.offsets.push(this.length)}reorder(t,e){let s=e.map(c=>c(this.length)),n=new Uint32Array(t.length+1),r=0;return t.forEach((c,a)=>{for(let u=this.offsets[c];u<this.offsets[c+1];u++,r++)s.forEach((i,d)=>{i[r]=this.values[d][u]});n[a+1]=r}),{offsets:n,columns:s}}},pt=t=>new Float64Array(t);function Fe(t,e,s={}){let{heads:n,from:r,to:c,steps:a,lengths:u,names:i,reference:d}="walks"in t?ce(t,s.referencePath,s.visits):fe(t),b=u.length,I=f=>f<b?u[f]:0,U=f=>dt(n[f].name).name,B=f=>n[f].start,st=(f,m)=>(B(f)??0)-(B(m)??0),X=n.map((f,m)=>m),F=X.filter(f=>U(f)===d).sort(st),nt=F[0],rt=X.filter(f=>U(f)!==d);if(nt===void 0||rt.length===0)return;let St=Lt(rt,U),O=new Float64Array(i.length).fill(Number.NaN),M=new Float64Array(i.length),K=!1,ot=!1,$=e&&e.end>e.start?e:void 0;for(let f of F){let m=B(f)??0;for(let S=r[f];S<c[f];S++){let P=a[S],k=I(P);Number.isNaN(O[P])&&(O[P]=m,M[P]=m+k,$&&(K||=m+k<=$.start,ot||=m>=$.end)),m+=k}}let gt=$===void 0||K&&ot,_=new Xt(5),q=[],z=new Map,vt=f=>{if(!f)return-1;let m=`${f.refName}:${f.start}`,S=z.get(m);return S===void 0&&(S=q.length,z.set(m,S),q.push(f)),S},G=s.visits?new Xt(4):void 0,R=new Int32Array(0),x=new Float64Array(0),at=f=>{let m=n[f[0]],S=[...f].sort(st),P=S.length-1;for(let l of S)P+=c[l]-r[l];P>R.length&&(R=new Int32Array(P*2),x=new Float64Array(P*2));let k=0,lt=!1,T,D,C=new Map;for(let l of S){let p=B(l);(p===void 0||T!==void 0&&p<T)&&(lt=!0),T!==void 0&&p!==void 0&&p>T&&(C.set(k,Zt(p-T,D===void 0?void 0:n[D].after,n[l].before)),R[k]=-(p-T),x[k++]=T);let v=p??0;for(let w=r[l];w<c[l];w++){let N=a[w];R[k]=N,x[k++]=v,v+=I(N)}T=p===void 0?void 0:v,D=l}let E=l=>l>=0&&l<k&&R[l]>=0&&!Number.isNaN(O[R[l]]),V=l=>{let p=R[l];return x[l]+(p<0?-p:I(p))},J=(l,p)=>{let v=l>=0?l:p,w=-1;for(let H=1;w<0&&H<k;H++)w=[v-H,v+H].find(E)??-1;if(w<0)return;let N=w>v==O[R[w]]>O[R[v]],A=N===(v===l),tt=A?k-1:0,et=E(tt);return{stop:{contigEnds:et,shortBp:et?Math.max(0,$.start-M[R[tt]],O[R[tt]]-$.end):0},lo:A?v+1:0,hi:A?k:v,backward:!N,axisStart:N?A?V(v):x[0]:A?V(k-1):x[v]}},Nt=0,kt=k,ht=!1,Pt=!0,ut=-1,Mt=-1,Rt,_t;if($&>){let l=-1,p=-1,v=-1/0,w=1/0;for(let N=0;N<k;N++)if(E(N)){let A=R[N];M[A]<=$.start&&M[A]>v&&(v=M[A],l=N),O[A]>=$.end&&O[A]<w&&(w=O[A],p=N)}if(l<0||p<0){Pt=!1;let N=l>=0||p>=0?J(l,p):void 0;N&&({lo:Nt,hi:kt,backward:ht,axisStart:Rt,stop:_t}=N)}else Nt=Math.min(l,p)+1,kt=Math.max(l,p),ht=l>p,ut=l,Mt=p}let bt=m.contig??Ct(dt(m.name).name),Yt=lt||k===0?void 0:Rt!==void 0?{contig:bt,start:Rt,reversed:ht}:ut<0?{contig:bt,start:x[0],reversed:!1}:ut<Mt?{contig:bt,start:V(ut),reversed:!1}:{contig:bt,start:x[ut],reversed:!0},Q=0,Ut=0,jt=0,mt=0,j=-1,Y=0,Z=Number.NaN,[,$t,te,ee]=_.values;for(let l=Nt;l<kt;l++){let p=ht?kt-1-(l-Nt):l,v=R[p];if(v<0){let H=C.get(p);for(let yt of ht?[...H].reverse():H)_.push(Q,yt.bp,Number.NaN,wt|ae(yt.cause),vt(yt.chunk)),Q+=yt.bp;j=_.length-1,Y=wt,Z=Number.NaN,jt+=-v,mt=0;continue}let w=I(v),N=O[v],A=!Number.isNaN(N);G?.push(v,Q,w,A?1:0);let tt=j>=0&&!Number.isNaN(Z),et=A&&tt&&mt>=0&&Z+($t[j]??0)===N,Bt=A&&tt&&mt<=0&&M[v]===Z;!A&&j>=0&&!(Y&(Tt|wt))?$t[j]+=w:et||Bt?($t[j]+=w,et||(Z=N,te[j]=N,Y|=Kt,ee[j]=Y),mt=et?1:-1):(Y=A?Tt:0,Z=A?N:Number.NaN,_.push(Q,w,Z,Y,-1),j=_.length-1,mt=0),Q+=w,A||(Ut+=w)}return _.endRow(),G?.endRow(),{name:m.name,label:le(m),sample:Qt(m),haplotype:m.haplotype,bp:Q,offReferenceBp:Ut,gapBp:jt,complete:Pt,stop:_t,axis:Yt}},L=at(F),W=[...St.values()].map((f,m)=>({meta:at(f),row:m+1}));W.sort((f,m)=>Number(m.meta.complete)-Number(f.meta.complete)||m.meta.bp-f.meta.bp||f.meta.label.localeCompare(m.meta.label));let it=[0,...W.map(f=>f.row)],{offsets:ft,columns:It}=_.reorder(it,[pt,pt,pt,f=>new Uint8Array(f),f=>new Int32Array(f)]),[ct,o,g,h,y]=It;return{origin:$?$.start:B(nt)??0,rows:[L,...W.map(f=>f.meta)],offsets:ft,start:ct,bp:o,referenceStart:g,flags:h,chunk:y,chunks:q,visits:G&&he(G,it,i)}}function he(t,e,s){let{offsets:n,columns:r}=t.reorder(e,[d=>new Int32Array(d),pt,pt,d=>new Uint8Array(d)]),[c,a,u,i]=r;return{offsets:n,node:c,start:a,bp:u,onReference:i,names:s}}function ue(t,e){let s=t.flags[e],n={start:t.start[e],bp:t.bp[e],onReference:(s&Tt)!==0};if(s&wt){n.gap=!0;let r=qt[s>>zt&3];r&&(n.cause=r);let c=t.chunks[t.chunk[e]];c&&(n.chunk=c)}else n.referenceStart=Number.isNaN(t.referenceStart[e])?void 0:t.referenceStart[e];return s&Kt&&(n.reversed=!0),n}function me(t,e){let s=[];for(let n=t.offsets[e];n<t.offsets[e+1];n++)s.push({node:t.names[t.node[n]],start:t.start[n],bp:t.bp[n],onReference:t.onReference[n]===1});return s}function Te(t,e){let s=t.rows.map((n,r)=>{let c=[];for(let a=t.offsets[r];a<t.offsets[r+1];a++)c.push(ue(t,a));return{...n,runs:c,visits:t.visits&&me(t.visits,r)}});return{origin:t.origin,unit:e,reference:s[0],rows:s.slice(1)}}export{se as a,ne as b,de as c,pe as d,Zt as e,Vt as f,ge as g,ve as h,xt as i,Se as j,ie as k,Fe as l,Te as m};
|
|
2
|
-
//# sourceMappingURL=chunk-6KUHGRHU.js.map
|