Add the mate mapping quality tag
This commit is contained in:
parent
d587a9561b
commit
f5e83701ab
1
bwamem.c
1
bwamem.c
|
|
@ -932,6 +932,7 @@ void mem_aln2sam(const mem_opt_t *opt, const bntseq_t *bns, kstring_t *str, bseq
|
|||
kputsn("\tMD:Z:", 6, str); kputs((char*)(p->cigar + p->n_cigar), str);
|
||||
}
|
||||
if (m && m->n_cigar) { kputsn("\tMC:Z:", 6, str); add_cigar(opt, m, str, which); }
|
||||
if (m) { kputsn("\tMQ:i:", 6, str); kputw(m->mapq, str);}
|
||||
if (p->score >= 0) { kputsn("\tAS:i:", 6, str); kputw(p->score, str); }
|
||||
if (p->sub >= 0) { kputsn("\tXS:i:", 6, str); kputw(p->sub, str); }
|
||||
if (bwa_rg_id[0]) { kputsn("\tRG:Z:", 6, str); kputs(bwa_rg_id, str); }
|
||||
|
|
|
|||
Loading…
Reference in New Issue