Skip Menu |
 

This queue is for tickets about the Bio-SamTools CPAN distribution.

Report information
The Basics
Id: 53705
Status: new
Priority: 0/
Queue: Bio-SamTools

People
Owner: Nobody in particular
Requestors: dan.kortschak [...] adelaide.edu.au
Cc:
AdminCc:

Bug Information
Severity: (no value)
Broken in: (no value)
Fixed in: (no value)



MIME-Version: 1.0
X-Spam-Status: No, hits=-1.0 required=8.0 tests=RCVD_IN_DNSWL_MED,RDNS_NONE
X-Mailer: Evolution 2.26.1
X-Ironport-Anti-Spam-Result: ApoEAEkpT0sKADUU/2dsb2JhbADNPYUUhkuEMAQ
X-Listener: MTA_SMTPS
Message-ID: <1263508313.4379.30.camel [...] epistle>
content-type: text/plain; charset="utf-8"
X-Ironport-Anti-Spam-Filtered: true
X-Ironport-Av: E=Sophos;i="4.49,277,1262525400"; d="scan'208";a="116796124"
Received: from la.mx.develooper.com (x1.develooper.com [207.171.7.70]) by diesel.bestpractical.com (Postfix) with SMTP id 521E54D820B for <bug-Bio-SamTools [...] rt.cpan.org>; Thu, 14 Jan 2010 17:32:13 -0500 (EST)
Received: (qmail 11567 invoked by uid 103); 14 Jan 2010 22:32:12 -0000
Received: from x16.dev (10.0.100.26) by x1.dev with QMQP; 14 Jan 2010 22:32:12 -0000
Received: from Unknown (HELO mailguard-send.adelaide.edu.au) (192.43.227.21) by 16.mx.develooper.com (qpsmtpd/0.80) with ESMTP; Thu, 14 Jan 2010 14:32:10 -0800
Received: from chickenhawk.services.adelaide.edu.au (HELO [129.127.136.4]) ([10.0.53.20]) by arrow-private.services.adelaide.edu.au with ESMTP; 15 Jan 2010 09:01:58 +1030
Delivered-To: cpan-bug+Bio-SamTools [...] diesel.bestpractical.com
Subject: $index undefined prior to coverage call
Return-Path: <dan.kortschak [...] adelaide.edu.au>
X-Original-To: bug-Bio-SamTools [...] rt.cpan.org
X-Spam-Check-BY: 16.mx.develooper.com
X-SBRS: None
Date: Fri, 15 Jan 2010 09:01:53 +1030
X-Spam-Level: *
To: bug-Bio-SamTools [...] rt.cpan.org
Content-Transfer-Encoding: 7bit
From: Dan Kortschak <dan.kortschak [...] adelaide.edu.au>
X-RT-Original-Encoding: ascii
Content-Length: 4010
Download (untitled) / with headers
text/plain 3.9k
Hi Lincoln, Preamble: I'm writing a bioperl wrapper for bowtie that makes use of the Bio::DB::Sam module to parse the bowtie output (via Bio::Assembly::IO::sam). There is a weird dependency issue that I don't understand between the bioperl WrapperBase/CommandExts module the bowtie module uses that has cause B:DB:Sam to fail (error below) - this was not the case in my first version that did not use CommandExts nor in a kludge version that modifies the inheritance of CommandExts to be a derived class from WrapperBase rather than MAJ's designed case where it is an extension to WrapperBase. Details: $ uname -a # debian lenny Linux spawn 2.6.26-2-amd64 #1 SMP Wed Aug 19 22:33:18 UTC 2009 x86_64 GNU/Linux $ perl -v This is perl, v5.10.0 built for x86_64-linux-gnu-thread-multi Details for 'Bio::DB::Sam' Author Lincoln D. Stein (lincoln.stein@gmail.com) Description None given Development Stage Unknown Installed File /usr/local/lib/perl/5.10.0/Bio/DB/Sam.pm Interface Style Unknown Language Used Unknown Package Bio-SamTools-1.09.tar.gz Public License Unknown Support Level Unknown Version Installed 1.09 Version on CPAN 1.09 Non-working bioperl version: r16693 (bioperl-run) and r16692 (bioperl-live) Working bioperl version: r16689 (bioperl-run) and r16690 (bioperl-live) Test script - requires bowtie to be install. All necessary data files are in the bioperl-run t/ filesystem or the bowtie src package: #!/usr/bin/perl use strict; use warnings; use Bio::Tools::Run::Bowtie; # these will depend on your src location for bioperl-live and bioperl-run my $rdq = '/usr/local/src/bioperl-run/t/data/bowtie/reads/e_coli_1000.fq'; my $refseq = '/usr/local/src/bioperl-run/t/data/bowtie/indexes/e_coli'; my $buildfac = Bio::Tools::Run::Bowtie->new( -command => 'build' ); # located in bowtie src directory my $index = $buildfac->run('/usr/local/src/bowtie-0.11.3/genomes/NC_008253.fna','e_coli'); my $bowtiefac = Bio::Tools::Run::Bowtie->new( -command => 'single', -max_seed_mismatches => 2, -seed_length => 28, -max_qual_mismatch => 70, -sam_format => 0 ); my $align = $bowtiefac->run($rdq,$index); print "$align\n"; open (IN,$align); my $lines =(my @lines)= <IN>; print @lines; print "\n\n$lines\n"; close IN; $bowtiefac = Bio::Tools::Run::Bowtie->new( -command => 'single', -max_seed_mismatches => 2, -seed_length => 28, -max_qual_mismatch => 70 ); $align = $bowtiefac->run($rdq,$index); print $align,"\n"; print "Contigs: ", $align->get_nof_contigs, "\tSinglets: ", $align->get_nof_singlets, "\n"; print "In index: $index\n"; print "Contigs:\n"; for ($align->all_contigs) { print $_->get_consensus_sequence->seq,"\n"; } Returned error: [fai_load] build FASTA index. [bam_index_build] creating index for /tmp/zJZ0YTTQ0F/JcG1TNxxJI.srt.bam [bam_index_load] fail to load BAM index. Can't call method "coverage" on an undefined value at /usr/local/lib/perl/5.10.0/Bio/DB/Sam.pm line 1778. I have looked into offending code and the seems to be the result of $index not be returned with a defined value at 1777 in `Bio::DB::Sam.pm'. I'm afraid even after an hour of following it with the debugger I have no idea why or why it depends on the inheritance of an indirectly calling object. Sorry. It may be an issue with bioperl itself (the Bio::Assembly::IO module most likely, though the files appear to be present after B:A:IO has passed control to B:DB:Sam, so I'm unclear on this). thanks for your time Dan Kortchak
MIME-Version: 1.0
X-Spam-Status: No, hits=-1.0 required=8.0 tests=RCVD_IN_DNSWL_MED,RDNS_NONE
In-Reply-To: <rt-3.8.HEAD-10117-1263508363-752.53705-3-0 [...] rt.cpan.org>
X-Ironport-Anti-Spam-Result: ApoEAGOaUEsKADUU/2dsb2JhbADOQYUmhlyEMQQ
X-Mailer: Evolution 2.26.1
X-Listener: MTA_SMTPS
References: <RT-Ticket-53705 [...] rt.cpan.org> <1263508313.4379.30.camel [...] epistle> <rt-3.8.HEAD-10117-1263508363-752.53705-3-0 [...] rt.cpan.org>
X-Ironport-Anti-Spam-Filtered: true
Content-Type: text/plain; charset="utf-8"
Message-ID: <1263603193.15729.5.camel [...] epistle>
X-RT-Original-Encoding: utf-8
X-Ironport-Av: E=Sophos;i="4.49,285,1262525400"; d="scan'208";a="116913166"
Received: from la.mx.develooper.com (x1.develooper.com [207.171.7.70]) by diesel.bestpractical.com (Postfix) with SMTP id 52C824D8002 for <bug-Bio-SamTools [...] rt.cpan.org>; Fri, 15 Jan 2010 19:53:33 -0500 (EST)
Received: (qmail 25846 invoked by uid 103); 16 Jan 2010 00:53:32 -0000
Received: from x16.dev (10.0.100.26) by x1.dev with QMQP; 16 Jan 2010 00:53:32 -0000
Received: from Unknown (HELO mailguard-send.adelaide.edu.au) (192.43.227.21) by 16.mx.develooper.com (qpsmtpd/0.80) with ESMTP; Fri, 15 Jan 2010 16:53:31 -0800
Received: from chickenhawk.services.adelaide.edu.au (HELO [129.127.136.5]) ([10.0.53.20]) by arrow-private.services.adelaide.edu.au with ESMTP; 16 Jan 2010 11:23:16 +1030
Delivered-To: cpan-bug+Bio-SamTools [...] diesel.bestpractical.com
Subject: Re: [rt.cpan.org #53705] AutoReply: $index undefined prior to coverage call
Return-Path: <dan.kortschak [...] adelaide.edu.au>
X-Spam-Check-BY: 16.mx.develooper.com
X-Original-To: bug-Bio-SamTools [...] rt.cpan.org
X-SBRS: None
Date: Sat, 16 Jan 2010 11:23:13 +1030
X-Spam-Level: *
To: bug-Bio-SamTools [...] rt.cpan.org
Content-Transfer-Encoding: 7bit
From: Dan Kortschak <dan.kortschak [...] adelaide.edu.au>
RT-Message-ID: <rt-3.8.HEAD-10124-1263603223-1470.53705-0-0 [...] rt.cpan.org>
Content-Length: 66
Resolved. Results from tempfile cleaning zealousness in bioperl.


This service is sponsored and maintained by Best Practical Solutions and runs on Perl.org infrastructure.

Please report any issues with rt.cpan.org to rt-cpan-admin@bestpractical.com.