[Bioperl-l] Bio::Align::DNAStatistics.pm

Dominic Rose dominic at bioinf.uni-leipzig.de
Fri Aug 22 13:35:21 UTC 2008


just a short suggestion to improve the code:

in function

sub _build_nt_matrix()

one finds the following lines:

         my $ti_index = $NucleotideIndexes{$ti};
         my $tj_index = $NucleotideIndexes{$tj};

         if( ! defined $ti_index ) {
             print "ti_index not defined for $ti\n";

However, it should be possible to stop/silence the printing of that 
error message. Many alignments contain N's what causes many many
"ti_index not defined for N" messages. That should be avoidable.


Dominic Rose
Professur für Bioinformatik
Institut für Informatik
Universität Leipzig
Härtelstr. 16-18
D-04107 Leipzig
WWW    http://www.bioinf.uni-leipzig.de
Phone: +49 341 97-16698
Fax:   +49 341 97-16679

More information about the Bioperl-l mailing list