mapped = {
let x = bioseq.cigar2gaps (ref_spike,qry_seq,aligned.position,aligned.CIGAR)
.map((sequence, i) => { return { header: i == 0 ? 'reference' : 'query', seq: sequence }});
x.number_of_sequences = 2;
x.number_of_sites = x[0].seq.length;
return x
}