Home
Reading
Searching
Subscribe
Sponsors
Statistics
Posting
Contact
Spam
Lists
Links
About
Hosting
Filtering
Features Download
Marketing
Archives
FAQ
Blog
 
Gmane
From: Hannes Hettling <hannes.hettling <at> gmail.com>
Subject: DNAStatistics distance between non-overlapping sequences
Newsgroups: gmane.comp.lang.perl.bio.general
Date: Thursday 27th March 2014 12:03:25 UTC (over 2 years ago)
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
 
CD: 3ms