This is a test version of Biostars. For the public version, visit https://www.biostars.org.
report all possible alignments

I am not interested in read mapping per se.

I am interested in reporting ALL possible end-to-end alignments of a read against a reference that have <= k mismatches (e.g. k = 5). I do not care about the "best" alignment.

Standard DNA & RNA read aligners (bwa, bowtie, STAR, etc.) that try to find the "best" alignment for a read place limits on number of mapping locations for short seeds, or something similar.

Does any read aligner allow FULL sensitivity?

For example, I can use MUMmer3 to report ALL exact matches of a read in a genome (of a given substring length). I am looking for something similar - i.e. full sensitivity - but allowing a limited number of mismatches.

Anybody know of a program that does this? I have been looking but am not convinced of the full sensitivity of the aligners.

alignment multimapping

1 answer

There are several "all" mappers, segemehl is one I see recommended frequently. Also, see this (long but even so incomplete) list, search for "A" under "alignment reported" for all-mappers. Finally, if you are wary of results, create a small test data set and benchmark the tools.

Log in to answer this question.