@genome-spy/app 1.2.0 → 1.3.0

package/dist/esm-DfMIJ_N2.js

@@ -0,0 +1,1522 @@
+import { n as e, r as t } from "./esm-C6erfjPM.js";
+import { t as n } from "./esm-C3A-AuRi.js";
+import { n as r, r as i } from "./esm-BUkWqSr2.js";
+//#region ../../node_modules/@gmod/bam/esm/chunk.js
+var a = class {
+	minv;
+	maxv;
+	bin;
+	endPosition;
+	constructor(e, t, n, r = t.blockPosition + 65536) {
+		this.minv = e, this.maxv = t, this.bin = n, this.endPosition = r;
+	}
+	toString() {
+		return `${this.minv.toString()}..${this.maxv.toString()} (bin ${this.bin}, fetchedSize ${this.fetchedSize()})`;
+	}
+	fetchedSize() {
+		return this.endPosition - this.minv.blockPosition;
+	}
+}, o = class {
+	blockPosition;
+	dataPosition;
+	constructor(e, t) {
+		this.blockPosition = e, this.dataPosition = t;
+	}
+	toString() {
+		return `${this.blockPosition}:${this.dataPosition}`;
+	}
+};
+function s(e, t = 0) {
+	return new o(e[t + 7] * 1099511627776 + e[t + 6] * 4294967296 + e[t + 5] * 16777216 + e[t + 4] * 65536 + e[t + 3] * 256 + e[t + 2], e[t + 1] << 8 | e[t]);
+}
+//#endregion
+//#region ../../node_modules/@gmod/bam/esm/util.js
+function c(e) {
+	if (e?.aborted) {
+		let { reason: t } = e;
+		throw t === void 0 ? new DOMException("This operation was aborted", "AbortError") : t;
+	}
+}
+function l(e, t) {
+	let n = e.length;
+	if (n === 0) return e;
+	let r;
+	if (t) {
+		let i = t.blockPosition, a = t.dataPosition;
+		r = [];
+		for (let t = 0; t < n; t++) {
+			let n = e[t], o = n.maxv;
+			(o.blockPosition - i || o.dataPosition - a) > 0 && r.push(n);
+		}
+		if (r.length === 0) return r;
+	} else r = e;
+	r.sort((e, t) => {
+		let n = e.minv.blockPosition - t.minv.blockPosition;
+		return n === 0 ? e.minv.dataPosition - t.minv.dataPosition : n;
+	});
+	let i = [r[0]], o = r[0].minv.blockPosition, s = r[0].maxv.blockPosition;
+	for (let e = 1; e < r.length; e++) {
+		let t = r[e], n = t.minv.blockPosition, c = t.maxv.blockPosition;
+		if (n - s < 65e3 && c - o < 5e6) {
+			let e = i[i.length - 1];
+			(c - s || t.maxv.dataPosition - e.maxv.dataPosition) > 0 && (i[i.length - 1] = new a(e.minv, t.maxv, e.bin, t.endPosition), s = c);
+		} else i.push(t), o = n, s = c;
+	}
+	return u(i);
+}
+function u(e) {
+	let t = [e[0]];
+	for (let n = 1; n < e.length; n++) {
+		let r = e[n], i = t[t.length - 1].maxv;
+		(r.minv.blockPosition - i.blockPosition || r.minv.dataPosition - i.dataPosition) >= 0 ? t.push(r) : (r.maxv.blockPosition - i.blockPosition || r.maxv.dataPosition - i.dataPosition) > 0 && t.push(new a(i, r.maxv, r.bin, r.endPosition));
+	}
+	return t;
+}
+function d(e, t) {
+	let n = e[t] | e[t + 1] << 8 | e[t + 2] << 16 | e[t + 3] << 24;
+	return { lineCount: ((e[t + 4] | e[t + 5] << 8 | e[t + 6] << 16 | e[t + 7] << 24) >>> 0) * 2 ** 32 + (n >>> 0) };
+}
+function f(e, t = []) {
+	if (e.length === 0) return;
+	let n = Array(t.length + e.length * 2), r = 0;
+	for (let e = 0, i = t.length; e < i; e++) n[r++] = t[e];
+	for (let t of e) n[r++] = t.minv.blockPosition, n[r++] = t.maxv.blockPosition;
+	n.sort((e, t) => e - t);
+	for (let t of e) {
+		let e = t.maxv.blockPosition, r = 0, i = n.length;
+		for (; r < i;) {
+			let t = r + i >>> 1;
+			n[t] > e ? i = t : r = t + 1;
+		}
+		r < n.length && (t.endPosition = Math.min(t.endPosition, n[r]));
+	}
+}
+function p(e, t, n) {
+	if (t + 4 > e.length) return;
+	let r = new DataView(e.buffer, e.byteOffset, e.byteLength), i = r.getInt32(t, !0), a = Object.create(null), o = [], s = new TextDecoder("utf8"), c = t + 4;
+	for (let t = 0; t < i; t++) {
+		if (c + 8 > e.length) return;
+		let i = r.getInt32(c, !0);
+		if (c + 8 + i > e.length) return;
+		let l = n(s.decode(e.subarray(c + 4, c + 4 + i - 1))), u = r.getInt32(c + i + 4, !0);
+		a[l] = t, o.push({
+			refName: l,
+			length: u
+		}), c += 8 + i;
+	}
+	return {
+		chrToIndex: a,
+		indexToChr: o,
+		end: c
+	};
+}
+function m(e, t, n, r) {
+	let i = r ? r.blockPosition : Infinity, a = r ? r.dataPosition : 0, s = !1;
+	for (let r = 0; r < n; r++) {
+		let n = t + r * 8, o = e[n + 7] * 1099511627776 + e[n + 6] * 4294967296 + e[n + 5] * 16777216 + e[n + 4] * 65536 + e[n + 3] * 256 + e[n + 2], c = e[n + 1] << 8 | e[n];
+		(o !== 0 || c !== 0) && (o < i || o === i && c < a) && (i = o, a = c, s = !0);
+	}
+	return s ? new o(i, a) : r;
+}
+function h(e, t = (e) => e) {
+	let n = new TextDecoder(), r = 0, i = 0, a = [], o = Object.create(null);
+	for (let s = 0; s < e.length; s++) if (!e[s]) {
+		if (i < s) {
+			let c = t(n.decode(e.subarray(i, s)));
+			a[r] = c, o[c] = r;
+		}
+		i = s + 1, r++;
+	}
+	return {
+		refNameToId: o,
+		refIdToName: a
+	};
+}
+function g(e) {
+	let t = 0;
+	for (let n of e) t += n.length;
+	let n = new Uint8Array(t), r = 0;
+	for (let t of e) n.set(t, r), r += t.length;
+	return n;
+}
+function _(e) {
+	return e.end > e.start ? e.end : e.start + 1;
+}
+function v(e, t, n, r, i = []) {
+	for (let a = 0, o = e.length; a < o; a++) {
+		let o = e[a];
+		if (o.ref_id === t) {
+			if (o.start >= r) break;
+			_(o) > n && i.push(o);
+		} else if (o.ref_id > t) break;
+	}
+	return i;
+}
+//#endregion
+//#region ../../node_modules/@gmod/bam/esm/indexFile.js
+function y(e, t = 5) {
+	let r = new n({ maxSize: t });
+	return (t) => {
+		let n = r.get(t);
+		if (n !== void 0) return n;
+		let i = e(t);
+		return i && r.set(t, i), i;
+	};
+}
+var b = class {
+	filehandle;
+	renameRefSeq;
+	setupP;
+	setupSignal;
+	constructor({ filehandle: e, renameRefSeq: t = (e) => e }) {
+		this.filehandle = e, this.renameRefSeq = t;
+	}
+	async blocksForRange(e, t, n, r = {}) {
+		t < 0 && (t = 0);
+		let i = (await this.parse(r)).indices(e);
+		if (!i) return [];
+		let a = this.reg2bins(t, n);
+		if (a.length === 0) return [];
+		let o = [], { binIndex: s } = i;
+		for (let [e, t] of a) for (let n = e; n <= t; n++) {
+			let e = s[n];
+			if (e) for (let t = 0, n = e.length; t < n; t++) o.push(e[t]);
+		}
+		return l(o, this.getLowestChunk(i, t));
+	}
+	async parse(e = {}, t = !1) {
+		c(e.signal);
+		let n = this.setupP;
+		if (!n) return this.startParse(e);
+		let r = this.setupSignal;
+		try {
+			return await n;
+		} catch (n) {
+			if (t || !r?.aborted || e.signal?.aborted) throw n;
+			return this.parse(e, !0);
+		}
+	}
+	startParse(e) {
+		let t = this._parse(e);
+		return this.setupP = t, this.setupSignal = e.signal, t.then(() => {
+			this.setupP === t && (this.setupSignal = void 0);
+		}, () => {
+			this.setupP === t && (this.setupP = void 0, this.setupSignal = void 0);
+		}), t;
+	}
+	async lineCount(e, t) {
+		return (await this.parse(t)).indices(e)?.stats?.lineCount ?? 0;
+	}
+	async hasRefSeq(e, t) {
+		return !!(await this.parse(t)).indices(e);
+	}
+	async estimatedBytesForRegions(e, t) {
+		let n = l((await Promise.all(e.map((e) => this.blocksForRange(e.refId, e.start, e.end, t)))).flat()), r = 0;
+		for (let e of n) r += e.fetchedSize();
+		return r;
+	}
+}, x = /* @__PURE__ */ new Float64Array(), ee = 21578050, S = 14, te = 16384, C = 2 ** 29;
+function ne(e, t) {
+	return e - e % t;
+}
+function w(e, t) {
+	let n = e % t;
+	return n === 0 ? e : e - n + t;
+}
+function re(e, t) {
+	return e > C && (e = C), t > C && (t = C), --t, [
+		[0, 0],
+		[1 + (e >> 26), 1 + (t >> 26)],
+		[9 + (e >> 23), 9 + (t >> 23)],
+		[73 + (e >> 20), 73 + (t >> 20)],
+		[585 + (e >> 17), 585 + (t >> 17)],
+		[4681 + (e >> S), 4681 + (t >> S)]
+	];
+}
+var T = class extends b {
+	async _parse(e) {
+		let t = await this.filehandle.readFile(e), n = new DataView(t.buffer, t.byteOffset, t.byteLength);
+		if (n.getUint32(0, !0) !== ee) throw Error("Not a BAI file");
+		let r = n.getInt32(4, !0), i = 8, o, c = [];
+		for (let e = 0; e < r; e++) {
+			c.push(i);
+			let e = n.getInt32(i, !0);
+			i += 4;
+			for (let t = 0; t < e; t += 1) {
+				let e = n.getUint32(i, !0);
+				if (i += 4, e === 37450) i += 4, i += 32;
+				else if (e > 37450) throw Error("bai index contains too many bins, please use CSI");
+				else {
+					let e = n.getInt32(i, !0);
+					i += 4 + 16 * e;
+				}
+			}
+			let r = n.getInt32(i, !0);
+			i += 4, o = m(t, i, r, o), i += 8 * r;
+		}
+		function l(e) {
+			let r = c[e];
+			if (r === void 0) return;
+			let i = n.getInt32(r, !0), o;
+			r += 4;
+			let l = {};
+			for (let e = 0; e < i; e += 1) {
+				let e = n.getUint32(r, !0);
+				if (r += 4, e === 37450) r += 4, o = d(t, r + 16), r += 32;
+				else if (e > 37450) throw Error("bai index contains too many bins, please use CSI");
+				else {
+					let i = n.getInt32(r, !0);
+					r += 4;
+					let o = Array(i);
+					for (let n = 0; n < i; n++) {
+						let i = s(t, r);
+						r += 8;
+						let c = s(t, r);
+						r += 8, o[n] = new a(i, c, e);
+					}
+					l[e] = o;
+				}
+			}
+			let u = n.getInt32(r, !0);
+			r += 4;
+			let p = u === 0 ? x : new Float64Array(u), m = u === 0 ? x : new Float64Array(u);
+			for (let e = 0; e < u; e++) p[e] = t[r + 7] * 1099511627776 + t[r + 6] * 4294967296 + t[r + 5] * 16777216 + t[r + 4] * 65536 + t[r + 3] * 256 + t[r + 2], m[e] = t[r + 1] << 8 | t[r], r += 8;
+			return f(Object.values(l).flat(), p), {
+				binIndex: l,
+				linearBlockPositions: p,
+				linearDataPositions: m,
+				stats: o
+			};
+		}
+		return {
+			bai: !0,
+			firstDataLine: o,
+			indices: y(l),
+			refCount: r
+		};
+	}
+	async indexCov(e, t, n, r) {
+		let i = te, a = (await this.parse(r)).indices(e);
+		if (!a) return [];
+		let { linearBlockPositions: o, stats: s } = a, c = o.length;
+		if (c === 0) return [];
+		let l = (c - 1) * i, u = t === void 0 ? 0 : Math.min(Math.max(ne(t, i), 0), l), d = n === void 0 ? l : Math.min(w(n, i), l);
+		if (d <= u) return [];
+		let f = Array((d - u) / i), p = o[c - 1], m = s?.lineCount ?? 0, h = o[u / i];
+		for (let e = u / i, t = 0; e < d / i; e++, t++) {
+			let n = o[e + 1];
+			f[t] = {
+				score: (n - h) * m / p,
+				start: e * i,
+				end: e * i + i
+			}, h = n;
+		}
+		return f;
+	}
+	reg2bins(e, t) {
+		return re(e, t);
+	}
+	getLowestChunk(e, t) {
+		let { linearBlockPositions: n, linearDataPositions: r } = e, i = Math.min(t >> S, n.length - 1);
+		return i < 0 ? void 0 : {
+			blockPosition: n[i],
+			dataPosition: r[i]
+		};
+	}
+}, E = [
+	0,
+	1996959894,
+	3993919788,
+	2567524794,
+	124634137,
+	1886057615,
+	3915621685,
+	2657392035,
+	249268274,
+	2044508324,
+	3772115230,
+	2547177864,
+	162941995,
+	2125561021,
+	3887607047,
+	2428444049,
+	498536548,
+	1789927666,
+	4089016648,
+	2227061214,
+	450548861,
+	1843258603,
+	4107580753,
+	2211677639,
+	325883990,
+	1684777152,
+	4251122042,
+	2321926636,
+	335633487,
+	1661365465,
+	4195302755,
+	2366115317,
+	997073096,
+	1281953886,
+	3579855332,
+	2724688242,
+	1006888145,
+	1258607687,
+	3524101629,
+	2768942443,
+	901097722,
+	1119000684,
+	3686517206,
+	2898065728,
+	853044451,
+	1172266101,
+	3705015759,
+	2882616665,
+	651767980,
+	1373503546,
+	3369554304,
+	3218104598,
+	565507253,
+	1454621731,
+	3485111705,
+	3099436303,
+	671266974,
+	1594198024,
+	3322730930,
+	2970347812,
+	795835527,
+	1483230225,
+	3244367275,
+	3060149565,
+	1994146192,
+	31158534,
+	2563907772,
+	4023717930,
+	1907459465,
+	112637215,
+	2680153253,
+	3904427059,
+	2013776290,
+	251722036,
+	2517215374,
+	3775830040,
+	2137656763,
+	141376813,
+	2439277719,
+	3865271297,
+	1802195444,
+	476864866,
+	2238001368,
+	4066508878,
+	1812370925,
+	453092731,
+	2181625025,
+	4111451223,
+	1706088902,
+	314042704,
+	2344532202,
+	4240017532,
+	1658658271,
+	366619977,
+	2362670323,
+	4224994405,
+	1303535960,
+	984961486,
+	2747007092,
+	3569037538,
+	1256170817,
+	1037604311,
+	2765210733,
+	3554079995,
+	1131014506,
+	879679996,
+	2909243462,
+	3663771856,
+	1141124467,
+	855842277,
+	2852801631,
+	3708648649,
+	1342533948,
+	654459306,
+	3188396048,
+	3373015174,
+	1466479909,
+	544179635,
+	3110523913,
+	3462522015,
+	1591671054,
+	702138776,
+	2966460450,
+	3352799412,
+	1504918807,
+	783551873,
+	3082640443,
+	3233442989,
+	3988292384,
+	2596254646,
+	62317068,
+	1957810842,
+	3939845945,
+	2647816111,
+	81470997,
+	1943803523,
+	3814918930,
+	2489596804,
+	225274430,
+	2053790376,
+	3826175755,
+	2466906013,
+	167816743,
+	2097651377,
+	4027552580,
+	2265490386,
+	503444072,
+	1762050814,
+	4150417245,
+	2154129355,
+	426522225,
+	1852507879,
+	4275313526,
+	2312317920,
+	282753626,
+	1742555852,
+	4189708143,
+	2394877945,
+	397917763,
+	1622183637,
+	3604390888,
+	2714866558,
+	953729732,
+	1340076626,
+	3518719985,
+	2797360999,
+	1068828381,
+	1219638859,
+	3624741850,
+	2936675148,
+	906185462,
+	1090812512,
+	3747672003,
+	2825379669,
+	829329135,
+	1181335161,
+	3412177804,
+	3160834842,
+	628085408,
+	1382605366,
+	3423369109,
+	3138078467,
+	570562233,
+	1426400815,
+	3317316542,
+	2998733608,
+	733239954,
+	1555261956,
+	3268935591,
+	3050360625,
+	752459403,
+	1541320221,
+	2607071920,
+	3965973030,
+	1969922972,
+	40735498,
+	2617837225,
+	3943577151,
+	1913087877,
+	83908371,
+	2512341634,
+	3803740692,
+	2075208622,
+	213261112,
+	2463272603,
+	3855990285,
+	2094854071,
+	198958881,
+	2262029012,
+	4057260610,
+	1759359992,
+	534414190,
+	2176718541,
+	4139329115,
+	1873836001,
+	414664567,
+	2282248934,
+	4279200368,
+	1711684554,
+	285281116,
+	2405801727,
+	4167216745,
+	1634467795,
+	376229701,
+	2685067896,
+	3608007406,
+	1308918612,
+	956543938,
+	2808555105,
+	3495958263,
+	1231636301,
+	1047427035,
+	2932959818,
+	3654703836,
+	1088359270,
+	936918e3,
+	2847714899,
+	3736837829,
+	1202900863,
+	817233897,
+	3183342108,
+	3401237130,
+	1404277552,
+	615818150,
+	3134207493,
+	3453421203,
+	1423857449,
+	601450431,
+	3009837614,
+	3294710456,
+	1567103746,
+	711928724,
+	3020668471,
+	3272380065,
+	1510334235,
+	755167117
+];
+typeof Int32Array < "u" && (E = new Int32Array(E));
+var ie = (e, t) => {
+	let n = t === 0 ? 0 : ~~t ^ -1;
+	for (let t = 0; t < e.length; t++) n = E[(n ^ e[t]) & 255] ^ n >>> 8;
+	return n ^ -1;
+}, ae = 21582659, oe = 38359875, se = new o(0, 0);
+function D(e, t) {
+	return e * 2 ** t;
+}
+function O(e, t) {
+	return Math.floor(e / 2 ** t);
+}
+var k = class extends b {
+	maxBinNumber = 0;
+	depth = 0;
+	minShift = 0;
+	async indexCov() {
+		return [];
+	}
+	parseAuxData(e, t) {
+		let n = new DataView(e.buffer, e.byteOffset, e.byteLength), r = n.getUint32(t, !0), i = r & 65536 ? "zero-based-half-open" : "1-based-closed", a = [
+			"generic",
+			"SAM",
+			"VCF"
+		][r & 15];
+		if (!a) throw Error(`invalid Tabix preset format flags ${r}`);
+		let o = {
+			ref: n.getInt32(t + 4, !0),
+			start: n.getInt32(t + 8, !0),
+			end: n.getInt32(t + 12, !0)
+		}, s = n.getInt32(t + 16, !0), c = s ? String.fromCharCode(s) : "", l = n.getInt32(t + 20, !0), u = n.getInt32(t + 24, !0);
+		return {
+			columnNumbers: o,
+			coordinateType: i,
+			metaValue: s,
+			metaChar: c,
+			skipLines: l,
+			format: a,
+			formatFlags: r,
+			...h(e.subarray(t + 28, t + 28 + u), this.renameRefSeq)
+		};
+	}
+	async _parse(e) {
+		let t = await this.filehandle.readFile(e), n = await r(t), i = new DataView(n.buffer, n.byteOffset, n.byteLength), o = i.getUint32(0, !0);
+		if (o !== ae && o !== oe) throw Error(`Not a CSI file ${o}`);
+		this.minShift = i.getInt32(4, !0), this.depth = i.getInt32(8, !0), this.maxBinNumber = ((1 << (this.depth + 1) * 3) - 1) / 7;
+		let c = this.maxBinNumber, l = i.getInt32(12, !0), u = l >= 30 ? this.parseAuxData(n, 16) : void 0, p = i.getInt32(16 + l, !0), h = 16 + l + 4, g, _ = [];
+		for (let e = 0; e < p; e++) {
+			_.push(h);
+			let e = i.getInt32(h, !0);
+			h += 4;
+			for (let t = 0; t < e; t++) {
+				let e = i.getUint32(h, !0);
+				if (h += 4, e > this.maxBinNumber) h += 44;
+				else {
+					g = m(n, h, 1, g), h += 8;
+					let e = i.getInt32(h, !0);
+					h += 4 + 16 * e;
+				}
+			}
+		}
+		function v(e) {
+			let t = _[e];
+			if (t === void 0) return;
+			let r = i.getInt32(t, !0);
+			t += 4;
+			let o = {}, l;
+			for (let e = 0; e < r; e++) {
+				let e = i.getUint32(t, !0);
+				if (t += 4, e > c) l = d(n, t + 28), t += 44;
+				else {
+					t += 8;
+					let r = i.getInt32(t, !0);
+					t += 4;
+					let c = Array(r);
+					for (let i = 0; i < r; i += 1) {
+						let r = s(n, t);
+						t += 8;
+						let o = s(n, t);
+						t += 8, c[i] = new a(r, o, e);
+					}
+					o[e] = c;
+				}
+			}
+			return f(Object.values(o).flat()), {
+				binIndex: o,
+				stats: l
+			};
+		}
+		return {
+			firstDataLine: g,
+			indices: y(v),
+			refCount: p,
+			csi: !0,
+			...u
+		};
+	}
+	getLowestChunk() {
+		return se;
+	}
+	reg2bins(e, t) {
+		let n = 2 ** (this.minShift + this.depth * 3);
+		t > n && (t = n), --t;
+		let r = 0, i = 0, a = this.minShift + this.depth * 3, o = [];
+		for (; r <= this.depth; a -= 3, i += D(1, r * 3), r += 1) {
+			let n = i + O(e, a), r = i + O(t, a);
+			if (r - n + o.length > this.maxBinNumber) throw Error(`query ${e}-${t} is too large for current binning scheme (shift ${this.minShift}, depth ${this.depth}), try a smaller query or a coarser index binning scheme`);
+			o.push([n, r]);
+		}
+		return o;
+	}
+}, ce = class {
+	read() {
+		throw Error("never called");
+	}
+	stat() {
+		throw Error("never called");
+	}
+	readFile() {
+		throw Error("never called");
+	}
+	close() {
+		throw Error("never called");
+	}
+}, A = {
+	BAM_FPAIRED: 1,
+	BAM_FPROPER_PAIR: 2,
+	BAM_FUNMAP: 4,
+	BAM_FMUNMAP: 8,
+	BAM_FREVERSE: 16,
+	BAM_FMREVERSE: 32,
+	BAM_FREAD1: 64,
+	BAM_FREAD2: 128,
+	BAM_FSECONDARY: 256,
+	BAM_FQCFAIL: 512,
+	BAM_FDUP: 1024,
+	BAM_FSUPPLEMENTARY: 2048
+}, j = "=ACMGRSVTWYHKDBN", M = j.split(""), N = Uint8Array.from(j, (e) => e.charCodeAt(0)), le = new Uint16Array(Uint8Array.of(1, 0).buffer)[0] === 1, P = /* @__PURE__ */ new Uint16Array(256), F = Array(256);
+for (let e = 0; e < 16; e++) for (let t = 0; t < 16; t++) {
+	let n = N[e], r = N[t];
+	P[e << 4 | t] = le ? n | r << 8 : n << 8 | r, F[e << 4 | t] = M[e] + M[t];
+}
+var ue = 300, de = 1024, fe = 0, I;
+function L() {
+	if (I === void 0) {
+		if (++fe < de) return;
+		let e = Array(65536);
+		for (let t = 0; t < 256; t++) {
+			let n = F[t], r = t << 8;
+			for (let t = 0; t < 256; t++) e[r | t] = n + F[t];
+		}
+		I = e;
+	}
+	return I;
+}
+var R = [
+	"F1F2",
+	"F1R2",
+	"R1F2",
+	"R1R2",
+	"F2F1",
+	"F2R1",
+	"R2F1",
+	"R2R1",
+	"F2F1",
+	"R2F1",
+	"F2R1",
+	"R2R1",
+	"F1F2",
+	"R1F2",
+	"F1R2",
+	"R1R2"
+], z = [
+	77,
+	73,
+	68,
+	78,
+	83,
+	72,
+	80,
+	61,
+	88,
+	63,
+	63,
+	63,
+	63,
+	63,
+	63,
+	63
+], B = new TextDecoder(), V = /* @__PURE__ */ new Map();
+function pe(e, t) {
+	let n = e << 8 | t, r = V.get(n);
+	return r === void 0 && (r = String.fromCharCode(e, t), V.set(n, r)), r;
+}
+var me = 32;
+function he(e, t, n) {
+	let r = n - t;
+	if (r < me) {
+		let i = Array(r);
+		for (let a = 0; a < r; a++) {
+			let r = e[t + a];
+			if (r > 127) return B.decode(e.subarray(t, n));
+			i[a] = r;
+		}
+		return String.fromCharCode(...i);
+	}
+	return B.decode(e.subarray(t, n));
+}
+var H = 397;
+function ge(e) {
+	return e instanceof Uint32Array || e instanceof Int32Array || Array.isArray(e);
+}
+function _e(e, t, n, r, i) {
+	let a = e.byteOffset + n;
+	switch (r) {
+		case 105: {
+			if (a % 4 == 0) return new Int32Array(e.buffer, a, i);
+			let r = Array(i);
+			for (let e = 0; e < i; e++) r[e] = t.getInt32(n + e * 4, !0);
+			return r;
+		}
+		case 73: {
+			if (a % 4 == 0) return new Uint32Array(e.buffer, a, i);
+			let r = Array(i);
+			for (let e = 0; e < i; e++) r[e] = t.getUint32(n + e * 4, !0);
+			return r;
+		}
+		case 115: {
+			if (a % 2 == 0) return new Int16Array(e.buffer, a, i);
+			let r = Array(i);
+			for (let e = 0; e < i; e++) r[e] = t.getInt16(n + e * 2, !0);
+			return r;
+		}
+		case 83: {
+			if (a % 2 == 0) return new Uint16Array(e.buffer, a, i);
+			let r = Array(i);
+			for (let e = 0; e < i; e++) r[e] = t.getUint16(n + e * 2, !0);
+			return r;
+		}
+		case 99: return new Int8Array(e.buffer, a, i);
+		case 67: return new Uint8Array(e.buffer, a, i);
+		case 102: {
+			if (a % 4 == 0) return new Float32Array(e.buffer, a, i);
+			let r = Array(i);
+			for (let e = 0; e < i; e++) r[e] = t.getFloat32(n + e * 4, !0);
+			return r;
+		}
+		default: return;
+	}
+}
+function ve(e, t) {
+	return e === 105 || e === 73 || e === 102 ? t << 2 : e === 115 || e === 83 ? t << 1 : e === 99 || e === 67 ? t : -1;
+}
+function U(e, t, n, r, i) {
+	switch (n) {
+		case 65:
+		case 99:
+		case 67: return r + 1;
+		case 115:
+		case 83: return r + 2;
+		case 105:
+		case 73:
+		case 102: return r + 4;
+		case 90:
+		case 72: {
+			let t = r;
+			for (; t < i && e[t] !== 0;) t++;
+			return t + 1;
+		}
+		case 66: {
+			let n = e[r], i = ve(n, t.getInt32(r + 1, !0));
+			return i < 0 ? (console.error("Unknown BAM B tag subtype", n), 0) : r + 5 + i;
+		}
+		default: return console.error("Unknown BAM tag type", n), 0;
+	}
+}
+function W(e, t, n, r, i, a) {
+	switch (n) {
+		case 65: return String.fromCharCode(e[r]);
+		case 105: return t.getInt32(r, !0);
+		case 73: return t.getUint32(r, !0);
+		case 99: return t.getInt8(r);
+		case 67: return t.getUint8(r);
+		case 115: return t.getInt16(r, !0);
+		case 83: return t.getUint16(r, !0);
+		case 102: return t.getFloat32(r, !0);
+		case 90:
+		case 72: return a ? e.subarray(r, i - 1) : he(e, r, i - 1);
+		default: {
+			let n = e[r], i = t.getInt32(r + 1, !0);
+			return _e(e, t, r + 5, n, i);
+		}
+	}
+}
+var G = class {
+	fileOffset;
+	_byteArray;
+	_start;
+	_end;
+	_dataView;
+	_cachedEnd;
+	_cachedTags;
+	_cachedLengthOnRef;
+	_cachedNumericCigar;
+	_cachedNUMERIC_MD;
+	_cachedSeqStart;
+	constructor(e, t, n, r, i) {
+		this._byteArray = e, this._start = t, this._end = n, this.fileOffset = r, this._dataView = i;
+	}
+	get byteArray() {
+		return this._byteArray;
+	}
+	get flags() {
+		return this._dataView.getUint16(this._start + 18, !0);
+	}
+	get ref_id() {
+		return this._dataView.getInt32(this._start + 4, !0);
+	}
+	get start() {
+		return this._dataView.getInt32(this._start + 8, !0);
+	}
+	get end() {
+		return this._cachedEnd === void 0 && (this._cachedEnd = this.start + this.length_on_ref), this._cachedEnd;
+	}
+	get mq() {
+		let e = (this.bin_mq_nl & 65280) >> 8;
+		return e === 255 ? void 0 : e;
+	}
+	get score() {
+		return this.mq;
+	}
+	get qual() {
+		let e = this.seq_length;
+		if (e === 0) return null;
+		{
+			let t = this.seqStart + (e + 1 >> 1);
+			return this._byteArray.subarray(t, t + e);
+		}
+	}
+	get strand() {
+		return this.isReverseComplemented() ? -1 : 1;
+	}
+	get b0() {
+		return this._start + 36;
+	}
+	get seqStart() {
+		return this._cachedSeqStart === void 0 && (this._cachedSeqStart = this.b0 + this.read_name_length + this.num_cigar_bytes), this._cachedSeqStart;
+	}
+	get tagsStart() {
+		let e = this.seq_length;
+		return this.seqStart + (e + 1 >> 1) + e;
+	}
+	get name() {
+		let e = this.read_name_length - 1, t = this.b0, n = this._byteArray, r = Array(e);
+		for (let i = 0; i < e; i++) r[i] = n[t + i];
+		return String.fromCharCode(...r);
+	}
+	get NUMERIC_MD() {
+		if (this._cachedNUMERIC_MD === void 0) {
+			let e = this.getTagRaw("MD");
+			this._cachedNUMERIC_MD = e instanceof Uint8Array ? e : null;
+		}
+		return this._cachedNUMERIC_MD === null ? void 0 : this._cachedNUMERIC_MD;
+	}
+	get tags() {
+		return this._cachedTags === void 0 && (this._cachedTags = this._computeTags()), this._cachedTags;
+	}
+	getTag(e) {
+		return this._cachedTags === void 0 ? this._findTag(e, !1) : this._cachedTags[e];
+	}
+	getTagRaw(e) {
+		return this._findTag(e, !0);
+	}
+	getTagAlt(e, t) {
+		if (this._cachedTags !== void 0) return this._cachedTags[e] ?? this._cachedTags[t];
+		let n = e.charCodeAt(0), r = e.charCodeAt(1), i = t.charCodeAt(0), a = t.charCodeAt(1), o = this._end, s = this._byteArray, c = this.tagsStart, l = -1, u = 0, d = 0;
+		for (; c < o;) {
+			let e = s[c], t = s[c + 1], f = s[c + 2], p = c + 3, m = U(s, this._dataView, f, p, o);
+			if (m === 0) break;
+			if (e === n && t === r) return W(s, this._dataView, f, p, m, !1);
+			l < 0 && e === i && t === a && (l = f, u = p, d = m), c = m;
+		}
+		return l < 0 ? void 0 : W(s, this._dataView, l, u, d, !1);
+	}
+	_findTag(e, t) {
+		let n = e.charCodeAt(0), r = e.charCodeAt(1), i = this._end, a = this._byteArray, o = this.tagsStart;
+		for (; o < i;) {
+			let e = a[o] === n && a[o + 1] === r, s = a[o + 2], c = o + 3, l = U(a, this._dataView, s, c, i);
+			if (l === 0) break;
+			if (e) return W(a, this._dataView, s, c, l, t);
+			o = l;
+		}
+	}
+	_computeTags() {
+		let e = this._end, t = this._byteArray, n = Object.create(null), r = this.tagsStart;
+		for (; r < e;) {
+			let i = pe(t[r], t[r + 1]), a = t[r + 2], o = r + 3, s = U(t, this._dataView, a, o, e);
+			if (s === 0) break;
+			n[i] = W(t, this._dataView, a, o, s, !1), r = s;
+		}
+		return n;
+	}
+	isPaired() {
+		return !!(this.flags & A.BAM_FPAIRED);
+	}
+	isProperlyPaired() {
+		return !!(this.flags & A.BAM_FPROPER_PAIR);
+	}
+	isSegmentUnmapped() {
+		return !!(this.flags & A.BAM_FUNMAP);
+	}
+	isMateUnmapped() {
+		return !!(this.flags & A.BAM_FMUNMAP);
+	}
+	isReverseComplemented() {
+		return !!(this.flags & A.BAM_FREVERSE);
+	}
+	isMateReverseComplemented() {
+		return !!(this.flags & A.BAM_FMREVERSE);
+	}
+	isRead1() {
+		return !!(this.flags & A.BAM_FREAD1);
+	}
+	isRead2() {
+		return !!(this.flags & A.BAM_FREAD2);
+	}
+	isSecondary() {
+		return !!(this.flags & A.BAM_FSECONDARY);
+	}
+	isFailedQc() {
+		return !!(this.flags & A.BAM_FQCFAIL);
+	}
+	isDuplicate() {
+		return !!(this.flags & A.BAM_FDUP);
+	}
+	isSupplementary() {
+		return !!(this.flags & A.BAM_FSUPPLEMENTARY);
+	}
+	_isCGTagPattern(e, t) {
+		if (t === 2) {
+			let t = this._dataView.getInt32(e, !0);
+			return (t & 15) == 4 && t >> 4 === this.seq_length;
+		}
+		return !1;
+	}
+	_computeLengthOnRef() {
+		let e = this._dataView.getInt32(this._start + 16, !0);
+		if (e & A.BAM_FUNMAP << 16) return 0;
+		let t = e & 65535, n = this.b0 + this.read_name_length;
+		if (this._isCGTagPattern(n, t)) {
+			let e = this._dataView.getInt32(n + 4, !0);
+			return (e & 15) != 3 && console.warn("CG tag with no N tag"), e >> 4;
+		}
+		let r = this._byteArray.byteOffset + n;
+		if (r % 4 == 0 && t > 50) {
+			let e = new Uint32Array(this._byteArray.buffer, r, t);
+			this._cachedNumericCigar = e;
+			let n = 0;
+			for (let r = 0; r < t; ++r) {
+				let t = e[r];
+				n += (t >> 4) * (H >> (t & 15) & 1);
+			}
+			return n;
+		}
+		let i = 0;
+		for (let e = 0; e < t; ++e) {
+			let t = this._dataView.getInt32(n + e * 4, !0);
+			i += (t >> 4) * (H >> (t & 15) & 1);
+		}
+		return i;
+	}
+	_computeNumericCigar() {
+		let e = this._dataView.getInt32(this._start + 16, !0) & 65535, t = this.b0 + this.read_name_length;
+		if (this._isCGTagPattern(t, e)) {
+			let e = this.getTag("CG");
+			return ge(e) ? e : /* @__PURE__ */ new Uint32Array();
+		}
+		let n = this._byteArray.byteOffset + t;
+		if (n % 4 == 0 && e > 50) return new Uint32Array(this._byteArray.buffer, n, e);
+		let r = Array(e);
+		for (let n = 0; n < e; ++n) r[n] = this._dataView.getInt32(t + n * 4, !0) | 0;
+		return r;
+	}
+	get length_on_ref() {
+		return this._cachedLengthOnRef === void 0 && (this._cachedLengthOnRef = this._computeLengthOnRef()), this._cachedLengthOnRef;
+	}
+	get NUMERIC_CIGAR() {
+		return this._cachedNumericCigar === void 0 && (this._cachedNumericCigar = this._computeNumericCigar()), this._cachedNumericCigar;
+	}
+	get CIGAR() {
+		let e = this.NUMERIC_CIGAR, t = "";
+		for (let n = 0, r = e.length; n < r; n++) {
+			let r = e[n];
+			t += r >> 4, t += String.fromCharCode(z[r & 15]);
+		}
+		return t;
+	}
+	get num_cigar_ops() {
+		return this.flag_nc & 65535;
+	}
+	get num_cigar_bytes() {
+		return this.num_cigar_ops << 2;
+	}
+	get read_name_length() {
+		return this.bin_mq_nl & 255;
+	}
+	get num_seq_bytes() {
+		return this.seq_length + 1 >> 1;
+	}
+	get NUMERIC_SEQ() {
+		let e = this.seqStart;
+		return this._byteArray.subarray(e, e + this.num_seq_bytes);
+	}
+	get seq() {
+		let e = this.seq_length, t = this._byteArray, n = this.seqStart, r = e >> 1, i;
+		if (e < ue) {
+			let a = L();
+			if (i = "", a === void 0) for (let e = 0; e < r; e++) i += F[t[n + e]];
+			else {
+				let e = r >> 1;
+				for (let r = 0; r < e; r++) {
+					let e = n + 2 * r;
+					i += a[t[e] << 8 | t[e + 1]];
+				}
+				r & 1 && (i += F[t[n + r - 1]]);
+			}
+			e & 1 && (i += M[(t[n + r] & 240) >> 4]);
+		} else {
+			let a = new Uint8Array(e + 1 & -2), o = new Uint16Array(a.buffer);
+			for (let e = 0; e < r; e++) o[e] = P[t[n + e]];
+			e & 1 && (a[e - 1] = N[(t[n + r] & 240) >> 4]), i = B.decode(a.subarray(0, e));
+		}
+		return i;
+	}
+	get pair_orientation() {
+		let e = this.flags;
+		if (!(e & A.BAM_FPAIRED)) return;
+		let t = this.ref_id, n = this.next_refid, r = this.start, i = this.next_pos, a = t === n ? r === i ? !!(e & A.BAM_FREAD1) : r < i : t < n;
+		return R[e >> 4 & 7 | (a ? 8 : 0)];
+	}
+	get bin_mq_nl() {
+		return this._dataView.getInt32(this._start + 12, !0);
+	}
+	get flag_nc() {
+		return this._dataView.getInt32(this._start + 16, !0);
+	}
+	get seq_length() {
+		return this._dataView.getInt32(this._start + 20, !0);
+	}
+	get next_refid() {
+		return this._dataView.getInt32(this._start + 24, !0);
+	}
+	get next_pos() {
+		return this._dataView.getInt32(this._start + 28, !0);
+	}
+	get template_length() {
+		return this._dataView.getInt32(this._start + 32, !0);
+	}
+	seqAt(e) {
+		if (e < this.seq_length) {
+			let t = this._byteArray[this.seqStart + (e >> 1)];
+			return e & 1 ? M[t & 15] : M[(t & 240) >> 4];
+		}
+	}
+	toJSON() {
+		return {
+			fileOffset: this.fileOffset,
+			ref_id: this.ref_id,
+			start: this.start,
+			end: this.end,
+			name: this.name,
+			flags: this.flags,
+			mq: this.mq,
+			CIGAR: this.CIGAR,
+			seq: this.seq,
+			next_refid: this.next_refid,
+			next_pos: this.next_pos,
+			template_length: this.template_length,
+			tags: this.tags
+		};
+	}
+};
+//#endregion
+//#region ../../node_modules/@gmod/bam/esm/sam.js
+function ye(e) {
+	let t = e.split(/\r?\n/), n = [];
+	for (let e of t) {
+		let [t, ...r] = e.split(/\t/);
+		t && n.push({
+			tag: t.slice(1),
+			data: r.map((e) => {
+				let t = e.indexOf(":");
+				return t === -1 ? {
+					tag: e,
+					value: ""
+				} : {
+					tag: e.slice(0, t),
+					value: e.slice(t + 1)
+				};
+			})
+		});
+	}
+	return n;
+}
+var K = 65536, be = 33554432;
+function q(n, r, i) {
+	return n ?? (r ? new e(r) : i ? new t(i) : void 0);
+}
+function J(e, t, n) {
+	let r = e?.[0];
+	return !!r && (r.ref_id > t || r.ref_id === t && r.start >= n);
+}
+function Y(e) {
+	let { minv: t, maxv: n } = e;
+	return `${t.blockPosition}:${t.dataPosition}-${n.blockPosition}:${n.dataPosition}`;
+}
+var X = 104857600, Z = 6, xe = class {
+	_maxBytes;
+	entries = /* @__PURE__ */ new Map();
+	bytes = 0;
+	constructor(e) {
+		this._maxBytes = e;
+	}
+	get maxBytes() {
+		return this._maxBytes;
+	}
+	set maxBytes(e) {
+		this._maxBytes = e, this.evict();
+	}
+	get size() {
+		return this.entries.size;
+	}
+	get byteSize() {
+		return this.bytes;
+	}
+	get(e) {
+		let t = this.entries.get(e);
+		return t && (this.entries.delete(e), this.entries.set(e, t)), t;
+	}
+	set(e, t) {
+		this.delete(e), this.entries.set(e, t), this.bytes += t.bytes, this.evict();
+	}
+	evict() {
+		let e = this.entries.keys();
+		for (; this.bytes > this._maxBytes && this.entries.size > 1;) this.delete(e.next().value);
+	}
+	delete(e) {
+		let t = this.entries.get(e);
+		t && (this.entries.delete(e), this.bytes -= t.bytes);
+	}
+	clear() {
+		this.entries.clear(), this.bytes = 0;
+	}
+	[Symbol.iterator]() {
+		return this.entries[Symbol.iterator]();
+	}
+}, Q = class {
+	renameRefSeq;
+	bam;
+	header;
+	chrToIndex;
+	indexToChr;
+	index;
+	htsget = !1;
+	headerP;
+	chunkFeatureCache;
+	inFlightChunks = /* @__PURE__ */ new Map();
+	RecordClass;
+	constructor({ bamFilehandle: e, bamPath: t, bamUrl: n, baiPath: r, baiFilehandle: i, baiUrl: a, csiPath: o, csiFilehandle: s, csiUrl: c, htsget: l, renameRefSeqs: u = (e) => e, recordClass: d, maxCacheBytes: f = X }) {
+		this.renameRefSeq = u, this.RecordClass = d ?? G, this.chunkFeatureCache = new xe(f);
+		let p = q(e, t, n);
+		if (p) this.bam = p;
+		else if (l) this.htsget = !0, this.bam = new ce();
+		else throw Error("no bam source: pass bamFilehandle, bamPath, bamUrl, or htsget: true");
+		let m = q(s, o, c), h = q(i, r, a) ?? q(void 0, t ? `${t}.bai` : void 0, n ? `${n}.bai` : void 0);
+		if (m) this.index = new k({ filehandle: m });
+		else if (h) this.index = new T({ filehandle: h });
+		else if (!l) throw Error("no index source: pass csi*/bai* options or a bamPath/bamUrl so the .bai sibling can be inferred");
+	}
+	async getHeaderPre(e = {}) {
+		if (!this.index) throw Error("no index to read a header from: use HtsgetFile for htsget sources");
+		let t = await this.index.parse(e), n = t.firstDataLine === void 0 ? K : t.firstDataLine.blockPosition + K, i, a;
+		do {
+			let t = await this.bam.read(n, 0, { signal: e.signal });
+			a = t.length < n, i = this.applyHeader(await r(t)), n *= 2;
+		} while (i === void 0 && !a && n <= be);
+		if (i === void 0) throw Error("Insufficient data for reference sequences");
+		return i;
+	}
+	applyHeader(e) {
+		let t = new DataView(e.buffer, e.byteOffset, e.byteLength);
+		if (t.getInt32(0, !0) !== 21840194) throw Error("Not a BAM file");
+		let n = t.getInt32(4, !0), r = p(e, n + 8, this.renameRefSeq), i;
+		if (r) {
+			let t = new TextDecoder("utf8").decode(e.subarray(8, 8 + n));
+			this.header = t, this.chrToIndex = r.chrToIndex, this.indexToChr = r.indexToChr, i = ye(t);
+		}
+		return i;
+	}
+	getHeader(e) {
+		return this.headerP ||= this.getHeaderPre(e).catch((e) => {
+			throw this.headerP = void 0, e;
+		}), this.headerP;
+	}
+	async getHeaderText(e = {}) {
+		return await this.getHeader(e), this.header;
+	}
+	async getSeqId(e, t) {
+		return await this.getHeader(t), this.chrToIndex?.[e];
+	}
+	async getRecordsForRange(e, t, n, r) {
+		c(r?.signal);
+		let i = await this.getSeqId(e, r);
+		if (i === void 0 || !this.index) return [];
+		let a = await this.index.blocksForRange(i, t - 1, n, r);
+		return this._fetchChunkFeatures(a, i, t, n, r);
+	}
+	_startChunkRead(e, t) {
+		let n = new AbortController(), r = this._readChunkFeatures(t, { signal: n.signal }).then((t) => (this.chunkFeatureCache.set(e, t), t)), i = {
+			promise: r,
+			signals: /* @__PURE__ */ new Set(),
+			pinned: !1,
+			controller: n,
+			dispose: new AbortController(),
+			settled: !1
+		};
+		this.inFlightChunks.set(e, i);
+		let a = () => {
+			i.settled = !0, i.dispose.abort(), i.signals.clear(), this.inFlightChunks.get(e) === i && this.inFlightChunks.delete(e);
+		};
+		return r.then(a, a), i;
+	}
+	joinChunkRead(e, t) {
+		t === void 0 ? e.pinned = !0 : t.aborted ? !e.pinned && e.signals.size === 0 && e.controller.abort(t.reason) : e.signals.has(t) || (e.signals.add(t), t.addEventListener("abort", () => {
+			e.signals.delete(t), !e.pinned && e.signals.size === 0 && e.controller.abort(t.reason);
+		}, {
+			once: !0,
+			signal: e.dispose.signal
+		}));
+	}
+	async _cachedChunkFeatures(e, t) {
+		c(t.signal);
+		let n = Y(e), r = this.chunkFeatureCache.get(n);
+		if (r) return r.features;
+		let i = this.inFlightChunks.get(n);
+		i?.controller.signal.aborted && !i.settled && (this.inFlightChunks.get(n) === i && this.inFlightChunks.delete(n), i = void 0), i ??= this._startChunkRead(n, e), i.settled || this.joinChunkRead(i, t.signal);
+		try {
+			let e = await i.promise;
+			return c(t.signal), e.features;
+		} catch (e) {
+			throw c(t.signal), e;
+		}
+	}
+	async _fetchChunkFeatures(e, t, n, r, i = {}) {
+		let { viewAsPairs: a, onProgress: o } = i, s = [], c = 0;
+		for (let t = 0, n = e.length; t < n; t++) c += e[t].fetchedSize();
+		let l = 0;
+		o?.(0, c);
+		let u = Array(e.length);
+		if (e.length === 1) u[0] = await this._cachedChunkFeatures(e[0], i), o?.(c, c);
+		else {
+			let n = async (t) => {
+				let n = e[t];
+				u[t] = await this._cachedChunkFeatures(n, i), l += n.fetchedSize(), o?.(l, c);
+			}, a = Math.min(Z, e.length);
+			await Promise.all(Array.from({ length: a }, (e, t) => n(t)));
+			let s = !1;
+			for (let e = 0; e < a; e++) if (J(u[e], t, r)) {
+				s = !0;
+				break;
+			}
+			if (!s && e.length > a) {
+				let t = a, r = async () => {
+					for (; t < e.length;) await n(t++);
+				}, i = Math.min(Z, e.length - a);
+				await Promise.all(Array.from({ length: i }, () => r()));
+			}
+			o?.(c, c);
+		}
+		for (let i = 0, a = e.length; i < a; i++) {
+			let e = u[i];
+			e && v(e, t, n, r, s);
+		}
+		if (a) {
+			let e = await this.fetchPairs(t, s, i);
+			for (let t = 0, n = e.length; t < n; t++) s.push(e[t]);
+		}
+		return s;
+	}
+	async fetchPairs(e, t, n) {
+		let { pairAcrossChr: r, maxInsertSize: i = 2e5 } = n, a = /* @__PURE__ */ new Map(), o = /* @__PURE__ */ new Set();
+		for (let e = 0, n = t.length; e < n; e++) {
+			let n = t[e], r = n.name;
+			a.set(r, (a.get(r) ?? 0) + 1), o.add(n.fileOffset);
+		}
+		let s = [];
+		for (let o = 0, c = t.length; o < c; o++) {
+			let c = t[o], l = c.name;
+			this.index && a.get(l) === 1 && (r || c.next_refid === e && Math.abs(c.start - c.next_pos) < i) && s.push(this.index.blocksForRange(c.next_refid, c.next_pos, c.next_pos + 1, n));
+		}
+		let c = /* @__PURE__ */ new Map(), l = await Promise.all(s);
+		for (let e = 0, t = l.length; e < t; e++) {
+			let t = l[e];
+			for (let e = 0, n = t.length; e < n; e++) {
+				let n = t[e];
+				c.set(Y(n), n);
+			}
+		}
+		let u = [...c.values()], d = Array(u.length), f = 0, p = async () => {
+			for (; f < u.length;) {
+				let e = f++, t = await this._cachedChunkFeatures(u[e], n), r = [];
+				for (let e = 0, n = t.length; e < n; e++) {
+					let n = t[e];
+					!o.has(n.fileOffset) && a.get(n.name) === 1 && r.push(n);
+				}
+				d[e] = r;
+			}
+		}, m = Math.min(Z, u.length);
+		return await Promise.all(Array.from({ length: m }, () => p())), d.flat();
+	}
+	async _readChunkFeatures(e, t) {
+		let n = await this.bam.read(e.fetchedSize(), e.minv.blockPosition, { signal: t.signal });
+		c(t.signal);
+		let { buffer: r, cpositions: a, dpositions: o } = await i(n, e);
+		return {
+			features: this.readBamFeatures(r, a, o, e),
+			bytes: r.byteLength
+		};
+	}
+	readBamFeatures(e, t, n, r) {
+		let i = 0, a = [], o = 0, s = new DataView(e.buffer, e.byteOffset, e.byteLength), c = n.length > 0, l = t.length > 0;
+		for (; i + 4 < e.length;) {
+			let u = s.getInt32(i, !0), d = i + 4 + u - 1;
+			if (c) {
+				let e = i + r.minv.dataPosition;
+				for (; o < n.length && e >= n[o];) o++;
+			}
+			if (d < e.length) {
+				let c = new this.RecordClass(e, i, d, l ? t[o] * 256 + (i - n[o]) + r.minv.dataPosition + 1 : ie(e.subarray(i, d)) >>> 0, s);
+				a.push(c);
+			}
+			i = d + 1;
+		}
+		return a;
+	}
+	async hasRefSeq(e, t) {
+		let n = await this.getSeqId(e, t);
+		return !this.index || n === void 0 ? !1 : this.index.hasRefSeq(n);
+	}
+	async lineCount(e, t) {
+		let n = await this.getSeqId(e, t);
+		return !this.index || n === void 0 ? 0 : this.index.lineCount(n);
+	}
+	async indexCov(e, t, n) {
+		let r = await this.getSeqId(e);
+		return !this.index || r === void 0 ? [] : this.index.indexCov(r, t, n);
+	}
+	async blocksForRange(e, t, n, r) {
+		let i = await this.getSeqId(e, r);
+		return !this.index || i === void 0 ? [] : this.index.blocksForRange(i, t, n, r);
+	}
+	clearFeatureCache() {
+		this.chunkFeatureCache.clear();
+	}
+	async estimatedBytesForRegions(e, t) {
+		if (!this.index) return 0;
+		await this.getHeader(t);
+		let n = e.flatMap((e) => {
+			let t = this.chrToIndex?.[e.refName];
+			return t === void 0 ? [] : [{
+				refId: t,
+				start: e.start,
+				end: e.end
+			}];
+		});
+		return this.index.estimatedBytesForRegions(n, t);
+	}
+};
+//#endregion
+//#region ../../node_modules/@gmod/bam/esm/htsget.js
+function Se(e) {
+	let t;
+	try {
+		t = JSON.parse(e);
+	} catch {}
+	let n = t?.htsget;
+	return n?.error === void 0 ? e : n.message === void 0 ? n.error : `${n.error}: ${n.message}`;
+}
+function Ce(e) {
+	let t = new DataView(e.buffer, e.byteOffset, e.byteLength);
+	if (e.byteLength < 8 || t.getInt32(0, !0) !== 21840194) return 0;
+	let n = p(e, 8 + t.getInt32(4, !0), (e) => e);
+	if (!n) throw Error("truncated BAM header in htsget response");
+	return n.end;
+}
+async function $(e, t, n) {
+	let r = await e(t, n);
+	if (!r.ok) throw Error(`HTTP ${r.status} fetching ${t}: ${Se(await r.text())}`);
+	return r;
+}
+async function we(e, { url: t, headers: n }, r) {
+	let { referer: i, ...a } = n ?? {}, o = t.startsWith("data:") ? await $(e, t) : await $(e, t, {
+		...r,
+		headers: a
+	});
+	return new Uint8Array(await o.arrayBuffer());
+}
+async function Te(e, t, n) {
+	return g(await Promise.all(t.map(async (t) => r(await we(e, t, n)))));
+}
+var Ee = class extends Q {
+	baseUrl;
+	trackId;
+	fetchFn;
+	constructor(e) {
+		super({
+			htsget: !0,
+			recordClass: e.recordClass
+		}), this.baseUrl = e.baseUrl, this.trackId = e.trackId, this.fetchFn = e.fetch ?? ((e, t) => fetch(e, t));
+	}
+	async fetchTicket(e, t) {
+		let n = `${this.baseUrl}/${this.trackId}?${e}`, r = await (await $(this.fetchFn, n, { signal: t?.signal })).json();
+		return Te(this.fetchFn, r.htsget.urls, { signal: t?.signal });
+	}
+	async getRecordsForRange(e, t, n, r) {
+		await this.getHeader(r);
+		let i = this.chrToIndex?.[e];
+		if (i === void 0) return [];
+		let s = await this.fetchTicket(`referenceName=${e}&start=${t}&end=${n}&format=BAM`, r), c = new o(0, 0);
+		return v(this.readBamFeatures(s.subarray(Ce(s)), [], [], new a(c, c, 0)), i, t, n);
+	}
+	async getHeaderPre(e = {}) {
+		let t = await this.fetchTicket("class=header&format=BAM", e), n = this.applyHeader(t);
+		if (!n) throw Error("Insufficient data for reference sequences");
+		return n;
+	}
+};
+//#endregion
+export { T as BAI, Q as BamFile, G as BamRecord, k as CSI, X as DEFAULT_MAX_CACHE_BYTES, Ee as HtsgetFile };