@gmod/cram 4.0.8 → 4.0.10

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (48) hide show
  1. package/dist/cram-bundle.js +1 -1
  2. package/dist/cramFile/file.js +2 -5
  3. package/dist/cramFile/file.js.map +1 -1
  4. package/dist/htscodecs/arith_gen.js +2 -2
  5. package/dist/htscodecs/arith_gen.js.map +1 -1
  6. package/dist/seek-bzip/bitreader.d.ts +27 -0
  7. package/dist/seek-bzip/bitreader.js +101 -0
  8. package/dist/seek-bzip/bitreader.js.map +1 -0
  9. package/dist/seek-bzip/crc32.d.ts +0 -0
  10. package/dist/seek-bzip/crc32.js +113 -0
  11. package/dist/seek-bzip/crc32.js.map +1 -0
  12. package/dist/seek-bzip/index.d.ts +9 -0
  13. package/dist/seek-bzip/index.js +585 -0
  14. package/dist/seek-bzip/index.js.map +1 -0
  15. package/dist/seek-bzip/stream.d.ts +1 -0
  16. package/dist/seek-bzip/stream.js +41 -0
  17. package/dist/seek-bzip/stream.js.map +1 -0
  18. package/dist/seek-bzip/toHex.d.ts +1 -0
  19. package/dist/seek-bzip/toHex.js +36 -0
  20. package/dist/seek-bzip/toHex.js.map +1 -0
  21. package/esm/cramFile/file.js +2 -5
  22. package/esm/cramFile/file.js.map +1 -1
  23. package/esm/htscodecs/arith_gen.js +2 -2
  24. package/esm/htscodecs/arith_gen.js.map +1 -1
  25. package/esm/seek-bzip/bitreader.d.ts +27 -0
  26. package/esm/seek-bzip/bitreader.js +98 -0
  27. package/esm/seek-bzip/bitreader.js.map +1 -0
  28. package/esm/seek-bzip/crc32.d.ts +0 -0
  29. package/esm/seek-bzip/crc32.js +113 -0
  30. package/esm/seek-bzip/crc32.js.map +1 -0
  31. package/esm/seek-bzip/index.d.ts +9 -0
  32. package/esm/seek-bzip/index.js +579 -0
  33. package/esm/seek-bzip/index.js.map +1 -0
  34. package/esm/seek-bzip/stream.d.ts +1 -0
  35. package/esm/seek-bzip/stream.js +41 -0
  36. package/esm/seek-bzip/stream.js.map +1 -0
  37. package/esm/seek-bzip/toHex.d.ts +1 -0
  38. package/esm/seek-bzip/toHex.js +33 -0
  39. package/esm/seek-bzip/toHex.js.map +1 -0
  40. package/package.json +1 -6
  41. package/src/cramFile/file.ts +2 -7
  42. package/src/htscodecs/arith_gen.ts +2 -2
  43. package/src/seek-bzip/README.md +6 -0
  44. package/src/seek-bzip/bitreader.ts +113 -0
  45. package/src/seek-bzip/crc32.ts +116 -0
  46. package/src/seek-bzip/index.ts +615 -0
  47. package/src/seek-bzip/stream.ts +42 -0
  48. package/src/seek-bzip/toHex.ts +33 -0
@@ -0,0 +1,33 @@
1
+ // from https://www.xaymar.com/articles/2020/12/08/fastest-uint8array-to-hex-string-conversion-in-javascript/
2
+ // Pre-Init
3
+ const LUT_HEX_4b = [
4
+ '0',
5
+ '1',
6
+ '2',
7
+ '3',
8
+ '4',
9
+ '5',
10
+ '6',
11
+ '7',
12
+ '8',
13
+ '9',
14
+ 'A',
15
+ 'B',
16
+ 'C',
17
+ 'D',
18
+ 'E',
19
+ 'F',
20
+ ];
21
+ const LUT_HEX_8b = new Array(0x100);
22
+ for (let n = 0; n < 0x100; n++) {
23
+ LUT_HEX_8b[n] = `${LUT_HEX_4b[(n >>> 4) & 0xf]}${LUT_HEX_4b[n & 0xf]}`;
24
+ }
25
+ // End Pre-Init
26
+ export function toHex(buffer) {
27
+ let out = '';
28
+ for (let idx = 0, edx = buffer.length; idx < edx; idx++) {
29
+ out += LUT_HEX_8b[buffer[idx]];
30
+ }
31
+ return out;
32
+ }
33
+ //# sourceMappingURL=toHex.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"toHex.js","sourceRoot":"","sources":["../../src/seek-bzip/toHex.ts"],"names":[],"mappings":"AAAA,6GAA6G;AAE7G,WAAW;AACX,MAAM,UAAU,GAAG;IACjB,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;CACJ,CAAA;AACD,MAAM,UAAU,GAAG,IAAI,KAAK,CAAC,KAAK,CAAC,CAAA;AACnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;IAC/B,UAAU,CAAC,CAAC,CAAC,GAAG,GAAG,UAAU,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,UAAU,CAAC,CAAC,GAAG,GAAG,CAAC,EAAE,CAAA;AACxE,CAAC;AACD,eAAe;AACf,MAAM,UAAU,KAAK,CAAC,MAAkB;IACtC,IAAI,GAAG,GAAG,EAAE,CAAA;IACZ,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,MAAM,CAAC,MAAM,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC;QACxD,GAAG,IAAI,UAAU,CAAC,MAAM,CAAC,GAAG,CAAE,CAAC,CAAA;IACjC,CAAC;IACD,OAAO,GAAG,CAAA;AACZ,CAAC"}
package/package.json CHANGED
@@ -1,6 +1,6 @@
1
1
  {
2
2
  "name": "@gmod/cram",
3
- "version": "4.0.8",
3
+ "version": "4.0.10",
4
4
  "description": "read CRAM files with pure Javascript",
5
5
  "license": "MIT",
6
6
  "repository": "GMOD/cram-js",
@@ -40,16 +40,11 @@
40
40
  "biojs"
41
41
  ],
42
42
  "dependencies": {
43
- "@foxglove/wasm-bz2": "^0.1.1",
44
- "bz2": "^1.0.1",
45
- "bzip2": "^0.1.1",
46
- "bzip2-wasm": "^1.0.1",
47
43
  "crc": "^4.3.2",
48
44
  "generic-filehandle2": "^1.0.0",
49
45
  "md5": "^2.2.1",
50
46
  "pako": "^1.0.4",
51
47
  "quick-lru": "^4.0.1",
52
- "seek-bzip": "^2.0.0",
53
48
  "xz-decompress": "^0.2.1"
54
49
  },
55
50
  "devDependencies": {
@@ -1,11 +1,5 @@
1
- // import bzip2 from 'bzip2'
2
- // import BZip2 from 'bzip2-wasm'
3
- // import { decompress } from 'bz2'
4
-
5
-
6
1
  import crc32 from 'crc/calculators/crc32'
7
2
  import QuickLRU from 'quick-lru'
8
- import Bunzip from 'seek-bzip'
9
3
  import { XzReadableStream } from 'xz-decompress'
10
4
 
11
5
  import { CramMalformedError, CramUnimplementedError } from '../errors'
@@ -13,6 +7,7 @@ import * as htscodecs from '../htscodecs'
13
7
  import { open } from '../io'
14
8
  import ransuncompress from '../rans'
15
9
  import { parseHeaderText } from '../sam'
10
+ import { decode } from '../seek-bzip'
16
11
  import { unzip } from '../unzip'
17
12
  import CramContainer from './container'
18
13
  import CramRecord from './record'
@@ -293,7 +288,7 @@ export default class CramFile {
293
288
  }
294
289
  return ret
295
290
  } else if (compressionMethod === 'bzip2') {
296
- return Bunzip.decode(inputBuffer)
291
+ return decode(inputBuffer)
297
292
  } else if (compressionMethod === 'lzma') {
298
293
  const decompressedResponse = new Response(
299
294
  new XzReadableStream(bufferToStream(inputBuffer)),
@@ -34,7 +34,7 @@
34
34
  * OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
35
35
  */
36
36
 
37
- import Bunzip from 'seek-bzip'
37
+ import { decode } from '../seek-bzip'
38
38
 
39
39
  import RangeCoder from './arith_sh'
40
40
  import ByteModel from './byte_model'
@@ -154,7 +154,7 @@ export default class RangeCoderGen {
154
154
  // ----------------------------------------------------------------------
155
155
  // External codec
156
156
  decodeExt(stream, n_out) {
157
- return Bunzip.decode(stream.buf.slice(stream.pos))
157
+ return decode(stream.buf.slice(stream.pos))
158
158
  }
159
159
 
160
160
  // ----------------------------------------------------------------------
@@ -0,0 +1,6 @@
1
+ ## seek-bzip
2
+
3
+ Vendored copy of seek-bzip, with code modifications to remove Buffer usage by
4
+ Colin Diesh (c) 2025
5
+
6
+ Original LICENSE MIT https://github.com/cscott/seek-bzip/blob/master/LICENSE
@@ -0,0 +1,113 @@
1
+ /*
2
+ node-bzip - a pure-javascript Node.JS module for decoding bzip2 data
3
+
4
+ Copyright (C) 2012 Eli Skeggs
5
+
6
+ Permission is hereby granted, free of charge, to any person obtaining
7
+ a copy of this software and associated documentation files (the
8
+ "Software"), to deal in the Software without restriction, including
9
+ without limitation the rights to use, copy, modify, merge, publish,
10
+ distribute, sublicense, and/or sell copies of the Software, and to
11
+ permit persons to whom the Software is furnished to do so, subject to
12
+ the following conditions:
13
+
14
+ The above copyright notice and this permission notice shall be
15
+ included in all copies or substantial portions of the Software.
16
+
17
+ THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
18
+ EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
19
+ MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
20
+ NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
21
+ LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
22
+ OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
23
+ WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
24
+
25
+ Adapted from bzip2.js, copyright 2011 antimatter15 (antimatter15@gmail.com).
26
+
27
+ Based on micro-bunzip by Rob Landley (rob@landley.net).
28
+
29
+ Based on bzip2 decompression code by Julian R Seward (jseward@acm.org),
30
+ which also acknowledges contributions by Mike Burrows, David Wheeler,
31
+ Peter Fenwick, Alistair Moffat, Radford Neal, Ian H. Witten,
32
+ Robert Sedgewick, and Jon L. Bentley.
33
+ */
34
+
35
+ import { toHex } from './toHex'
36
+
37
+ const BITMASK = [0x00, 0x01, 0x03, 0x07, 0x0f, 0x1f, 0x3f, 0x7f, 0xff]
38
+
39
+ interface Stream {
40
+ readByte(): number
41
+ seek(pos: number): void
42
+ }
43
+
44
+ export default class BitReader {
45
+ private stream: Stream
46
+ private bitOffset: number
47
+ private curByte: number
48
+ private hasByte: boolean
49
+
50
+ constructor(stream: Stream) {
51
+ this.stream = stream
52
+ this.bitOffset = 0
53
+ this.curByte = 0
54
+ this.hasByte = false
55
+ }
56
+
57
+ private _ensureByte(): void {
58
+ if (!this.hasByte) {
59
+ this.curByte = this.stream.readByte()
60
+ this.hasByte = true
61
+ }
62
+ }
63
+
64
+ /**
65
+ * Reads bits from the buffer
66
+ * @param bits Number of bits to read
67
+ */
68
+ read(bits: number): number {
69
+ let result = 0
70
+ while (bits > 0) {
71
+ this._ensureByte()
72
+ const remaining = 8 - this.bitOffset
73
+ // if we're in a byte
74
+ if (bits >= remaining) {
75
+ result <<= remaining
76
+ result |= BITMASK[remaining]! & this.curByte
77
+ this.hasByte = false
78
+ this.bitOffset = 0
79
+ bits -= remaining
80
+ } else {
81
+ result <<= bits
82
+ const shift = remaining - bits
83
+ result |= (this.curByte & (BITMASK[bits]! << shift)) >> shift
84
+ this.bitOffset += bits
85
+ bits = 0
86
+ }
87
+ }
88
+ return result
89
+ }
90
+
91
+ /**
92
+ * Seek to an arbitrary point in the buffer (expressed in bits)
93
+ * @param pos Position in bits
94
+ */
95
+ seek(pos: number): void {
96
+ const n_bit = pos % 8
97
+ const n_byte = (pos - n_bit) / 8
98
+ this.bitOffset = n_bit
99
+ this.stream.seek(n_byte)
100
+ this.hasByte = false
101
+ }
102
+
103
+ /**
104
+ * Reads 6 bytes worth of data using the read method
105
+ */
106
+ pi(): string {
107
+ const buf = new Uint8Array(6)
108
+ for (let i = 0; i < buf.length; i++) {
109
+ buf[i] = this.read(8)
110
+ }
111
+ return toHex(buf)
112
+ }
113
+ }
@@ -0,0 +1,116 @@
1
+ // @ts-nocheck
2
+
3
+ /* CRC32, used in Bzip2 implementation.
4
+ * This is a port of CRC32.java from the jbzip2 implementation at
5
+ * https://code.google.com/p/jbzip2
6
+ * which is:
7
+ * Copyright (c) 2011 Matthew Francis
8
+ *
9
+ * Permission is hereby granted, free of charge, to any person
10
+ * obtaining a copy of this software and associated documentation
11
+ * files (the "Software"), to deal in the Software without
12
+ * restriction, including without limitation the rights to use,
13
+ * copy, modify, merge, publish, distribute, sublicense, and/or sell
14
+ * copies of the Software, and to permit persons to whom the
15
+ * Software is furnished to do so, subject to the following
16
+ * conditions:
17
+ *
18
+ * The above copyright notice and this permission notice shall be
19
+ * included in all copies or substantial portions of the Software.
20
+ *
21
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
22
+ * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES
23
+ * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
24
+ * NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
25
+ * HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY,
26
+ * WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
27
+ * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR
28
+ * OTHER DEALINGS IN THE SOFTWARE.
29
+ * This JavaScript implementation is:
30
+ * Copyright (c) 2013 C. Scott Ananian
31
+ * with the same licensing terms as Matthew Francis' original implementation.
32
+ */
33
+ module.exports = (function () {
34
+ /**
35
+ * A static CRC lookup table
36
+ */
37
+ const crc32Lookup = new Uint32Array([
38
+ 0x00000000, 0x04c11db7, 0x09823b6e, 0x0d4326d9, 0x130476dc, 0x17c56b6b,
39
+ 0x1a864db2, 0x1e475005, 0x2608edb8, 0x22c9f00f, 0x2f8ad6d6, 0x2b4bcb61,
40
+ 0x350c9b64, 0x31cd86d3, 0x3c8ea00a, 0x384fbdbd, 0x4c11db70, 0x48d0c6c7,
41
+ 0x4593e01e, 0x4152fda9, 0x5f15adac, 0x5bd4b01b, 0x569796c2, 0x52568b75,
42
+ 0x6a1936c8, 0x6ed82b7f, 0x639b0da6, 0x675a1011, 0x791d4014, 0x7ddc5da3,
43
+ 0x709f7b7a, 0x745e66cd, 0x9823b6e0, 0x9ce2ab57, 0x91a18d8e, 0x95609039,
44
+ 0x8b27c03c, 0x8fe6dd8b, 0x82a5fb52, 0x8664e6e5, 0xbe2b5b58, 0xbaea46ef,
45
+ 0xb7a96036, 0xb3687d81, 0xad2f2d84, 0xa9ee3033, 0xa4ad16ea, 0xa06c0b5d,
46
+ 0xd4326d90, 0xd0f37027, 0xddb056fe, 0xd9714b49, 0xc7361b4c, 0xc3f706fb,
47
+ 0xceb42022, 0xca753d95, 0xf23a8028, 0xf6fb9d9f, 0xfbb8bb46, 0xff79a6f1,
48
+ 0xe13ef6f4, 0xe5ffeb43, 0xe8bccd9a, 0xec7dd02d, 0x34867077, 0x30476dc0,
49
+ 0x3d044b19, 0x39c556ae, 0x278206ab, 0x23431b1c, 0x2e003dc5, 0x2ac12072,
50
+ 0x128e9dcf, 0x164f8078, 0x1b0ca6a1, 0x1fcdbb16, 0x018aeb13, 0x054bf6a4,
51
+ 0x0808d07d, 0x0cc9cdca, 0x7897ab07, 0x7c56b6b0, 0x71159069, 0x75d48dde,
52
+ 0x6b93dddb, 0x6f52c06c, 0x6211e6b5, 0x66d0fb02, 0x5e9f46bf, 0x5a5e5b08,
53
+ 0x571d7dd1, 0x53dc6066, 0x4d9b3063, 0x495a2dd4, 0x44190b0d, 0x40d816ba,
54
+ 0xaca5c697, 0xa864db20, 0xa527fdf9, 0xa1e6e04e, 0xbfa1b04b, 0xbb60adfc,
55
+ 0xb6238b25, 0xb2e29692, 0x8aad2b2f, 0x8e6c3698, 0x832f1041, 0x87ee0df6,
56
+ 0x99a95df3, 0x9d684044, 0x902b669d, 0x94ea7b2a, 0xe0b41de7, 0xe4750050,
57
+ 0xe9362689, 0xedf73b3e, 0xf3b06b3b, 0xf771768c, 0xfa325055, 0xfef34de2,
58
+ 0xc6bcf05f, 0xc27dede8, 0xcf3ecb31, 0xcbffd686, 0xd5b88683, 0xd1799b34,
59
+ 0xdc3abded, 0xd8fba05a, 0x690ce0ee, 0x6dcdfd59, 0x608edb80, 0x644fc637,
60
+ 0x7a089632, 0x7ec98b85, 0x738aad5c, 0x774bb0eb, 0x4f040d56, 0x4bc510e1,
61
+ 0x46863638, 0x42472b8f, 0x5c007b8a, 0x58c1663d, 0x558240e4, 0x51435d53,
62
+ 0x251d3b9e, 0x21dc2629, 0x2c9f00f0, 0x285e1d47, 0x36194d42, 0x32d850f5,
63
+ 0x3f9b762c, 0x3b5a6b9b, 0x0315d626, 0x07d4cb91, 0x0a97ed48, 0x0e56f0ff,
64
+ 0x1011a0fa, 0x14d0bd4d, 0x19939b94, 0x1d528623, 0xf12f560e, 0xf5ee4bb9,
65
+ 0xf8ad6d60, 0xfc6c70d7, 0xe22b20d2, 0xe6ea3d65, 0xeba91bbc, 0xef68060b,
66
+ 0xd727bbb6, 0xd3e6a601, 0xdea580d8, 0xda649d6f, 0xc423cd6a, 0xc0e2d0dd,
67
+ 0xcda1f604, 0xc960ebb3, 0xbd3e8d7e, 0xb9ff90c9, 0xb4bcb610, 0xb07daba7,
68
+ 0xae3afba2, 0xaafbe615, 0xa7b8c0cc, 0xa379dd7b, 0x9b3660c6, 0x9ff77d71,
69
+ 0x92b45ba8, 0x9675461f, 0x8832161a, 0x8cf30bad, 0x81b02d74, 0x857130c3,
70
+ 0x5d8a9099, 0x594b8d2e, 0x5408abf7, 0x50c9b640, 0x4e8ee645, 0x4a4ffbf2,
71
+ 0x470cdd2b, 0x43cdc09c, 0x7b827d21, 0x7f436096, 0x7200464f, 0x76c15bf8,
72
+ 0x68860bfd, 0x6c47164a, 0x61043093, 0x65c52d24, 0x119b4be9, 0x155a565e,
73
+ 0x18197087, 0x1cd86d30, 0x029f3d35, 0x065e2082, 0x0b1d065b, 0x0fdc1bec,
74
+ 0x3793a651, 0x3352bbe6, 0x3e119d3f, 0x3ad08088, 0x2497d08d, 0x2056cd3a,
75
+ 0x2d15ebe3, 0x29d4f654, 0xc5a92679, 0xc1683bce, 0xcc2b1d17, 0xc8ea00a0,
76
+ 0xd6ad50a5, 0xd26c4d12, 0xdf2f6bcb, 0xdbee767c, 0xe3a1cbc1, 0xe760d676,
77
+ 0xea23f0af, 0xeee2ed18, 0xf0a5bd1d, 0xf464a0aa, 0xf9278673, 0xfde69bc4,
78
+ 0x89b8fd09, 0x8d79e0be, 0x803ac667, 0x84fbdbd0, 0x9abc8bd5, 0x9e7d9662,
79
+ 0x933eb0bb, 0x97ffad0c, 0xafb010b1, 0xab710d06, 0xa6322bdf, 0xa2f33668,
80
+ 0xbcb4666d, 0xb8757bda, 0xb5365d03, 0xb1f740b4,
81
+ ])
82
+
83
+ const CRC32 = function () {
84
+ /**
85
+ * The current CRC
86
+ */
87
+ let crc = 0xffffffff
88
+
89
+ /**
90
+ * @return The current CRC
91
+ */
92
+ this.getCRC = function () {
93
+ return ~crc >>> 0 // return an unsigned value
94
+ }
95
+
96
+ /**
97
+ * Update the CRC with a single byte
98
+ * @param value The value to update the CRC with
99
+ */
100
+ this.updateCRC = function (value) {
101
+ crc = (crc << 8) ^ crc32Lookup[((crc >>> 24) ^ value) & 0xff]
102
+ }
103
+
104
+ /**
105
+ * Update the CRC with a sequence of identical bytes
106
+ * @param value The value to update the CRC with
107
+ * @param count The number of bytes
108
+ */
109
+ this.updateCRCRun = function (value, count) {
110
+ while (count-- > 0) {
111
+ crc = (crc << 8) ^ crc32Lookup[((crc >>> 24) ^ value) & 0xff]
112
+ }
113
+ }
114
+ }
115
+ return CRC32
116
+ })()