[Bioperl-l] a problem when using the Bioperl modules

Jason Stajich jason.stajich at ucr.edu
Mon Aug 23 11:58:07 EDT 2010


You haven't defined this variable $db - you need to not skip the part 
that initializes the Bio::DB::Fasta object that you had previous asked 
about.
Please send all your future queries to the mailing list.


Guifeng Wei wrote, On 8/23/10 8:14 AM:
> Dear professor,
> after that, i revised my scripts, which is that i divide the genomic 
> sequences into 7 single file, every file contains the sequence from a 
> chromosome.
> however, when i try to run the scripts, the following error was coming.
> Can't call method "seq" on an undefined value at bed_to_fasta.pl 
> <http://bed_to_fasta.pl> line 29, <IN> line 1.
> while(<IN>){
>         chomp $_;
>         my @bed=split(/\s+/, $_ );
>     #print length($db->seq('chrI'));
>         my $chr_id=$bed[0];
>         my $start=$bed[1];
>         my $end=$bed[2];
>         my $seq_name=$bed[3];
>         my $strand=$bed[5];
> my $segment =  $db ->seq($chr_id,$start=>$end);
>         print ">",$seq_name,"_",$chr_id,":",$start=>$end;
>         print "$segment\n";
> }
> the blue line is .
> why?

-- 
Jason E. Stajich, PhD
Assistant Professor
Department of Plant Pathology & Microbiology
University of California
Riverside, CA 92521
jason.stajich at ucr.edu
office: 951.827.2363

http://lab.stajich.org/
http://twitter.com/stajichlab
http://fungalgenomes.org/blog/

http://plantpathology.ucr.edu/
http://genomics.ucr.edu/
http://cepceb.ucr.edu/



More information about the Bioperl-l mailing list