Question: how to output the best alignment rather than multiple alignments for each read in BWA
0
gravatar for zhangdengwei
9 months ago by
zhangdengwei50
zhangdengwei50 wrote:

Dear all, I would like to use BWA mem to output the best alignment for each read instead of all alignments, just like the "-k" argument below in Bowtie2, however, I haven't found the same function in BWA mem, which my version is 0.7.17-r1188, did I miss something or how can I tackle it?

  -k <int>           report up to <int> alns per read; MAPQ not meaningful
   OR
  -a/--all           report all alignments; very slow, MAPQ not meaningful
alignment • 334 views
ADD COMMENTlink modified 9 months ago by d-cameron2.1k • written 9 months ago by zhangdengwei50
0
gravatar for d-cameron
9 months ago by
d-cameron2.1k
Australia
d-cameron2.1k wrote:

did I miss something

Reporting only the best alignment is the default in bwa mem (and bowtie2). By default, bwa reports some alternate alignments in the XA tag but you can ignore that. Note that for some reads the best alignment is chimeric so multiple alignment record will be reported for that read, linked by the SA tag. Note that these record collectively form a single alignment and are not considered to be alternate alignments.

ADD COMMENTlink written 9 months ago by d-cameron2.1k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 709 users visited in the last hour