[Bioperl-l] Quickest Codon Based MSA?

Fairley, Derek derek.fairley at belfasttrust.hscni.net
Fri Jan 25 03:31:28 EST 2008


Johan,

There is currently no Bioperl-run wrapper for this program, but you
might want to have a look at Codon Align 2.0 as well:
http://homepage.mac.com/barryghall/CodonAlign.html

Derek

-----Original Message-----
From: bioperl-l-bounces at lists.open-bio.org
[mailto:bioperl-l-bounces at lists.open-bio.org] On Behalf Of Johan Nilsson
Sent: 24 January 2008 22:34
To: bioperl-l at lists.open-bio.org
Subject: [Bioperl-l] Quickest Codon Based MSA?

Hello,

I have a question which might not necessarily be related to Bioperl, 
although I do believe the expertise is available here. I have a couple 
of thousand FASTA files, each containing 20 CDS sequence orthologues of 
rather high sequence similarity. I would like to create a codon-based 
multiple sequence alignment for each of these FASTA files (i.e. a 
nucleotide sequence alignment inferred from alignment of the translated 
peptide sequences, to assure that no frame shifts will occur). I first 
tried running Dialign2, which can perform the 
translation/back-translation in one go, but this turned out to be far 
too slow. I next tried to build protein alignments using ClustalW and 
subsequently built the coding region alignment using EMBOSS 'tranalign',

but this also was too slow.

Is there any method available which significantly speeds up the 
codon-preserving alignment??? As I mentioned, the sequences to be 
aligned are in general very conserved, so any heuristic taking advantage

of the low divergence would be very helpful! Also, is there any 
adjustable parameter in dialign2/dialign-T that might speed up the 
program when looking at highly similar sequences?

Best regards
/Johan Nilsson
_______________________________________________
Bioperl-l mailing list
Bioperl-l at lists.open-bio.org
http://lists.open-bio.org/mailman/listinfo/bioperl-l



More information about the Bioperl-l mailing list