[Bioperl-l] Problem with spliced_seq in BioPerl 1.5.2
Gathman, Allen
agathman at semo.edu
Tue May 1 23:10:20 UTC 2007
Hi, all --
I've been using BioPerl 1.4 for a while; recently I installed 1.5.2, and
found that scripts that had been using spliced_seq are now broken. Any
thoughts on what might be going on?
Here's a sample script:
*********************************************
#!/usr/bin/perl -w
use strict;
use Bio::DB::GFF;
my $db = Bio::DB::GFF-> new ( - adaptor => 'dbi::mysql',
-dsn =>
'dbi:mysql:database=cc;host=localhost',
-fasta => '/gbrowse/databases/cc'
);
$db->add_aggregator('transcript{CDS/mRNA}');
my $seg=$db->segment('ccin_Contig120');
my @genes=$seg->features(-types=>('transcript:GLEAN_alt'));
for my $gene (@genes) {
my $gid = $gene->display_id;
print STDERR "Gene is $gid\n";
my $splgene = $gene->spliced_seq();
}
********************************************
The line with "spliced_seq" in it crashes the program. Here's the
STDERR output:
Gene is Jan06m400_GLEAN_11487
-------------------- WARNING ---------------------
MSG: Calling spliced_seq with a Bio::Das::SegmentI which does have
absolute set to 1 -- be warned you may not be getting things on the
correct strand
---------------------------------------------------
-------------------- WARNING ---------------------
MSG: seq doesn't validate, mismatch is
::,(0,88211,0),::,(0,8821170),::,(0,8821260),::,(0,8821308),::,(0,881935
,),::,(0,881,468),::,(0,881,4,4),::,(0,8818,0),::,(0,8819098)
---------------------------------------------------
------------- EXCEPTION: Bio::Root::Exception -------------
MSG: Attempting to set the sequence to
[Bio::PrimarySeq=HASH(0x88211d0)Bio::PrimarySeq=HASH(0x8821170)Bio::Prim
arySeq=HASH(0x8821260)Bio::PrimarySeq=HASH(0x8821308)Bio::PrimarySeq=HAS
H(0x881935c)Bio::PrimarySeq=HASH(0x881f468)Bio::PrimarySeq=HASH(0x881f4a
4)Bio::PrimarySeq=HASH(0x8818ff0)Bio::PrimarySeq=HASH(0x8819098)] which
does not look healthy
STACK: Error::throw
STACK: Bio::Root::Root::throw
/usr/lib/perl5/site_perl/5.8.6/Bio/Root/Root.pm:359
STACK: Bio::PrimarySeq::seq
/usr/lib/perl5/site_perl/5.8.6/Bio/PrimarySeq.pm:258
STACK: Bio::PrimarySeq::new
/usr/lib/perl5/site_perl/5.8.6/Bio/PrimarySeq.pm:210
STACK: Bio::Seq::new /usr/lib/perl5/site_perl/5.8.6/Bio/Seq.pm:484
STACK: Bio::SeqFeatureI::spliced_seq
/usr/lib/perl5/site_perl/5.8.6/Bio/SeqFeatureI.pm:498
STACK: /transfer/testsplice.pl:20
-----------------------------------------------------------
Allen Gathman
http://cstl-csm.semo.edu/gathman
More information about the Bioperl-l
mailing list