LAST for global pairwise alignment of DNA sequences
2
0
Entering edit mode
4.9 years ago

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 • 2.8k views
ADD COMMENT
2
Entering edit mode
4.9 years ago
h.mon 34k

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.

ADD COMMENT
0
Entering edit mode
4.9 years ago

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.

ADD COMMENT
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 1485 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6