[Bioperl-l] Getting pairwise alignment scores for existing multiple alignment
Fields, Christopher J
cjfields at illinois.edu
Tue Sep 30 16:54:33 UTC 2014
On Sep 30, 2014, at 1:20 AM, Alexey Morozov <alexeymorozov1991 at gmail.com<mailto:alexeymorozov1991 at gmail.com>> wrote:
Is there a method in bioperl that will calculate pairwise alignment scores for any given pair of genes in MSA (according to a given matrix and gap opening/extension cost)? It seems that Bio::SimpleAlign methods only work with score if it has been described in MSA file and can only hold a general multiple sequence alignment score.
LIN SB RAS, bioinformatics group.
Bio::Align::PairwiseStatistics appears to deal with this, see if it fits your needs. You may need to extract the pairwise alignment of the two genes if they are in a multiple sequence alignment.
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Bioperl-l