Hannes Hettling | 27 Mar 13:03 2014
Picon

DNAStatistics distance between non-overlapping sequences

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

Gmane