@gmod/gbz-base 0.0.1 → 1.0.0

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 (56) hide show
  1. package/README.md +3 -3
  2. package/bin/query.js +3 -1
  3. package/dist/cli.js +26 -7
  4. package/dist/cli.js.map +1 -0
  5. package/dist/db.d.ts +5 -5
  6. package/dist/db.js +35 -14
  7. package/dist/db.js.map +1 -0
  8. package/dist/filehandle.js +1 -0
  9. package/dist/filehandle.js.map +1 -0
  10. package/dist/gbwt/bytecode.d.ts +1 -1
  11. package/dist/gbwt/bytecode.js +7 -3
  12. package/dist/gbwt/bytecode.js.map +1 -0
  13. package/dist/gbwt/node.js +7 -2
  14. package/dist/gbwt/node.js.map +1 -0
  15. package/dist/gbwt/record.js +8 -2
  16. package/dist/gbwt/record.js.map +1 -0
  17. package/dist/gbwt/sequence.js +8 -1
  18. package/dist/gbwt/sequence.js.map +1 -0
  19. package/dist/index.d.ts +2 -2
  20. package/dist/index.js +3 -2
  21. package/dist/index.js.map +1 -0
  22. package/dist/lcs.js +32 -10
  23. package/dist/lcs.js.map +1 -0
  24. package/dist/query.d.ts +1 -1
  25. package/dist/query.js +1 -0
  26. package/dist/query.js.map +1 -0
  27. package/dist/sqlite/btree.d.ts +1 -1
  28. package/dist/sqlite/btree.js +17 -5
  29. package/dist/sqlite/btree.js.map +1 -0
  30. package/dist/sqlite/database.d.ts +1 -1
  31. package/dist/sqlite/database.js +12 -2
  32. package/dist/sqlite/database.js.map +1 -0
  33. package/dist/sqlite/pager.d.ts +3 -3
  34. package/dist/sqlite/pager.js +6 -4
  35. package/dist/sqlite/pager.js.map +1 -0
  36. package/dist/sqlite/record.js +3 -1
  37. package/dist/sqlite/record.js.map +1 -0
  38. package/dist/subgraph.d.ts +2 -2
  39. package/dist/subgraph.js +116 -29
  40. package/dist/subgraph.js.map +1 -0
  41. package/package.json +36 -11
  42. package/src/cli.ts +203 -0
  43. package/src/db.ts +322 -0
  44. package/src/filehandle.ts +4 -0
  45. package/src/gbwt/bytecode.ts +87 -0
  46. package/src/gbwt/node.ts +64 -0
  47. package/src/gbwt/record.ts +154 -0
  48. package/src/gbwt/sequence.ts +50 -0
  49. package/src/index.ts +31 -0
  50. package/src/lcs.ts +282 -0
  51. package/src/query.ts +79 -0
  52. package/src/sqlite/btree.ts +255 -0
  53. package/src/sqlite/database.ts +117 -0
  54. package/src/sqlite/pager.ts +78 -0
  55. package/src/sqlite/record.ts +100 -0
  56. package/src/subgraph.ts +1066 -0
package/README.md CHANGED
@@ -62,9 +62,9 @@ cd tools/haplotype-index && cargo build --release
62
62
  ./target/release/gbz-haplotype-index --interval 4096 --from-db graph.gbz.db
63
63
  ```
64
64
 
65
- The second form walks the paths through the database's own node records, so
66
- a database whose GBZ is no longer at hand can still be augmented; the two
67
- forms write identical tables.
65
+ The second form walks the paths through the database's own node records, so a
66
+ database whose GBZ is no longer at hand can still be augmented; the two forms
67
+ write identical tables.
68
68
 
69
69
  `HaplotypeSamples` holds one GBWT position every `--interval` bp along every
70
70
  path in both orientations, with the path handle and the forward coordinate of
package/bin/query.js CHANGED
@@ -2,6 +2,8 @@
2
2
  import { main } from '../dist/cli.js'
3
3
 
4
4
  main(process.argv.slice(2)).catch(error => {
5
- process.stderr.write(`${error instanceof Error ? error.message : String(error)}\n`)
5
+ process.stderr.write(
6
+ `${error instanceof Error ? error.message : String(error)}\n`,
7
+ )
6
8
  process.exit(1)
7
9
  })
package/dist/cli.js CHANGED
@@ -1,6 +1,6 @@
1
1
  import { LocalFile, RemoteFile } from 'generic-filehandle2';
2
2
  import { GBZBase, formatPathName } from "./db.js";
3
- import { subgraphAroundNodes, subgraphAtOffset, subgraphInInterval } from "./query.js";
3
+ import { subgraphAroundNodes, subgraphAtOffset, subgraphInInterval, } from "./query.js";
4
4
  const USAGE = `Usage: gbz-base-query [options] graph.gbz.db
5
5
 
6
6
  --sample STR sample name (default: generic path)
@@ -111,10 +111,20 @@ function parseArgs(argv) {
111
111
  }
112
112
  export async function main(argv) {
113
113
  const args = parseArgs(argv);
114
- const source = /^https?:\/\//.test(args.file) ? new RemoteFile(args.file) : new LocalFile(args.file);
114
+ const source = /^https?:\/\//.test(args.file)
115
+ ? new RemoteFile(args.file)
116
+ : new LocalFile(args.file);
115
117
  const db = await GBZBase.open(source, { blockSize: args.blockSize });
116
- const opts = { context: args.context, haplotypes: args.haplotypes, ...(args.limit === undefined ? {} : { limit: args.limit }) };
117
- const query = { contig: args.contig ?? '', haplotype: args.haplotype, ...(args.sample === undefined ? {} : { sample: args.sample }) };
118
+ const opts = {
119
+ context: args.context,
120
+ haplotypes: args.haplotypes,
121
+ ...(args.limit === undefined ? {} : { limit: args.limit }),
122
+ };
123
+ const query = {
124
+ contig: args.contig ?? '',
125
+ haplotype: args.haplotype,
126
+ ...(args.sample === undefined ? {} : { sample: args.sample }),
127
+ };
118
128
  const subgraph = args.nodes.length > 0
119
129
  ? await subgraphAroundNodes(db, args.nodes, opts)
120
130
  : args.interval
@@ -129,12 +139,21 @@ export async function main(argv) {
129
139
  await subgraph.identifyPaths();
130
140
  }
131
141
  const output = args.alignments
132
- ? subgraph.alignments().map(a => ({ ...a, name: a.name ? formatPathName(a.name, a.name.fragment) : undefined, start: undefined }))
133
- : subgraph.toJSON(args.cigar, { names: args.resolve ? 'resolved' : 'anonymous' });
142
+ ? subgraph.alignments().map(a => ({
143
+ ...a,
144
+ name: a.name && a.hapStart !== undefined && a.hapEnd !== undefined
145
+ ? formatPathName({ ...a.name, fragment: a.hapStart }, a.hapEnd)
146
+ : undefined,
147
+ start: undefined,
148
+ }))
149
+ : subgraph.toJSON(args.cigar, {
150
+ names: args.resolve ? 'resolved' : 'anonymous',
151
+ });
134
152
  process.stdout.write(`${JSON.stringify(output)}\n`);
135
153
  if (args.stats) {
136
154
  const { fetches, bytesFetched } = db.sqlite.pager;
137
- const { orderedAlignments, lcsAlignments, identificationSteps, identificationFetches } = subgraph.stats;
155
+ const { orderedAlignments, lcsAlignments, identificationSteps, identificationFetches, } = subgraph.stats;
138
156
  process.stderr.write(`Subgraph contains ${subgraph.nodeCount} nodes and ${subgraph.pathCount} paths; ${fetches} fetches, ${bytesFetched} bytes; ${orderedAlignments} ordered + ${lcsAlignments} lcs alignments; identification ${identificationSteps} steps, ${identificationFetches} lookups\n`);
139
157
  }
140
158
  }
159
+ //# sourceMappingURL=cli.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"cli.js","sourceRoot":"","sources":["../src/cli.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,SAAS,EAAE,UAAU,EAAE,MAAM,qBAAqB,CAAA;AAE3D,OAAO,EAAE,OAAO,EAAE,cAAc,EAAE,MAAM,SAAS,CAAA;AACjD,OAAO,EACL,mBAAmB,EACnB,gBAAgB,EAChB,kBAAkB,GACnB,MAAM,YAAY,CAAA;AAInB,MAAM,KAAK,GAAG;;;;;;;;;;;;;;;;CAgBb,CAAA;AAoBD,SAAS,SAAS,CAAC,IAAc;IAC/B,MAAM,IAAI,GAAS;QACjB,IAAI,EAAE,EAAE;QACR,SAAS,EAAE,CAAC;QACZ,KAAK,EAAE,EAAE;QACT,OAAO,EAAE,GAAG;QACZ,UAAU,EAAE,KAAK;QACjB,KAAK,EAAE,KAAK;QACZ,OAAO,EAAE,KAAK;QACd,UAAU,EAAE,KAAK;QACjB,SAAS,EAAE,KAAK;QAChB,KAAK,EAAE,KAAK;KACb,CAAA;IACD,MAAM,IAAI,GAAG,CAAC,CAAS,EAAE,EAAE;QACzB,MAAM,KAAK,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;QACzB,IAAI,KAAK,KAAK,SAAS,EAAE,CAAC;YACxB,MAAM,IAAI,KAAK,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,gBAAgB,CAAC,CAAA;QAC7C,CAAC;QACD,OAAO,KAAK,CAAA;IACd,CAAC,CAAA;IACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACrC,MAAM,GAAG,GAAG,IAAI,CAAC,CAAC,CAAE,CAAA;QACpB,QAAQ,GAAG,EAAE,CAAC;YACZ,KAAK,UAAU;gBACb,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,CAAA;gBACvB,MAAK;YACP,KAAK,UAAU;gBACb,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,CAAA;gBACvB,MAAK;YACP,KAAK,aAAa;gBAChB,IAAI,CAAC,SAAS,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAA;gBAClC,MAAK;YACP,KAAK,IAAI,CAAC;YACV,KAAK,UAAU;gBACb,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAA;gBAC/B,MAAK;YACP,KAAK,IAAI,CAAC;YACV,KAAK,YAAY,CAAC,CAAC,CAAC;gBAClB,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAA;gBACpC,IAAI,CAAC,QAAQ,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAA;gBACtC,MAAK;YACP,CAAC;YACD,KAAK,IAAI,CAAC;YACV,KAAK,QAAQ;gBACX,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAA;gBAClC,MAAK;YACP,KAAK,WAAW;gBACd,IAAI,CAAC,OAAO,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAA;gBAChC,MAAK;YACP,KAAK,SAAS;gBACZ,IAAI,CAAC,KAAK,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAA;gBAC9B,MAAK;YACP,KAAK,cAAc;gBACjB,IAAI,CAAC,UAAU,GAAG,IAAI,CAAC,CAAC,EAAE,CAAoB,CAAA;gBAC9C,MAAK;YACP,KAAK,SAAS;gBACZ,IAAI,CAAC,KAAK,GAAG,IAAI,CAAA;gBACjB,MAAK;YACP,KAAK,WAAW;gBACd,IAAI,CAAC,OAAO,GAAG,IAAI,CAAA;gBACnB,MAAK;YACP,KAAK,cAAc;gBACjB,IAAI,CAAC,UAAU,GAAG,IAAI,CAAA;gBACtB,IAAI,CAAC,OAAO,GAAG,IAAI,CAAA;gBACnB,MAAK;YACP,KAAK,cAAc;gBACjB,IAAI,CAAC,SAAS,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAA;gBAClC,MAAK;YACP,KAAK,SAAS;gBACZ,IAAI,CAAC,KAAK,GAAG,IAAI,CAAA;gBACjB,MAAK;YACP,KAAK,UAAU;gBACb,IAAI,CAAC,CAAC,EAAE,CAAC,CAAA;gBACT,MAAK;YACP,KAAK,IAAI,CAAC;YACV,KAAK,QAAQ;gBACX,OAAO,CAAC,MAAM,CAAC,KAAK,CAAC,KAAK,CAAC,CAAA;gBAC3B,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAA;gBACf,MAAK;YACP;gBACE,IAAI,GAAG,CAAC,UAAU,CAAC,GAAG,CAAC,EAAE,CAAC;oBACxB,MAAM,IAAI,KAAK,CAAC,kBAAkB,GAAG,EAAE,CAAC,CAAA;gBAC1C,CAAC;gBACD,IAAI,CAAC,IAAI,GAAG,GAAG,CAAA;QACnB,CAAC;IACH,CAAC;IACD,IAAI,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC;QACf,MAAM,IAAI,KAAK,CAAC,KAAK,CAAC,CAAA;IACxB,CAAC;IACD,OAAO,IAAI,CAAA;AACb,CAAC;AAED,MAAM,CAAC,KAAK,UAAU,IAAI,CAAC,IAAc;IACvC,MAAM,IAAI,GAAG,SAAS,CAAC,IAAI,CAAC,CAAA;IAC5B,MAAM,MAAM,GAAG,cAAc,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC;QAC3C,CAAC,CAAC,IAAI,UAAU,CAAC,IAAI,CAAC,IAAI,CAAC;QAC3B,CAAC,CAAC,IAAI,SAAS,CAAC,IAAI,CAAC,IAAI,CAAC,CAAA;IAC5B,MAAM,EAAE,GAAG,MAAM,OAAO,CAAC,IAAI,CAAC,MAAM,EAAE,EAAE,SAAS,EAAE,IAAI,CAAC,SAAS,EAAE,CAAC,CAAA;IACpE,MAAM,IAAI,GAAG;QACX,OAAO,EAAE,IAAI,CAAC,OAAO;QACrB,UAAU,EAAE,IAAI,CAAC,UAAU;QAC3B,GAAG,CAAC,IAAI,CAAC,KAAK,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,IAAI,CAAC,KAAK,EAAE,CAAC;KAC3D,CAAA;IACD,MAAM,KAAK,GAAG;QACZ,MAAM,EAAE,IAAI,CAAC,MAAM,IAAI,EAAE;QACzB,SAAS,EAAE,IAAI,CAAC,SAAS;QACzB,GAAG,CAAC,IAAI,CAAC,MAAM,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,IAAI,CAAC,MAAM,EAAE,CAAC;KAC9D,CAAA;IACD,MAAM,QAAQ,GACZ,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC;QACnB,CAAC,CAAC,MAAM,mBAAmB,CAAC,EAAE,EAAE,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC;QACjD,CAAC,CAAC,IAAI,CAAC,QAAQ;YACb,CAAC,CAAC,MAAM,kBAAkB,CACtB,EAAE,EACF,KAAK,EACL,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,EAChB,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,EAChB,IAAI,CACL;YACH,CAAC,CAAC,IAAI,CAAC,MAAM,KAAK,SAAS;gBACzB,CAAC,CAAC,MAAM,gBAAgB,CAAC,EAAE,EAAE,KAAK,EAAE,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC;gBACtD,CAAC,CAAC,SAAS,CAAA;IACnB,IAAI,CAAC,QAAQ,EAAE,CAAC;QACd,MAAM,IAAI,KAAK,CACb,mEAAmE,CACpE,CAAA;IACH,CAAC;IACD,IAAI,IAAI,CAAC,OAAO,EAAE,CAAC;QACjB,MAAM,QAAQ,CAAC,aAAa,EAAE,CAAA;IAChC,CAAC;IACD,MAAM,MAAM,GAAG,IAAI,CAAC,UAAU;QAC5B,CAAC,CAAC,QAAQ,CAAC,UAAU,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YAC9B,GAAG,CAAC;YACJ,IAAI,EACF,CAAC,CAAC,IAAI,IAAI,CAAC,CAAC,QAAQ,KAAK,SAAS,IAAI,CAAC,CAAC,MAAM,KAAK,SAAS;gBAC1D,CAAC,CAAC,cAAc,CAAC,EAAE,GAAG,CAAC,CAAC,IAAI,EAAE,QAAQ,EAAE,CAAC,CAAC,QAAQ,EAAE,EAAE,CAAC,CAAC,MAAM,CAAC;gBAC/D,CAAC,CAAC,SAAS;YACf,KAAK,EAAE,SAAS;SACjB,CAAC,CAAC;QACL,CAAC,CAAC,QAAQ,CAAC,MAAM,CAAC,IAAI,CAAC,KAAK,EAAE;YAC1B,KAAK,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,WAAW;SAC/C,CAAC,CAAA;IACN,OAAO,CAAC,MAAM,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC,MAAM,CAAC,IAAI,CAAC,CAAA;IACnD,IAAI,IAAI,CAAC,KAAK,EAAE,CAAC;QACf,MAAM,EAAE,OAAO,EAAE,YAAY,EAAE,GAAG,EAAE,CAAC,MAAM,CAAC,KAAK,CAAA;QACjD,MAAM,EACJ,iBAAiB,EACjB,aAAa,EACb,mBAAmB,EACnB,qBAAqB,GACtB,GAAG,QAAQ,CAAC,KAAK,CAAA;QAClB,OAAO,CAAC,MAAM,CAAC,KAAK,CAClB,qBAAqB,QAAQ,CAAC,SAAS,cAAc,QAAQ,CAAC,SAAS,WAAW,OAAO,aAAa,YAAY,WAAW,iBAAiB,cAAc,aAAa,mCAAmC,mBAAmB,WAAW,qBAAqB,YAAY,CAC5Q,CAAA;IACH,CAAC;AACH,CAAC"}
package/dist/db.d.ts CHANGED
@@ -1,7 +1,7 @@
1
- import type { ByteSource } from './filehandle.ts';
2
1
  import { GbwtRecord } from './gbwt/record.ts';
3
- import type { Pos } from './gbwt/record.ts';
4
2
  import { SqliteDatabase } from './sqlite/database.ts';
3
+ import type { ByteSource } from './filehandle.ts';
4
+ import type { Pos } from './gbwt/record.ts';
5
5
  import type { PagerOptions } from './sqlite/pager.ts';
6
6
  export interface PathName {
7
7
  sample: string;
@@ -12,8 +12,8 @@ export interface PathName {
12
12
  export declare const GENERIC_SAMPLE = "_gbwt_ref";
13
13
  export declare const SCHEMA_VERSION = "GBZ-base version 4";
14
14
  export declare class SchemaVersionError extends Error {
15
- readonly found: string | undefined;
16
15
  name: string;
16
+ readonly found: string | undefined;
17
17
  constructor(found: string | undefined);
18
18
  }
19
19
  export declare function formatPathName(name: PathName, end: number): string;
@@ -32,12 +32,12 @@ export interface GbzPath {
32
32
  isIndexed: boolean;
33
33
  }
34
34
  export declare class GbzRecord {
35
+ readonly sequenceLen: number;
35
36
  readonly handle: number;
36
37
  readonly edges: Pos[];
37
38
  readonly bwt: Uint8Array;
38
39
  readonly encodedSequence: Uint8Array;
39
40
  readonly next: number | undefined;
40
- readonly sequenceLen: number;
41
41
  private decoded;
42
42
  constructor(handle: number, edges: Pos[], bwt: Uint8Array, encodedSequence: Uint8Array, next: number | undefined);
43
43
  get id(): number;
@@ -47,9 +47,9 @@ export declare class GbzRecord {
47
47
  gbwt(): GbwtRecord;
48
48
  }
49
49
  export declare class GBZBase {
50
- readonly sqlite: SqliteDatabase;
51
50
  private tagCache;
52
51
  private pathCache;
52
+ readonly sqlite: SqliteDatabase;
53
53
  private constructor();
54
54
  static open(source: ByteSource, opts?: PagerOptions): Promise<GBZBase>;
55
55
  tags(): Promise<Map<string, string>>;
package/dist/db.js CHANGED
@@ -1,12 +1,12 @@
1
- import { GbwtRecord, decompressEdges } from "./gbwt/record.js";
2
1
  import { ENDMARKER, nodeId, nodeOrientation } from "./gbwt/node.js";
2
+ import { GbwtRecord, decompressEdges } from "./gbwt/record.js";
3
3
  import { decodeSequence, encodedSequenceLength } from "./gbwt/sequence.js";
4
4
  import { SqliteDatabase } from "./sqlite/database.js";
5
5
  export const GENERIC_SAMPLE = '_gbwt_ref';
6
6
  export const SCHEMA_VERSION = 'GBZ-base version 4';
7
7
  export class SchemaVersionError extends Error {
8
- found;
9
8
  name = 'SchemaVersionError';
9
+ found;
10
10
  constructor(found) {
11
11
  super(found === undefined
12
12
  ? `not a gbz-base database: its Tags table has no version`
@@ -18,12 +18,12 @@ export function formatPathName(name, end) {
18
18
  return `${name.sample}#${name.haplotype}#${name.contig}[${name.fragment}-${end}]`;
19
19
  }
20
20
  export class GbzRecord {
21
+ sequenceLen;
21
22
  handle;
22
23
  edges;
23
24
  bwt;
24
25
  encodedSequence;
25
26
  next;
26
- sequenceLen;
27
27
  decoded;
28
28
  constructor(handle, edges, bwt, encodedSequence, next) {
29
29
  this.handle = handle;
@@ -74,8 +74,14 @@ function blob(value, what) {
74
74
  function rowToPath(rowid, values) {
75
75
  return {
76
76
  handle: rowid,
77
- fwStart: { node: num(values[1], 'Paths.fw_node'), offset: num(values[2], 'Paths.fw_offset') },
78
- revStart: { node: num(values[3], 'Paths.rev_node'), offset: num(values[4], 'Paths.rev_offset') },
77
+ fwStart: {
78
+ node: num(values[1], 'Paths.fw_node'),
79
+ offset: num(values[2], 'Paths.fw_offset'),
80
+ },
81
+ revStart: {
82
+ node: num(values[3], 'Paths.rev_node'),
83
+ offset: num(values[4], 'Paths.rev_offset'),
84
+ },
79
85
  name: {
80
86
  sample: str(values[5], 'Paths.sample'),
81
87
  contig: str(values[6], 'Paths.contig'),
@@ -86,9 +92,9 @@ function rowToPath(rowid, values) {
86
92
  };
87
93
  }
88
94
  export class GBZBase {
89
- sqlite;
90
95
  tagCache;
91
96
  pathCache;
97
+ sqlite;
92
98
  constructor(sqlite) {
93
99
  this.sqlite = sqlite;
94
100
  }
@@ -150,7 +156,7 @@ export class GBZBase {
150
156
  return (await this.paths()).filter(p => p.name.sample === sample);
151
157
  }
152
158
  get hasHaplotypeIndex() {
153
- return this.sqlite.has('HaplotypeSamples') && this.sqlite.has('HaplotypeLengths');
159
+ return (this.sqlite.has('HaplotypeSamples') && this.sqlite.has('HaplotypeLengths'));
154
160
  }
155
161
  async haplotypeSampleInterval() {
156
162
  const value = await this.tag('haplotype_index_interval');
@@ -161,13 +167,18 @@ export class GBZBase {
161
167
  node: num(values[0], 'HaplotypeSamples.node_handle'),
162
168
  offset: num(values[1], 'HaplotypeSamples.node_offset'),
163
169
  pathHandle: num(values[2], 'HaplotypeSamples.path_handle'),
164
- orientation: num(values[3], 'HaplotypeSamples.orientation') === 0 ? 'forward' : 'reverse',
170
+ orientation: num(values[3], 'HaplotypeSamples.orientation') === 0
171
+ ? 'forward'
172
+ : 'reverse',
165
173
  pathOffset: num(values[4], 'HaplotypeSamples.path_offset'),
166
174
  };
167
175
  }
168
176
  async haplotypeSamplesInRange(minHandle, maxHandle) {
169
177
  const samples = [];
170
- for await (const key of this.sqlite.indexScanFrom('HaplotypeSamples', [minHandle, 0])) {
178
+ for await (const key of this.sqlite.indexScanFrom('HaplotypeSamples', [
179
+ minHandle,
180
+ 0,
181
+ ])) {
171
182
  const node = num(key[0], 'HaplotypeSamples.node_handle');
172
183
  if (node > maxHandle) {
173
184
  break;
@@ -180,8 +191,11 @@ export class GBZBase {
180
191
  return samples;
181
192
  }
182
193
  async haplotypeSampleAt(node, offset) {
183
- const key = await this.sqlite.indexSeekLE('HaplotypeSamples', [node, offset]);
184
- if (!key || key[0] !== node || key[1] !== offset) {
194
+ const key = await this.sqlite.indexSeekLE('HaplotypeSamples', [
195
+ node,
196
+ offset,
197
+ ]);
198
+ if (key?.[0] !== node || key[1] !== offset) {
185
199
  return undefined;
186
200
  }
187
201
  const row = await this.sqlite.byRowid('HaplotypeSamples', num(key[2], 'HaplotypeSamples rowid'));
@@ -192,8 +206,11 @@ export class GBZBase {
192
206
  return row ? num(row[1], 'HaplotypeLengths.length') : undefined;
193
207
  }
194
208
  async indexedPosition(pathHandle, pathOffset) {
195
- const key = await this.sqlite.indexSeekLE('ReferenceIndex', [pathHandle, pathOffset]);
196
- if (!key || key[0] !== pathHandle) {
209
+ const key = await this.sqlite.indexSeekLE('ReferenceIndex', [
210
+ pathHandle,
211
+ pathOffset,
212
+ ]);
213
+ if (key?.[0] !== pathHandle) {
197
214
  return undefined;
198
215
  }
199
216
  const row = await this.sqlite.byRowid('ReferenceIndex', num(key[2], 'ReferenceIndex rowid'));
@@ -202,7 +219,11 @@ export class GBZBase {
202
219
  }
203
220
  return {
204
221
  pathOffset: num(row[1], 'ReferenceIndex.path_offset'),
205
- pos: { node: num(row[2], 'ReferenceIndex.node_handle'), offset: num(row[3], 'ReferenceIndex.node_offset') },
222
+ pos: {
223
+ node: num(row[2], 'ReferenceIndex.node_handle'),
224
+ offset: num(row[3], 'ReferenceIndex.node_offset'),
225
+ },
206
226
  };
207
227
  }
208
228
  }
229
+ //# sourceMappingURL=db.js.map
package/dist/db.js.map ADDED
@@ -0,0 +1 @@
1
+ {"version":3,"file":"db.js","sourceRoot":"","sources":["../src/db.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,SAAS,EAAE,MAAM,EAAE,eAAe,EAAE,MAAM,gBAAgB,CAAA;AACnE,OAAO,EAAE,UAAU,EAAE,eAAe,EAAE,MAAM,kBAAkB,CAAA;AAC9D,OAAO,EAAE,cAAc,EAAE,qBAAqB,EAAE,MAAM,oBAAoB,CAAA;AAC1E,OAAO,EAAE,cAAc,EAAE,MAAM,sBAAsB,CAAA;AAcrD,MAAM,CAAC,MAAM,cAAc,GAAG,WAAW,CAAA;AACzC,MAAM,CAAC,MAAM,cAAc,GAAG,oBAAoB,CAAA;AAElD,MAAM,OAAO,kBAAmB,SAAQ,KAAK;IAClC,IAAI,GAAG,oBAAoB,CAAA;IAE3B,KAAK,CAAoB;IAElC,YAAY,KAAyB;QACnC,KAAK,CACH,KAAK,KAAK,SAAS;YACjB,CAAC,CAAC,wDAAwD;YAC1D,CAAC,CAAC,gCAAgC,KAAK,+BAA+B,cAAc,GAAG,CAC1F,CAAA;QACD,IAAI,CAAC,KAAK,GAAG,KAAK,CAAA;IACpB,CAAC;CACF;AAED,MAAM,UAAU,cAAc,CAAC,IAAc,EAAE,GAAW;IACxD,OAAO,GAAG,IAAI,CAAC,MAAM,IAAI,IAAI,CAAC,SAAS,IAAI,IAAI,CAAC,MAAM,IAAI,IAAI,CAAC,QAAQ,IAAI,GAAG,GAAG,CAAA;AACnF,CAAC;AAkBD,MAAM,OAAO,SAAS;IACX,WAAW,CAAQ;IACnB,MAAM,CAAQ;IACd,KAAK,CAAO;IACZ,GAAG,CAAY;IACf,eAAe,CAAY;IAC3B,IAAI,CAAoB;IACzB,OAAO,CAAoB;IAEnC,YACE,MAAc,EACd,KAAY,EACZ,GAAe,EACf,eAA2B,EAC3B,IAAwB;QAExB,IAAI,CAAC,MAAM,GAAG,MAAM,CAAA;QACpB,IAAI,CAAC,KAAK,GAAG,KAAK,CAAA;QAClB,IAAI,CAAC,GAAG,GAAG,GAAG,CAAA;QACd,IAAI,CAAC,eAAe,GAAG,eAAe,CAAA;QACtC,IAAI,CAAC,IAAI,GAAG,IAAI,CAAA;QAChB,IAAI,CAAC,WAAW,GAAG,qBAAqB,CAAC,eAAe,CAAC,CAAA;IAC3D,CAAC;IAED,IAAI,EAAE;QACJ,OAAO,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,CAAA;IAC5B,CAAC;IAED,IAAI,WAAW;QACb,OAAO,eAAe,CAAC,IAAI,CAAC,MAAM,CAAC,CAAA;IACrC,CAAC;IAED,IAAI,QAAQ;QACV,IAAI,CAAC,OAAO,KAAK,cAAc,CAAC,IAAI,CAAC,eAAe,CAAC,CAAA;QACrD,OAAO,IAAI,CAAC,OAAO,CAAA;IACrB,CAAC;IAED,UAAU;QACR,OAAO,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,KAAK,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAA;IACtE,CAAC;IAED,IAAI;QACF,IAAI,IAAI,CAAC,KAAK,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;YAC5B,MAAM,IAAI,KAAK,CAAC,0BAA0B,IAAI,CAAC,MAAM,WAAW,CAAC,CAAA;QACnE,CAAC;QACD,OAAO,IAAI,UAAU,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,CAAA;IAC7C,CAAC;CACF;AAED,SAAS,GAAG,CAAC,KAA2B,EAAE,IAAY;IACpD,IAAI,OAAO,KAAK,KAAK,QAAQ,EAAE,CAAC;QAC9B,MAAM,IAAI,KAAK,CAAC,GAAG,IAAI,kCAAkC,CAAC,CAAA;IAC5D,CAAC;IACD,OAAO,KAAK,CAAA;AACd,CAAC;AAED,SAAS,GAAG,CAAC,KAA2B,EAAE,IAAY;IACpD,IAAI,OAAO,KAAK,KAAK,QAAQ,EAAE,CAAC;QAC9B,MAAM,IAAI,KAAK,CAAC,GAAG,IAAI,8BAA8B,CAAC,CAAA;IACxD,CAAC;IACD,OAAO,KAAK,CAAA;AACd,CAAC;AAED,SAAS,IAAI,CAAC,KAA2B,EAAE,IAAY;IACrD,IAAI,CAAC,CAAC,KAAK,YAAY,UAAU,CAAC,EAAE,CAAC;QACnC,MAAM,IAAI,KAAK,CAAC,GAAG,IAAI,gCAAgC,CAAC,CAAA;IAC1D,CAAC;IACD,OAAO,KAAK,CAAA;AACd,CAAC;AAED,SAAS,SAAS,CAAC,KAAa,EAAE,MAAkB;IAClD,OAAO;QACL,MAAM,EAAE,KAAK;QACb,OAAO,EAAE;YACP,IAAI,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,eAAe,CAAC;YACrC,MAAM,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,iBAAiB,CAAC;SAC1C;QACD,QAAQ,EAAE;YACR,IAAI,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,gBAAgB,CAAC;YACtC,MAAM,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,kBAAkB,CAAC;SAC3C;QACD,IAAI,EAAE;YACJ,MAAM,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,cAAc,CAAC;YACtC,MAAM,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,cAAc,CAAC;YACtC,SAAS,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,iBAAiB,CAAC;YAC5C,QAAQ,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,gBAAgB,CAAC;SAC3C;QACD,SAAS,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,kBAAkB,CAAC,KAAK,CAAC;KACpD,CAAA;AACH,CAAC;AAED,MAAM,OAAO,OAAO;IACV,QAAQ,CAA0C;IAClD,SAAS,CAAgC;IAExC,MAAM,CAAgB;IAE/B,YAAoB,MAAsB;QACxC,IAAI,CAAC,MAAM,GAAG,MAAM,CAAA;IACtB,CAAC;IAED,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,MAAkB,EAAE,OAAqB,EAAE;QAC3D,MAAM,MAAM,GAAG,MAAM,cAAc,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,CAAA;QACtD,KAAK,MAAM,KAAK,IAAI,CAAC,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,gBAAgB,CAAC,EAAE,CAAC;YACjE,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAA;QACxB,CAAC;QACD,MAAM,EAAE,GAAG,IAAI,OAAO,CAAC,MAAM,CAAC,CAAA;QAC9B,MAAM,OAAO,GAAG,MAAM,EAAE,CAAC,GAAG,CAAC,SAAS,CAAC,CAAA;QACvC,IAAI,OAAO,KAAK,cAAc,EAAE,CAAC;YAC/B,MAAM,IAAI,kBAAkB,CAAC,OAAO,CAAC,CAAA;QACvC,CAAC;QACD,OAAO,EAAE,CAAA;IACX,CAAC;IAED,IAAI;QACF,IAAI,CAAC,QAAQ,KAAK,CAAC,KAAK,IAAI,EAAE;YAC5B,MAAM,IAAI,GAAG,IAAI,GAAG,EAAkB,CAAA;YACtC,IAAI,KAAK,EAAE,MAAM,EAAE,MAAM,EAAE,IAAI,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,EAAE,CAAC;gBACxD,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,UAAU,CAAC,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,YAAY,CAAC,CAAC,CAAA;YACpE,CAAC;YACD,OAAO,IAAI,CAAA;QACb,CAAC,CAAC,EAAE,CAAA;QACJ,OAAO,IAAI,CAAC,QAAQ,CAAA;IACtB,CAAC;IAED,KAAK,CAAC,GAAG,CAAC,GAAW;QACnB,OAAO,CAAC,MAAM,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;IACrC,CAAC;IAED,KAAK,CAAC,SAAS,CAAC,MAAc;QAC5B,MAAM,GAAG,GAAG,MAAM,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC,OAAO,EAAE,MAAM,CAAC,CAAA;QACtD,IAAI,CAAC,GAAG,EAAE,CAAC;YACT,OAAO,SAAS,CAAA;QAClB,CAAC;QACD,MAAM,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAA;QACnB,OAAO,IAAI,SAAS,CAClB,MAAM,EACN,eAAe,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,aAAa,CAAC,CAAC,EAC5C,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,WAAW,CAAC,EACzB,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,gBAAgB,CAAC,EAC9B,OAAO,IAAI,KAAK,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,SAAS,CAC5C,CAAA;IACH,CAAC;IAED,KAAK;QACH,IAAI,CAAC,SAAS,KAAK,CAAC,KAAK,IAAI,EAAE;YAC7B,MAAM,KAAK,GAAc,EAAE,CAAA;YAC3B,IAAI,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,IAAI,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,OAAO,CAAC,EAAE,CAAC;gBAChE,KAAK,CAAC,IAAI,CAAC,SAAS,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,CAAA;YACtC,CAAC;YACD,OAAO,KAAK,CAAA;QACd,CAAC,CAAC,EAAE,CAAA;QACJ,OAAO,IAAI,CAAC,SAAS,CAAA;IACvB,CAAC;IAED,KAAK,CAAC,OAAO,CAAC,MAAc;QAC1B,MAAM,GAAG,GAAG,MAAM,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC,OAAO,EAAE,MAAM,CAAC,CAAA;QACtD,OAAO,GAAG,CAAC,CAAC,CAAC,SAAS,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,SAAS,CAAA;IACjD,CAAC;IAED,KAAK,CAAC,QAAQ,CAAC,IAAc;QAC3B,MAAM,UAAU,GAAG,CAAC,MAAM,IAAI,CAAC,KAAK,EAAE,CAAC,CAAC,MAAM,CAC5C,CAAC,CAAC,EAAE,CACF,CAAC,CAAC,IAAI,CAAC,MAAM,KAAK,IAAI,CAAC,MAAM;YAC7B,CAAC,CAAC,IAAI,CAAC,MAAM,KAAK,IAAI,CAAC,MAAM;YAC7B,CAAC,CAAC,IAAI,CAAC,SAAS,KAAK,IAAI,CAAC,SAAS;YACnC,CAAC,CAAC,IAAI,CAAC,QAAQ,IAAI,IAAI,CAAC,QAAQ,CACnC,CAAA;QACD,OAAO,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,GAAG,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAA;IACxE,CAAC;IAED,KAAK,CAAC,cAAc,CAAC,MAAc;QACjC,OAAO,CAAC,MAAM,IAAI,CAAC,KAAK,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,KAAK,MAAM,CAAC,CAAA;IACnE,CAAC;IAED,IAAI,iBAAiB;QACnB,OAAO,CACL,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,kBAAkB,CAAC,IAAI,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,kBAAkB,CAAC,CAC3E,CAAA;IACH,CAAC;IAED,KAAK,CAAC,uBAAuB;QAC3B,MAAM,KAAK,GAAG,MAAM,IAAI,CAAC,GAAG,CAAC,0BAA0B,CAAC,CAAA;QACxD,OAAO,KAAK,KAAK,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,MAAM,CAAC,KAAK,CAAC,CAAA;IACxD,CAAC;IAEO,aAAa,CAAC,MAAkB;QACtC,OAAO;YACL,IAAI,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,8BAA8B,CAAC;YACpD,MAAM,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,8BAA8B,CAAC;YACtD,UAAU,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,8BAA8B,CAAC;YAC1D,WAAW,EACT,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,8BAA8B,CAAC,KAAK,CAAC;gBAClD,CAAC,CAAC,SAAS;gBACX,CAAC,CAAC,SAAS;YACf,UAAU,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,8BAA8B,CAAC;SAC3D,CAAA;IACH,CAAC;IAED,KAAK,CAAC,uBAAuB,CAAC,SAAiB,EAAE,SAAiB;QAChE,MAAM,OAAO,GAAsB,EAAE,CAAA;QACrC,IAAI,KAAK,EAAE,MAAM,GAAG,IAAI,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,kBAAkB,EAAE;YACpE,SAAS;YACT,CAAC;SACF,CAAC,EAAE,CAAC;YACH,MAAM,IAAI,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,8BAA8B,CAAC,CAAA;YACxD,IAAI,IAAI,GAAG,SAAS,EAAE,CAAC;gBACrB,MAAK;YACP,CAAC;YACD,MAAM,GAAG,GAAG,MAAM,IAAI,CAAC,MAAM,CAAC,OAAO,CACnC,kBAAkB,EAClB,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,wBAAwB,CAAC,CACtC,CAAA;YACD,IAAI,GAAG,EAAE,CAAC;gBACR,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,aAAa,CAAC,GAAG,CAAC,CAAC,CAAA;YACvC,CAAC;QACH,CAAC;QACD,OAAO,OAAO,CAAA;IAChB,CAAC;IAED,KAAK,CAAC,iBAAiB,CAAC,IAAY,EAAE,MAAc;QAClD,MAAM,GAAG,GAAG,MAAM,IAAI,CAAC,MAAM,CAAC,WAAW,CAAC,kBAAkB,EAAE;YAC5D,IAAI;YACJ,MAAM;SACP,CAAC,CAAA;QACF,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,KAAK,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC,KAAK,MAAM,EAAE,CAAC;YAC3C,OAAO,SAAS,CAAA;QAClB,CAAC;QACD,MAAM,GAAG,GAAG,MAAM,IAAI,CAAC,MAAM,CAAC,OAAO,CACnC,kBAAkB,EAClB,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,wBAAwB,CAAC,CACtC,CAAA;QACD,OAAO,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,SAAS,CAAA;IAClD,CAAC;IAED,KAAK,CAAC,eAAe,CAAC,UAAkB;QACtC,MAAM,GAAG,GAAG,MAAM,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC,kBAAkB,EAAE,UAAU,CAAC,CAAA;QACrE,OAAO,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,yBAAyB,CAAC,CAAC,CAAC,CAAC,SAAS,CAAA;IACjE,CAAC;IAED,KAAK,CAAC,eAAe,CACnB,UAAkB,EAClB,UAAkB;QAElB,MAAM,GAAG,GAAG,MAAM,IAAI,CAAC,MAAM,CAAC,WAAW,CAAC,gBAAgB,EAAE;YAC1D,UAAU;YACV,UAAU;SACX,CAAC,CAAA;QACF,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,KAAK,UAAU,EAAE,CAAC;YAC5B,OAAO,SAAS,CAAA;QAClB,CAAC;QACD,MAAM,GAAG,GAAG,MAAM,IAAI,CAAC,MAAM,CAAC,OAAO,CACnC,gBAAgB,EAChB,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,sBAAsB,CAAC,CACpC,CAAA;QACD,IAAI,CAAC,GAAG,EAAE,CAAC;YACT,MAAM,IAAI,KAAK,CAAC,uDAAuD,CAAC,CAAA;QAC1E,CAAC;QACD,OAAO;YACL,UAAU,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,4BAA4B,CAAC;YACrD,GAAG,EAAE;gBACH,IAAI,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,4BAA4B,CAAC;gBAC/C,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,4BAA4B,CAAC;aAClD;SACF,CAAA;IACH,CAAC;CACF"}
@@ -1 +1,2 @@
1
1
  export {};
2
+ //# sourceMappingURL=filehandle.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"filehandle.js","sourceRoot":"","sources":["../src/filehandle.ts"],"names":[],"mappings":""}
@@ -1,6 +1,6 @@
1
1
  export declare class ByteCodeReader {
2
- private bytes;
3
2
  offset: number;
3
+ private bytes;
4
4
  constructor(bytes: Uint8Array);
5
5
  get done(): boolean;
6
6
  byte(): number | undefined;
@@ -1,6 +1,6 @@
1
1
  export class ByteCodeReader {
2
- bytes;
3
2
  offset = 0;
3
+ bytes;
4
4
  constructor(bytes) {
5
5
  this.bytes = bytes;
6
6
  }
@@ -39,13 +39,16 @@ export class RunReader {
39
39
  constructor(bytes, sigma) {
40
40
  this.source = new ByteCodeReader(bytes);
41
41
  this.sigma = sigma === 0 ? Number.MAX_SAFE_INTEGER : sigma;
42
- this.threshold = this.sigma < RLE_THRESHOLD ? Math.floor(RLE_UNIVERSE / this.sigma) : 0;
42
+ this.threshold =
43
+ this.sigma < RLE_THRESHOLD ? Math.floor(RLE_UNIVERSE / this.sigma) : 0;
43
44
  }
44
45
  next() {
45
46
  if (this.sigma >= RLE_THRESHOLD) {
46
47
  const value = this.source.int();
47
48
  const len = this.source.int();
48
- return value === undefined || len === undefined ? undefined : { value, len: len + 1 };
49
+ return value === undefined || len === undefined
50
+ ? undefined
51
+ : { value, len: len + 1 };
49
52
  }
50
53
  const byte = this.source.byte();
51
54
  if (byte === undefined) {
@@ -68,3 +71,4 @@ export class RunReader {
68
71
  }
69
72
  }
70
73
  }
74
+ //# sourceMappingURL=bytecode.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"bytecode.js","sourceRoot":"","sources":["../../src/gbwt/bytecode.ts"],"names":[],"mappings":"AAAA,MAAM,OAAO,cAAc;IACzB,MAAM,GAAG,CAAC,CAAA;IACF,KAAK,CAAY;IAEzB,YAAY,KAAiB;QAC3B,IAAI,CAAC,KAAK,GAAG,KAAK,CAAA;IACpB,CAAC;IAED,IAAI,IAAI;QACN,OAAO,IAAI,CAAC,MAAM,IAAI,IAAI,CAAC,KAAK,CAAC,MAAM,CAAA;IACzC,CAAC;IAED,IAAI;QACF,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,MAAM,CAAC,CAAA;QACrC,IAAI,KAAK,KAAK,SAAS,EAAE,CAAC;YACxB,OAAO,SAAS,CAAA;QAClB,CAAC;QACD,IAAI,CAAC,MAAM,IAAI,CAAC,CAAA;QAChB,OAAO,KAAK,CAAA;IACd,CAAC;IAED,GAAG;QACD,IAAI,KAAK,GAAG,CAAC,CAAA;QACb,IAAI,MAAM,GAAG,CAAC,CAAA;QACd,OAAO,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC;YACvC,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,MAAM,CAAE,CAAA;YACtC,IAAI,CAAC,MAAM,IAAI,CAAC,CAAA;YAChB,MAAM,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,KAAK,CAAA;YAChC,KAAK,IAAI,GAAG,CAAA;YACZ,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;gBACzB,OAAO,MAAM,CAAA;YACf,CAAC;QACH,CAAC;QACD,OAAO,SAAS,CAAA;IAClB,CAAC;CACF;AAOD,MAAM,aAAa,GAAG,GAAG,CAAA;AACzB,MAAM,YAAY,GAAG,GAAG,CAAA;AAExB,MAAM,OAAO,SAAS;IACZ,MAAM,CAAgB;IACtB,KAAK,CAAQ;IACb,SAAS,CAAQ;IAEzB,YAAY,KAAiB,EAAE,KAAa;QAC1C,IAAI,CAAC,MAAM,GAAG,IAAI,cAAc,CAAC,KAAK,CAAC,CAAA;QACvC,IAAI,CAAC,KAAK,GAAG,KAAK,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,gBAAgB,CAAC,CAAC,CAAC,KAAK,CAAA;QAC1D,IAAI,CAAC,SAAS;YACZ,IAAI,CAAC,KAAK,GAAG,aAAa,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,YAAY,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;IAC1E,CAAC;IAED,IAAI;QACF,IAAI,IAAI,CAAC,KAAK,IAAI,aAAa,EAAE,CAAC;YAChC,MAAM,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,GAAG,EAAE,CAAA;YAC/B,MAAM,GAAG,GAAG,IAAI,CAAC,MAAM,CAAC,GAAG,EAAE,CAAA;YAC7B,OAAO,KAAK,KAAK,SAAS,IAAI,GAAG,KAAK,SAAS;gBAC7C,CAAC,CAAC,SAAS;gBACX,CAAC,CAAC,EAAE,KAAK,EAAE,GAAG,EAAE,GAAG,GAAG,CAAC,EAAE,CAAA;QAC7B,CAAC;QACD,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,CAAC,IAAI,EAAE,CAAA;QAC/B,IAAI,IAAI,KAAK,SAAS,EAAE,CAAC;YACvB,OAAO,SAAS,CAAA;QAClB,CAAC;QACD,MAAM,KAAK,GAAG,IAAI,GAAG,IAAI,CAAC,KAAK,CAAA;QAC/B,IAAI,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAA;QAC3C,IAAI,GAAG,KAAK,IAAI,CAAC,SAAS,EAAE,CAAC;YAC3B,MAAM,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,GAAG,EAAE,CAAA;YAC/B,IAAI,KAAK,KAAK,SAAS,EAAE,CAAC;gBACxB,OAAO,SAAS,CAAA;YAClB,CAAC;YACD,GAAG,IAAI,KAAK,CAAA;QACd,CAAC;QACD,OAAO,EAAE,KAAK,EAAE,GAAG,EAAE,CAAA;IACvB,CAAC;IAED,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC;QAChB,KAAK,IAAI,GAAG,GAAG,IAAI,CAAC,IAAI,EAAE,EAAE,GAAG,KAAK,SAAS,EAAE,GAAG,GAAG,IAAI,CAAC,IAAI,EAAE,EAAE,CAAC;YACjE,MAAM,GAAG,CAAA;QACX,CAAC;IACH,CAAC;CACF"}
package/dist/gbwt/node.js CHANGED
@@ -32,7 +32,9 @@ export function exitOrientation(side) {
32
32
  export function edgeIsCanonical(from, to) {
33
33
  const fromId = nodeId(from);
34
34
  const toId = nodeId(to);
35
- return isReverse(from) ? toId > fromId || (toId === fromId && !isReverse(to)) : toId >= fromId;
35
+ return isReverse(from)
36
+ ? toId > fromId || (toId === fromId && !isReverse(to))
37
+ : toId >= fromId;
36
38
  }
37
39
  export function pathIsCanonical(path) {
38
40
  const first = path[0];
@@ -40,5 +42,8 @@ export function pathIsCanonical(path) {
40
42
  if (first === undefined || last === undefined) {
41
43
  return true;
42
44
  }
43
- return isReverse(first) === isReverse(last) ? !isReverse(first) : edgeIsCanonical(first, last);
45
+ return isReverse(first) === isReverse(last)
46
+ ? !isReverse(first)
47
+ : edgeIsCanonical(first, last);
44
48
  }
49
+ //# sourceMappingURL=node.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"node.js","sourceRoot":"","sources":["../../src/gbwt/node.ts"],"names":[],"mappings":"AAAA,MAAM,CAAC,MAAM,SAAS,GAAG,CAAC,CAAA;AAI1B,MAAM,UAAU,UAAU,CAAC,EAAU,EAAE,WAAwB;IAC7D,OAAO,CAAC,GAAG,EAAE,GAAG,CAAC,WAAW,KAAK,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;AACrD,CAAC;AAED,MAAM,UAAU,MAAM,CAAC,MAAc;IACnC,OAAO,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,CAAA;AAC/B,CAAC;AAED,MAAM,UAAU,eAAe,CAAC,MAAc;IAC5C,OAAO,MAAM,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,SAAS,CAAA;AACjD,CAAC;AAED,MAAM,UAAU,SAAS,CAAC,MAAc;IACtC,OAAO,MAAM,GAAG,CAAC,KAAK,CAAC,CAAA;AACzB,CAAC;AAED,MAAM,UAAU,QAAQ,CAAC,MAAc;IACrC,OAAO,MAAM,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAA;AACnD,CAAC;AAID,MAAM,UAAU,QAAQ,CAAC,IAAc;IACrC,OAAO,IAAI,KAAK,MAAM,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,CAAA;AAC3C,CAAC;AAED,MAAM,UAAU,SAAS,CAAC,WAAwB;IAChD,OAAO,WAAW,KAAK,SAAS,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,OAAO,CAAA;AACrD,CAAC;AAED,MAAM,UAAU,QAAQ,CAAC,WAAwB;IAC/C,OAAO,WAAW,KAAK,SAAS,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,CAAA;AACrD,CAAC;AAED,MAAM,UAAU,gBAAgB,CAAC,IAAc;IAC7C,OAAO,IAAI,KAAK,MAAM,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,SAAS,CAAA;AAChD,CAAC;AAED,MAAM,UAAU,eAAe,CAAC,IAAc;IAC5C,OAAO,IAAI,KAAK,OAAO,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,SAAS,CAAA;AACjD,CAAC;AAED,MAAM,UAAU,eAAe,CAAC,IAAY,EAAE,EAAU;IACtD,MAAM,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC,CAAA;IAC3B,MAAM,IAAI,GAAG,MAAM,CAAC,EAAE,CAAC,CAAA;IACvB,OAAO,SAAS,CAAC,IAAI,CAAC;QACpB,CAAC,CAAC,IAAI,GAAG,MAAM,IAAI,CAAC,IAAI,KAAK,MAAM,IAAI,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC;QACtD,CAAC,CAAC,IAAI,IAAI,MAAM,CAAA;AACpB,CAAC;AAED,MAAM,UAAU,eAAe,CAAC,IAAc;IAC5C,MAAM,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC,CAAA;IACrB,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,CAAA;IAClC,IAAI,KAAK,KAAK,SAAS,IAAI,IAAI,KAAK,SAAS,EAAE,CAAC;QAC9C,OAAO,IAAI,CAAA;IACb,CAAC;IACD,OAAO,SAAS,CAAC,KAAK,CAAC,KAAK,SAAS,CAAC,IAAI,CAAC;QACzC,CAAC,CAAC,CAAC,SAAS,CAAC,KAAK,CAAC;QACnB,CAAC,CAAC,eAAe,CAAC,KAAK,EAAE,IAAI,CAAC,CAAA;AAClC,CAAC"}
@@ -43,7 +43,9 @@ export class GbwtRecord {
43
43
  const edge = this.edge(run.value);
44
44
  const soFar = offsets[run.value];
45
45
  if (offset + run.len > i) {
46
- return edge.node === ENDMARKER ? undefined : { node: edge.node, offset: soFar + (i - offset) };
46
+ return edge.node === ENDMARKER
47
+ ? undefined
48
+ : { node: edge.node, offset: soFar + (i - offset) };
47
49
  }
48
50
  offsets[run.value] = soFar + run.len;
49
51
  offset += run.len;
@@ -94,7 +96,10 @@ export class GbwtRecord {
94
96
  return undefined;
95
97
  }
96
98
  predecessorAt(i) {
97
- const counts = this.edges.map(e => ({ node: e.node === ENDMARKER ? ENDMARKER : flipNode(e.node), count: 0 }));
99
+ const counts = this.edges.map(e => ({
100
+ node: e.node === ENDMARKER ? ENDMARKER : flipNode(e.node),
101
+ count: 0,
102
+ }));
98
103
  for (const run of this.runs()) {
99
104
  const entry = counts[run.value];
100
105
  if (entry) {
@@ -131,3 +136,4 @@ export class GbwtRecord {
131
136
  return result;
132
137
  }
133
138
  }
139
+ //# sourceMappingURL=record.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"record.js","sourceRoot":"","sources":["../../src/gbwt/record.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,cAAc,EAAE,SAAS,EAAE,MAAM,eAAe,CAAA;AACzD,OAAO,EAAE,SAAS,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,WAAW,CAAA;AAOvD,MAAM,UAAU,eAAe,CAAC,KAAiB;IAC/C,MAAM,MAAM,GAAG,IAAI,cAAc,CAAC,KAAK,CAAC,CAAA;IACxC,MAAM,KAAK,GAAG,MAAM,CAAC,GAAG,EAAE,CAAA;IAC1B,IAAI,KAAK,KAAK,SAAS,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;QACvC,OAAO,EAAE,CAAA;IACX,CAAC;IACD,MAAM,KAAK,GAAU,EAAE,CAAA;IACvB,IAAI,IAAI,GAAG,CAAC,CAAA;IACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,MAAM,KAAK,GAAG,MAAM,CAAC,GAAG,EAAE,CAAA;QAC1B,MAAM,MAAM,GAAG,MAAM,CAAC,GAAG,EAAE,CAAA;QAC3B,IAAI,KAAK,KAAK,SAAS,IAAI,MAAM,KAAK,SAAS,EAAE,CAAC;YAChD,MAAM,IAAI,KAAK,CAAC,2BAA2B,CAAC,CAAA;QAC9C,CAAC;QACD,IAAI,IAAI,KAAK,CAAA;QACb,KAAK,CAAC,IAAI,CAAC,EAAE,IAAI,EAAE,IAAI,EAAE,MAAM,EAAE,CAAC,CAAA;IACpC,CAAC;IACD,OAAO,KAAK,CAAA;AACd,CAAC;AAED,MAAM,OAAO,UAAU;IACZ,KAAK,CAAO;IACZ,GAAG,CAAY;IAExB,YAAY,KAAY,EAAE,GAAe;QACvC,IAAI,CAAC,KAAK,GAAG,KAAK,CAAA;QAClB,IAAI,CAAC,GAAG,GAAG,GAAG,CAAA;IAChB,CAAC;IAED,IAAI;QACF,OAAO,IAAI,SAAS,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,CAAA;IACnD,CAAC;IAEO,IAAI,CAAC,IAAY;QACvB,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAA;QAC7B,IAAI,IAAI,KAAK,SAAS,EAAE,CAAC;YACvB,MAAM,IAAI,KAAK,CACb,gCAAgC,IAAI,OAAO,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,CAC/D,CAAA;QACH,CAAC;QACD,OAAO,IAAI,CAAA;IACb,CAAC;IAED,EAAE,CAAC,CAAS;QACV,MAAM,OAAO,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,CAAA;QAC7C,IAAI,MAAM,GAAG,CAAC,CAAA;QACd,KAAK,MAAM,GAAG,IAAI,IAAI,CAAC,IAAI,EAAE,EAAE,CAAC;YAC9B,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAA;YACjC,MAAM,KAAK,GAAG,OAAO,CAAC,GAAG,CAAC,KAAK,CAAE,CAAA;YACjC,IAAI,MAAM,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,EAAE,CAAC;gBACzB,OAAO,IAAI,CAAC,IAAI,KAAK,SAAS;oBAC5B,CAAC,CAAC,SAAS;oBACX,CAAC,CAAC,EAAE,IAAI,EAAE,IAAI,CAAC,IAAI,EAAE,MAAM,EAAE,KAAK,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,EAAE,CAAA;YACvD,CAAC;YACD,OAAO,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,KAAK,GAAG,GAAG,CAAC,GAAG,CAAA;YACpC,MAAM,IAAI,GAAG,CAAC,GAAG,CAAA;QACnB,CAAC;QACD,OAAO,SAAS,CAAA;IAClB,CAAC;IAEO,MAAM,CAAC,IAAY;QACzB,IAAI,GAAG,GAAG,CAAC,CAAA;QACX,IAAI,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAA;QAC5B,OAAO,GAAG,GAAG,IAAI,EAAE,CAAC;YAClB,MAAM,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAA;YAC7B,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAE,CAAA;YAC7B,IAAI,IAAI,CAAC,IAAI,KAAK,IAAI,EAAE,CAAC;gBACvB,OAAO,GAAG,CAAA;YACZ,CAAC;YACD,IAAI,IAAI,CAAC,IAAI,GAAG,IAAI,EAAE,CAAC;gBACrB,GAAG,GAAG,GAAG,GAAG,CAAC,CAAA;YACf,CAAC;iBAAM,CAAC;gBACN,IAAI,GAAG,GAAG,CAAA;YACZ,CAAC;QACH,CAAC;QACD,OAAO,SAAS,CAAA;IAClB,CAAC;IAED,QAAQ,CAAC,GAAQ;QACf,IAAI,GAAG,CAAC,IAAI,KAAK,SAAS,EAAE,CAAC;YAC3B,OAAO,SAAS,CAAA;QAClB,CAAC;QACD,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,CAAA;QAClC,IAAI,IAAI,KAAK,SAAS,EAAE,CAAC;YACvB,OAAO,SAAS,CAAA;QAClB,CAAC;QACD,IAAI,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAE,CAAC,MAAM,CAAA;QACvC,IAAI,QAAQ,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC;YAC1B,OAAO,SAAS,CAAA;QAClB,CAAC;QACD,IAAI,MAAM,GAAG,CAAC,CAAA;QACd,KAAK,MAAM,GAAG,IAAI,IAAI,CAAC,IAAI,EAAE,EAAE,CAAC;YAC9B,MAAM,IAAI,GAAG,CAAC,GAAG,CAAA;YACjB,IAAI,GAAG,CAAC,KAAK,KAAK,IAAI,EAAE,CAAC;gBACvB,SAAQ;YACV,CAAC;YACD,QAAQ,IAAI,GAAG,CAAC,GAAG,CAAA;YACnB,IAAI,QAAQ,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC;gBAC1B,OAAO,MAAM,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC,MAAM,CAAC,CAAA;YACzC,CAAC;QACH,CAAC;QACD,OAAO,SAAS,CAAA;IAClB,CAAC;IAED,aAAa,CAAC,CAAS;QACrB,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YAClC,IAAI,EAAE,CAAC,CAAC,IAAI,KAAK,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC;YACzD,KAAK,EAAE,CAAC;SACT,CAAC,CAAC,CAAA;QACH,KAAK,MAAM,GAAG,IAAI,IAAI,CAAC,IAAI,EAAE,EAAE,CAAC;YAC9B,MAAM,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,CAAA;YAC/B,IAAI,KAAK,EAAE,CAAC;gBACV,KAAK,CAAC,KAAK,IAAI,GAAG,CAAC,GAAG,CAAA;YACxB,CAAC;QACH,CAAC;QACD,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,MAAM,CAAC,MAAM,EAAE,IAAI,EAAE,EAAE,CAAC;YAChD,MAAM,IAAI,GAAG,MAAM,CAAC,IAAI,GAAG,CAAC,CAAoC,CAAA;YAChE,MAAM,IAAI,GAAG,MAAM,CAAC,IAAI,CAAoC,CAAA;YAC5D,IAAI,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;gBAC5C,MAAM,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,IAAI,CAAA;gBACvB,MAAM,CAAC,IAAI,CAAC,GAAG,IAAI,CAAA;YACrB,CAAC;QACH,CAAC;QACD,IAAI,MAAM,GAAG,CAAC,CAAA;QACd,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;YAC3B,MAAM,IAAI,KAAK,CAAC,KAAK,CAAA;YACrB,IAAI,MAAM,GAAG,CAAC,EAAE,CAAC;gBACf,OAAO,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAA;YAC1D,CAAC;QACH,CAAC;QACD,OAAO,SAAS,CAAA;IAClB,CAAC;IAED,UAAU;QACR,MAAM,OAAO,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,CAAA;QAC7C,MAAM,MAAM,GAAU,EAAE,CAAA;QACxB,KAAK,MAAM,GAAG,IAAI,IAAI,CAAC,IAAI,EAAE,EAAE,CAAC;YAC9B,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAA;YACjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;gBACjC,MAAM,CAAC,IAAI,CAAC,EAAE,IAAI,EAAE,IAAI,CAAC,IAAI,EAAE,MAAM,EAAE,OAAO,CAAC,GAAG,CAAC,KAAK,CAAE,EAAE,CAAC,CAAA;gBAC7D,OAAO,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,KAAK,CAAE,GAAG,CAAC,CAAA;YAC9C,CAAC;QACH,CAAC;QACD,OAAO,MAAM,CAAA;IACf,CAAC;CACF"}
@@ -30,7 +30,13 @@ export function decodeSequence(encoded) {
30
30
  }
31
31
  return result;
32
32
  }
33
- const COMPLEMENT = { A: 'T', C: 'G', G: 'C', T: 'A', N: 'N' };
33
+ const COMPLEMENT = {
34
+ A: 'T',
35
+ C: 'G',
36
+ G: 'C',
37
+ T: 'A',
38
+ N: 'N',
39
+ };
34
40
  export function reverseComplement(sequence) {
35
41
  let result = '';
36
42
  for (let i = sequence.length - 1; i >= 0; i--) {
@@ -38,3 +44,4 @@ export function reverseComplement(sequence) {
38
44
  }
39
45
  return result;
40
46
  }
47
+ //# sourceMappingURL=sequence.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"sequence.js","sourceRoot":"","sources":["../../src/gbwt/sequence.ts"],"names":[],"mappings":"AAAA,MAAM,MAAM,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAA;AAE5C,MAAM,UAAU,qBAAqB,CAAC,OAAmB;IACvD,MAAM,IAAI,GAAG,OAAO,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC,CAAC,CAAA;IACxC,IAAI,IAAI,KAAK,SAAS,EAAE,CAAC;QACvB,OAAO,CAAC,CAAA;IACV,CAAC;IACD,IAAI,KAAK,GAAG,IAAI,CAAA;IAChB,IAAI,MAAM,GAAG,CAAC,CAAA;IACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,IAAI,KAAK,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC;YACpB,MAAK;QACP,CAAC;QACD,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,CAAA;QAC7B,MAAM,IAAI,CAAC,CAAA;IACb,CAAC;IACD,OAAO,CAAC,GAAG,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,MAAM,CAAA;AAC1C,CAAC;AAED,MAAM,UAAU,cAAc,CAAC,OAAmB;IAChD,IAAI,MAAM,GAAG,EAAE,CAAA;IACf,KAAK,MAAM,IAAI,IAAI,OAAO,EAAE,CAAC;QAC3B,IAAI,KAAK,GAAG,IAAI,CAAA;QAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,IAAI,GAAG,MAAM,CAAC,KAAK,GAAG,CAAC,CAAE,CAAA;YAC/B,IAAI,IAAI,KAAK,EAAE,EAAE,CAAC;gBAChB,OAAO,MAAM,CAAA;YACf,CAAC;YACD,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,CAAA;YAC7B,MAAM,IAAI,IAAI,CAAA;QAChB,CAAC;IACH,CAAC;IACD,OAAO,MAAM,CAAA;AACf,CAAC;AAED,MAAM,UAAU,GAA2B;IACzC,CAAC,EAAE,GAAG;IACN,CAAC,EAAE,GAAG;IACN,CAAC,EAAE,GAAG;IACN,CAAC,EAAE,GAAG;IACN,CAAC,EAAE,GAAG;CACP,CAAA;AAED,MAAM,UAAU,iBAAiB,CAAC,QAAgB;IAChD,IAAI,MAAM,GAAG,EAAE,CAAA;IACf,KAAK,IAAI,CAAC,GAAG,QAAQ,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC9C,MAAM,IAAI,UAAU,CAAC,QAAQ,CAAC,CAAC,CAAE,CAAC,IAAI,GAAG,CAAA;IAC3C,CAAC;IACD,OAAO,MAAM,CAAA;AACf,CAAC"}
package/dist/index.d.ts CHANGED
@@ -1,10 +1,10 @@
1
- export { GBZBase, GbzRecord, SchemaVersionError, formatPathName, GENERIC_SAMPLE, SCHEMA_VERSION } from './db.ts';
1
+ export { GBZBase, GENERIC_SAMPLE, GbzRecord, SCHEMA_VERSION, SchemaVersionError, formatPathName, } from './db.ts';
2
2
  export type { GbzPath, HaplotypeSample, PathName } from './db.ts';
3
3
  export type { ByteSource } from './filehandle.ts';
4
4
  export type { Pos } from './gbwt/record.ts';
5
5
  export { Subgraph } from './subgraph.ts';
6
6
  export type { HaplotypeAlignment, HaplotypeOutput, PathIdentity, PathPosition, ReferencePath, SubgraphJson, SubgraphPath, ToJsonOptions, } from './subgraph.ts';
7
- export { subgraphAtOffset, subgraphInInterval, subgraphAroundNodes } from './query.ts';
7
+ export { subgraphAroundNodes, subgraphAtOffset, subgraphInInterval, } from './query.ts';
8
8
  export type { PathQuery, QueryOptions } from './query.ts';
9
9
  export { SqliteDatabase } from './sqlite/database.ts';
10
10
  export { weightedLcs } from './lcs.ts';
package/dist/index.js CHANGED
@@ -1,6 +1,7 @@
1
- export { GBZBase, GbzRecord, SchemaVersionError, formatPathName, GENERIC_SAMPLE, SCHEMA_VERSION } from "./db.js";
1
+ export { GBZBase, GENERIC_SAMPLE, GbzRecord, SCHEMA_VERSION, SchemaVersionError, formatPathName, } from "./db.js";
2
2
  export { Subgraph } from "./subgraph.js";
3
- export { subgraphAtOffset, subgraphInInterval, subgraphAroundNodes } from "./query.js";
3
+ export { subgraphAroundNodes, subgraphAtOffset, subgraphInInterval, } from "./query.js";
4
4
  export { SqliteDatabase } from "./sqlite/database.js";
5
5
  export { weightedLcs } from "./lcs.js";
6
6
  export * as nodes from "./gbwt/node.js";
7
+ //# sourceMappingURL=index.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA,OAAO,EACL,OAAO,EACP,cAAc,EACd,SAAS,EACT,cAAc,EACd,kBAAkB,EAClB,cAAc,GACf,MAAM,SAAS,CAAA;AAIhB,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAA;AAWxC,OAAO,EACL,mBAAmB,EACnB,gBAAgB,EAChB,kBAAkB,GACnB,MAAM,YAAY,CAAA;AAEnB,OAAO,EAAE,cAAc,EAAE,MAAM,sBAAsB,CAAA;AACrD,OAAO,EAAE,WAAW,EAAE,MAAM,UAAU,CAAA;AACtC,OAAO,KAAK,KAAK,MAAM,gBAAgB,CAAA"}
package/dist/lcs.js CHANGED
@@ -54,10 +54,10 @@ class MinHeap {
54
54
  }
55
55
  }
56
56
  class Matrix {
57
- a;
58
- b;
59
57
  aSums;
60
58
  bSums;
59
+ a;
60
+ b;
61
61
  points = new Map();
62
62
  pendingEdits = new MinHeap();
63
63
  constructor(a, b, weight) {
@@ -103,7 +103,9 @@ class Matrix {
103
103
  continue;
104
104
  }
105
105
  const point = { ...found };
106
- while (point.a < this.a.length && point.b < this.b.length && this.a[point.a] === this.b[point.b]) {
106
+ while (point.a < this.a.length &&
107
+ point.b < this.b.length &&
108
+ this.a[point.a] === this.b[point.b]) {
107
109
  point.weight += 2 * this.aWeight(point.a);
108
110
  point.a += 1;
109
111
  point.b += 1;
@@ -117,27 +119,44 @@ class Matrix {
117
119
  }
118
120
  if (point.a < this.a.length) {
119
121
  const w = this.aWeight(point.a);
120
- this.tryInsert(edits + w, diagonal + w, { weight: point.weight, a: point.a + 1, b: point.b, matches: 0 });
122
+ this.tryInsert(edits + w, diagonal + w, {
123
+ weight: point.weight,
124
+ a: point.a + 1,
125
+ b: point.b,
126
+ matches: 0,
127
+ });
121
128
  }
122
129
  if (point.b < this.b.length) {
123
130
  const w = this.bWeight(point.b);
124
- this.tryInsert(edits + w, diagonal - w, { weight: point.weight, a: point.a, b: point.b + 1, matches: 0 });
131
+ this.tryInsert(edits + w, diagonal - w, {
132
+ weight: point.weight,
133
+ a: point.a,
134
+ b: point.b + 1,
135
+ matches: 0,
136
+ });
125
137
  }
126
138
  }
127
139
  return undefined;
128
140
  }
129
141
  nextEdits(edits) {
130
- while (this.pendingEdits.peek() !== undefined && this.pendingEdits.peek() <= edits) {
142
+ while (this.pendingEdits.peek() !== undefined &&
143
+ this.pendingEdits.peek() <= edits) {
131
144
  this.pendingEdits.pop();
132
145
  }
133
146
  return this.pendingEdits.peek();
134
147
  }
135
148
  predecessor(a, b, edits) {
136
149
  const diagonal = this.aSums[a] - this.bSums[b];
137
- const prev = a > 0 && this.aWeight(a - 1) <= edits ? this.get(edits - this.aWeight(a - 1), diagonal - this.aWeight(a - 1)) : undefined;
138
- const next = b > 0 && this.bWeight(b - 1) <= edits ? this.get(edits - this.bWeight(b - 1), diagonal + this.bWeight(b - 1)) : undefined;
150
+ const prev = a > 0 && this.aWeight(a - 1) <= edits
151
+ ? this.get(edits - this.aWeight(a - 1), diagonal - this.aWeight(a - 1))
152
+ : undefined;
153
+ const next = b > 0 && this.bWeight(b - 1) <= edits
154
+ ? this.get(edits - this.bWeight(b - 1), diagonal + this.bWeight(b - 1))
155
+ : undefined;
139
156
  if (prev && next) {
140
- return prev.weight > next.weight ? [prev, edits - this.aWeight(a - 1)] : [next, edits - this.bWeight(b - 1)];
157
+ return prev.weight > next.weight
158
+ ? [prev, edits - this.aWeight(a - 1)]
159
+ : [next, edits - this.bWeight(b - 1)];
141
160
  }
142
161
  if (prev) {
143
162
  return [prev, edits - this.aWeight(a - 1)];
@@ -154,7 +173,9 @@ export function weightedLcs(a, b, weight) {
154
173
  prefix += 1;
155
174
  }
156
175
  let suffix = 0;
157
- while (suffix < a.length - prefix && suffix < b.length - prefix && a[a.length - 1 - suffix] === b[b.length - 1 - suffix]) {
176
+ while (suffix < a.length - prefix &&
177
+ suffix < b.length - prefix &&
178
+ a[a.length - 1 - suffix] === b[b.length - 1 - suffix]) {
158
179
  suffix += 1;
159
180
  }
160
181
  const pairs = [];
@@ -212,3 +233,4 @@ function weightedLcsCore(a, b, weight) {
212
233
  result.reverse();
213
234
  return [result, finalWeight];
214
235
  }
236
+ //# sourceMappingURL=lcs.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"lcs.js","sourceRoot":"","sources":["../src/lcs.ts"],"names":[],"mappings":"AAOA,SAAS,UAAU,CAAC,QAAkB,EAAE,MAA6B;IACnE,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC,CAAA;IAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACzC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAE,GAAG,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAE,CAAC,CAAC,CAAA;IAC5C,CAAC;IACD,OAAO,IAAI,CAAA;AACb,CAAC;AAED,MAAM,OAAO;IACH,KAAK,GAAa,EAAE,CAAA;IAE5B,IAAI,CAAC,KAAa;QAChB,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAA;QACxB,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAA;QACjB,IAAI,CAAC,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,CAAA;QACxB,OAAO,CAAC,GAAG,CAAC,EAAE,CAAC;YACb,MAAM,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAA;YAC3B,IAAI,KAAK,CAAC,MAAM,CAAE,IAAI,KAAK,EAAE,CAAC;gBAC5B,MAAK;YACP,CAAC;YACD,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAE,CAAA;YACzB,CAAC,GAAG,MAAM,CAAA;QACZ,CAAC;QACD,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAA;IAClB,CAAC;IAED,IAAI;QACF,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAA;IACtB,CAAC;IAED,GAAG;QACD,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAA;QACxB,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAA;QACpB,MAAM,IAAI,GAAG,KAAK,CAAC,GAAG,EAAE,CAAA;QACxB,IAAI,KAAK,CAAC,MAAM,GAAG,CAAC,IAAI,IAAI,KAAK,SAAS,EAAE,CAAC;YAC3C,IAAI,CAAC,GAAG,CAAC,CAAA;YACT,SAAS,CAAC;gBACR,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;gBACtB,MAAM,KAAK,GAAG,IAAI,GAAG,CAAC,CAAA;gBACtB,IAAI,QAAQ,GAAG,CAAC,CAAA;gBAChB,IAAI,KAAK,GAAG,IAAI,CAAA;gBAChB,IAAI,IAAI,GAAG,KAAK,CAAC,MAAM,IAAI,KAAK,CAAC,IAAI,CAAE,GAAG,KAAK,EAAE,CAAC;oBAChD,QAAQ,GAAG,IAAI,CAAA;oBACf,KAAK,GAAG,KAAK,CAAC,IAAI,CAAE,CAAA;gBACtB,CAAC;gBACD,IAAI,KAAK,GAAG,KAAK,CAAC,MAAM,IAAI,KAAK,CAAC,KAAK,CAAE,GAAG,KAAK,EAAE,CAAC;oBAClD,QAAQ,GAAG,KAAK,CAAA;gBAClB,CAAC;gBACD,IAAI,QAAQ,KAAK,CAAC,EAAE,CAAC;oBACnB,MAAK;gBACP,CAAC;gBACD,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAE,CAAA;gBAC3B,CAAC,GAAG,QAAQ,CAAA;YACd,CAAC;YACD,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAA;QACjB,CAAC;QACD,OAAO,GAAG,CAAA;IACZ,CAAC;CACF;AAED,MAAM,MAAM;IACD,KAAK,CAAU;IACf,KAAK,CAAU;IACf,CAAC,CAAU;IACX,CAAC,CAAU;IACX,MAAM,GAAG,IAAI,GAAG,EAA8B,CAAA;IAC/C,YAAY,GAAG,IAAI,OAAO,EAAE,CAAA;IAEpC,YAAY,CAAW,EAAE,CAAW,EAAE,MAA6B;QACjE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAA;QACV,IAAI,CAAC,CAAC,GAAG,CAAC,CAAA;QACV,IAAI,CAAC,KAAK,GAAG,UAAU,CAAC,CAAC,EAAE,MAAM,CAAC,CAAA;QAClC,IAAI,CAAC,KAAK,GAAG,UAAU,CAAC,CAAC,EAAE,MAAM,CAAC,CAAA;QAClC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,CAAC,CAAA;IACvD,CAAC;IAEO,GAAG,CAAC,KAAa,EAAE,QAAgB,EAAE,KAAY;QACvD,IAAI,GAAG,GAAG,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,CAAA;QAChC,IAAI,CAAC,GAAG,EAAE,CAAC;YACT,GAAG,GAAG,IAAI,GAAG,EAAE,CAAA;YACf,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,KAAK,EAAE,GAAG,CAAC,CAAA;YAC3B,IAAI,CAAC,YAAY,CAAC,IAAI,CAAC,KAAK,CAAC,CAAA;QAC/B,CAAC;QACD,GAAG,CAAC,GAAG,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAA;IAC1B,CAAC;IAED,GAAG,CAAC,KAAa,EAAE,QAAgB;QACjC,OAAO,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAA;IAC9C,CAAC;IAEO,SAAS,CAAC,KAAa,EAAE,QAAgB,EAAE,KAAY;QAC7D,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,QAAQ,CAAC,CAAA;QAC1C,IAAI,CAAC,QAAQ,IAAI,KAAK,CAAC,MAAM,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC;YAChD,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,QAAQ,EAAE,KAAK,CAAC,CAAA;QAClC,CAAC;IACH,CAAC;IAED,OAAO,CAAC,MAAc;QACpB,OAAO,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAE,CAAA;IACtD,CAAC;IAED,OAAO,CAAC,MAAc;QACpB,OAAO,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAE,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAE,CAAA;IACtD,CAAC;IAED,MAAM,CAAC,KAAa;QAClB,MAAM,GAAG,GAAG,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,CAAA;QAClC,IAAI,CAAC,GAAG,EAAE,CAAC;YACT,OAAO,SAAS,CAAA;QAClB,CAAC;QACD,MAAM,SAAS,GAAG,CAAC,GAAG,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;QACvD,KAAK,MAAM,QAAQ,IAAI,SAAS,EAAE,CAAC;YACjC,MAAM,KAAK,GAAG,GAAG,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAA;YAC/B,IAAI,CAAC,KAAK,EAAE,CAAC;gBACX,SAAQ;YACV,CAAC;YACD,MAAM,KAAK,GAAG,EAAE,GAAG,KAAK,EAAE,CAAA;YAC1B,OACE,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,MAAM;gBACvB,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,MAAM;gBACvB,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EACnC,CAAC;gBACD,KAAK,CAAC,MAAM,IAAI,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAA;gBACzC,KAAK,CAAC,CAAC,IAAI,CAAC,CAAA;gBACZ,KAAK,CAAC,CAAC,IAAI,CAAC,CAAA;gBACZ,KAAK,CAAC,OAAO,IAAI,CAAC,CAAA;YACpB,CAAC;YACD,IAAI,KAAK,CAAC,OAAO,GAAG,CAAC,EAAE,CAAC;gBACtB,GAAG,CAAC,GAAG,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAA;YAC1B,CAAC;YACD,IAAI,KAAK,CAAC,CAAC,KAAK,IAAI,CAAC,CAAC,CAAC,MAAM,IAAI,KAAK,CAAC,CAAC,KAAK,IAAI,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC;gBAC3D,OAAO,KAAK,CAAA;YACd,CAAC;YACD,IAAI,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC;gBAC5B,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAA;gBAC/B,IAAI,CAAC,SAAS,CAAC,KAAK,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC,EAAE;oBACtC,MAAM,EAAE,KAAK,CAAC,MAAM;oBACpB,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC;oBACd,CAAC,EAAE,KAAK,CAAC,CAAC;oBACV,OAAO,EAAE,CAAC;iBACX,CAAC,CAAA;YACJ,CAAC;YACD,IAAI,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC;gBAC5B,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAA;gBAC/B,IAAI,CAAC,SAAS,CAAC,KAAK,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC,EAAE;oBACtC,MAAM,EAAE,KAAK,CAAC,MAAM;oBACpB,CAAC,EAAE,KAAK,CAAC,CAAC;oBACV,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC;oBACd,OAAO,EAAE,CAAC;iBACX,CAAC,CAAA;YACJ,CAAC;QACH,CAAC;QACD,OAAO,SAAS,CAAA;IAClB,CAAC;IAED,SAAS,CAAC,KAAa;QACrB,OACE,IAAI,CAAC,YAAY,CAAC,IAAI,EAAE,KAAK,SAAS;YACtC,IAAI,CAAC,YAAY,CAAC,IAAI,EAAG,IAAI,KAAK,EAClC,CAAC;YACD,IAAI,CAAC,YAAY,CAAC,GAAG,EAAE,CAAA;QACzB,CAAC;QACD,OAAO,IAAI,CAAC,YAAY,CAAC,IAAI,EAAE,CAAA;IACjC,CAAC;IAED,WAAW,CACT,CAAS,EACT,CAAS,EACT,KAAa;QAEb,MAAM,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAE,CAAA;QAChD,MAAM,IAAI,GACR,CAAC,GAAG,CAAC,IAAI,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,KAAK;YACnC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,QAAQ,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACvE,CAAC,CAAC,SAAS,CAAA;QACf,MAAM,IAAI,GACR,CAAC,GAAG,CAAC,IAAI,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,KAAK;YACnC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,QAAQ,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACvE,CAAC,CAAC,SAAS,CAAA;QACf,IAAI,IAAI,IAAI,IAAI,EAAE,CAAC;YACjB,OAAO,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,MAAM;gBAC9B,CAAC,CAAC,CAAC,IAAI,EAAE,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBACrC,CAAC,CAAC,CAAC,IAAI,EAAE,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;QACzC,CAAC;QACD,IAAI,IAAI,EAAE,CAAC;YACT,OAAO,CAAC,IAAI,EAAE,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;QAC5C,CAAC;QACD,IAAI,IAAI,EAAE,CAAC;YACT,OAAO,CAAC,IAAI,EAAE,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;QAC5C,CAAC;QACD,OAAO,SAAS,CAAA;IAClB,CAAC;CACF;AAED,MAAM,UAAU,WAAW,CACzB,CAAW,EACX,CAAW,EACX,MAA6B;IAE7B,IAAI,MAAM,GAAG,CAAC,CAAA;IACd,OAAO,MAAM,GAAG,CAAC,CAAC,MAAM,IAAI,MAAM,GAAG,CAAC,CAAC,MAAM,IAAI,CAAC,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,MAAM,CAAC,EAAE,CAAC;QACzE,MAAM,IAAI,CAAC,CAAA;IACb,CAAC;IACD,IAAI,MAAM,GAAG,CAAC,CAAA;IACd,OACE,MAAM,GAAG,CAAC,CAAC,MAAM,GAAG,MAAM;QAC1B,MAAM,GAAG,CAAC,CAAC,MAAM,GAAG,MAAM;QAC1B,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,GAAG,MAAM,CAAC,EACrD,CAAC;QACD,MAAM,IAAI,CAAC,CAAA;IACb,CAAC;IACD,MAAM,KAAK,GAAuB,EAAE,CAAA;IACpC,IAAI,KAAK,GAAG,CAAC,CAAA;IACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAA;QAClB,KAAK,IAAI,MAAM,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC,CAAA;IACxB,CAAC;IACD,MAAM,CAAC,MAAM,EAAE,YAAY,CAAC,GAAG,eAAe,CAC5C,CAAC,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC,CAAC,MAAM,GAAG,MAAM,CAAC,EAClC,CAAC,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC,CAAC,MAAM,GAAG,MAAM,CAAC,EAClC,MAAM,CACP,CAAA;IACD,KAAK,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,MAAM,EAAE,CAAC;QAC5B,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC,CAAA;IACtC,CAAC;IACD,KAAK,IAAI,YAAY,CAAA;IACrB,KAAK,IAAI,CAAC,GAAG,MAAM,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAA;QACxC,KAAK,IAAI,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAE,CAAC,CAAA;IACnC,CAAC;IACD,OAAO,CAAC,KAAK,EAAE,KAAK,CAAC,CAAA;AACvB,CAAC;AAED,SAAS,eAAe,CACtB,CAAW,EACX,CAAW,EACX,MAA6B;IAE7B,IAAI,CAAC,CAAC,MAAM,KAAK,CAAC,IAAI,CAAC,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QACrC,OAAO,CAAC,EAAE,EAAE,CAAC,CAAC,CAAA;IAChB,CAAC;IACD,MAAM,MAAM,GAAG,IAAI,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAA;IACvC,IAAI,KAAK,GAAG,CAAC,CAAA;IACb,IAAI,KAAK,GAAU,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,CAAA;IACxD,SAAS,CAAC;QACR,MAAM,GAAG,GAAG,MAAM,CAAC,MAAM,CAAC,KAAK,CAAC,CAAA;QAChC,IAAI,GAAG,EAAE,CAAC;YACR,KAAK,GAAG,GAAG,CAAA;YACX,MAAK;QACP,CAAC;QACD,MAAM,IAAI,GAAG,MAAM,CAAC,SAAS,CAAC,KAAK,CAAC,CAAA;QACpC,IAAI,IAAI,KAAK,SAAS,EAAE,CAAC;YACvB,MAAK;QACP,CAAC;QACD,KAAK,GAAG,IAAI,CAAA;IACd,CAAC;IACD,MAAM,MAAM,GAAuB,EAAE,CAAA;IACrC,MAAM,WAAW,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,CAAA;IACpC,KAAK,GAAG,EAAE,GAAG,KAAK,EAAE,CAAA;IACpB,SAAS,CAAC;QACR,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;YACvC,KAAK,CAAC,CAAC,IAAI,CAAC,CAAA;YACZ,KAAK,CAAC,CAAC,IAAI,CAAC,CAAA;YACZ,MAAM,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAA;QACjC,CAAC;QACD,MAAM,IAAI,GAAG,MAAM,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAA;QACxD,IAAI,CAAC,IAAI,EAAE,CAAC;YACV,MAAK;QACP,CAAC;QACD,KAAK,GAAG,EAAE,GAAG,IAAI,CAAC,CAAC,CAAC,EAAE,CAAA;QACtB,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC,CAAA;IACjB,CAAC;IACD,MAAM,CAAC,OAAO,EAAE,CAAA;IAChB,OAAO,CAAC,MAAM,EAAE,WAAW,CAAC,CAAA;AAC9B,CAAC"}
package/dist/query.d.ts CHANGED
@@ -1,5 +1,5 @@
1
- import type { GBZBase } from './db.ts';
2
1
  import { Subgraph } from './subgraph.ts';
2
+ import type { GBZBase } from './db.ts';
3
3
  import type { HaplotypeOutput } from './subgraph.ts';
4
4
  export interface QueryOptions {
5
5
  context?: number;
package/dist/query.js CHANGED
@@ -35,3 +35,4 @@ export async function subgraphAroundNodes(db, nodes, opts = {}) {
35
35
  subgraph.extractPaths(undefined, haplotypes);
36
36
  return subgraph;
37
37
  }
38
+ //# sourceMappingURL=query.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"query.js","sourceRoot":"","sources":["../src/query.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,cAAc,EAAE,MAAM,SAAS,CAAA;AACxC,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAA;AAiBxC,SAAS,QAAQ,CAAC,KAAgB,EAAE,QAAgB;IAClD,OAAO;QACL,MAAM,EAAE,KAAK,CAAC,MAAM,IAAI,cAAc;QACtC,MAAM,EAAE,KAAK,CAAC,MAAM;QACpB,SAAS,EAAE,KAAK,CAAC,SAAS,IAAI,CAAC;QAC/B,QAAQ;KACT,CAAA;AACH,CAAC;AAED,MAAM,CAAC,KAAK,UAAU,gBAAgB,CACpC,EAAW,EACX,KAAgB,EAChB,MAAc,EACd,OAAqB,EAAE;IAEvB,MAAM,QAAQ,GAAG,IAAI,QAAQ,CAAC,EAAE,CAAC,CAAA;IACjC,QAAQ,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,CAAA;IAC3B,MAAM,SAAS,GAAG,MAAM,QAAQ,CAAC,YAAY,CAAC,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,CAAA;IACtE,MAAM,QAAQ,CAAC,cAAc,CAC3B,SAAS,CAAC,QAAQ,CAAC,MAAM,EACzB,SAAS,CAAC,QAAQ,CAAC,UAAU,EAC7B,IAAI,CAAC,OAAO,IAAI,GAAG,CACpB,CAAA;IACD,QAAQ,CAAC,YAAY,CAAC,SAAS,EAAE,IAAI,CAAC,UAAU,IAAI,KAAK,CAAC,CAAA;IAC1D,OAAO,QAAQ,CAAA;AACjB,CAAC;AAED,MAAM,CAAC,KAAK,UAAU,kBAAkB,CACtC,EAAW,EACX,KAAgB,EAChB,KAAa,EACb,GAAW,EACX,OAAqB,EAAE;IAEvB,MAAM,QAAQ,GAAG,IAAI,QAAQ,CAAC,EAAE,CAAC,CAAA;IACjC,QAAQ,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,CAAA;IAC3B,MAAM,SAAS,GAAG,MAAM,QAAQ,CAAC,YAAY,CAAC,QAAQ,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC,CAAA;IACrE,MAAM,QAAQ,CAAC,cAAc,CAC3B,SAAS,CAAC,QAAQ,EAClB,GAAG,GAAG,KAAK,EACX,IAAI,CAAC,OAAO,IAAI,GAAG,CACpB,CAAA;IACD,QAAQ,CAAC,YAAY,CAAC,SAAS,EAAE,IAAI,CAAC,UAAU,IAAI,KAAK,CAAC,CAAA;IAC1D,OAAO,QAAQ,CAAA;AACjB,CAAC;AAED,MAAM,CAAC,KAAK,UAAU,mBAAmB,CACvC,EAAW,EACX,KAAe,EACf,OAAqB,EAAE;IAEvB,MAAM,UAAU,GAAG,IAAI,CAAC,UAAU,IAAI,KAAK,CAAA;IAC3C,IAAI,UAAU,KAAK,gBAAgB,EAAE,CAAC;QACpC,MAAM,IAAI,KAAK,CAAC,sDAAsD,CAAC,CAAA;IACzE,CAAC;IACD,MAAM,QAAQ,GAAG,IAAI,QAAQ,CAAC,EAAE,CAAC,CAAA;IACjC,QAAQ,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,CAAA;IAC3B,MAAM,QAAQ,CAAC,WAAW,CAAC,KAAK,EAAE,IAAI,CAAC,OAAO,IAAI,GAAG,CAAC,CAAA;IACtD,QAAQ,CAAC,YAAY,CAAC,SAAS,EAAE,UAAU,CAAC,CAAA;IAC5C,OAAO,QAAQ,CAAA;AACjB,CAAC"}
@@ -1,8 +1,8 @@
1
1
  import type { Pager } from './pager.ts';
2
2
  import type { SqlValue } from './record.ts';
3
3
  export declare class BTree {
4
- private pager;
5
4
  private readonly usable;
5
+ private pager;
6
6
  constructor(pager: Pager, reservedBytes: number);
7
7
  private localPayloadSize;
8
8
  private payload;