"use strict"; var __awaiter = (this && this.__awaiter) || function (thisArg, _arguments, P, generator) { function adopt(value) { return value instanceof P ? value : new P(function (resolve) { resolve(value); }); } return new (P || (P = Promise))(function (resolve, reject) { function fulfilled(value) { try { step(generator.next(value)); } catch (e) { reject(e); } } function rejected(value) { try { step(generator["throw"](value)); } catch (e) { reject(e); } } function step(result) { result.done ? resolve(result.value) : adopt(result.value).then(fulfilled, rejected); } step((generator = generator.apply(thisArg, _arguments || [])).next()); }); }; var __generator = (this && this.__generator) || function (thisArg, body) { var _ = { label: 0, sent: function() { if (t[0] & 1) throw t[1]; return t[1]; }, trys: [], ops: [] }, f, y, t, g; return g = { next: verb(0), "throw": verb(1), "return": verb(2) }, typeof Symbol === "function" && (g[Symbol.iterator] = function() { return this; }), g; function verb(n) { return function (v) { return step([n, v]); }; } function step(op) { if (f) throw new TypeError("Generator is already executing."); while (_) try { if (f = 1, y && (t = op[0] & 2 ? y["return"] : op[0] ? y["throw"] || ((t = y["return"]) && t.call(y), 0) : y.next) && !(t = t.call(y, op[1])).done) return t; if (y = 0, t) op = [op[0] & 2, t.value]; switch (op[0]) { case 0: case 1: t = op; break; case 4: _.label++; return { value: op[1], done: false }; case 5: _.label++; y = op[1]; op = [0]; continue; case 7: op = _.ops.pop(); _.trys.pop(); continue; default: if (!(t = _.trys, t = t.length > 0 && t[t.length - 1]) && (op[0] === 6 || op[0] === 2)) { _ = 0; continue; } if (op[0] === 3 && (!t || (op[1] > t[0] && op[1] < t[3]))) { _.label = op[1]; break; } if (op[0] === 6 && _.label < t[1]) { _.label = t[1]; t = op; break; } if (t && _.label < t[2]) { _.label = t[2]; _.ops.push(op); break; } if (t[2]) _.ops.pop(); _.trys.pop(); continue; } op = body.call(thisArg, _); } catch (e) { op = [6, e]; y = 0; } finally { f = t = 0; } if (op[0] & 5) throw op[1]; return { value: op[0] ? op[1] : void 0, done: true }; } }; var __importDefault = (this && this.__importDefault) || function (mod) { return (mod && mod.__esModule) ? mod : { "default": mod }; }; Object.defineProperty(exports, "__esModule", { value: true }); var child_process_1 = require("child_process"); var fs_1 = __importDefault(require("fs")); var path_1 = __importDefault(require("path")); var async_exec = function (prog, args, onData) { return new Promise(function (resolve, reject) { var child = (0, child_process_1.spawn)(prog, args, { shell: true }); child.stdout.on('data', function (data) { return onData(data.toString().trim()); }); child.stderr.on('data', function (data) { return onData(data.toString().trim()); }); child.on('error', function (err) { return reject(err); }); child.on('exit', function (code) { return resolve(code); }); }); }; var asyncBwaMem = function (refPath, R1, R2, runName, libName, output_dir, onData, options) { if (typeof options === 'undefined') options = {}; var refName = path_1.default.parse(refPath).name; var bwa = 'bwa'; var samblaster = 'samblaster'; var samtools = 'samtools'; var sambamba = 'sambamba'; var R1_arr = Array.isArray(R1) ? R1.join(' ') : R1; var R2_arr = Array.isArray(R2) ? R2.join(' ') : R2; var R1_kitty = R1_arr.slice(-2) === 'gz' ? 'zcat' : 'cat'; var R2_kitty = R2_arr.slice(-2) === 'gz' ? 'zcat' : 'cat'; var R1_in = "'< ".concat(R1_kitty, " ").concat(R1_arr, "'"); var R2_in = "'< ".concat(R2_kitty, " ").concat(R2_arr, "'"); var discordantFile = path_1.default.join(output_dir, "bwa_mem_discordants_on_".concat(refName, ".sam")); var splitterFile = path_1.default.join(output_dir, "bwa_mem_splitters_on_".concat(refName, ".sam")); var unmappedFile = path_1.default.join(output_dir, "bwa_mem_unmapped_on_".concat(refName, ".fq")); var bam = path_1.default.join(output_dir, "bwa_mem_properly_on_".concat(refName, ".bam")); var bamSorted = path_1.default.join(output_dir, "bwa_mem_properly_on_".concat(refName, ".sorted.bam")); return new Promise(function (resolve, reject) { return __awaiter(void 0, void 0, void 0, function () { var code, code_sort, err_1; return __generator(this, function (_a) { switch (_a.label) { case 0: _a.trys.push([0, 3, , 4]); return [4 /*yield*/, async_exec(bwa, ['mem', '-t', String(require('os').cpus().length), '-R', "\"@RG\\tPL:Illumina\\tID:".concat(+(new Date), "\\tSM:").concat(runName, "\\tLB:").concat(libName, "\""), refPath, R1_in, R2_in, '|', samblaster, '--addMateTags', '-a', '-e', '-d', discordantFile, '-s', splitterFile, '-u', unmappedFile, '|', samtools, 'view', '-b', '-', '>', bam], function (message) { return onData('[BWA-MEM] ' + message); })]; case 1: code = _a.sent(); onData('[BWA-MEM][EXIT CODE] ' + code); return [4 /*yield*/, async_exec(sambamba, ['sort', '-t', String(require('os').cpus().length), bam ], function (message) { return onData('[SAMBAMBA-SORT] ' + message); })]; case 2: code_sort = _a.sent(); onData('[SAMBAMBA-SORT][EXIT CODE] ' + code_sort); fs_1.default.unlinkSync(bam); resolve([bamSorted, discordantFile, splitterFile, unmappedFile]); return [3 /*break*/, 4]; case 3: err_1 = _a.sent(); reject(err_1); return [3 /*break*/, 4]; case 4: return [2 /*return*/]; } }); }); }); }; exports.default = asyncBwaMem;