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?