This is a test version of Biostars. For the public version, visit https://www.biostars.org.
LAST for global pairwise alignment of DNA sequences

Hello,

I am trying to use LAST to perform a pairwise global alignment between two DNA sequences. Although LAST very quickly produces a MAF file containing local alignments, I am wondering if there is a way to automatically stitch the local alignments together to produce a global alignment.

Any guidance would be greatly appreciated! Thank you.

(LAST - http://last.cbrc.jp/)

sequence alignment last
• 4,471 views • link
written by dustin.farkas.yo • dustin.farkas.yo

2 answers

You are using the wrong software:

This program finds local alignments between query sequences, and reference sequences that have been prepared using lastdb.

However, you may try -T 1:

-T NUMBER

Type of alignment: 0 means "local alignment" and 1 means "overlap alignment". Local alignments can end anywhere in the middle or at the ends of the sequences. Overlap alignments must extend to the left until they hit the end of a sequence (either query or reference), and to the right until they hit the end of a sequence.

Warning: it's often a bad idea to use -T1. This setting does not change the maximum score drop allowed inside alignments, so if an alignment cannot be extended to the end of a sequence without exceeding this drop, it will be discarded.

• 0 views • link
written by h.mon h.mon

Thank you very much for your reply and suggestion! Unfortunately software such as LAST/MUMmer is the only kind I can find that can handle the 2mbp sequences I am trying to align. I was hoping that they could also perform global alignment as well. I will try your suggestion now.

• 0 views • link
written by dustin.farkas.yo • dustin.farkas.yo

Do not use SUBMIT ANSWER to add comments, use the ADD COMMENT at the appropriated place to, well, add comments.

• 0 views • link
written by h.mon h.mon

Log in to answer this question.

More posts like this