[Bioperl-l] DNAStatistics distance between non-overlapping sequences
Hannes Hettling
hannes.hettling at gmail.com
Thu Mar 27 12:03:25 UTC 2014
A few years ago, a bug was reported (bug 2901) and apparently fixed
(http://bioperl.org/pipermail/bioperl-l/2010-January/032026.html) in
DNAStatistics when calculating distances of non overlapping sequences.
However, I am still getting the error
'Illegal division by zero
at /usr/local/src/bioperl-live/Bio/Align/DNAStatistics.pm line 527'
when calculating distances for the following sequences:
>seq1
GGTACCAGCAGGTGGGCCGCCTACTGCGCACGCGCGGGTTTGCGGGCAGCCGC
>seq2
GGTACCAGCTGGTGGGCCGCCTACTGCGCACGCGCGGGTTTGCGGGCAGCCGC
>seq3
GGTACCAGCAGGTGGTCCGCCTA------------------------------
>seq4
--------------------------CGCACGCGCGTGTTTGCGGGCAGCCGC
I am using the latest bioperl-live from github.
Could anyone suggest a way around this?
Many thanks in advance,
Hannes
More information about the Bioperl-l
mailing list