fixed secondary alignment bug

walaj · May 18, 2024 · 740d895 · 740d895
1 parent 8402085
commit 740d895
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/bwamem_extra.c b/bwamem_extra.c
@@ -109,8 +109,7 @@ mem_alnreg_v mem_align1(const mem_opt_t *opt, const bwt_t *bwt, const bntseq_t *
 	seq = malloc(l_seq);
 	memcpy(seq, seq_, l_seq); // makes a copy of seq_
 	ar = mem_align1_core(opt, bwt, bns, pac, l_seq, seq, 0);
-	//mem_mark_primary_se(opt, ar.n, ar.a, lrand48());
-	//uint32_t k = __ac_Wang_hash(__ac_X31_hash_string(seq) ^ 1337); //Jeremiah
+	mem_mark_primary_se(opt, ar.n, ar.a, lrand48());
 	free(seq);
 	return ar;
 }