r341: bugfix - wrong mate position

when one end is mapped with a score less than -T. Caused by the -T option.
This commit is contained in:
Heng Li 2013-03-07 21:35:57 -05:00
parent b0a76884e8
commit b5b50ac8da
2 changed files with 2 additions and 2 deletions

View File

@ -303,7 +303,7 @@ int mem_sam_pe(const mem_opt_t *opt, const bntseq_t *bns, const uint8_t *pac, co
no_pairing: no_pairing:
for (i = 0; i < 2; ++i) { for (i = 0; i < 2; ++i) {
if (a[i].n) { if (a[i].n && a[i].a[0].score >= opt->T) {
mem_alnreg2hit(&a[i].a[0], &h[i]); mem_alnreg2hit(&a[i].a[0], &h[i]);
bwa_fix_xref(opt->mat, opt->q, opt->r, opt->w, bns, pac, (uint8_t*)s[i].seq, &h[i].qb, &h[i].qe, &h[i].rb, &h[i].re); bwa_fix_xref(opt->mat, opt->q, opt->r, opt->w, bns, pac, (uint8_t*)s[i].seq, &h[i].qb, &h[i].qe, &h[i].rb, &h[i].re);
} else h[i].rb = h[i].re = -1; } else h[i].rb = h[i].re = -1;

2
main.c
View File

@ -3,7 +3,7 @@
#include "utils.h" #include "utils.h"
#ifndef PACKAGE_VERSION #ifndef PACKAGE_VERSION
#define PACKAGE_VERSION "0.7.0-r340-beta" #define PACKAGE_VERSION "0.7.0-r341-beta"
#endif #endif
int bwa_fa2pac(int argc, char *argv[]); int bwa_fa2pac(int argc, char *argv[]);