[Bioperl-l] splice alignment

Filipe Garrett fgarret at ub.edu
Fri Aug 28 12:37:54 EDT 2009


Hi all,

I need to analyse the 1st, 2nd and 3rd positions of an alignment separately. I've been through 
BioPerl pages but couldn't find no direct way to do it. The closest I fond was "slice" (AlignI) but 
it just extracts a contiguous subsequence.

Is there any subroutine that does the job?
Or maybe a more generic one, so we can select the columns to be extracted; eg:

@aln_pos = qw/1,4,7,10,13,14,17,20/;

$aln_1 = $aln->get_pos(@aln_pos);


thanks in adv,
FG


-- 
Filipe G. Vieira
Departament de Genetica
Universitat de Barcelona
Av. Diagonal, 645
08028 Barcelona
SPAIN
Phone: +34 934 035 306
Fax: +34 934 034 420
fgarret at ub.edu
http://www.ub.edu/molevol/


More information about the Bioperl-l mailing list