clustal-js 2.0.0 → 2.0.1

package/CHANGELOG.md

@@ -1,3 +1,7 @@
+# v2.0.1
+
+- Fix preversion build script
+
 # v2.0.0
 
 - Parse pairwise format, outputted by e.g. EMBOSS needle, etc

package/dist/index.d.ts

@@ -0,0 +1,36 @@
+export declare function parseClustalIter(arr: Iterator<string>): {
+    consensus: string;
+    alns: {
+        id: any;
+        seq: any;
+    }[];
+    header: {
+        info: string;
+        version: string;
+    };
+};
+export declare function parsePairwiseIter(arr: string): {
+    consensus: string;
+    alns: {
+        id: any;
+        seq: any;
+    }[];
+};
+export declare function parse(contents: string): {
+    consensus: string;
+    alns: {
+        id: any;
+        seq: any;
+    }[];
+    header: {
+        info: string;
+        version: string;
+    };
+};
+export declare function parsePairwise(contents: string): {
+    consensus: string;
+    alns: {
+        id: any;
+        seq: any;
+    }[];
+};

package/dist/index.js

@@ -0,0 +1,49 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.parsePairwise = exports.parse = exports.parsePairwiseIter = exports.parseClustalIter = void 0;
+const pairwise_1 = require("./pairwise");
+const util_1 = require("./util");
+function parseClustalIter(arr) {
+    const line = (0, util_1.getFirstNonEmptyLine)(arr);
+    if (!line) {
+        throw new Error('Empty file received');
+    }
+    const header = (0, util_1.parseHeader)(line);
+    const res = (0, util_1.parseBlocks)(arr);
+    if (res === undefined) {
+        throw new Error('No blocks parsed');
+    }
+    const alns = res.seqs.map((n, index) => ({ id: res.ids[index], seq: n }));
+    const { consensus } = res;
+    if (consensus.length != alns[0].seq.length) {
+        throw new Error(`Consensus length != sequence length. Con ${consensus.length} seq ${alns[0].seq.length}`);
+    }
+    return { consensus, alns, header };
+}
+exports.parseClustalIter = parseClustalIter;
+function parsePairwiseIter(arr) {
+    const res = (0, pairwise_1.parsePairwiseBlocks)(arr.split('\n')[Symbol.iterator]());
+    if (res === undefined) {
+        throw new Error('No blocks parsed');
+    }
+    const alns = res.seqs.map((n, index) => ({ id: res.ids[index], seq: n }));
+    const { consensus } = res;
+    if (consensus.length != alns[0].seq.length) {
+        throw new Error(`Consensus length != sequence length. Con ${consensus.length} seq ${alns[0].seq.length}`);
+    }
+    return { consensus, alns };
+}
+exports.parsePairwiseIter = parsePairwiseIter;
+function parse(contents) {
+    const iter = contents.split('\n')[Symbol.iterator]();
+    return parseClustalIter(iter);
+}
+exports.parse = parse;
+function parsePairwise(contents) {
+    const res = contents
+        .split('\n')
+        .filter(f => !f.startsWith('#'))
+        .join('\n');
+    return parsePairwiseIter(res);
+}
+exports.parsePairwise = parsePairwise;

package/dist/pairwise.d.ts

@@ -0,0 +1,11 @@
+export declare function getSeqBounds(line: string): readonly [number, number];
+export declare function parsePairwiseBlock(arr: Iterator<string>): {
+    ids: any[];
+    seqs: any[];
+    consensus: string;
+} | undefined;
+export declare function parsePairwiseBlocks(arr: Iterator<string>): {
+    ids: any[];
+    seqs: any[];
+    consensus: string;
+} | undefined;

package/dist/pairwise.js

@@ -0,0 +1,57 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.parsePairwiseBlocks = exports.parsePairwiseBlock = exports.getSeqBounds = void 0;
+const util_1 = require("./util");
+function getSeqBounds(line) {
+    const fields = line.split(/\s+/);
+    const k1 = fields[0].length + fields[1].length;
+    const temp = line.slice(k1);
+    const s = k1 + temp.indexOf(fields[2]);
+    const e = s + fields[2].length;
+    return [s, e];
+}
+exports.getSeqBounds = getSeqBounds;
+// Use the first block to get the sequence identifiers
+function parsePairwiseBlock(arr) {
+    let line = (0, util_1.getFirstNonEmptyLine)(arr);
+    const block = [];
+    let consensusLine = '';
+    if (!line) {
+        return undefined;
+    }
+    while (line) {
+        if (line[0] !== ' ') {
+            block.push(line);
+        }
+        else {
+            consensusLine = line;
+        }
+        line = arr.next().value;
+    }
+    const [start, end] = getSeqBounds(block[0]);
+    const fields = block.map(s => s.split(/\s+/));
+    const ids = fields.map(s => s[0]);
+    const seqs = fields.map(s => s[2]);
+    let consensus = consensusLine.slice(start, end);
+    // handle if the consensus trailing whitespace got trimmed
+    const remainder = seqs[0].length - consensus.length;
+    if (remainder) {
+        consensus += ' '.repeat(remainder);
+    }
+    return { ids, seqs, consensus };
+}
+exports.parsePairwiseBlock = parsePairwiseBlock;
+function parsePairwiseBlocks(arr) {
+    let block;
+    const res = parsePairwiseBlock(arr);
+    if (res !== undefined) {
+        while ((block = parsePairwiseBlock(arr))) {
+            for (let i = 0; i < block.seqs.length; i++) {
+                res.seqs[i] += block.seqs[i];
+            }
+            res.consensus += block.consensus;
+        }
+    }
+    return res;
+}
+exports.parsePairwiseBlocks = parsePairwiseBlocks;

package/dist/util.d.ts

@@ -0,0 +1,17 @@
+export declare function parseVersion(line: string): string;
+export declare function parseHeader(info: string): {
+    info: string;
+    version: string;
+};
+export declare function getFirstNonEmptyLine(arr: Iterator<string>): any;
+export declare function getSeqBounds(line: string): readonly [number, number];
+export declare function parseBlock(arr: Iterator<string>): {
+    ids: any[];
+    seqs: any[];
+    consensus: string;
+} | undefined;
+export declare function parseBlocks(arr: Iterator<string>): {
+    ids: any[];
+    seqs: any[];
+    consensus: string;
+} | undefined;

package/dist/util.js

@@ -0,0 +1,78 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.parseBlocks = exports.parseBlock = exports.getSeqBounds = exports.getFirstNonEmptyLine = exports.parseHeader = exports.parseVersion = void 0;
+function parseVersion(line) {
+    const res = line.match(/\(?(\d+(\.\d+)+)\)?/);
+    return res && res.length > 1 ? res[1] : '';
+}
+exports.parseVersion = parseVersion;
+function parseHeader(info) {
+    const knownHeaders = ['CLUSTAL', 'PROBCONS', 'MUSCLE', 'MSAPROBS', 'Kalign'];
+    if (!knownHeaders.some(l => info.startsWith(l))) {
+        console.warn(`${info} is not a known CLUSTAL header: ${knownHeaders.join(',')}, proceeding but could indicate an issue`);
+    }
+    const version = parseVersion(info);
+    return { info, version };
+}
+exports.parseHeader = parseHeader;
+function getFirstNonEmptyLine(arr) {
+    // There should be two blank lines after the header line
+    let line = arr.next();
+    while (!line.done && line.value.trim() === '') {
+        line = arr.next();
+    }
+    return line.value;
+}
+exports.getFirstNonEmptyLine = getFirstNonEmptyLine;
+function getSeqBounds(line) {
+    const fields = line.split(/\s+/);
+    const temp = line.slice(fields[0].length);
+    const s = fields[0].length + temp.indexOf(fields[1]);
+    const e = s + fields[1].length;
+    return [s, e];
+}
+exports.getSeqBounds = getSeqBounds;
+// Use the first block to get the sequence identifiers
+function parseBlock(arr) {
+    let line = getFirstNonEmptyLine(arr);
+    const block = [];
+    let consensusLine = '';
+    if (!line) {
+        return undefined;
+    }
+    while (line) {
+        if (line[0] !== ' ') {
+            block.push(line);
+        }
+        else {
+            consensusLine = line;
+        }
+        line = arr.next().value;
+    }
+    const [start, end] = getSeqBounds(block[0]);
+    const fields = block.map(s => s.split(/\s+/));
+    const ids = fields.map(s => s[0]);
+    const seqs = block.map(s => s.slice(start, end));
+    let consensus = consensusLine.slice(start, end);
+    // handle if the consensus trailing whitespace got trimmed
+    const remainder = seqs[0].length - consensus.length;
+    if (remainder) {
+        consensus += ' '.repeat(remainder);
+    }
+    return { ids, seqs, consensus };
+}
+exports.parseBlock = parseBlock;
+function parseBlocks(arr) {
+    let block;
+    const res = parseBlock(arr);
+    if (res !== undefined) {
+        while ((block = parseBlock(arr))) {
+            for (let i = 0; i < block.seqs.length; i++) {
+                res.seqs[i] += block.seqs[i];
+            }
+            res.consensus += block.consensus;
+        }
+    }
+    return res;
+}
+exports.parseBlocks = parseBlocks;

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "clustal-js",
-  "version": "2.0.0",
+  "version": "2.0.1",
   "main": "dist/index.js",
   "repository": "cmdcolin/clustal-js",
   "files": [
@@ -22,9 +22,12 @@
     "typescript": "^5.3.3"
   },
   "scripts": {
+    "clean": "rimraf dist",
     "test": "jest",
     "lint": "eslint --ext .ts src test",
     "build": "tsc",
+    "prebuild": "npm run clean",
+    "preversion": "npm run lint && npm test && npm run build",
     "postversion": "git push --follow-tags"
   }
 }