[Bioperl-l] is the Bio::Align::AlignI doc right?

Ben Bimber bimber at wisc.edu
Thu Feb 25 12:35:01 UTC 2010


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?



More information about the Bioperl-l mailing list