Multiple Sequence Alignment Editor With Option To Generate Consensus Sequence
2
6
Entering edit mode
10.9 years ago
Pals ★ 1.3k

Hello, could anyone suggest me an alignment editor program that generates a reliable consensus sequence out the MSA? I tried a number of such programs (Bioedit, Jalview, Seaview etc.). Problem: there is option to change the parameters (threshold values) however, the length of consensus sequence is different while changing the parameters. It can be best explained by the image below. For me, the program should generate the consensus with equal length by ignoring obviously underrepresented positions.

The figure shows part of an alignment having hundreds of sequences. The number after consensus (consensxx) denotes the threshold value (in Seaview program).

consensus alignment multiple • 7.5k views
0
Entering edit mode

You could write a python script to calculate the consensus you had in mind and write it as a .fa for importing into the MSA editor.

6
Entering edit mode
10.7 years ago

I've used CodonCode Aligner and Geneious previously and both seem to work well at this, however they are commercial apps and I'm not sure if the trial versions will work with any more than a few sequences (although Geneious seems to have quite a lot of free functionality).

You could use BioPerl's Bio::SimpleAlign to generate a consensus from an MSA, or BioPython's Bio.AlignIO module.

Perl:

use strict;
use warnings;
use Bio::SimpleAlign;

my alignment = Bio::SimpleAlign->new('alignment.fas'); myconsensus = \$aln->consensus_iupac();


Python:

from Bio import AlignIO
from Bio.Align import AlignInfo

summary_align = AlignInfo.SummaryInfo(alignment)
consensus = summary_align.dumb_consensus()

1
Entering edit mode

Thanks a lot !!

0
Entering edit mode

No problem!! :)

1
Entering edit mode
10.7 years ago
Hari ▴ 280

Try MEGA..http://www.megasoftware.net/