[Bioperl-l] problem with Bio::Align::Utilities - aa_to_dna_aln
Sofia
sofia at bioinfo.hr
Tue Jun 28 15:29:12 UTC 2011
Hi,
I'm using the aa_to_dna_aln function from the BioPerl module
Bio::Align::Utilities.
I've used it some time ago and now i'm just re-running an old script and
i get this error repeatedly:
substr outside of string at
/common/software/API/bioperl-1.6.1/lib/perl5/site_perl/5.8.8//Bio/Align/Utilities.pm
line 160.
Use of uninitialized value in string eq at
/common/software/API/bioperl-1.6.1/lib/perl5/site_perl/5.8.8//Bio/Align/Utilities.pm
line 161.
--------------------- WARNING ---------------------
MSG: In sequence ENSP00000372224 residue count gives end value 1161.
Overriding value [1107] with value 1161 for Bio::LocatableSeq::end().
ATGGGGCGCTGGGCCTGGGTCCCCAGCCCCTGGCCCCCACCGGGGCTGGGCCCCTTCCTCCTCCTCCTCCTGCTGCTGCTGCTGCTGCCACGGGGGTTCCAGCCCCAGCCTGGCGGGAACCGTACGGAGTCCCCAGAACCTAATGCCACAGCGACCCCTGCGATCCCCACTATCCTGGTGACCTCTGTGACCTCTGAGACCCCAGCAACAAGTGCTCCAGAGGCAGAGGGACCCCAAAGTGGGGGGCTCCCGCCCCCGCCCAGGGCAGTTCCCTCGAGCAGTAGCCCCCAGGCCCAAGCACTCACCGAGGAC
---------------------------------------------------
I had a look at the code of Utilities.pm and the line 160 is:
my $char = substr($aa_seqstr,$i + $start_offset,1);
I was wondering why is the $start_offset applied to the amino acid
sequence $aa_seqstr and not to the dna sequence $nt_seqstr(line 164)?
When i remove the $start_offset from the line 160 and add it to the line
164 i don't get any errors.
I don't know if that is a problem or if the problem are the arguments
sending to the function.
The arguments i'm using are the same i used before: an alignment of
protein sequences and a set of dna sequences. Did something change
regarding the arguments?
I would appreciate any help you could provide.
Best regards,
Sofia Pinto
More information about the Bioperl-l
mailing list