Question: How to get all possible alignments from MAFFT or another global aligner
0
gravatar for Elephantdumbo
23 days ago by
Elephantdumbo20 wrote:

Hey guys I have two sequences

> sequence1
 TAGTAGTAG
> Sequence2 
TAG

MAFFT aligns it as

TAGTAGTAG

------TAG

I want it to also output

TAGTAGTAG

TAG------

and

TAGTAGTAG

---TAG---

Does anyone knows how to do that ?

Thank you

sequence alignment • 89 views
ADD COMMENTlink modified 23 days ago by genomax91k • written 23 days ago by Elephantdumbo20
1
gravatar for Mensur Dlakic
23 days ago by
Mensur Dlakic6.9k
USA
Mensur Dlakic6.9k wrote:

The first two alignments are equivalent in terms of score, but the third one isn't because the gap is opened twice.

Alignment programs are geared towards finding an optimal alignment rather than finding all alignments. It really is wasteful to output all the possible sub-optimal alignments. It may not seem so for your toy example, but for 500 sequences of average length 300 and with some repeats in them, this would be a huge output of suboptimal alignments. It is difficult enough (and sometimes slow enough) to output the best alignment only, let alone to do so for thousands of suboptimal results.

ADD COMMENTlink written 23 days ago by Mensur Dlakic6.9k

To actually answer your question, if you are willing to work a bit: there is a toy global alignment program written by Sean Eddy in the supplement of this paper:

https://www.nature.com/articles/nbt0704-909

You should be able to modify it such that it outputs alignments that fit your need.

ADD REPLYlink written 23 days ago by Mensur Dlakic6.9k

Hey Thanks for your response !

I am only interested in pairwise alignment. I will not be aligning 500 sequences. Thanks for the paper. I think it will be helpful

I am interested in sub-optimal alignment because I am trying to model repeat clusters misaligning. If you have any advice regarding that it will be really appreciated.

Thank you so much

ADD REPLYlink modified 23 days ago • written 23 days ago by Elephantdumbo20
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: 1010 users visited in the last hour