Best tools to search for and report all alignments

Question

[hersh]

I am looking to align a fastq file of 50mer sequences against a reference genome. I am looking to return ALL possible alignments, not just the best alignments. I looked at using Bowtie2 with the -a flag, however in the documenation it says "Some tools are designed with this reporting mode in mind. Bowtie 2 is not! For very large genomes, this mode is very slow.".

Which tools are designed with this in mind? Ideally the output would return alignments in descending order by alignment score, equivalent to Bowtie2. If there is not a better option for this workflow, I will stick to using Bowtie2. Thanks in advance.

Answer 1

[h.mon]

There are several "all" mappers, segemehl is one I see recommended frequently. Also, see the HTS Mappers (long but even so incomplete) list, search for "A" under "alignment reported" for all-mappers.

Another suggestion at Best aligner for reporting all exact matches of multi-mapping short reads? .

[hersh]

Thanks for the response. I will take a look at segemehl for my analysis, and I appreciate the additional links and info.

Answer 2

[Biogeek]

Subread align or STAR with custom settings to your desire? . These two are pretty fast and accurate.subread slightly pipping STAR according to literature