Multiple Sequence Alignment Editor With Option To Generate Consensus Sequence
2
6
Entering edit mode
9.7 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).

Thanks in advance

consensus alignment multiple • 7.0k views
ADD COMMENT
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.

ADD REPLY
6
Entering edit mode
9.5 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');
my $consensus = $aln->consensus_iupac();

Python:

from Bio import AlignIO
from Bio.Align import AlignInfo

alignment = AlignIO.read(open("alignment.fas"), "fasta")
summary_align = AlignInfo.SummaryInfo(alignment)
consensus = summary_align.dumb_consensus()
ADD COMMENT
1
Entering edit mode

Thanks a lot !!

ADD REPLY
0
Entering edit mode

No problem!! :)

ADD REPLY
1
Entering edit mode
9.5 years ago
Hari ▴ 280

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

ADD COMMENT

Login before adding your answer.

Traffic: 1481 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6