Subject: | bug in Bio::DB::BigWig |
Date: | Sat, 15 Aug 2015 04:53:25 +0000 |
To: | "bug-Bio-BigFile [...] rt.cpan.org" <bug-Bio-BigFile [...] rt.cpan.org> |
From: | Andre Minoche <a.minoche [...] garvan.org.au> |
Dear Prof Lincoln Stein,
Thank you very much for writing and sharing Bio::DB::BigWig, it is very handy module!
I found a bug which skips one base after every bin. It is problematic if the query interval size is small. E.g.
$bwObj->features(-seq_id=>”chr1”, -start=> 500101,-end=> 500400,-type=>”bin:100”);
Here how to fix it:
in Bio-BigFile-1.07/lib/Bio/DB/BigFile/Iterators.pm
change line 145
$self->{start} += $self->{binsize} + 1;
to: $self->{start} += $self->{binsize};
and line: 157
my $end = int($self->{start}+$self->{binsize});
to: my $end = int($self->{start}+$self->{binsize})-1;
Best regards,
André
Dr Andre Minoche | Research Officer
Cancer Division
Garvan Institute of Medical Research
The Kinghorn Cancer Centre, 370 Victoria Street, Darlinghurst, NSW 2010
NOTICE
Please consider the environment before printing this email. This message and any attachments are intended for the addressee named and may contain legally privileged/confidential/copyright information. If you are not the intended recipient, you should not read, use, disclose, copy or distribute this communication. If you have received this message in error please notify us at once by return email and then delete both messages. We accept no liability for the distribution of viruses or similar in electronic communications. This notice should not be removed.