Question: ClustalO computing suboptimal alignment
0
gravatar for mschmid
11 months ago by
mschmid150
Switzerland
mschmid150 wrote:

Clustal Omega - 1.2.4

I align about 600 (n) highly identical sequences. n-1 sequences are identical. There is one sequence which has some differences. For some regions I get an alignment which is obviously suboptiomal (gap is one base to early).

EDIT: This is just one example. Obviously I could just correct this manually. But I automate a process where such alignments should not happen.

ClustalO call is the following:

clustalo -i inp.fasta -o out.clo -t DNA --outfmt=clu --full --dealign --auto

I get the following suboptimal alignment (just a small part not showing all the n-1 identical sequences):

CTCAAGTAACTTCTTTT
CTCAAGTAACTTCTTTT
CTCAAGTAACTTCTTTT
CTCAAGTAACTTCTTTT
CTCAAGTA-ACTTCTTT
********   *  ***

How can I:

  1. Optimize the aligment with ClustalO?
  2. Do a post-optimization?
  3. Or do I have to use a different aligner? For mostly highly identical sequences.
clustalomega alignment • 346 views
ADD COMMENTlink modified 11 months ago • written 11 months ago by mschmid150
1

Any aligner should have spotted that, but its not uncommon to see things like this.

I would try with a different aligner (Clustal typically tends to be better with protein sequences).

Why are you using --dealign?

ADD REPLYlink modified 11 months ago • written 11 months ago by Joe16k

Yes, MUSCLE seems to do a much better job for my application. Strange though....

ADD REPLYlink written 11 months ago by mschmid150
1
gravatar for mschmid
11 months ago by
mschmid150
Switzerland
mschmid150 wrote:

I tried MUSCLE and it seems to do a much better job.

ADD COMMENTlink modified 11 months ago • written 11 months ago by mschmid150
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: 1739 users visited in the last hour