The mapper supports two different modes of alignment: local (default) and global.
The alignment mode is set using the
In local mode,
-a local, alignments are chosen to be locally optimal, i.e. the mapper looks for the highest scoring alignment of part of the read. That is, unaligned ends are not penalized in any way.
In global mode4.2,
-a global, the mapper is forced to look for the highest scoring alignment of the entire read. That is, unaligned ends are no longer free.
Generally, we recommend that you stick to local alignment, as there are several reasons, why you may not want to force the ends of a raw read sequence to align:
- Read quality deteriorates towards the end of a read, so you may end up aligning noise.
- If present, untrimmed adapter sequence cannot be aligned in any meaningful way.
- Structural variation may be present in the sample, so the reference may not have corresponding sequence for all of the read.
- ... mode4.2
- This is sometimes called semi-global alignment, because the alignment is global with respect to the read, but remains local with respect to the reference.