Question: Are secondary alignments used by any tool?
gravatar for e.benn
6 months ago by
e.benn70 wrote:

I am trying to work out what tools use the secondary alignment feature of the sam format. My understanding is that there will always be exactly 1 primary alignment for each read, and secondary alignments could arise because the aligner cannot decide the exact best alignment.

Does BWA / other aligners emit secondary alignments?

Which downstream tools use secondary alignments? Haplotypecaller?

alignment • 327 views
ADD COMMENTlink modified 6 months ago by h.mon24k • written 6 months ago by e.benn70
gravatar for h.mon
6 months ago by
h.mon24k wrote:

For RNAseq quantification, RSEM, eXpress and Salmon (when quantifying reads from BAMs) include secondary alignments in their quantification. Probably all mappers output secondary alignments, but implementation details may vary, and often there are flags to alter this behaviour. For example, bwa mem -M will mark shorter split hits as secondary instead of suplementary, or STAR --outSAMprimaryFlag AllBestScore will mark all multiple alignments with the same best score as primary.

I am not aware of other tools that use secondary alignments, but this doesn't mean they don't exist.

ADD COMMENTlink modified 6 months ago • written 6 months ago by h.mon24k
Please log in to add an answer.


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