From 603d1bd1a92b28ea0df3139ba4cf60e1350ea9cc Mon Sep 17 00:00:00 2001 From: Nils Homer Date: Mon, 16 Dec 2024 11:01:07 -0700 Subject: [PATCH] debug --- bwapy/libbwapy.pyx | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/bwapy/libbwapy.pyx b/bwapy/libbwapy.pyx index 621726a..2a1289c 100644 --- a/bwapy/libbwapy.pyx +++ b/bwapy/libbwapy.pyx @@ -301,10 +301,10 @@ cdef class BwaAln: # # tags if seq.type != BWA_TYPE_NO_MATCH: attrs = dict() - attrs["MD"] = bwa_cal_md1(seq.n_cigar, seq.cigar, seq.len, seq.pos, - seq.rseq if seq.strand else seq.seq, self._index._bns.l_pac, - self._pacseq, kstr, &nm) - # rec.set_tags(list(attrs.items())) + # attrs["MD"] = bwa_cal_md1(seq.n_cigar, seq.cigar, seq.len, seq.pos, + # seq.rseq if seq.strand else seq.seq, self._index._bns.l_pac, + # self._pacseq, kstr, &nm) + rec.set_tags(list(attrs.items())) # # TODO:the custom bwa tags: XT, NM, XN, SM, AM, X0, X1, XM, XO, XG, XA, HN return rec