Question: Reliable alignment tool for similar sequences of similar length other than Mummer?
1
gravatar for nchuang
4.6 years ago by
nchuang230
United States
nchuang230 wrote:

Still new to all this, but I am looking at a collection of similar sequences that are of the same length. I had been using Mummer to align the sequences to a reference. Unfortunately I have been doing it in a cumbersome way by extracting the aligned sequence with Bioperl based on the coord file from Mummer. I did not realize that by doing it this way I am not accounting for large deletions and possibly insertions (meaning my extracted sequences has no gaps when there should be).

My understanding is I should be looking for global alignment but with these indels would I want glocal alignment? What would be best for this if Mummer cannot do this? I feel mummer's algorithm would be perfect for glocal alignment but I don't know how to make it just give me a fasta file as oppose to delta/coord files.

Also I have been testing MAFFT and Tcoffee for this and they work well as far as I can tell. I am just wondering if that is overkill since they are Multiple sequence aligners.

alignment mummer • 1.5k views
ADD COMMENTlink modified 3.4 years ago by novice990 • written 4.6 years ago by nchuang230

Take a look at LASTZ. You may want to look at the link that says (newer versions, unreleased).

ADD REPLYlink written 4.6 years ago by genomax92k

Note that at least one option to get the sequence outputted is to use the delta2maf program. The program is not part of mummer but I uploaded it here https://github.com/cmdcolin/whole_genome_alignment_pipeline/tree/master/util

ADD REPLYlink written 3.4 years ago by cmdcolin1.4k
0
gravatar for novice
3.4 years ago by
novice990
United States
novice990 wrote:

From what I understand, you're looking for an aligner optimized for split-alignment. If that is the case, then standard BWA-SW does that. I would also recommend YAHA, which I found more sensitive than BWA. Both of these support relatively large queries.

ADD COMMENTlink written 3.4 years ago by novice990
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1261 users visited in the last hour