Question: How to get all possible alignments from MAFFT or another global aligner
0
gravatar for Elephantdumbo
5 months 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 • 162 views
ADD COMMENTlink modified 5 months ago by GenoMax96k • written 5 months ago by Elephantdumbo20
1
gravatar for Mensur Dlakic
5 months ago by
Mensur Dlakic9.0k
USA
Mensur Dlakic9.0k 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 5 months ago by Mensur Dlakic9.0k

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 5 months ago by Mensur Dlakic9.0k

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 5 months ago • written 5 months 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: 1658 users visited in the last hour
_