@gmod/cram 1.5.9 → 1.6.2

package/esm/cramFile/slice/index.js

@@ -0,0 +1,364 @@
+import { CramMalformedError, CramBufferOverrunError, CramArgumentError, } from '../../errors';
+import { parseItem, tinyMemoize, sequenceMD5 } from '../util';
+import Constants from '../constants';
+import decodeRecord from './decodeRecord';
+/**
+ * @private
+ * Try to estimate the template length from a bunch of interrelated multi-segment reads.
+ * @param {Array[CramRecord]} allRecords
+ * @param {number} currentRecordNumber
+ * @param {CramRecord} thisRecord
+ */
+function calculateMultiSegmentMatedTemplateLength(allRecords, currentRecordNumber, thisRecord) {
+    function getAllMatedRecords(startRecord) {
+        const records = [startRecord];
+        if (startRecord.mateRecordNumber >= 0) {
+            const mateRecord = allRecords[startRecord.mateRecordNumber];
+            if (!mateRecord) {
+                throw new CramMalformedError('intra-slice mate record not found, this file seems malformed');
+            }
+            records.push(...getAllMatedRecords(mateRecord));
+        }
+        return records;
+    }
+    const matedRecords = getAllMatedRecords(thisRecord);
+    const starts = matedRecords.map(r => r.alignmentStart);
+    const ends = matedRecords.map(r => r.alignmentStart + r.readLength - 1);
+    const estimatedTemplateLength = Math.max(...ends) - Math.min(...starts) + 1;
+    if (estimatedTemplateLength >= 0) {
+        matedRecords.forEach(r => {
+            if (r.templateLength !== undefined) {
+                throw new CramMalformedError('mate pair group has some members that have template lengths already, this file seems malformed');
+            }
+            r.templateLength = estimatedTemplateLength;
+        });
+    }
+}
+/**
+ * @private
+ * Attempt to calculate the `templateLength` for a pair of intra-slice paired reads.
+ * Ported from htslib. Algorithm is imperfect.
+ * @param {CramRecord} thisRecord
+ * @param {CramRecord} mateRecord
+ */
+function calculateIntraSliceMatePairTemplateLength(thisRecord, mateRecord) {
+    // this just estimates the template length by using the simple (non-gapped) end coordinate of each
+    // read, because gapping in the alignment doesn't mean the template is longer or shorter
+    const start = Math.min(thisRecord.alignmentStart, mateRecord.alignmentStart);
+    const end = Math.max(thisRecord.alignmentStart + thisRecord.readLength - 1, mateRecord.alignmentStart + mateRecord.readLength - 1);
+    const lengthEstimate = end - start + 1;
+    thisRecord.templateLength = lengthEstimate;
+    mateRecord.templateLength = lengthEstimate;
+}
+/**
+ * @private establishes a mate-pair relationship between two records in the same slice.
+ * CRAM compresses mate-pair relationships between records in the same slice down into
+ * just one record having the index in the slice of its mate
+ */
+function associateIntraSliceMate(allRecords, currentRecordNumber, thisRecord, mateRecord) {
+    if (!mateRecord) {
+        throw new CramMalformedError('could not resolve intra-slice mate pairs, file seems truncated or malformed');
+    }
+    const complicatedMultiSegment = !!(mateRecord.mate ||
+        (mateRecord.mateRecordNumber !== undefined &&
+            mateRecord.mateRecordNumber !== currentRecordNumber));
+    // Deal with lossy read names
+    if (!thisRecord.readName) {
+        thisRecord.readName = String(thisRecord.uniqueId);
+        mateRecord.readName = thisRecord.readName;
+    }
+    thisRecord.mate = {
+        sequenceId: mateRecord.sequenceId,
+        alignmentStart: mateRecord.alignmentStart,
+        uniqueId: mateRecord.uniqueId,
+    };
+    if (mateRecord.readName) {
+        thisRecord.mate.readName = mateRecord.readName;
+    }
+    // the mate record might have its own mate pointer, if this is some kind of
+    // multi-segment (more than paired) scheme, so only relate that one back to this one
+    // if it does not have any other relationship
+    if (!mateRecord.mate && mateRecord.mateRecordNumber === undefined) {
+        mateRecord.mate = {
+            sequenceId: thisRecord.sequenceId,
+            alignmentStart: thisRecord.alignmentStart,
+            uniqueId: thisRecord.uniqueId,
+        };
+        if (thisRecord.readName) {
+            mateRecord.mate.readName = thisRecord.readName;
+        }
+    }
+    // make sure the proper flags and cramFlags are set on both records
+    // paired
+    thisRecord.flags |= Constants.BAM_FPAIRED;
+    // set mate unmapped if needed
+    if (mateRecord.flags & Constants.BAM_FUNMAP) {
+        thisRecord.flags |= Constants.BAM_FMUNMAP;
+        // thisRecord.templateLength = 0
+    }
+    if (thisRecord.flags & Constants.BAM_FUNMAP) {
+        // thisRecord.templateLength = 0
+        mateRecord.flags |= Constants.BAM_FMUNMAP;
+    }
+    // set mate reversed if needed
+    if (mateRecord.flags & Constants.BAM_FREVERSE) {
+        thisRecord.flags |= Constants.BAM_FMREVERSE;
+    }
+    if (thisRecord.flags & Constants.BAM_FREVERSE) {
+        mateRecord.flags |= Constants.BAM_FMREVERSE;
+    }
+    if (thisRecord.templateLength === undefined) {
+        if (complicatedMultiSegment) {
+            calculateMultiSegmentMatedTemplateLength(allRecords, currentRecordNumber, thisRecord);
+        }
+        else {
+            calculateIntraSliceMatePairTemplateLength(thisRecord, mateRecord);
+        }
+    }
+    // delete this last because it's used by the
+    // complicated template length estimation
+    delete thisRecord.mateRecordNumber;
+}
+export default class CramSlice {
+    constructor(container, position) {
+        this.container = container;
+        this.file = container.file;
+        this.containerPosition = position;
+    }
+    // memoize
+    async getHeader() {
+        // fetch and parse the slice header
+        const sectionParsers = await this.file.getSectionParsers();
+        const containerHeader = await this.container.getHeader();
+        const header = await this.file.readBlock(containerHeader._endPosition + this.containerPosition);
+        if (header.contentType === 'MAPPED_SLICE_HEADER') {
+            header.content = parseItem(header.content, sectionParsers.cramMappedSliceHeader.parser, 0, containerHeader._endPosition);
+        }
+        else if (header.contentType === 'UNMAPPED_SLICE_HEADER') {
+            header.content = parseItem(header.content, sectionParsers.cramUnmappedSliceHeader.parser, 0, containerHeader._endPosition);
+        }
+        else {
+            throw new CramMalformedError(`error reading slice header block, invalid content type ${header._contentType}`);
+        }
+        return header;
+    }
+    // memoize
+    async getBlocks() {
+        const header = await this.getHeader();
+        // read all the blocks into memory and store them
+        let blockPosition = header._endPosition;
+        const blocks = new Array(header.content.numBlocks);
+        for (let i = 0; i < blocks.length; i += 1) {
+            blocks[i] = await this.file.readBlock(blockPosition);
+            blockPosition = blocks[i]._endPosition;
+        }
+        return blocks;
+    }
+    // no memoize
+    async getCoreDataBlock() {
+        const blocks = await this.getBlocks();
+        // the core data block is always the first block in the slice
+        return blocks[0];
+    }
+    // memoize
+    async _getBlocksContentIdIndex() {
+        const blocks = await this.getBlocks();
+        const blocksByContentId = {};
+        blocks.forEach(block => {
+            if (block.contentType === 'EXTERNAL_DATA') {
+                blocksByContentId[block.contentId] = block;
+            }
+        });
+        return blocksByContentId;
+    }
+    async getBlockByContentId(id) {
+        const blocksByContentId = await this._getBlocksContentIdIndex();
+        return blocksByContentId[id];
+    }
+    async getReferenceRegion() {
+        // read the slice header
+        const sliceHeader = (await this.getHeader()).content;
+        if (sliceHeader.refSeqId < 0) {
+            return undefined;
+        }
+        const compressionScheme = await this.container.getCompressionScheme();
+        // console.log(JSON.stringify(sliceHeader, null, '  '))
+        if (sliceHeader.refBaseBlockId >= 0) {
+            const refBlock = this.getBlockByContentId(sliceHeader.refBaseBlockId);
+            if (!refBlock) {
+                throw new CramMalformedError('embedded reference specified, but reference block does not exist');
+            }
+            if (sliceHeader.span > refBlock.uncompressedSize) {
+                throw new CramMalformedError('Embedded reference is too small');
+            }
+            return {
+                seq: refBlock.data.toString('utf8'),
+                start: sliceHeader.refSeqStart,
+                end: sliceHeader.refSeqStart + sliceHeader.refSeqSpan - 1,
+                span: sliceHeader.refSeqSpan,
+            };
+        }
+        if (compressionScheme.referenceRequired ||
+            this.file.fetchReferenceSequenceCallback) {
+            if (!this.file.fetchReferenceSequenceCallback) {
+                throw new Error('reference sequence not embedded, and seqFetch callback not provided, cannot fetch reference sequence');
+            }
+            const seq = await this.file.fetchReferenceSequenceCallback(sliceHeader.refSeqId, sliceHeader.refSeqStart, sliceHeader.refSeqStart + sliceHeader.refSeqSpan - 1);
+            if (seq.length !== sliceHeader.refSeqSpan) {
+                throw new CramArgumentError('seqFetch callback returned a reference sequence of the wrong length');
+            }
+            return {
+                seq,
+                start: sliceHeader.refSeqStart,
+                end: sliceHeader.refSeqStart + sliceHeader.refSeqSpan - 1,
+                span: sliceHeader.refSeqSpan,
+            };
+        }
+        return undefined;
+    }
+    getAllRecords() {
+        return this.getRecords(() => true);
+    }
+    async _fetchRecords() {
+        const { majorVersion } = await this.file.getDefinition();
+        const compressionScheme = await this.container.getCompressionScheme();
+        const sliceHeader = await this.getHeader();
+        const blocksByContentId = await this._getBlocksContentIdIndex();
+        // check MD5 of reference if available
+        if (majorVersion > 1 &&
+            this.file.options.checkSequenceMD5 &&
+            sliceHeader.content.refSeqId >= 0 &&
+            sliceHeader.content.md5.join('') !== '0000000000000000') {
+            const refRegion = await this.getReferenceRegion();
+            if (refRegion) {
+                const { seq, start, end } = refRegion;
+                const seqMd5 = sequenceMD5(seq);
+                const storedMd5 = sliceHeader.content.md5
+                    .map(byte => (byte < 16 ? '0' : '') + byte.toString(16))
+                    .join('');
+                if (seqMd5 !== storedMd5) {
+                    throw new CramMalformedError(`MD5 checksum reference mismatch for ref ${sliceHeader.content.refSeqId} pos ${start}..${end}. recorded MD5: ${storedMd5}, calculated MD5: ${seqMd5}`);
+                }
+            }
+        }
+        // tracks the read position within the block. codec.decode() methods
+        // advance the byte and bit positions in the cursor as they decode data
+        // note that we are only decoding a single block here, the core data block
+        const coreDataBlock = await this.getCoreDataBlock();
+        const cursors = {
+            lastAlignmentStart: sliceHeader.content.refSeqStart || 0,
+            coreBlock: { bitPosition: 7, bytePosition: 0 },
+            externalBlocks: {
+                getCursor(contentId) {
+                    if (!this[contentId]) {
+                        this[contentId] = { bitPosition: 7, bytePosition: 0 };
+                    }
+                    return this[contentId];
+                },
+            },
+        };
+        const decodeDataSeries = dataSeriesName => {
+            const codec = compressionScheme.getCodecForDataSeries(dataSeriesName);
+            if (!codec) {
+                throw new CramMalformedError(`no codec defined for ${dataSeriesName} data series`);
+            }
+            // console.log(dataSeriesName, Object.getPrototypeOf(codec))
+            return codec.decode(this, coreDataBlock, blocksByContentId, cursors);
+        };
+        let records = new Array(sliceHeader.content.numRecords);
+        for (let i = 0; i < records.length; i += 1) {
+            try {
+                records[i] = decodeRecord(this, decodeDataSeries, compressionScheme, sliceHeader, coreDataBlock, blocksByContentId, cursors, majorVersion, i);
+                records[i].uniqueId =
+                    sliceHeader.contentPosition +
+                        sliceHeader.content.recordCounter +
+                        i +
+                        1;
+            }
+            catch (e) {
+                if (e instanceof CramBufferOverrunError) {
+                    console.warn('read attempted beyond end of buffer, file seems truncated.');
+                    records = records.filter(r => !!r);
+                    break;
+                }
+                else {
+                    throw e;
+                }
+            }
+        }
+        // interpret `recordsToNextFragment` attributes to make standard `mate` objects
+        // Resolve mate pair cross-references between records in this slice
+        for (let i = 0; i < records.length; i += 1) {
+            const { mateRecordNumber } = records[i];
+            if (mateRecordNumber >= 0) {
+                associateIntraSliceMate(records, i, records[i], records[mateRecordNumber]);
+            }
+        }
+        return records;
+    }
+    async getRecords(filterFunction) {
+        // fetch the features if necessary, using the file-level feature cache
+        const cacheKey = this.container.filePosition + this.containerPosition;
+        let recordsPromise = this.file.featureCache.get(cacheKey);
+        if (!recordsPromise) {
+            recordsPromise = this._fetchRecords();
+            this.file.featureCache.set(cacheKey, recordsPromise);
+        }
+        const records = (await recordsPromise).filter(filterFunction);
+        // if we can fetch reference sequence, add the reference sequence to the records
+        if (records.length && this.file.fetchReferenceSequenceCallback) {
+            const sliceHeader = await this.getHeader();
+            if (sliceHeader.content.refSeqId >= 0 || // single-ref slice
+                sliceHeader.content.refSeqId === -2 // multi-ref slice
+            ) {
+                const singleRefId = sliceHeader.content.refSeqId >= 0
+                    ? sliceHeader.content.refSeqId
+                    : undefined;
+                const compressionScheme = await this.container.getCompressionScheme();
+                const refRegions = {}; // seqId => { start, end, seq }
+                // iterate over the records to find the spans of the reference sequences we need to fetch
+                for (let i = 0; i < records.length; i += 1) {
+                    const seqId = singleRefId !== undefined ? singleRefId : records[i].sequenceId;
+                    let refRegion = refRegions[seqId];
+                    if (!refRegion) {
+                        refRegion = {
+                            id: seqId,
+                            start: records[i].alignmentStart,
+                            end: -Infinity,
+                        };
+                        refRegions[seqId] = refRegion;
+                    }
+                    const end = records[i].alignmentStart +
+                        (records[i].lengthOnRef || records[i].readLength) -
+                        1;
+                    if (end > refRegion.end) {
+                        refRegion.end = end;
+                    }
+                    if (records[i].alignmentStart < refRegion.start) {
+                        refRegion.start = records[i].alignmentStart;
+                    }
+                }
+                // fetch the `seq` for all of the ref regions
+                await Promise.all(Object.values(refRegions).map(async (refRegion) => {
+                    if (refRegion.id !== -1 && refRegion.start <= refRegion.end) {
+                        refRegion.seq = await this.file.fetchReferenceSequenceCallback(refRegion.id, refRegion.start, refRegion.end);
+                    }
+                }));
+                // now decorate all the records with them
+                for (let i = 0; i < records.length; i += 1) {
+                    const seqId = singleRefId !== undefined ? singleRefId : records[i].sequenceId;
+                    const refRegion = refRegions[seqId];
+                    if (refRegion && refRegion.seq) {
+                        records[i].addReferenceSequence(refRegion, compressionScheme);
+                    }
+                }
+            }
+        }
+        return records;
+    }
+}
+// memoize several methods in the class for performance
+'getHeader getBlocks _getBlocksContentIdIndex'
+    .split(' ')
+    .forEach(method => tinyMemoize(CramSlice, method));
+//# sourceMappingURL=index.js.map

package/esm/cramFile/slice/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../src/cramFile/slice/index.js"],"names":[],"mappings":"AAAA,OAAO,EACL,kBAAkB,EAClB,sBAAsB,EACtB,iBAAiB,GAClB,MAAM,cAAc,CAAA;AACrB,OAAO,EAAE,SAAS,EAAE,WAAW,EAAE,WAAW,EAAE,MAAM,SAAS,CAAA;AAE7D,OAAO,SAAS,MAAM,cAAc,CAAA;AACpC,OAAO,YAAY,MAAM,gBAAgB,CAAA;AAEzC;;;;;;GAMG;AACH,SAAS,wCAAwC,CAC/C,UAAU,EACV,mBAAmB,EACnB,UAAU;IAEV,SAAS,kBAAkB,CAAC,WAAW;QACrC,MAAM,OAAO,GAAG,CAAC,WAAW,CAAC,CAAA;QAC7B,IAAI,WAAW,CAAC,gBAAgB,IAAI,CAAC,EAAE;YACrC,MAAM,UAAU,GAAG,UAAU,CAAC,WAAW,CAAC,gBAAgB,CAAC,CAAA;YAC3D,IAAI,CAAC,UAAU,EAAE;gBACf,MAAM,IAAI,kBAAkB,CAC1B,8DAA8D,CAC/D,CAAA;aACF;YACD,OAAO,CAAC,IAAI,CAAC,GAAG,kBAAkB,CAAC,UAAU,CAAC,CAAC,CAAA;SAChD;QACD,OAAO,OAAO,CAAA;IAChB,CAAC;IAED,MAAM,YAAY,GAAG,kBAAkB,CAAC,UAAU,CAAC,CAAA;IACnD,MAAM,MAAM,GAAG,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,cAAc,CAAC,CAAA;IACtD,MAAM,IAAI,GAAG,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,cAAc,GAAG,CAAC,CAAC,UAAU,GAAG,CAAC,CAAC,CAAA;IACvE,MAAM,uBAAuB,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAA;IAC3E,IAAI,uBAAuB,IAAI,CAAC,EAAE;QAChC,YAAY,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE;YACvB,IAAI,CAAC,CAAC,cAAc,KAAK,SAAS,EAAE;gBAClC,MAAM,IAAI,kBAAkB,CAC1B,gGAAgG,CACjG,CAAA;aACF;YACD,CAAC,CAAC,cAAc,GAAG,uBAAuB,CAAA;QAC5C,CAAC,CAAC,CAAA;KACH;AACH,CAAC;AAED;;;;;;GAMG;AACH,SAAS,yCAAyC,CAAC,UAAU,EAAE,UAAU;IACvE,kGAAkG;IAClG,wFAAwF;IACxF,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,cAAc,EAAE,UAAU,CAAC,cAAc,CAAC,CAAA;IAC5E,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAClB,UAAU,CAAC,cAAc,GAAG,UAAU,CAAC,UAAU,GAAG,CAAC,EACrD,UAAU,CAAC,cAAc,GAAG,UAAU,CAAC,UAAU,GAAG,CAAC,CACtD,CAAA;IACD,MAAM,cAAc,GAAG,GAAG,GAAG,KAAK,GAAG,CAAC,CAAA;IACtC,UAAU,CAAC,cAAc,GAAG,cAAc,CAAA;IAC1C,UAAU,CAAC,cAAc,GAAG,cAAc,CAAA;AAC5C,CAAC;AAED;;;;GAIG;AACH,SAAS,uBAAuB,CAC9B,UAAU,EACV,mBAAmB,EACnB,UAAU,EACV,UAAU;IAEV,IAAI,CAAC,UAAU,EAAE;QACf,MAAM,IAAI,kBAAkB,CAC1B,6EAA6E,CAC9E,CAAA;KACF;IAED,MAAM,uBAAuB,GAAG,CAAC,CAAC,CAChC,UAAU,CAAC,IAAI;QACf,CAAC,UAAU,CAAC,gBAAgB,KAAK,SAAS;YACxC,UAAU,CAAC,gBAAgB,KAAK,mBAAmB,CAAC,CACvD,CAAA;IAED,6BAA6B;IAC7B,IAAI,CAAC,UAAU,CAAC,QAAQ,EAAE;QACxB,UAAU,CAAC,QAAQ,GAAG,MAAM,CAAC,UAAU,CAAC,QAAQ,CAAC,CAAA;QACjD,UAAU,CAAC,QAAQ,GAAG,UAAU,CAAC,QAAQ,CAAA;KAC1C;IAED,UAAU,CAAC,IAAI,GAAG;QAChB,UAAU,EAAE,UAAU,CAAC,UAAU;QACjC,cAAc,EAAE,UAAU,CAAC,cAAc;QACzC,QAAQ,EAAE,UAAU,CAAC,QAAQ;KAC9B,CAAA;IACD,IAAI,UAAU,CAAC,QAAQ,EAAE;QACvB,UAAU,CAAC,IAAI,CAAC,QAAQ,GAAG,UAAU,CAAC,QAAQ,CAAA;KAC/C;IAED,2EAA2E;IAC3E,oFAAoF;IACpF,6CAA6C;IAC7C,IAAI,CAAC,UAAU,CAAC,IAAI,IAAI,UAAU,CAAC,gBAAgB,KAAK,SAAS,EAAE;QACjE,UAAU,CAAC,IAAI,GAAG;YAChB,UAAU,EAAE,UAAU,CAAC,UAAU;YACjC,cAAc,EAAE,UAAU,CAAC,cAAc;YACzC,QAAQ,EAAE,UAAU,CAAC,QAAQ;SAC9B,CAAA;QACD,IAAI,UAAU,CAAC,QAAQ,EAAE;YACvB,UAAU,CAAC,IAAI,CAAC,QAAQ,GAAG,UAAU,CAAC,QAAQ,CAAA;SAC/C;KACF;IAED,mEAAmE;IACnE,SAAS;IACT,UAAU,CAAC,KAAK,IAAI,SAAS,CAAC,WAAW,CAAA;IAEzC,8BAA8B;IAC9B,IAAI,UAAU,CAAC,KAAK,GAAG,SAAS,CAAC,UAAU,EAAE;QAC3C,UAAU,CAAC,KAAK,IAAI,SAAS,CAAC,WAAW,CAAA;QACzC,gCAAgC;KACjC;IACD,IAAI,UAAU,CAAC,KAAK,GAAG,SAAS,CAAC,UAAU,EAAE;QAC3C,gCAAgC;QAChC,UAAU,CAAC,KAAK,IAAI,SAAS,CAAC,WAAW,CAAA;KAC1C;IAED,8BAA8B;IAC9B,IAAI,UAAU,CAAC,KAAK,GAAG,SAAS,CAAC,YAAY,EAAE;QAC7C,UAAU,CAAC,KAAK,IAAI,SAAS,CAAC,aAAa,CAAA;KAC5C;IACD,IAAI,UAAU,CAAC,KAAK,GAAG,SAAS,CAAC,YAAY,EAAE;QAC7C,UAAU,CAAC,KAAK,IAAI,SAAS,CAAC,aAAa,CAAA;KAC5C;IAED,IAAI,UAAU,CAAC,cAAc,KAAK,SAAS,EAAE;QAC3C,IAAI,uBAAuB,EAAE;YAC3B,wCAAwC,CACtC,UAAU,EACV,mBAAmB,EACnB,UAAU,CACX,CAAA;SACF;aAAM;YACL,yCAAyC,CAAC,UAAU,EAAE,UAAU,CAAC,CAAA;SAClE;KACF;IAED,4CAA4C;IAC5C,yCAAyC;IACzC,OAAO,UAAU,CAAC,gBAAgB,CAAA;AACpC,CAAC;AAED,MAAM,CAAC,OAAO,OAAO,SAAS;IAC5B,YAAY,SAAS,EAAE,QAAQ;QAC7B,IAAI,CAAC,SAAS,GAAG,SAAS,CAAA;QAC1B,IAAI,CAAC,IAAI,GAAG,SAAS,CAAC,IAAI,CAAA;QAC1B,IAAI,CAAC,iBAAiB,GAAG,QAAQ,CAAA;IACnC,CAAC;IAED,UAAU;IACV,KAAK,CAAC,SAAS;QACb,mCAAmC;QACnC,MAAM,cAAc,GAAG,MAAM,IAAI,CAAC,IAAI,CAAC,iBAAiB,EAAE,CAAA;QAC1D,MAAM,eAAe,GAAG,MAAM,IAAI,CAAC,SAAS,CAAC,SAAS,EAAE,CAAA;QACxD,MAAM,MAAM,GAAG,MAAM,IAAI,CAAC,IAAI,CAAC,SAAS,CACtC,eAAe,CAAC,YAAY,GAAG,IAAI,CAAC,iBAAiB,CACtD,CAAA;QACD,IAAI,MAAM,CAAC,WAAW,KAAK,qBAAqB,EAAE;YAChD,MAAM,CAAC,OAAO,GAAG,SAAS,CACxB,MAAM,CAAC,OAAO,EACd,cAAc,CAAC,qBAAqB,CAAC,MAAM,EAC3C,CAAC,EACD,eAAe,CAAC,YAAY,CAC7B,CAAA;SACF;aAAM,IAAI,MAAM,CAAC,WAAW,KAAK,uBAAuB,EAAE;YACzD,MAAM,CAAC,OAAO,GAAG,SAAS,CACxB,MAAM,CAAC,OAAO,EACd,cAAc,CAAC,uBAAuB,CAAC,MAAM,EAC7C,CAAC,EACD,eAAe,CAAC,YAAY,CAC7B,CAAA;SACF;aAAM;YACL,MAAM,IAAI,kBAAkB,CAC1B,0DAA0D,MAAM,CAAC,YAAY,EAAE,CAChF,CAAA;SACF;QACD,OAAO,MAAM,CAAA;IACf,CAAC;IAED,UAAU;IACV,KAAK,CAAC,SAAS;QACb,MAAM,MAAM,GAAG,MAAM,IAAI,CAAC,SAAS,EAAE,CAAA;QACrC,iDAAiD;QACjD,IAAI,aAAa,GAAG,MAAM,CAAC,YAAY,CAAA;QACvC,MAAM,MAAM,GAAG,IAAI,KAAK,CAAC,MAAM,CAAC,OAAO,CAAC,SAAS,CAAC,CAAA;QAClD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,EAAE;YACzC,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,aAAa,CAAC,CAAA;YACpD,aAAa,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,YAAY,CAAA;SACvC;QAED,OAAO,MAAM,CAAA;IACf,CAAC;IAED,aAAa;IACb,KAAK,CAAC,gBAAgB;QACpB,MAAM,MAAM,GAAG,MAAM,IAAI,CAAC,SAAS,EAAE,CAAA;QACrC,6DAA6D;QAC7D,OAAO,MAAM,CAAC,CAAC,CAAC,CAAA;IAClB,CAAC;IAED,UAAU;IACV,KAAK,CAAC,wBAAwB;QAC5B,MAAM,MAAM,GAAG,MAAM,IAAI,CAAC,SAAS,EAAE,CAAA;QACrC,MAAM,iBAAiB,GAAG,EAAE,CAAA;QAC5B,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE;YACrB,IAAI,KAAK,CAAC,WAAW,KAAK,eAAe,EAAE;gBACzC,iBAAiB,CAAC,KAAK,CAAC,SAAS,CAAC,GAAG,KAAK,CAAA;aAC3C;QACH,CAAC,CAAC,CAAA;QACF,OAAO,iBAAiB,CAAA;IAC1B,CAAC;IAED,KAAK,CAAC,mBAAmB,CAAC,EAAE;QAC1B,MAAM,iBAAiB,GAAG,MAAM,IAAI,CAAC,wBAAwB,EAAE,CAAA;QAC/D,OAAO,iBAAiB,CAAC,EAAE,CAAC,CAAA;IAC9B,CAAC;IAED,KAAK,CAAC,kBAAkB;QACtB,wBAAwB;QACxB,MAAM,WAAW,GAAG,CAAC,MAAM,IAAI,CAAC,SAAS,EAAE,CAAC,CAAC,OAAO,CAAA;QAEpD,IAAI,WAAW,CAAC,QAAQ,GAAG,CAAC,EAAE;YAC5B,OAAO,SAAS,CAAA;SACjB;QAED,MAAM,iBAAiB,GAAG,MAAM,IAAI,CAAC,SAAS,CAAC,oBAAoB,EAAE,CAAA;QAErE,uDAAuD;QAEvD,IAAI,WAAW,CAAC,cAAc,IAAI,CAAC,EAAE;YACnC,MAAM,QAAQ,GAAG,IAAI,CAAC,mBAAmB,CAAC,WAAW,CAAC,cAAc,CAAC,CAAA;YACrE,IAAI,CAAC,QAAQ,EAAE;gBACb,MAAM,IAAI,kBAAkB,CAC1B,kEAAkE,CACnE,CAAA;aACF;YAED,IAAI,WAAW,CAAC,IAAI,GAAG,QAAQ,CAAC,gBAAgB,EAAE;gBAChD,MAAM,IAAI,kBAAkB,CAAC,iCAAiC,CAAC,CAAA;aAChE;YAED,OAAO;gBACL,GAAG,EAAE,QAAQ,CAAC,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC;gBACnC,KAAK,EAAE,WAAW,CAAC,WAAW;gBAC9B,GAAG,EAAE,WAAW,CAAC,WAAW,GAAG,WAAW,CAAC,UAAU,GAAG,CAAC;gBACzD,IAAI,EAAE,WAAW,CAAC,UAAU;aAC7B,CAAA;SACF;QACD,IACE,iBAAiB,CAAC,iBAAiB;YACnC,IAAI,CAAC,IAAI,CAAC,8BAA8B,EACxC;YACA,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,8BAA8B,EAAE;gBAC7C,MAAM,IAAI,KAAK,CACb,sGAAsG,CACvG,CAAA;aACF;YAED,MAAM,GAAG,GAAG,MAAM,IAAI,CAAC,IAAI,CAAC,8BAA8B,CACxD,WAAW,CAAC,QAAQ,EACpB,WAAW,CAAC,WAAW,EACvB,WAAW,CAAC,WAAW,GAAG,WAAW,CAAC,UAAU,GAAG,CAAC,CACrD,CAAA;YAED,IAAI,GAAG,CAAC,MAAM,KAAK,WAAW,CAAC,UAAU,EAAE;gBACzC,MAAM,IAAI,iBAAiB,CACzB,qEAAqE,CACtE,CAAA;aACF;YAED,OAAO;gBACL,GAAG;gBACH,KAAK,EAAE,WAAW,CAAC,WAAW;gBAC9B,GAAG,EAAE,WAAW,CAAC,WAAW,GAAG,WAAW,CAAC,UAAU,GAAG,CAAC;gBACzD,IAAI,EAAE,WAAW,CAAC,UAAU;aAC7B,CAAA;SACF;QAED,OAAO,SAAS,CAAA;IAClB,CAAC;IAED,aAAa;QACX,OAAO,IAAI,CAAC,UAAU,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,CAAA;IACpC,CAAC;IAED,KAAK,CAAC,aAAa;QACjB,MAAM,EAAE,YAAY,EAAE,GAAG,MAAM,IAAI,CAAC,IAAI,CAAC,aAAa,EAAE,CAAA;QAExD,MAAM,iBAAiB,GAAG,MAAM,IAAI,CAAC,SAAS,CAAC,oBAAoB,EAAE,CAAA;QAErE,MAAM,WAAW,GAAG,MAAM,IAAI,CAAC,SAAS,EAAE,CAAA;QAE1C,MAAM,iBAAiB,GAAG,MAAM,IAAI,CAAC,wBAAwB,EAAE,CAAA;QAE/D,sCAAsC;QACtC,IACE,YAAY,GAAG,CAAC;YAChB,IAAI,CAAC,IAAI,CAAC,OAAO,CAAC,gBAAgB;YAClC,WAAW,CAAC,OAAO,CAAC,QAAQ,IAAI,CAAC;YACjC,WAAW,CAAC,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,KAAK,kBAAkB,EACvD;YACA,MAAM,SAAS,GAAG,MAAM,IAAI,CAAC,kBAAkB,EAAE,CAAA;YACjD,IAAI,SAAS,EAAE;gBACb,MAAM,EAAE,GAAG,EAAE,KAAK,EAAE,GAAG,EAAE,GAAG,SAAS,CAAA;gBACrC,MAAM,MAAM,GAAG,WAAW,CAAC,GAAG,CAAC,CAAA;gBAC/B,MAAM,SAAS,GAAG,WAAW,CAAC,OAAO,CAAC,GAAG;qBACtC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;qBACvD,IAAI,CAAC,EAAE,CAAC,CAAA;gBACX,IAAI,MAAM,KAAK,SAAS,EAAE;oBACxB,MAAM,IAAI,kBAAkB,CAC1B,2CAA2C,WAAW,CAAC,OAAO,CAAC,QAAQ,QAAQ,KAAK,KAAK,GAAG,mBAAmB,SAAS,qBAAqB,MAAM,EAAE,CACtJ,CAAA;iBACF;aACF;SACF;QAED,oEAAoE;QACpE,uEAAuE;QACvE,0EAA0E;QAC1E,MAAM,aAAa,GAAG,MAAM,IAAI,CAAC,gBAAgB,EAAE,CAAA;QACnD,MAAM,OAAO,GAAG;YACd,kBAAkB,EAAE,WAAW,CAAC,OAAO,CAAC,WAAW,IAAI,CAAC;YACxD,SAAS,EAAE,EAAE,WAAW,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE;YAC9C,cAAc,EAAE;gBACd,SAAS,CAAC,SAAS;oBACjB,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,EAAE;wBACpB,IAAI,CAAC,SAAS,CAAC,GAAG,EAAE,WAAW,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,CAAA;qBACtD;oBACD,OAAO,IAAI,CAAC,SAAS,CAAC,CAAA;gBACxB,CAAC;aACF;SACF,CAAA;QAED,MAAM,gBAAgB,GAAG,cAAc,CAAC,EAAE;YACxC,MAAM,KAAK,GAAG,iBAAiB,CAAC,qBAAqB,CAAC,cAAc,CAAC,CAAA;YACrE,IAAI,CAAC,KAAK,EAAE;gBACV,MAAM,IAAI,kBAAkB,CAC1B,wBAAwB,cAAc,cAAc,CACrD,CAAA;aACF;YACD,4DAA4D;YAC5D,OAAO,KAAK,CAAC,MAAM,CAAC,IAAI,EAAE,aAAa,EAAE,iBAAiB,EAAE,OAAO,CAAC,CAAA;QACtE,CAAC,CAAA;QACD,IAAI,OAAO,GAAG,IAAI,KAAK,CAAC,WAAW,CAAC,OAAO,CAAC,UAAU,CAAC,CAAA;QACvD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,EAAE;YAC1C,IAAI;gBACF,OAAO,CAAC,CAAC,CAAC,GAAG,YAAY,CACvB,IAAI,EACJ,gBAAgB,EAChB,iBAAiB,EACjB,WAAW,EACX,aAAa,EACb,iBAAiB,EACjB,OAAO,EACP,YAAY,EACZ,CAAC,CACF,CAAA;gBACD,OAAO,CAAC,CAAC,CAAC,CAAC,QAAQ;oBACjB,WAAW,CAAC,eAAe;wBAC3B,WAAW,CAAC,OAAO,CAAC,aAAa;wBACjC,CAAC;wBACD,CAAC,CAAA;aACJ;YAAC,OAAO,CAAC,EAAE;gBACV,IAAI,CAAC,YAAY,sBAAsB,EAAE;oBACvC,OAAO,CAAC,IAAI,CACV,4DAA4D,CAC7D,CAAA;oBACD,OAAO,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;oBAClC,MAAK;iBACN;qBAAM;oBACL,MAAM,CAAC,CAAA;iBACR;aACF;SACF;QAED,+EAA+E;QAC/E,mEAAmE;QACnE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,EAAE;YAC1C,MAAM,EAAE,gBAAgB,EAAE,GAAG,OAAO,CAAC,CAAC,CAAC,CAAA;YACvC,IAAI,gBAAgB,IAAI,CAAC,EAAE;gBACzB,uBAAuB,CACrB,OAAO,EACP,CAAC,EACD,OAAO,CAAC,CAAC,CAAC,EACV,OAAO,CAAC,gBAAgB,CAAC,CAC1B,CAAA;aACF;SACF;QAED,OAAO,OAAO,CAAA;IAChB,CAAC;IAED,KAAK,CAAC,UAAU,CAAC,cAAc;QAC7B,sEAAsE;QACtE,MAAM,QAAQ,GAAG,IAAI,CAAC,SAAS,CAAC,YAAY,GAAG,IAAI,CAAC,iBAAiB,CAAA;QACrE,IAAI,cAAc,GAAG,IAAI,CAAC,IAAI,CAAC,YAAY,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAA;QACzD,IAAI,CAAC,cAAc,EAAE;YACnB,cAAc,GAAG,IAAI,CAAC,aAAa,EAAE,CAAA;YACrC,IAAI,CAAC,IAAI,CAAC,YAAY,CAAC,GAAG,CAAC,QAAQ,EAAE,cAAc,CAAC,CAAA;SACrD;QAED,MAAM,OAAO,GAAG,CAAC,MAAM,cAAc,CAAC,CAAC,MAAM,CAAC,cAAc,CAAC,CAAA;QAE7D,gFAAgF;QAChF,IAAI,OAAO,CAAC,MAAM,IAAI,IAAI,CAAC,IAAI,CAAC,8BAA8B,EAAE;YAC9D,MAAM,WAAW,GAAG,MAAM,IAAI,CAAC,SAAS,EAAE,CAAA;YAC1C,IACE,WAAW,CAAC,OAAO,CAAC,QAAQ,IAAI,CAAC,IAAI,mBAAmB;gBACxD,WAAW,CAAC,OAAO,CAAC,QAAQ,KAAK,CAAC,CAAC,CAAC,kBAAkB;cACtD;gBACA,MAAM,WAAW,GACf,WAAW,CAAC,OAAO,CAAC,QAAQ,IAAI,CAAC;oBAC/B,CAAC,CAAC,WAAW,CAAC,OAAO,CAAC,QAAQ;oBAC9B,CAAC,CAAC,SAAS,CAAA;gBACf,MAAM,iBAAiB,GAAG,MAAM,IAAI,CAAC,SAAS,CAAC,oBAAoB,EAAE,CAAA;gBACrE,MAAM,UAAU,GAAG,EAAE,CAAA,CAAC,+BAA+B;gBAErD,yFAAyF;gBACzF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,EAAE;oBAC1C,MAAM,KAAK,GACT,WAAW,KAAK,SAAS,CAAC,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,UAAU,CAAA;oBACjE,IAAI,SAAS,GAAG,UAAU,CAAC,KAAK,CAAC,CAAA;oBACjC,IAAI,CAAC,SAAS,EAAE;wBACd,SAAS,GAAG;4BACV,EAAE,EAAE,KAAK;4BACT,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,cAAc;4BAChC,GAAG,EAAE,CAAC,QAAQ;yBACf,CAAA;wBACD,UAAU,CAAC,KAAK,CAAC,GAAG,SAAS,CAAA;qBAC9B;oBAED,MAAM,GAAG,GACP,OAAO,CAAC,CAAC,CAAC,CAAC,cAAc;wBACzB,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,WAAW,IAAI,OAAO,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC;wBACjD,CAAC,CAAA;oBACH,IAAI,GAAG,GAAG,SAAS,CAAC,GAAG,EAAE;wBACvB,SAAS,CAAC,GAAG,GAAG,GAAG,CAAA;qBACpB;oBACD,IAAI,OAAO,CAAC,CAAC,CAAC,CAAC,cAAc,GAAG,SAAS,CAAC,KAAK,EAAE;wBAC/C,SAAS,CAAC,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,cAAc,CAAA;qBAC5C;iBACF;gBAED,6CAA6C;gBAC7C,MAAM,OAAO,CAAC,GAAG,CACf,MAAM,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,GAAG,CAAC,KAAK,EAAC,SAAS,EAAC,EAAE;oBAC9C,IAAI,SAAS,CAAC,EAAE,KAAK,CAAC,CAAC,IAAI,SAAS,CAAC,KAAK,IAAI,SAAS,CAAC,GAAG,EAAE;wBAC3D,SAAS,CAAC,GAAG,GAAG,MAAM,IAAI,CAAC,IAAI,CAAC,8BAA8B,CAC5D,SAAS,CAAC,EAAE,EACZ,SAAS,CAAC,KAAK,EACf,SAAS,CAAC,GAAG,CACd,CAAA;qBACF;gBACH,CAAC,CAAC,CACH,CAAA;gBAED,yCAAyC;gBACzC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,EAAE;oBAC1C,MAAM,KAAK,GACT,WAAW,KAAK,SAAS,CAAC,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,UAAU,CAAA;oBACjE,MAAM,SAAS,GAAG,UAAU,CAAC,KAAK,CAAC,CAAA;oBACnC,IAAI,SAAS,IAAI,SAAS,CAAC,GAAG,EAAE;wBAC9B,OAAO,CAAC,CAAC,CAAC,CAAC,oBAAoB,CAAC,SAAS,EAAE,iBAAiB,CAAC,CAAA;qBAC9D;iBACF;aACF;SACF;QAED,OAAO,OAAO,CAAA;IAChB,CAAC;CACF;AAED,uDAAuD;AACvD,8CAA8C;KAC3C,KAAK,CAAC,GAAG,CAAC;KACV,OAAO,CAAC,MAAM,CAAC,EAAE,CAAC,WAAW,CAAC,SAAS,EAAE,MAAM,CAAC,CAAC,CAAA"}

package/esm/cramFile/util.d.ts

@@ -0,0 +1,5 @@
+export function itf8Size(v: any): 1 | 2 | 3 | 4 | 5;
+export function parseItf8(buffer: any, initialOffset: any): any[];
+export function parseItem(buffer: any, parser: any, startBufferPosition?: number, startFilePosition?: number): any;
+export function tinyMemoize(_class: any, methodName: any): void;
+export function sequenceMD5(seq: any): any;

package/esm/cramFile/util.js

@@ -0,0 +1,161 @@
+import md5 from 'md5';
+import { CramBufferOverrunError } from '../errors';
+export function itf8Size(v) {
+    if (!(v & ~0x7f)) {
+        return 1;
+    }
+    if (!(v & ~0x3fff)) {
+        return 2;
+    }
+    if (!(v & ~0x1fffff)) {
+        return 3;
+    }
+    if (!(v & ~0xfffffff)) {
+        return 4;
+    }
+    return 5;
+}
+export function parseItf8(buffer, initialOffset) {
+    let offset = initialOffset;
+    const countFlags = buffer[offset];
+    let result;
+    if (countFlags < 0x80) {
+        result = countFlags;
+        offset += 1;
+    }
+    else if (countFlags < 0xc0) {
+        result = ((countFlags << 8) | buffer[offset + 1]) & 0x3fff;
+        offset += 2;
+    }
+    else if (countFlags < 0xe0) {
+        result =
+            ((countFlags << 16) | (buffer[offset + 1] << 8) | buffer[offset + 2]) &
+                0x1fffff;
+        offset += 3;
+    }
+    else if (countFlags < 0xf0) {
+        result =
+            ((countFlags << 24) |
+                (buffer[offset + 1] << 16) |
+                (buffer[offset + 2] << 8) |
+                buffer[offset + 3]) &
+                0x0fffffff;
+        offset += 4;
+    }
+    else {
+        result =
+            ((countFlags & 0x0f) << 28) |
+                (buffer[offset + 1] << 20) |
+                (buffer[offset + 2] << 12) |
+                (buffer[offset + 3] << 4) |
+                (buffer[offset + 4] & 0x0f);
+        // x=((0xff & 0x0f)<<28) | (0xff<<20) | (0xff<<12) | (0xff<<4) | (0x0f & 0x0f);
+        // TODO *val_p = uv < 0x80000000UL ? uv : -((int32_t) (0xffffffffUL - uv)) - 1;
+        offset += 5;
+    }
+    if (offset > buffer.length) {
+        throw new CramBufferOverrunError('Attempted to read beyond end of buffer; this file seems truncated.');
+    }
+    return [result, offset - initialOffset];
+}
+// parseLtf8(buffer, initialOffset) {
+//   let offset = initialOffset
+//   const countFlags = buffer[offset]
+//   let result
+//   if (countFlags < 0x80) {
+//     result = countFlags
+//     offset += 1
+//   } else if (countFlags < 0xc0) {
+//     result = ((buffer[offset] << 8) | buffer[offset + 1]) & 0x3fff
+//     offset += 2
+//   } else if (countFlags < 0xe0) {
+//     result =
+//       ((buffer[offset] << 16) |
+//         (buffer[offset + 1] << 8) |
+//         buffer[offset + 2]) &
+//       0x1fffff
+//     offset += 3
+//   } else if (countFlags < 0xf0) {
+//     result =
+//       ((buffer[offset] << 24) |
+//         (buffer[offset + 1] << 16) |
+//         (buffer[offset + 2] << 8) |
+//         buffer[offset + 3]) &
+//       0x0fffffff
+//     offset += 4
+//   } else if (countFlags < 0xf8) {
+//     result =
+//       ((buffer[offset] & 15) * Math.pow(2,32) + (buffer[offset + 1] << 24)) |
+//       ((buffer[offset + 2] << 16) |
+//         (buffer[offset + 3] << 8) |
+//         buffer[offset + 4])
+//     // TODO *val_p = uv < 0x80000000UL ? uv : -((int32_t) (0xffffffffUL - uv)) - 1;
+//     offset += 5
+//   } else if (countFlags < 0xfc) {
+//     result =
+//       ((((buffer[offset] & 7) << 8) | buffer[offset + 1]) * Math.pow(2,32) +
+//         (buffer[offset + 2] << 24)) |
+//       ((buffer[offset + 3] << 16) |
+//         (buffer[offset + 4] << 8) |
+//         buffer[offset + 5])
+//     offset += 6
+//   } else if (countFlags < 0xfe) {
+//     result =
+//       ((((buffer[offset] & 3) << 16) |
+//         (buffer[offset + 1] << 8) |
+//         buffer[offset + 2]) *
+//         Math.pow(2,32) +
+//         (buffer[offset + 3] << 24)) |
+//       ((buffer[offset + 4] << 16) |
+//         (buffer[offset + 5] << 8) |
+//         buffer[offset + 6])
+//     offset += 7
+//   } else if (countFlags < 0xff) {
+//     result = Long.fromBytesBE(buffer.slice(offset + 1, offset + 8))
+//     if (
+//       result.greaterThan(Number.MAX_SAFE_INTEGER) ||
+//       result.lessThan(Number.MIN_SAFE_INTEGER)
+//     )
+//       throw new CramUnimplementedError('integer overflow')
+//     result = result.toNumber()
+//     offset += 8
+//   } else {
+//     result = Long.fromBytesBE(buffer.slice(offset + 1, offset + 9))
+//     if (
+//       result.greaterThan(Number.MAX_SAFE_INTEGER) ||
+//       result.lessThan(Number.MIN_SAFE_INTEGER)
+//     )
+//       throw new CramUnimplementedError('integer overflow')
+//     result = result.toNumber()
+//     offset += 9
+//   }
+//   return [result, offset - initialOffset]
+// },
+export function parseItem(buffer, parser, startBufferPosition = 0, startFilePosition = 0) {
+    const { offset, result } = parser.parse(buffer);
+    result._endPosition = offset + startFilePosition;
+    result._size = offset - startBufferPosition;
+    return result;
+}
+// this would be nice as a decorator, but i'm a little worried about
+// babel support for it going away or changing.
+// memoizes a method in the stupidest possible way, with no regard for the
+// arguments.  actually, this only works on methods that take no arguments
+export function tinyMemoize(_class, methodName) {
+    const method = _class.prototype[methodName];
+    const memoAttrName = `_memo_${methodName}`;
+    _class.prototype[methodName] = function _tinyMemoized() {
+        if (!(memoAttrName in this)) {
+            const res = method.call(this);
+            this[memoAttrName] = res;
+            Promise.resolve(res).catch(() => {
+                delete this[memoAttrName];
+            });
+        }
+        return this[memoAttrName];
+    };
+}
+export function sequenceMD5(seq) {
+    return md5(seq.toUpperCase().replace(/[^\x21-\x7e]/g, ''));
+}
+//# sourceMappingURL=util.js.map

package/esm/cramFile/util.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"util.js","sourceRoot":"","sources":["../../src/cramFile/util.js"],"names":[],"mappings":"AAAA,OAAO,GAAG,MAAM,KAAK,CAAA;AACrB,OAAO,EAAE,sBAAsB,EAAE,MAAM,WAAW,CAAA;AAElD,MAAM,UAAU,QAAQ,CAAC,CAAC;IACxB,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE;QAChB,OAAO,CAAC,CAAA;KACT;IACD,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE;QAClB,OAAO,CAAC,CAAA;KACT;IACD,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE;QACpB,OAAO,CAAC,CAAA;KACT;IACD,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,SAAS,CAAC,EAAE;QACrB,OAAO,CAAC,CAAA;KACT;IACD,OAAO,CAAC,CAAA;AACV,CAAC;AAED,MAAM,UAAU,SAAS,CAAC,MAAM,EAAE,aAAa;IAC7C,IAAI,MAAM,GAAG,aAAa,CAAA;IAC1B,MAAM,UAAU,GAAG,MAAM,CAAC,MAAM,CAAC,CAAA;IACjC,IAAI,MAAM,CAAA;IACV,IAAI,UAAU,GAAG,IAAI,EAAE;QACrB,MAAM,GAAG,UAAU,CAAA;QACnB,MAAM,IAAI,CAAC,CAAA;KACZ;SAAM,IAAI,UAAU,GAAG,IAAI,EAAE;QAC5B,MAAM,GAAG,CAAC,CAAC,UAAU,IAAI,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,GAAG,MAAM,CAAA;QAC1D,MAAM,IAAI,CAAC,CAAA;KACZ;SAAM,IAAI,UAAU,GAAG,IAAI,EAAE;QAC5B,MAAM;YACJ,CAAC,CAAC,UAAU,IAAI,EAAE,CAAC,GAAG,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;gBACrE,QAAQ,CAAA;QACV,MAAM,IAAI,CAAC,CAAA;KACZ;SAAM,IAAI,UAAU,GAAG,IAAI,EAAE;QAC5B,MAAM;YACJ,CAAC,CAAC,UAAU,IAAI,EAAE,CAAC;gBACjB,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC;gBAC1B,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;gBACzB,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;gBACrB,UAAU,CAAA;QACZ,MAAM,IAAI,CAAC,CAAA;KACZ;SAAM;QACL,MAAM;YACJ,CAAC,CAAC,UAAU,GAAG,IAAI,CAAC,IAAI,EAAE,CAAC;gBAC3B,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC;gBAC1B,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC;gBAC1B,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;gBACzB,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAA;QAC7B,+EAA+E;QAC/E,+EAA+E;QAC/E,MAAM,IAAI,CAAC,CAAA;KACZ;IACD,IAAI,MAAM,GAAG,MAAM,CAAC,MAAM,EAAE;QAC1B,MAAM,IAAI,sBAAsB,CAC9B,oEAAoE,CACrE,CAAA;KACF;IACD,OAAO,CAAC,MAAM,EAAE,MAAM,GAAG,aAAa,CAAC,CAAA;AACzC,CAAC;AAED,qCAAqC;AACrC,+BAA+B;AAC/B,sCAAsC;AACtC,eAAe;AACf,6BAA6B;AAC7B,0BAA0B;AAC1B,kBAAkB;AAClB,oCAAoC;AACpC,qEAAqE;AACrE,kBAAkB;AAClB,oCAAoC;AACpC,eAAe;AACf,kCAAkC;AAClC,sCAAsC;AACtC,gCAAgC;AAChC,iBAAiB;AACjB,kBAAkB;AAClB,oCAAoC;AACpC,eAAe;AACf,kCAAkC;AAClC,uCAAuC;AACvC,sCAAsC;AACtC,gCAAgC;AAChC,mBAAmB;AACnB,kBAAkB;AAClB,oCAAoC;AACpC,eAAe;AACf,gFAAgF;AAChF,sCAAsC;AACtC,sCAAsC;AACtC,8BAA8B;AAC9B,sFAAsF;AACtF,kBAAkB;AAClB,oCAAoC;AACpC,eAAe;AACf,+EAA+E;AAC/E,wCAAwC;AACxC,sCAAsC;AACtC,sCAAsC;AACtC,8BAA8B;AAC9B,kBAAkB;AAClB,oCAAoC;AACpC,eAAe;AACf,yCAAyC;AACzC,sCAAsC;AACtC,gCAAgC;AAChC,2BAA2B;AAC3B,wCAAwC;AACxC,sCAAsC;AACtC,sCAAsC;AACtC,8BAA8B;AAC9B,kBAAkB;AAClB,oCAAoC;AACpC,sEAAsE;AACtE,WAAW;AACX,uDAAuD;AACvD,iDAAiD;AACjD,QAAQ;AACR,6DAA6D;AAC7D,iCAAiC;AACjC,kBAAkB;AAClB,aAAa;AACb,sEAAsE;AACtE,WAAW;AACX,uDAAuD;AACvD,iDAAiD;AACjD,QAAQ;AACR,6DAA6D;AAC7D,iCAAiC;AACjC,kBAAkB;AAClB,MAAM;AACN,4CAA4C;AAC5C,KAAK;AAEL,MAAM,UAAU,SAAS,CACvB,MAAM,EACN,MAAM,EACN,mBAAmB,GAAG,CAAC,EACvB,iBAAiB,GAAG,CAAC;IAErB,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,CAAA;IAC/C,MAAM,CAAC,YAAY,GAAG,MAAM,GAAG,iBAAiB,CAAA;IAChD,MAAM,CAAC,KAAK,GAAG,MAAM,GAAG,mBAAmB,CAAA;IAC3C,OAAO,MAAM,CAAA;AACf,CAAC;AAED,oEAAoE;AACpE,+CAA+C;AAC/C,0EAA0E;AAC1E,0EAA0E;AAC1E,MAAM,UAAU,WAAW,CAAC,MAAM,EAAE,UAAU;IAC5C,MAAM,MAAM,GAAG,MAAM,CAAC,SAAS,CAAC,UAAU,CAAC,CAAA;IAC3C,MAAM,YAAY,GAAG,SAAS,UAAU,EAAE,CAAA;IAC1C,MAAM,CAAC,SAAS,CAAC,UAAU,CAAC,GAAG,SAAS,aAAa;QACnD,IAAI,CAAC,CAAC,YAAY,IAAI,IAAI,CAAC,EAAE;YAC3B,MAAM,GAAG,GAAG,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAA;YAC7B,IAAI,CAAC,YAAY,CAAC,GAAG,GAAG,CAAA;YACxB,OAAO,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,GAAG,EAAE;gBAC9B,OAAO,IAAI,CAAC,YAAY,CAAC,CAAA;YAC3B,CAAC,CAAC,CAAA;SACH;QACD,OAAO,IAAI,CAAC,YAAY,CAAC,CAAA;IAC3B,CAAC,CAAA;AACH,CAAC;AAED,MAAM,UAAU,WAAW,CAAC,GAAG;IAC7B,OAAO,GAAG,CAAC,GAAG,CAAC,WAAW,EAAE,CAAC,OAAO,CAAC,eAAe,EAAE,EAAE,CAAC,CAAC,CAAA;AAC5D,CAAC"}

package/esm/errors.d.ts

@@ -0,0 +1,23 @@
+export class CramError extends Error {
+}
+/** Error caused by encountering a part of the CRAM spec that has not yet been implemented */
+export class CramUnimplementedError extends Error {
+}
+/** An error caused by malformed data.  */
+export class CramMalformedError extends CramError {
+}
+/**
+ * An error caused by attempting to read beyond the end of the defined data.
+ */
+export class CramBufferOverrunError extends CramMalformedError {
+}
+/**
+ * An error caused by data being too big, exceeding a size limit.
+ */
+export class CramSizeLimitError extends CramError {
+}
+/**
+ * An invalid argument was supplied to a cram-js method or object.
+ */
+export class CramArgumentError extends CramError {
+}

package/esm/errors.js

@@ -0,0 +1,24 @@
+export class CramError extends Error {
+}
+/** Error caused by encountering a part of the CRAM spec that has not yet been implemented */
+export class CramUnimplementedError extends Error {
+}
+/** An error caused by malformed data.  */
+export class CramMalformedError extends CramError {
+}
+/**
+ * An error caused by attempting to read beyond the end of the defined data.
+ */
+export class CramBufferOverrunError extends CramMalformedError {
+}
+/**
+ * An error caused by data being too big, exceeding a size limit.
+ */
+export class CramSizeLimitError extends CramError {
+}
+/**
+ * An invalid argument was supplied to a cram-js method or object.
+ */
+export class CramArgumentError extends CramError {
+}
+//# sourceMappingURL=errors.js.map

package/esm/errors.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"errors.js","sourceRoot":"","sources":["../src/errors.js"],"names":[],"mappings":"AAAA,MAAM,OAAO,SAAU,SAAQ,KAAK;CAAG;AAEvC,6FAA6F;AAC7F,MAAM,OAAO,sBAAuB,SAAQ,KAAK;CAAG;AAEpD,0CAA0C;AAC1C,MAAM,OAAO,kBAAmB,SAAQ,SAAS;CAAG;AAEpD;;GAEG;AACH,MAAM,OAAO,sBAAuB,SAAQ,kBAAkB;CAAG;AAEjE;;GAEG;AACH,MAAM,OAAO,kBAAmB,SAAQ,SAAS;CAAG;AAEpD;;GAEG;AACH,MAAM,OAAO,iBAAkB,SAAQ,SAAS;CAAG"}

package/esm/index.d.ts

@@ -0,0 +1,4 @@
+import CramFile from "./cramFile";
+import IndexedCramFile from "./indexedCramFile";
+import CraiIndex from "./craiIndex";
+export { CramFile, IndexedCramFile, CraiIndex };

package/esm/index.js

@@ -0,0 +1,5 @@
+import CramFile from './cramFile';
+import IndexedCramFile from './indexedCramFile';
+import CraiIndex from './craiIndex';
+export { CramFile, IndexedCramFile, CraiIndex };
+//# sourceMappingURL=index.js.map

package/esm/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.js"],"names":[],"mappings":"AAAA,OAAO,QAAQ,MAAM,YAAY,CAAA;AACjC,OAAO,eAAe,MAAM,mBAAmB,CAAA;AAC/C,OAAO,SAAS,MAAM,aAAa,CAAA;AAEnC,OAAO,EAAE,QAAQ,EAAE,eAAe,EAAE,SAAS,EAAE,CAAA"}