[Bioperl-l] is the Bio::Align::AlignI doc right?
Chris Fields
cjfields at illinois.edu
Thu Feb 25 08:24:17 EST 2010
Yes, but this has been long corrected on main trunk and the latest CPAN code (1.6.1). v.1.2.3 is very old, very very very old. Dinosaurs and perl 5.6.1-old. You should update, even if you are using ensembl's perl API ;>
chris
On Feb 25, 2010, at 6:35 AM, Ben Bimber wrote:
> below is a Bio::Align::AlignI example found here:
>
> http://search.cpan.org/~birney/bioperl-1.2.3/Bio/Align/AlignI.pm#column_from_residue_number
>
>
> Title : column_from_residue_number
> Usage : $col = $ali->column_from_residue_number( $seqname, $resnumber)
> Function:
>
> This function gives the position in the alignment
> (i.e. column number) of the given residue number in the
> sequence with the given name. For example, for the
> alignment
>
> Seq1/91-97 AC..DEF.GH
> Seq2/24-30 ACGG.RTY..
> Seq3/43-51 AC.DDEFGHI
>
> column_from_residue_number( "Seq1", 94 ) returns 5.
> column_from_residue_number( "Seq2", 25 ) returns 2.
> column_from_residue_number( "Seq3", 50 ) returns 9.
>
>
> Perhaps I am simply misunderstanding the method, but shouldnt the first
> example return 6, not 5?
>
> E is the 94th residue. It is in the 6th column. Am I missing something?
> _______________________________________________
> 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