@gmod/bam 7.6.2 → 7.7.0

package/esm/util.d.ts

@@ -1,6 +1,6 @@
 import Chunk from './chunk.ts';
 import { VirtualOffset } from './virtualOffset.ts';
-import type { Offset } from './virtualOffset.ts';
+import type { OffsetCoords } from './virtualOffset.ts';
 export interface BamOpts {
     viewAsPairs?: boolean;
     pairAcrossChr?: boolean;
@@ -25,11 +25,11 @@ export interface BaseOpts {
      */
     onProgress?: (bytesDownloaded: number, totalBytes?: number) => void;
 }
-export declare function optimizeChunks(chunks: Chunk[], lowest?: Offset): Chunk[];
+export declare function optimizeChunks(chunks: Chunk[], lowest?: OffsetCoords): Chunk[];
 export declare function parsePseudoBin(bytes: Uint8Array, offset: number): {
     lineCount: number;
 };
-export declare function clampChunkEnds(chunks: Chunk[], extraBoundaries?: number[]): void;
+export declare function clampChunkEnds(chunks: Chunk[], extraBoundaries?: ArrayLike<number>): void;
 export declare function parseRefSeqs(uncba: Uint8Array, start: number, renameRefSeq: (s: string) => string): {
     chrToIndex: Record<string, number>;
     indexToChr: {

package/esm/util.js

@@ -1,5 +1,4 @@
 import Chunk from "./chunk.js";
-import { longFromBytesToUnsigned } from "./long.js";
 import { VirtualOffset } from "./virtualOffset.js";
 export function optimizeChunks(chunks, lowest) {
     const n = chunks.length;
@@ -60,9 +59,20 @@ export function optimizeChunks(chunks, lowest) {
     }
     return mergedChunks;
 }
+// The pseudo-bin's mapped-record count, a little-endian uint64. Read as two
+// 32-bit halves rather than via BigInt: exact up to Number.MAX_SAFE_INTEGER,
+// which is far past any real record count, and allocates nothing.
 export function parsePseudoBin(bytes, offset) {
+    const low = bytes[offset] |
+        (bytes[offset + 1] << 8) |
+        (bytes[offset + 2] << 16) |
+        (bytes[offset + 3] << 24);
+    const high = bytes[offset + 4] |
+        (bytes[offset + 5] << 8) |
+        (bytes[offset + 6] << 16) |
+        (bytes[offset + 7] << 24);
     return {
-        lineCount: longFromBytesToUnsigned(bytes, offset),
+        lineCount: (high >>> 0) * 2 ** 32 + (low >>> 0),
     };
 }
 // Tighten each chunk's endPosition (default: a full max-size BGZF block past
@@ -73,11 +83,17 @@ export function parsePseudoBin(bytes, offset) {
 // least the true block end, so the clamped fetch still contains the whole block.
 // Shrinks both the byte estimate and the actual fetch with no extra I/O.
 export function clampChunkEnds(chunks, extraBoundaries = []) {
-    const boundaries = [...extraBoundaries];
+    // Float64Array rather than number[]: `extraBoundaries` is the linear index,
+    // which runs to tens of thousands of entries on a human-sized reference, and
+    // this is the only place they get copied.
+    const boundaries = new Float64Array(extraBoundaries.length + chunks.length * 2);
+    boundaries.set(extraBoundaries);
+    let n = extraBoundaries.length;
     for (const c of chunks) {
-        boundaries.push(c.minv.blockPosition, c.maxv.blockPosition);
+        boundaries[n++] = c.minv.blockPosition;
+        boundaries[n++] = c.maxv.blockPosition;
     }
-    boundaries.sort((a, b) => a - b);
+    boundaries.sort();
     for (const c of chunks) {
         const max = c.maxv.blockPosition;
         // first boundary strictly greater than max

package/esm/util.js.map

@@ -1 +1 @@
-{"version":3,"file":"util.js","sourceRoot":"","sources":["../src/util.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,MAAM,YAAY,CAAA;AAC9B,OAAO,EAAE,uBAAuB,EAAE,MAAM,WAAW,CAAA;AACnD,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAA;AA8BlD,MAAM,UAAU,cAAc,CAAC,MAAe,EAAE,MAAe;IAC7D,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAA;IACvB,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACZ,OAAO,MAAM,CAAA;IACf,CAAC;IAED,0DAA0D;IAC1D,IAAI,QAAiB,CAAA;IACrB,IAAI,MAAM,EAAE,CAAC;QACX,MAAM,WAAW,GAAG,MAAM,CAAC,aAAa,CAAA;QACxC,MAAM,UAAU,GAAG,MAAM,CAAC,YAAY,CAAA;QACtC,QAAQ,GAAG,EAAE,CAAA;QACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,KAAK,GAAG,MAAM,CAAC,CAAC,CAAE,CAAA;YACxB,MAAM,IAAI,GAAG,KAAK,CAAC,IAAI,CAAA;YACvB,MAAM,GAAG,GACP,IAAI,CAAC,aAAa,GAAG,WAAW,IAAI,IAAI,CAAC,YAAY,GAAG,UAAU,CAAA;YACpE,IAAI,GAAG,GAAG,CAAC,EAAE,CAAC;gBACZ,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,CAAA;YACtB,CAAC;QACH,CAAC;QACD,IAAI,QAAQ,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;YAC1B,OAAO,QAAQ,CAAA;QACjB,CAAC;IACH,CAAC;SAAM,CAAC;QACN,QAAQ,GAAG,MAAM,CAAA;IACnB,CAAC;IAED,QAAQ,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE;QACvB,MAAM,GAAG,GAAG,EAAE,CAAC,IAAI,CAAC,aAAa,GAAG,EAAE,CAAC,IAAI,CAAC,aAAa,CAAA;QACzD,OAAO,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,YAAY,GAAG,EAAE,CAAC,IAAI,CAAC,YAAY,CAAA;IACtE,CAAC,CAAC,CAAA;IAEF,yEAAyE;IACzE,uEAAuE;IACvE,MAAM,YAAY,GAAY,CAAC,QAAQ,CAAC,CAAC,CAAE,CAAC,CAAA;IAC5C,IAAI,YAAY,GAAG,QAAQ,CAAC,CAAC,CAAE,CAAC,IAAI,CAAC,aAAa,CAAA;IAClD,IAAI,YAAY,GAAG,QAAQ,CAAC,CAAC,CAAE,CAAC,IAAI,CAAC,aAAa,CAAA;IAElD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACzC,MAAM,KAAK,GAAG,QAAQ,CAAC,CAAC,CAAE,CAAA;QAC1B,MAAM,aAAa,GAAG,KAAK,CAAC,IAAI,CAAC,aAAa,CAAA;QAC9C,MAAM,aAAa,GAAG,KAAK,CAAC,IAAI,CAAC,aAAa,CAAA;QAC9C,oEAAoE;QACpE,yEAAyE;QACzE,IACE,aAAa,GAAG,YAAY,GAAG,KAAK;YACpC,aAAa,GAAG,YAAY,GAAG,OAAO,EACtC,CAAC;YACD,MAAM,SAAS,GAAG,YAAY,CAAC,YAAY,CAAC,MAAM,GAAG,CAAC,CAAE,CAAA;YACxD,MAAM,GAAG,GACP,aAAa,GAAG,YAAY;gBAC5B,KAAK,CAAC,IAAI,CAAC,YAAY,GAAG,SAAS,CAAC,IAAI,CAAC,YAAY,CAAA;YACvD,IAAI,GAAG,GAAG,CAAC,EAAE,CAAC;gBACZ,YAAY,CAAC,YAAY,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,IAAI,KAAK,CAC/C,SAAS,CAAC,IAAI,EACd,KAAK,CAAC,IAAI,EACV,SAAS,CAAC,GAAG,EACb,KAAK,CAAC,WAAW,CAClB,CAAA;gBACD,YAAY,GAAG,aAAa,CAAA;YAC9B,CAAC;QACH,CAAC;aAAM,CAAC;YACN,YAAY,CAAC,IAAI,CAAC,KAAK,CAAC,CAAA;YACxB,YAAY,GAAG,aAAa,CAAA;YAC5B,YAAY,GAAG,aAAa,CAAA;QAC9B,CAAC;IACH,CAAC;IAED,OAAO,YAAY,CAAA;AACrB,CAAC;AAED,MAAM,UAAU,cAAc,CAAC,KAAiB,EAAE,MAAc;IAC9D,OAAO;QACL,SAAS,EAAE,uBAAuB,CAAC,KAAK,EAAE,MAAM,CAAC;KAClD,CAAA;AACH,CAAC;AAED,6EAA6E;AAC7E,kFAAkF;AAClF,kFAAkF;AAClF,iFAAiF;AACjF,kFAAkF;AAClF,iFAAiF;AACjF,yEAAyE;AACzE,MAAM,UAAU,cAAc,CAC5B,MAAe,EACf,kBAA4B,EAAE;IAE9B,MAAM,UAAU,GAAG,CAAC,GAAG,eAAe,CAAC,CAAA;IACvC,KAAK,MAAM,CAAC,IAAI,MAAM,EAAE,CAAC;QACvB,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,aAAa,EAAE,CAAC,CAAC,IAAI,CAAC,aAAa,CAAC,CAAA;IAC7D,CAAC;IACD,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;IAEhC,KAAK,MAAM,CAAC,IAAI,MAAM,EAAE,CAAC;QACvB,MAAM,GAAG,GAAG,CAAC,CAAC,IAAI,CAAC,aAAa,CAAA;QAChC,2CAA2C;QAC3C,IAAI,EAAE,GAAG,CAAC,CAAA;QACV,IAAI,EAAE,GAAG,UAAU,CAAC,MAAM,CAAA;QAC1B,OAAO,EAAE,GAAG,EAAE,EAAE,CAAC;YACf,MAAM,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,CAAA;YAC3B,IAAI,UAAU,CAAC,GAAG,CAAE,GAAG,GAAG,EAAE,CAAC;gBAC3B,EAAE,GAAG,GAAG,CAAA;YACV,CAAC;iBAAM,CAAC;gBACN,EAAE,GAAG,GAAG,GAAG,CAAC,CAAA;YACd,CAAC;QACH,CAAC;QACD,IAAI,EAAE,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC;YAC3B,CAAC,CAAC,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,WAAW,EAAE,UAAU,CAAC,EAAE,CAAE,CAAC,CAAA;QAC1D,CAAC;IACH,CAAC;AACH,CAAC;AAED,6EAA6E;AAC7E,4EAA4E;AAC5E,eAAe;AACf,MAAM,UAAU,YAAY,CAC1B,KAAiB,EACjB,KAAa,EACb,YAAmC;IAEnC,IAAI,KAAK,GAAG,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC;QAC7B,OAAO,SAAS,CAAA;IAClB,CAAC;IACD,MAAM,QAAQ,GAAG,IAAI,QAAQ,CAC3B,KAAK,CAAC,MAAM,EACZ,KAAK,CAAC,UAAU,EAChB,KAAK,CAAC,UAAU,CACjB,CAAA;IACD,MAAM,IAAI,GAAG,QAAQ,CAAC,QAAQ,CAAC,KAAK,EAAE,IAAI,CAAC,CAAA;IAC3C,kEAAkE;IAClE,uDAAuD;IACvD,MAAM,UAAU,GAA2B,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC,CAAA;IAC9D,MAAM,UAAU,GAA0C,EAAE,CAAA;IAC5D,MAAM,OAAO,GAAG,IAAI,WAAW,CAAC,MAAM,CAAC,CAAA;IAEvC,IAAI,CAAC,GAAG,KAAK,GAAG,CAAC,CAAA;IACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;QAC9B,IAAI,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC;YACzB,OAAO,SAAS,CAAA;QAClB,CAAC;QACD,MAAM,KAAK,GAAG,QAAQ,CAAC,QAAQ,CAAC,CAAC,EAAE,IAAI,CAAC,CAAA;QACxC,IAAI,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC;YACjC,OAAO,SAAS,CAAA;QAClB,CAAC;QACD,MAAM,OAAO,GAAG,YAAY,CAC1B,OAAO,CAAC,MAAM,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CACzD,CAAA;QACD,MAAM,IAAI,GAAG,QAAQ,CAAC,QAAQ,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,EAAE,IAAI,CAAC,CAAA;QACnD,UAAU,CAAC,OAAO,CAAC,GAAG,CAAC,CAAA;QACvB,UAAU,CAAC,IAAI,CAAC,EAAE,OAAO,EAAE,MAAM,EAAE,IAAI,EAAE,CAAC,CAAA;QAC1C,CAAC,IAAI,CAAC,GAAG,KAAK,CAAA;IAChB,CAAC;IACD,6EAA6E;IAC7E,OAAO,EAAE,UAAU,EAAE,UAAU,EAAE,GAAG,EAAE,CAAC,EAAE,CAAA;AAC3C,CAAC;AAED,2EAA2E;AAC3E,+EAA+E;AAC/E,4CAA4C;AAC5C;;;;;;;;;;;;;;GAcG;AACH,MAAM,UAAU,gBAAgB,CAC9B,KAAiB,EACjB,MAAc,EACd,KAAa,EACb,OAAkC;IAElC,IAAI,QAAQ,GAAG,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC,aAAa,CAAC,CAAC,CAAC,QAAQ,CAAA;IACzD,IAAI,OAAO,GAAG,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC,YAAY,CAAC,CAAC,CAAC,CAAC,CAAA;IAChD,IAAI,KAAK,GAAG,KAAK,CAAA;IACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,MAAM,CAAC,GAAG,MAAM,GAAG,CAAC,GAAG,CAAC,CAAA;QACxB,MAAM,KAAK,GACT,KAAK,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,aAAa;YAC7B,KAAK,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,WAAW;YAC3B,KAAK,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,SAAS;YACzB,KAAK,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,OAAO;YACvB,KAAK,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,KAAK;YACrB,KAAK,CAAC,CAAC,GAAG,CAAC,CAAE,CAAA;QACf,MAAM,IAAI,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAE,IAAI,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAE,CAAA;QAC7C,IACE,CAAC,KAAK,KAAK,CAAC,IAAI,IAAI,KAAK,CAAC,CAAC;YAC3B,CAAC,KAAK,GAAG,QAAQ,IAAI,CAAC,KAAK,KAAK,QAAQ,IAAI,IAAI,GAAG,OAAO,CAAC,CAAC,EAC5D,CAAC;YACD,QAAQ,GAAG,KAAK,CAAA;YAChB,OAAO,GAAG,IAAI,CAAA;YACd,KAAK,GAAG,IAAI,CAAA;QACd,CAAC;IACH,CAAC;IACD,OAAO,KAAK,CAAC,CAAC,CAAC,IAAI,aAAa,CAAC,QAAQ,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,OAAO,CAAA;AAC/D,CAAC;AAED,4FAA4F;AAC5F,MAAM,UAAU,cAAc,CAC5B,UAAsB,EACtB,eAAwC,CAAC,CAAC,EAAE,CAAC,CAAC;IAE9C,MAAM,OAAO,GAAG,IAAI,WAAW,EAAE,CAAA;IACjC,IAAI,SAAS,GAAG,CAAC,CAAA;IACjB,IAAI,aAAa,GAAG,CAAC,CAAA;IACrB,MAAM,WAAW,GAAa,EAAE,CAAA;IAChC,MAAM,WAAW,GAA2B,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC,CAAA;IAC/D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3C,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,CAAC;YACnB,IAAI,aAAa,GAAG,CAAC,EAAE,CAAC;gBACtB,MAAM,OAAO,GAAG,YAAY,CAC1B,OAAO,CAAC,MAAM,CAAC,UAAU,CAAC,QAAQ,CAAC,aAAa,EAAE,CAAC,CAAC,CAAC,CACtD,CAAA;gBACD,WAAW,CAAC,SAAS,CAAC,GAAG,OAAO,CAAA;gBAChC,WAAW,CAAC,OAAO,CAAC,GAAG,SAAS,CAAA;YAClC,CAAC;YACD,aAAa,GAAG,CAAC,GAAG,CAAC,CAAA;YACrB,SAAS,EAAE,CAAA;QACb,CAAC;IACH,CAAC;IACD,OAAO,EAAE,WAAW,EAAE,WAAW,EAAE,CAAA;AACrC,CAAC;AAED,MAAM,UAAU,gBAAgB,CAAC,IAAkB;IACjD,IAAI,WAAW,GAAG,CAAC,CAAA;IACnB,KAAK,MAAM,KAAK,IAAI,IAAI,EAAE,CAAC;QACzB,WAAW,IAAI,KAAK,CAAC,MAAM,CAAA;IAC7B,CAAC;IACD,MAAM,WAAW,GAAG,IAAI,UAAU,CAAC,WAAW,CAAC,CAAA;IAC/C,IAAI,MAAM,GAAG,CAAC,CAAA;IACd,KAAK,MAAM,KAAK,IAAI,IAAI,EAAE,CAAC;QACzB,WAAW,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,CAAA;QAC9B,MAAM,IAAI,KAAK,CAAC,MAAM,CAAA;IACxB,CAAC;IACD,OAAO,WAAW,CAAA;AACpB,CAAC;AAQD,0EAA0E;AAC1E,4EAA4E;AAC5E,6EAA6E;AAC7E,iDAAiD;AACjD,MAAM,UAAU,aAAa,CAC3B,OAAY,EACZ,KAAa,EACb,GAAW,EACX,GAAW,EACX,MAAW,EAAE;IAEb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/C,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAE,CAAA;QACrB,IAAI,CAAC,CAAC,MAAM,KAAK,KAAK,EAAE,CAAC;YACvB,IAAI,CAAC,CAAC,KAAK,IAAI,GAAG,EAAE,CAAC;gBACnB,MAAK;YACP,CAAC;iBAAM,IAAI,CAAC,CAAC,GAAG,IAAI,GAAG,EAAE,CAAC;gBACxB,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAA;YACb,CAAC;QACH,CAAC;aAAM,IAAI,CAAC,CAAC,MAAM,GAAG,KAAK,EAAE,CAAC;YAC5B,MAAK;QACP,CAAC;IACH,CAAC;IACD,OAAO,GAAG,CAAA;AACZ,CAAC"}
+{"version":3,"file":"util.js","sourceRoot":"","sources":["../src/util.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,MAAM,YAAY,CAAA;AAC9B,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAA;AA8BlD,MAAM,UAAU,cAAc,CAAC,MAAe,EAAE,MAAqB;IACnE,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAA;IACvB,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACZ,OAAO,MAAM,CAAA;IACf,CAAC;IAED,0DAA0D;IAC1D,IAAI,QAAiB,CAAA;IACrB,IAAI,MAAM,EAAE,CAAC;QACX,MAAM,WAAW,GAAG,MAAM,CAAC,aAAa,CAAA;QACxC,MAAM,UAAU,GAAG,MAAM,CAAC,YAAY,CAAA;QACtC,QAAQ,GAAG,EAAE,CAAA;QACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,KAAK,GAAG,MAAM,CAAC,CAAC,CAAE,CAAA;YACxB,MAAM,IAAI,GAAG,KAAK,CAAC,IAAI,CAAA;YACvB,MAAM,GAAG,GACP,IAAI,CAAC,aAAa,GAAG,WAAW,IAAI,IAAI,CAAC,YAAY,GAAG,UAAU,CAAA;YACpE,IAAI,GAAG,GAAG,CAAC,EAAE,CAAC;gBACZ,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,CAAA;YACtB,CAAC;QACH,CAAC;QACD,IAAI,QAAQ,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;YAC1B,OAAO,QAAQ,CAAA;QACjB,CAAC;IACH,CAAC;SAAM,CAAC;QACN,QAAQ,GAAG,MAAM,CAAA;IACnB,CAAC;IAED,QAAQ,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE;QACvB,MAAM,GAAG,GAAG,EAAE,CAAC,IAAI,CAAC,aAAa,GAAG,EAAE,CAAC,IAAI,CAAC,aAAa,CAAA;QACzD,OAAO,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,YAAY,GAAG,EAAE,CAAC,IAAI,CAAC,YAAY,CAAA;IACtE,CAAC,CAAC,CAAA;IAEF,yEAAyE;IACzE,uEAAuE;IACvE,MAAM,YAAY,GAAY,CAAC,QAAQ,CAAC,CAAC,CAAE,CAAC,CAAA;IAC5C,IAAI,YAAY,GAAG,QAAQ,CAAC,CAAC,CAAE,CAAC,IAAI,CAAC,aAAa,CAAA;IAClD,IAAI,YAAY,GAAG,QAAQ,CAAC,CAAC,CAAE,CAAC,IAAI,CAAC,aAAa,CAAA;IAElD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACzC,MAAM,KAAK,GAAG,QAAQ,CAAC,CAAC,CAAE,CAAA;QAC1B,MAAM,aAAa,GAAG,KAAK,CAAC,IAAI,CAAC,aAAa,CAAA;QAC9C,MAAM,aAAa,GAAG,KAAK,CAAC,IAAI,CAAC,aAAa,CAAA;QAC9C,oEAAoE;QACpE,yEAAyE;QACzE,IACE,aAAa,GAAG,YAAY,GAAG,KAAK;YACpC,aAAa,GAAG,YAAY,GAAG,OAAO,EACtC,CAAC;YACD,MAAM,SAAS,GAAG,YAAY,CAAC,YAAY,CAAC,MAAM,GAAG,CAAC,CAAE,CAAA;YACxD,MAAM,GAAG,GACP,aAAa,GAAG,YAAY;gBAC5B,KAAK,CAAC,IAAI,CAAC,YAAY,GAAG,SAAS,CAAC,IAAI,CAAC,YAAY,CAAA;YACvD,IAAI,GAAG,GAAG,CAAC,EAAE,CAAC;gBACZ,YAAY,CAAC,YAAY,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,IAAI,KAAK,CAC/C,SAAS,CAAC,IAAI,EACd,KAAK,CAAC,IAAI,EACV,SAAS,CAAC,GAAG,EACb,KAAK,CAAC,WAAW,CAClB,CAAA;gBACD,YAAY,GAAG,aAAa,CAAA;YAC9B,CAAC;QACH,CAAC;aAAM,CAAC;YACN,YAAY,CAAC,IAAI,CAAC,KAAK,CAAC,CAAA;YACxB,YAAY,GAAG,aAAa,CAAA;YAC5B,YAAY,GAAG,aAAa,CAAA;QAC9B,CAAC;IACH,CAAC;IAED,OAAO,YAAY,CAAA;AACrB,CAAC;AAED,4EAA4E;AAC5E,6EAA6E;AAC7E,kEAAkE;AAClE,MAAM,UAAU,cAAc,CAAC,KAAiB,EAAE,MAAc;IAC9D,MAAM,GAAG,GACP,KAAK,CAAC,MAAM,CAAE;QACd,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,IAAI,CAAC,CAAC;QACzB,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,IAAI,EAAE,CAAC;QAC1B,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,IAAI,EAAE,CAAC,CAAA;IAC5B,MAAM,IAAI,GACR,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE;QAClB,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,IAAI,CAAC,CAAC;QACzB,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,IAAI,EAAE,CAAC;QAC1B,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,IAAI,EAAE,CAAC,CAAA;IAC5B,OAAO;QACL,SAAS,EAAE,CAAC,IAAI,KAAK,CAAC,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC;KAChD,CAAA;AACH,CAAC;AAED,6EAA6E;AAC7E,kFAAkF;AAClF,kFAAkF;AAClF,iFAAiF;AACjF,kFAAkF;AAClF,iFAAiF;AACjF,yEAAyE;AACzE,MAAM,UAAU,cAAc,CAC5B,MAAe,EACf,kBAAqC,EAAE;IAEvC,4EAA4E;IAC5E,6EAA6E;IAC7E,0CAA0C;IAC1C,MAAM,UAAU,GAAG,IAAI,YAAY,CACjC,eAAe,CAAC,MAAM,GAAG,MAAM,CAAC,MAAM,GAAG,CAAC,CAC3C,CAAA;IACD,UAAU,CAAC,GAAG,CAAC,eAAe,CAAC,CAAA;IAC/B,IAAI,CAAC,GAAG,eAAe,CAAC,MAAM,CAAA;IAC9B,KAAK,MAAM,CAAC,IAAI,MAAM,EAAE,CAAC;QACvB,UAAU,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,aAAa,CAAA;QACtC,UAAU,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,aAAa,CAAA;IACxC,CAAC;IACD,UAAU,CAAC,IAAI,EAAE,CAAA;IAEjB,KAAK,MAAM,CAAC,IAAI,MAAM,EAAE,CAAC;QACvB,MAAM,GAAG,GAAG,CAAC,CAAC,IAAI,CAAC,aAAa,CAAA;QAChC,2CAA2C;QAC3C,IAAI,EAAE,GAAG,CAAC,CAAA;QACV,IAAI,EAAE,GAAG,UAAU,CAAC,MAAM,CAAA;QAC1B,OAAO,EAAE,GAAG,EAAE,EAAE,CAAC;YACf,MAAM,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,CAAA;YAC3B,IAAI,UAAU,CAAC,GAAG,CAAE,GAAG,GAAG,EAAE,CAAC;gBAC3B,EAAE,GAAG,GAAG,CAAA;YACV,CAAC;iBAAM,CAAC;gBACN,EAAE,GAAG,GAAG,GAAG,CAAC,CAAA;YACd,CAAC;QACH,CAAC;QACD,IAAI,EAAE,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC;YAC3B,CAAC,CAAC,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,WAAW,EAAE,UAAU,CAAC,EAAE,CAAE,CAAC,CAAA;QAC1D,CAAC;IACH,CAAC;AACH,CAAC;AAED,6EAA6E;AAC7E,4EAA4E;AAC5E,eAAe;AACf,MAAM,UAAU,YAAY,CAC1B,KAAiB,EACjB,KAAa,EACb,YAAmC;IAEnC,IAAI,KAAK,GAAG,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC;QAC7B,OAAO,SAAS,CAAA;IAClB,CAAC;IACD,MAAM,QAAQ,GAAG,IAAI,QAAQ,CAC3B,KAAK,CAAC,MAAM,EACZ,KAAK,CAAC,UAAU,EAChB,KAAK,CAAC,UAAU,CACjB,CAAA;IACD,MAAM,IAAI,GAAG,QAAQ,CAAC,QAAQ,CAAC,KAAK,EAAE,IAAI,CAAC,CAAA;IAC3C,kEAAkE;IAClE,uDAAuD;IACvD,MAAM,UAAU,GAA2B,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC,CAAA;IAC9D,MAAM,UAAU,GAA0C,EAAE,CAAA;IAC5D,MAAM,OAAO,GAAG,IAAI,WAAW,CAAC,MAAM,CAAC,CAAA;IAEvC,IAAI,CAAC,GAAG,KAAK,GAAG,CAAC,CAAA;IACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;QAC9B,IAAI,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC;YACzB,OAAO,SAAS,CAAA;QAClB,CAAC;QACD,MAAM,KAAK,GAAG,QAAQ,CAAC,QAAQ,CAAC,CAAC,EAAE,IAAI,CAAC,CAAA;QACxC,IAAI,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC;YACjC,OAAO,SAAS,CAAA;QAClB,CAAC;QACD,MAAM,OAAO,GAAG,YAAY,CAC1B,OAAO,CAAC,MAAM,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CACzD,CAAA;QACD,MAAM,IAAI,GAAG,QAAQ,CAAC,QAAQ,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,EAAE,IAAI,CAAC,CAAA;QACnD,UAAU,CAAC,OAAO,CAAC,GAAG,CAAC,CAAA;QACvB,UAAU,CAAC,IAAI,CAAC,EAAE,OAAO,EAAE,MAAM,EAAE,IAAI,EAAE,CAAC,CAAA;QAC1C,CAAC,IAAI,CAAC,GAAG,KAAK,CAAA;IAChB,CAAC;IACD,6EAA6E;IAC7E,OAAO,EAAE,UAAU,EAAE,UAAU,EAAE,GAAG,EAAE,CAAC,EAAE,CAAA;AAC3C,CAAC;AAED,2EAA2E;AAC3E,+EAA+E;AAC/E,4CAA4C;AAC5C;;;;;;;;;;;;;;GAcG;AACH,MAAM,UAAU,gBAAgB,CAC9B,KAAiB,EACjB,MAAc,EACd,KAAa,EACb,OAAkC;IAElC,IAAI,QAAQ,GAAG,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC,aAAa,CAAC,CAAC,CAAC,QAAQ,CAAA;IACzD,IAAI,OAAO,GAAG,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC,YAAY,CAAC,CAAC,CAAC,CAAC,CAAA;IAChD,IAAI,KAAK,GAAG,KAAK,CAAA;IACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,MAAM,CAAC,GAAG,MAAM,GAAG,CAAC,GAAG,CAAC,CAAA;QACxB,MAAM,KAAK,GACT,KAAK,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,aAAa;YAC7B,KAAK,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,WAAW;YAC3B,KAAK,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,SAAS;YACzB,KAAK,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,OAAO;YACvB,KAAK,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,KAAK;YACrB,KAAK,CAAC,CAAC,GAAG,CAAC,CAAE,CAAA;QACf,MAAM,IAAI,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAE,IAAI,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAE,CAAA;QAC7C,IACE,CAAC,KAAK,KAAK,CAAC,IAAI,IAAI,KAAK,CAAC,CAAC;YAC3B,CAAC,KAAK,GAAG,QAAQ,IAAI,CAAC,KAAK,KAAK,QAAQ,IAAI,IAAI,GAAG,OAAO,CAAC,CAAC,EAC5D,CAAC;YACD,QAAQ,GAAG,KAAK,CAAA;YAChB,OAAO,GAAG,IAAI,CAAA;YACd,KAAK,GAAG,IAAI,CAAA;QACd,CAAC;IACH,CAAC;IACD,OAAO,KAAK,CAAC,CAAC,CAAC,IAAI,aAAa,CAAC,QAAQ,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,OAAO,CAAA;AAC/D,CAAC;AAED,4FAA4F;AAC5F,MAAM,UAAU,cAAc,CAC5B,UAAsB,EACtB,eAAwC,CAAC,CAAC,EAAE,CAAC,CAAC;IAE9C,MAAM,OAAO,GAAG,IAAI,WAAW,EAAE,CAAA;IACjC,IAAI,SAAS,GAAG,CAAC,CAAA;IACjB,IAAI,aAAa,GAAG,CAAC,CAAA;IACrB,MAAM,WAAW,GAAa,EAAE,CAAA;IAChC,MAAM,WAAW,GAA2B,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC,CAAA;IAC/D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3C,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,CAAC;YACnB,IAAI,aAAa,GAAG,CAAC,EAAE,CAAC;gBACtB,MAAM,OAAO,GAAG,YAAY,CAC1B,OAAO,CAAC,MAAM,CAAC,UAAU,CAAC,QAAQ,CAAC,aAAa,EAAE,CAAC,CAAC,CAAC,CACtD,CAAA;gBACD,WAAW,CAAC,SAAS,CAAC,GAAG,OAAO,CAAA;gBAChC,WAAW,CAAC,OAAO,CAAC,GAAG,SAAS,CAAA;YAClC,CAAC;YACD,aAAa,GAAG,CAAC,GAAG,CAAC,CAAA;YACrB,SAAS,EAAE,CAAA;QACb,CAAC;IACH,CAAC;IACD,OAAO,EAAE,WAAW,EAAE,WAAW,EAAE,CAAA;AACrC,CAAC;AAED,MAAM,UAAU,gBAAgB,CAAC,IAAkB;IACjD,IAAI,WAAW,GAAG,CAAC,CAAA;IACnB,KAAK,MAAM,KAAK,IAAI,IAAI,EAAE,CAAC;QACzB,WAAW,IAAI,KAAK,CAAC,MAAM,CAAA;IAC7B,CAAC;IACD,MAAM,WAAW,GAAG,IAAI,UAAU,CAAC,WAAW,CAAC,CAAA;IAC/C,IAAI,MAAM,GAAG,CAAC,CAAA;IACd,KAAK,MAAM,KAAK,IAAI,IAAI,EAAE,CAAC;QACzB,WAAW,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,CAAA;QAC9B,MAAM,IAAI,KAAK,CAAC,MAAM,CAAA;IACxB,CAAC;IACD,OAAO,WAAW,CAAA;AACpB,CAAC;AAQD,0EAA0E;AAC1E,4EAA4E;AAC5E,6EAA6E;AAC7E,iDAAiD;AACjD,MAAM,UAAU,aAAa,CAC3B,OAAY,EACZ,KAAa,EACb,GAAW,EACX,GAAW,EACX,MAAW,EAAE;IAEb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/C,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAE,CAAA;QACrB,IAAI,CAAC,CAAC,MAAM,KAAK,KAAK,EAAE,CAAC;YACvB,IAAI,CAAC,CAAC,KAAK,IAAI,GAAG,EAAE,CAAC;gBACnB,MAAK;YACP,CAAC;iBAAM,IAAI,CAAC,CAAC,GAAG,IAAI,GAAG,EAAE,CAAC;gBACxB,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAA;YACb,CAAC;QACH,CAAC;aAAM,IAAI,CAAC,CAAC,MAAM,GAAG,KAAK,EAAE,CAAC;YAC5B,MAAK;QACP,CAAC;IACH,CAAC;IACD,OAAO,GAAG,CAAA;AACZ,CAAC"}

package/esm/virtualOffset.d.ts

@@ -1,14 +1,19 @@
-export interface Offset {
+/**
+ * The two coordinates a virtual offset packs, without the `blockPos:dataPos`
+ * string form. What every consumer that only compares positions needs — the
+ * BAI linear index stores these as raw numbers and materializes no object.
+ */
+export interface OffsetCoords {
     blockPosition: number;
     dataPosition: number;
+}
+export interface Offset extends OffsetCoords {
     toString(): string;
-    compareTo(arg: Offset): number;
 }
 export declare class VirtualOffset {
     blockPosition: number;
     dataPosition: number;
     constructor(blockPosition: number, dataPosition: number);
     toString(): string;
-    compareTo(b: VirtualOffset): number;
 }
 export declare function fromBytes(bytes: Uint8Array, offset?: number): VirtualOffset;

package/esm/virtualOffset.js

@@ -8,9 +8,6 @@ export class VirtualOffset {
     toString() {
         return `${this.blockPosition}:${this.dataPosition}`;
     }
-    compareTo(b) {
-        return (this.blockPosition - b.blockPosition || this.dataPosition - b.dataPosition);
-    }
 }
 export function fromBytes(bytes, offset = 0) {
     return new VirtualOffset(bytes[offset + 7] * 0x10000000000 +

package/esm/virtualOffset.js.map

@@ -1 +1 @@
-{"version":3,"file":"virtualOffset.js","sourceRoot":"","sources":["../src/virtualOffset.ts"],"names":[],"mappings":"AAOA,MAAM,OAAO,aAAa;IACjB,aAAa,CAAQ;IACrB,YAAY,CAAQ;IAC3B,YAAY,aAAqB,EAAE,YAAoB;QACrD,IAAI,CAAC,aAAa,GAAG,aAAa,CAAA,CAAC,wCAAwC;QAC3E,IAAI,CAAC,YAAY,GAAG,YAAY,CAAA,CAAC,sCAAsC;IACzE,CAAC;IAED,QAAQ;QACN,OAAO,GAAG,IAAI,CAAC,aAAa,IAAI,IAAI,CAAC,YAAY,EAAE,CAAA;IACrD,CAAC;IAED,SAAS,CAAC,CAAgB;QACxB,OAAO,CACL,IAAI,CAAC,aAAa,GAAG,CAAC,CAAC,aAAa,IAAI,IAAI,CAAC,YAAY,GAAG,CAAC,CAAC,YAAY,CAC3E,CAAA;IACH,CAAC;CACF;AACD,MAAM,UAAU,SAAS,CAAC,KAAiB,EAAE,MAAM,GAAG,CAAC;IACrD,OAAO,IAAI,aAAa,CACtB,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,GAAG,aAAa;QAChC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,GAAG,WAAW;QAChC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,GAAG,SAAS;QAC9B,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,GAAG,OAAO;QAC5B,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,GAAG,KAAK;QAC1B,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,EACpB,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,IAAI,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAE,CAC3C,CAAA;AACH,CAAC"}
+{"version":3,"file":"virtualOffset.js","sourceRoot":"","sources":["../src/virtualOffset.ts"],"names":[],"mappings":"AAcA,MAAM,OAAO,aAAa;IACjB,aAAa,CAAQ;IACrB,YAAY,CAAQ;IAC3B,YAAY,aAAqB,EAAE,YAAoB;QACrD,IAAI,CAAC,aAAa,GAAG,aAAa,CAAA,CAAC,wCAAwC;QAC3E,IAAI,CAAC,YAAY,GAAG,YAAY,CAAA,CAAC,sCAAsC;IACzE,CAAC;IAED,QAAQ;QACN,OAAO,GAAG,IAAI,CAAC,aAAa,IAAI,IAAI,CAAC,YAAY,EAAE,CAAA;IACrD,CAAC;CACF;AACD,MAAM,UAAU,SAAS,CAAC,KAAiB,EAAE,MAAM,GAAG,CAAC;IACrD,OAAO,IAAI,aAAa,CACtB,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,GAAG,aAAa;QAChC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,GAAG,WAAW;QAChC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,GAAG,SAAS;QAC9B,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,GAAG,OAAO;QAC5B,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,GAAG,KAAK;QAC1B,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,EACpB,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,IAAI,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAE,CAC3C,CAAA;AACH,CAAC"}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@gmod/bam",
-  "version": "7.6.2",
+  "version": "7.7.0",
   "packageManager": "pnpm@11.15.1",
   "description": "Parser for BAM and BAM index (bai) files",
   "license": "MIT",

package/src/bai.ts

@@ -7,8 +7,17 @@ import type { ParsedIndexBase, RefIndex } from './indexFile.ts'
 import type { BaseOpts } from './util.ts'
 import type { VirtualOffset } from './virtualOffset.ts'
 
+// The linear index as two parallel Float64Arrays rather than an array of
+// VirtualOffset objects. A human-sized reference has one entry per 16kb window
+// — ~15k for chr1 — and an object apiece costs roughly an order of magnitude
+// more memory than the packed form, retained for as long as memoizeByRefId
+// holds the reference. Both consumers want the raw numbers anyway: indexCov and
+// clampChunkEnds read blockPosition only, and getLowestChunk builds the one
+// VirtualOffset a query actually needs. Same shape as bgzf-filehandle's
+// GziIndex, for the same reason.
 interface BaiRefIndex extends RefIndex {
-  linearIndex: VirtualOffset[]
+  linearBlockPositions: Float64Array
+  linearDataPositions: Float64Array
 }
 
 interface BaiParsed extends ParsedIndexBase<BaiRefIndex> {
@@ -142,19 +151,27 @@ export default class BAI extends IndexFile<BaiParsed> {
 
       const linearCount = dataView.getInt32(curr, true)
       curr += 4
-      const linearIndex = new Array<VirtualOffset>(linearCount)
+      const linearBlockPositions = new Float64Array(linearCount)
+      const linearDataPositions = new Float64Array(linearCount)
       for (let j = 0; j < linearCount; j++) {
-        linearIndex[j] = fromBytes(bytes, curr)
+        // a virtual offset is a 48-bit block position in the high bytes and a
+        // 16-bit data position in the low two
+        linearBlockPositions[j] =
+          bytes[curr + 7]! * 0x10000000000 +
+          bytes[curr + 6]! * 0x100000000 +
+          bytes[curr + 5]! * 0x1000000 +
+          bytes[curr + 4]! * 0x10000 +
+          bytes[curr + 3]! * 0x100 +
+          bytes[curr + 2]!
+        linearDataPositions[j] = (bytes[curr + 1]! << 8) | bytes[curr]!
         curr += 8
       }
 
-      clampChunkEnds(
-        Object.values(binIndex).flat(),
-        linearIndex.map(v => v.blockPosition),
-      )
+      clampChunkEnds(Object.values(binIndex).flat(), linearBlockPositions)
       return {
         binIndex,
-        linearIndex,
+        linearBlockPositions,
+        linearDataPositions,
         stats,
       }
     }
@@ -162,7 +179,6 @@ export default class BAI extends IndexFile<BaiParsed> {
     return {
       bai: true,
       firstDataLine,
-      maxBlockSize: 1 << 16,
       indices: memoizeByRefId(getIndices),
       refCount,
     }
@@ -182,17 +198,18 @@ export default class BAI extends IndexFile<BaiParsed> {
     if (!seqIdx) {
       return []
     }
-    const { linearIndex, stats } = seqIdx
-    if (linearIndex.length === 0) {
+    const { linearBlockPositions, stats } = seqIdx
+    const nintv = linearBlockPositions.length
+    if (nintv === 0) {
       return []
     }
-    const e = end === undefined ? (linearIndex.length - 1) * v : roundUp(end, v)
+    const e = end === undefined ? (nintv - 1) * v : roundUp(end, v)
     const s = start === undefined ? 0 : roundDown(start, v)
     const depths: IndexCovEntry[] = range
       ? new Array((e - s) / v)
-      : new Array(linearIndex.length - 1)
-    const totalSize = linearIndex[linearIndex.length - 1]!.blockPosition
-    if (e > (linearIndex.length - 1) * v) {
+      : new Array(nintv - 1)
+    const totalSize = linearBlockPositions[nintv - 1]!
+    if (e > (nintv - 1) * v) {
       throw new Error('query outside of range of linear index')
     }
     // Scale the block-delta into a read count as we go, rather than building the
@@ -200,9 +217,9 @@ export default class BAI extends IndexFile<BaiParsed> {
     // multiply-then-divide order: hoisting lineCount/totalSize into a factor
     // reassociates the arithmetic and shifts scores by an ulp.
     const lineCount = stats?.lineCount ?? 0
-    let currentPos = linearIndex[s / v]!.blockPosition
+    let currentPos = linearBlockPositions[s / v]!
     for (let i = s / v, j = 0; i < e / v; i++, j++) {
-      const nextPos = linearIndex[i + 1]!.blockPosition
+      const nextPos = linearBlockPositions[i + 1]!
       depths[j] = {
         score: ((nextPos - currentPos) * lineCount) / totalSize,
         start: i * v,
@@ -221,8 +238,13 @@ export default class BAI extends IndexFile<BaiParsed> {
   // contain alignments in the region. Linear index entries are monotonically
   // non-decreasing, so the first entry at minLin is the minimum.
   protected getLowestChunk(refIndex: BaiRefIndex, min: number) {
-    const { linearIndex } = refIndex
-    const nintv = linearIndex.length
-    return linearIndex[Math.min(min >> BAI_LINEAR_SHIFT, nintv - 1)]
+    const { linearBlockPositions, linearDataPositions } = refIndex
+    const i = Math.min(min >> BAI_LINEAR_SHIFT, linearBlockPositions.length - 1)
+    return i < 0
+      ? undefined
+      : {
+          blockPosition: linearBlockPositions[i]!,
+          dataPosition: linearDataPositions[i]!,
+        }
   }
 }

package/src/bamFile.ts

@@ -232,7 +232,14 @@ export default class BamFile<T extends BamRecordLike = BAMFeature> {
 
   async getHeaderPre(opts: BaseOpts = {}) {
     if (!this.index) {
-      return undefined
+      // Only reachable via `new BamFile({htsget: true})` used directly rather
+      // than through HtsgetFile, which overrides this. Throwing rather than
+      // returning undefined keeps `| undefined` off getHeader()'s public type,
+      // and makes that misuse loud instead of silently answering every query
+      // with zero records (getSeqId would find no chrToIndex).
+      throw new Error(
+        'no index to read a header from: use HtsgetFile for htsget sources',
+      )
     }
     const indexData = await this.index.parse(opts)
 
@@ -512,13 +519,26 @@ export default class BamFile<T extends BamRecordLike = BAMFeature> {
       const chunks = res[i]!
       for (let j = 0, jl = chunks.length; j < jl; j++) {
         const m = chunks[j]!
-        map.set(m.toString(), m)
+        // Key on the virtual-offset span — the same key _cachedChunkFeatures
+        // uses. Chunk.toString() also folds in `bin` and fetchedSize(), which
+        // keeps two chunks covering an identical span apart here even though
+        // the cache below collapses them, so their records came back twice.
+        map.set(chunkCacheKey(m), m)
       }
     }
-
-    const mateFeatLists = await Promise.all(
-      [...map.values()].map(async c => {
-        const features = await this._cachedChunkFeatures(c, opts)
+    const mateChunks = [...map.values()]
+
+    // Bounded for the reason ADR 0008 bounds the main query path: a viewAsPairs
+    // query over a busy region resolves to many distinct mate chunks, and an
+    // unbounded fan-out inflates every one of them at once. Same inline pool
+    // shape as _fetchChunkFeatures — ADR 0009 measured extracting a shared
+    // mapConcurrent helper at ~11% on the hot path, so it stays inline.
+    const mateFeatLists = new Array<T[]>(mateChunks.length)
+    let next = 0
+    const readNext = async () => {
+      while (next < mateChunks.length) {
+        const ci = next++
+        const features = await this._cachedChunkFeatures(mateChunks[ci]!, opts)
         const mateRecs = [] as T[]
         for (let i = 0, l = features.length; i < l; i++) {
           const feature = features[i]!
@@ -529,9 +549,11 @@ export default class BamFile<T extends BamRecordLike = BAMFeature> {
             mateRecs.push(feature)
           }
         }
-        return mateRecs
-      }),
-    )
+        mateFeatLists[ci] = mateRecs
+      }
+    }
+    const workers = Math.min(MAX_CONCURRENT_CHUNK_READS, mateChunks.length)
+    await Promise.all(Array.from({ length: workers }, () => readNext()))
     return mateFeatLists.flat()
   }
 
@@ -559,10 +581,15 @@ export default class BamFile<T extends BamRecordLike = BAMFeature> {
     }
   }
 
+  // cpositions/dpositions are `ArrayLike`, not `number[]`, because the wasm
+  // decompressor produces them as Float64Arrays — @gmod/bgzf-filehandle spreads
+  // them into plain arrays on the way out, and this signature is what lets that
+  // copy be dropped there without a change here. Only indexed reads and
+  // `.length` are used below either way.
   readBamFeatures(
     ba: Uint8Array,
-    cpositions: number[],
-    dpositions: number[],
+    cpositions: ArrayLike<number>,
+    dpositions: ArrayLike<number>,
     chunk: Chunk,
   ) {
     let blockStart = 0

package/src/chunk.ts

@@ -29,14 +29,6 @@ export default class Chunk {
     }, fetchedSize ${this.fetchedSize()})`
   }
 
-  compareTo(b: Chunk) {
-    return (
-      this.minv.compareTo(b.minv) ||
-      this.maxv.compareTo(b.maxv) ||
-      this.bin - b.bin
-    )
-  }
-
   fetchedSize() {
     return this.endPosition - this.minv.blockPosition
   }

package/src/csi.ts

@@ -83,16 +83,12 @@ export default class CSI extends IndexFile {
       bytes.byteOffset,
       bytes.byteLength,
     )
-    let csiVersion
+    // CSI\1 and CSI\2 are parsed identically here — v2 only adds fields this
+    // reader doesn't consume — so the version is validated, not retained.
     const magic = dataView.getUint32(0, true)
-
-    if (magic === CSI1_MAGIC) {
-      csiVersion = 1
-    } else if (magic === CSI2_MAGIC) {
-      csiVersion = 2
-    } else {
-      throw new Error(`Not a CSI file ${magic}`)
+    if (magic !== CSI1_MAGIC && magic !== CSI2_MAGIC) {
       // TODO: do we need to support big-endian CSI files?
+      throw new Error(`Not a CSI file ${magic}`)
     }
 
     this.minShift = dataView.getInt32(4, true)
@@ -170,12 +166,10 @@ export default class CSI extends IndexFile {
     }
 
     return {
-      csiVersion,
       firstDataLine,
       indices: memoizeByRefId(getIndices),
       refCount,
       csi: true,
-      maxBlockSize: 1 << 16,
       ...aux,
     }
   }

package/src/indexFile.ts

@@ -4,7 +4,7 @@ import { optimizeChunks } from './util.ts'
 
 import type Chunk from './chunk.ts'
 import type { BaseOpts } from './util.ts'
-import type { Offset, VirtualOffset } from './virtualOffset.ts'
+import type { OffsetCoords, VirtualOffset } from './virtualOffset.ts'
 import type { GenericFilehandle } from 'generic-filehandle2'
 
 export interface Region {
@@ -21,7 +21,6 @@ export interface RefIndex {
 export interface ParsedIndexBase<R extends RefIndex = RefIndex> {
   firstDataLine: VirtualOffset | undefined
   refCount: number
-  maxBlockSize: number
   indices: (refId: number) => R | undefined
 }
 
@@ -84,7 +83,7 @@ export default abstract class IndexFile<
   protected abstract getLowestChunk(
     refIndex: RefIndex,
     min: number,
-  ): Offset | undefined
+  ): OffsetCoords | undefined
 
   async blocksForRange(
     refId: number,

package/src/sam.ts

@@ -8,9 +8,14 @@ export function parseHeaderText(text: string) {
         tag: tag.slice(1),
         data: fields.map(f => {
           const r = f.indexOf(':')
-          const fieldTag = f.slice(0, r)
-          const value = f.slice(r + 1)
-          return { tag: fieldTag, value }
+          // A field with no colon is not a TAG:VALUE pair — @CO lines are free
+          // text, one comment per line (`samtools view -H c2#pad.3.0.cram`).
+          // Without this branch `slice(0, -1)` silently drops the comment's last
+          // character into the tag. Same shape @gmod/cram's parseHeaderText
+          // returns, so both feed a consumer's header parser identically.
+          return r === -1
+            ? { tag: f, value: '' }
+            : { tag: f.slice(0, r), value: f.slice(r + 1) }
         }),
       })
     }

package/src/util.ts

@@ -1,8 +1,7 @@
 import Chunk from './chunk.ts'
-import { longFromBytesToUnsigned } from './long.ts'
 import { VirtualOffset } from './virtualOffset.ts'
 
-import type { Offset } from './virtualOffset.ts'
+import type { OffsetCoords } from './virtualOffset.ts'
 
 export interface BamOpts {
   viewAsPairs?: boolean
@@ -30,7 +29,7 @@ export interface BaseOpts {
   onProgress?: (bytesDownloaded: number, totalBytes?: number) => void
 }
 
-export function optimizeChunks(chunks: Chunk[], lowest?: Offset) {
+export function optimizeChunks(chunks: Chunk[], lowest?: OffsetCoords) {
   const n = chunks.length
   if (n === 0) {
     return chunks
@@ -102,9 +101,22 @@ export function optimizeChunks(chunks: Chunk[], lowest?: Offset) {
   return mergedChunks
 }
 
+// The pseudo-bin's mapped-record count, a little-endian uint64. Read as two
+// 32-bit halves rather than via BigInt: exact up to Number.MAX_SAFE_INTEGER,
+// which is far past any real record count, and allocates nothing.
 export function parsePseudoBin(bytes: Uint8Array, offset: number) {
+  const low =
+    bytes[offset]! |
+    (bytes[offset + 1]! << 8) |
+    (bytes[offset + 2]! << 16) |
+    (bytes[offset + 3]! << 24)
+  const high =
+    bytes[offset + 4]! |
+    (bytes[offset + 5]! << 8) |
+    (bytes[offset + 6]! << 16) |
+    (bytes[offset + 7]! << 24)
   return {
-    lineCount: longFromBytesToUnsigned(bytes, offset),
+    lineCount: (high >>> 0) * 2 ** 32 + (low >>> 0),
   }
 }
 
@@ -117,13 +129,21 @@ export function parsePseudoBin(bytes: Uint8Array, offset: number) {
 // Shrinks both the byte estimate and the actual fetch with no extra I/O.
 export function clampChunkEnds(
   chunks: Chunk[],
-  extraBoundaries: number[] = [],
+  extraBoundaries: ArrayLike<number> = [],
 ) {
-  const boundaries = [...extraBoundaries]
+  // Float64Array rather than number[]: `extraBoundaries` is the linear index,
+  // which runs to tens of thousands of entries on a human-sized reference, and
+  // this is the only place they get copied.
+  const boundaries = new Float64Array(
+    extraBoundaries.length + chunks.length * 2,
+  )
+  boundaries.set(extraBoundaries)
+  let n = extraBoundaries.length
   for (const c of chunks) {
-    boundaries.push(c.minv.blockPosition, c.maxv.blockPosition)
+    boundaries[n++] = c.minv.blockPosition
+    boundaries[n++] = c.maxv.blockPosition
   }
-  boundaries.sort((a, b) => a - b)
+  boundaries.sort()
 
   for (const c of chunks) {
     const max = c.maxv.blockPosition

package/src/virtualOffset.ts

@@ -1,8 +1,15 @@
-export interface Offset {
+/**
+ * The two coordinates a virtual offset packs, without the `blockPos:dataPos`
+ * string form. What every consumer that only compares positions needs — the
+ * BAI linear index stores these as raw numbers and materializes no object.
+ */
+export interface OffsetCoords {
   blockPosition: number
   dataPosition: number
+}
+
+export interface Offset extends OffsetCoords {
   toString(): string
-  compareTo(arg: Offset): number
 }
 
 export class VirtualOffset {
@@ -16,12 +23,6 @@ export class VirtualOffset {
   toString() {
     return `${this.blockPosition}:${this.dataPosition}`
   }
-
-  compareTo(b: VirtualOffset) {
-    return (
-      this.blockPosition - b.blockPosition || this.dataPosition - b.dataPosition
-    )
-  }
 }
 export function fromBytes(bytes: Uint8Array, offset = 0) {
   return new VirtualOffset(

package/dist/long.d.ts

@@ -1 +0,0 @@
-export declare function longFromBytesToUnsigned(source: Uint8Array, i?: number): number;

package/dist/long.js

@@ -1,17 +0,0 @@
-"use strict";
-Object.defineProperty(exports, "__esModule", { value: true });
-exports.longFromBytesToUnsigned = longFromBytesToUnsigned;
-const TWO_PWR_32_DBL = 2 ** 32;
-// avoids dependency on long.js
-function longFromBytesToUnsigned(source, i = 0) {
-    const low = source[i] |
-        (source[i + 1] << 8) |
-        (source[i + 2] << 16) |
-        (source[i + 3] << 24);
-    const high = source[i + 4] |
-        (source[i + 5] << 8) |
-        (source[i + 6] << 16) |
-        (source[i + 7] << 24);
-    return (high >>> 0) * TWO_PWR_32_DBL + (low >>> 0);
-}
-//# sourceMappingURL=long.js.map

package/dist/long.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"long.js","sourceRoot":"","sources":["../src/long.ts"],"names":[],"mappings":";;AAGA,0DAYC;AAfD,MAAM,cAAc,GAAG,CAAC,IAAI,EAAE,CAAA;AAE9B,+BAA+B;AAC/B,SAAgB,uBAAuB,CAAC,MAAkB,EAAE,CAAC,GAAG,CAAC;IAC/D,MAAM,GAAG,GACP,MAAM,CAAC,CAAC,CAAE;QACV,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAE,IAAI,CAAC,CAAC;QACrB,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAE,IAAI,EAAE,CAAC;QACtB,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAE,IAAI,EAAE,CAAC,CAAA;IACxB,MAAM,IAAI,GACR,MAAM,CAAC,CAAC,GAAG,CAAC,CAAE;QACd,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAE,IAAI,CAAC,CAAC;QACrB,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAE,IAAI,EAAE,CAAC;QACtB,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAE,IAAI,EAAE,CAAC,CAAA;IACxB,OAAO,CAAC,IAAI,KAAK,CAAC,CAAC,GAAG,cAAc,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAA;AACpD,CAAC"}

package/esm/long.d.ts

@@ -1 +0,0 @@
-export declare function longFromBytesToUnsigned(source: Uint8Array, i?: number): number;

package/esm/long.js

@@ -1,14 +0,0 @@
-const TWO_PWR_32_DBL = 2 ** 32;
-// avoids dependency on long.js
-export function longFromBytesToUnsigned(source, i = 0) {
-    const low = source[i] |
-        (source[i + 1] << 8) |
-        (source[i + 2] << 16) |
-        (source[i + 3] << 24);
-    const high = source[i + 4] |
-        (source[i + 5] << 8) |
-        (source[i + 6] << 16) |
-        (source[i + 7] << 24);
-    return (high >>> 0) * TWO_PWR_32_DBL + (low >>> 0);
-}
-//# sourceMappingURL=long.js.map

package/esm/long.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"long.js","sourceRoot":"","sources":["../src/long.ts"],"names":[],"mappings":"AAAA,MAAM,cAAc,GAAG,CAAC,IAAI,EAAE,CAAA;AAE9B,+BAA+B;AAC/B,MAAM,UAAU,uBAAuB,CAAC,MAAkB,EAAE,CAAC,GAAG,CAAC;IAC/D,MAAM,GAAG,GACP,MAAM,CAAC,CAAC,CAAE;QACV,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAE,IAAI,CAAC,CAAC;QACrB,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAE,IAAI,EAAE,CAAC;QACtB,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAE,IAAI,EAAE,CAAC,CAAA;IACxB,MAAM,IAAI,GACR,MAAM,CAAC,CAAC,GAAG,CAAC,CAAE;QACd,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAE,IAAI,CAAC,CAAC;QACrB,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAE,IAAI,EAAE,CAAC;QACtB,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAE,IAAI,EAAE,CAAC,CAAA;IACxB,OAAO,CAAC,IAAI,KAAK,CAAC,CAAC,GAAG,cAAc,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAA;AACpD,CAAC"}

package/src/long.ts

@@ -1,16 +0,0 @@
-const TWO_PWR_32_DBL = 2 ** 32
-
-// avoids dependency on long.js
-export function longFromBytesToUnsigned(source: Uint8Array, i = 0) {
-  const low =
-    source[i]! |
-    (source[i + 1]! << 8) |
-    (source[i + 2]! << 16) |
-    (source[i + 3]! << 24)
-  const high =
-    source[i + 4]! |
-    (source[i + 5]! << 8) |
-    (source[i + 6]! << 16) |
-    (source[i + 7]! << 24)
-  return (high >>> 0) * TWO_PWR_32_DBL + (low >>> 0)
-}