Question: ClustalO computing suboptimal alignment
0
gravatar for mschmid
22 months ago by
mschmid170
Switzerland
mschmid170 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 • 586 views
ADD COMMENTlink modified 22 months ago • written 22 months ago by mschmid170
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 22 months ago • written 22 months ago by Joe19k

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

ADD REPLYlink written 22 months ago by mschmid170
1
gravatar for mschmid
22 months ago by
mschmid170
Switzerland
mschmid170 wrote:

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

ADD COMMENTlink modified 22 months ago • written 22 months ago by mschmid170
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: 2256 users visited in the last hour
_