[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