This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Aligning with huge insertion variations

What would be the best way to align multiple gene sequences that all have various large insertions like SINEs and LINEs and transposable elements.. I cannot get a good alignment. What would be the best algorithm, parameters, and platform to accomplish this?

alignment

which software have you tried? I would expect bwa-mem to be able to handle those alignments with the 'SA' tags

1 answer

Here it depends on what you are looking for. If you want a single, semi global alignment (with at most non-overlapping alternatives) then bwa-mem will do the job.

If you want to look at all possible alignments then you will need to use a generic aligner such as lastz or last

consolidating the information out of all the alternatives these tools may produce may be more challenging.

Log in to answer this question.