r987: position on for strand wrongly outputted

This commit is contained in:
Heng Li
2020-04-22 10:31:25 -04:00
parent 9c0e2c67f8
commit cd2b19035b
2 changed files with 2 additions and 2 deletions
+1 -1
View File
@@ -323,7 +323,7 @@ void mm_write_paf3(kstring_t *s, const mm_idx_t *mi, const mm_bseq1_t *t, const
if (mi->seq[r->rid].name) mm_sprintf_lite(s, "%s", mi->seq[r->rid].name);
else mm_sprintf_lite(s, "%d", r->rid);
mm_sprintf_lite(s, "\t%d", mi->seq[r->rid].len);
if (opt_flag & MM_F_QSTRAND)
if ((opt_flag & MM_F_QSTRAND) && r->rev)
mm_sprintf_lite(s, "\t%d\t%d", mi->seq[r->rid].len - r->re, mi->seq[r->rid].len - r->rs);
else
mm_sprintf_lite(s, "\t%d\t%d", r->rs, r->re);