 +		kputw(p->pair_score, str);
 +		kputc(',', str);
 +		if(p->pair_sub >= 0) {
-+			//assert(p->pair_score >= p->pair_sub);
 +			kputw(p->pair_score - p->pair_sub, str);
 +		} else {
 +			kputsn("NA", 2, str);
 +		}
++		kputc(',', str);
++		kputw(p->n_comp, str);
++		kputc(',', str);
++		kputw(p->seedlen0, str);
 +		kputc(',', str);
 +		kputw(p->sub_n, str);
 +		kputc(',', str);

 	kputc('\n', str);
 }
 
@@ -1131,6 +1164,11 @@ mem_aln_t mem_reg2aln(const mem_opt_t *o
 	a.pos = pos - bns->anns[a.rid].offset;
 	a.score = ar->score; a.sub = ar->sub > ar->csub? ar->sub : ar->csub;
 	a.is_alt = ar->is_alt; a.alt_sc = ar->alt_sc;
+	a.seedcov = ((float)ar->seedcov) / l_query;
+	a.frac_rep = ar->frac_rep;
+	a.sub_n = ar->sub_n;

 		for (i = 0; i < n_aa[0]; ++i)
 			mem_aln2sam(opt, bns, &str, &s[0], n_aa[0], aa[0], i, &h[1]); // write read1 hits
 		s[0].sam = strdup(str.s); str.l = 0;
