[Bioperl-l] Problem with Bio::SeqIO opening gzipped files

Zayed Albertyn zayed at sanbi.ac.za
Fri Nov 14 05:34:53 EST 2003


Hi Andreas

Adding the -c switch still doesnt work. I still get the same error
message. Input is the full path to the file e.g.

/cip0/db/GENBANK/RELEASE137/gbest13.seq.gz  == $path/$file

I've written another script that does the normal
open(FILE,"/bin/gunzip -c file1 |")

and it works fine

Z

>
>     my $seq_in = Bio::SeqIO::new(
> 	'-file'   => "/bin/gunzip -c $path/$file|",
> 	'-format' => 'genbank'
>     );
>
>
>
> --
> |()()|      Andreas Kähäri                                |(==)|
> |)()(|      EMBL, European Bioinformatics Institute       |=)(=|
> |()()|      Wellcome Trust Genome Campus, Hinxton         |(==)|
> |)()(|      Cambridge, CB10 1SD                           |=)(=|
> |()()|      United Kingdom                                |(==)|
>

-----------------------------------------------
From: Zayed Albertyn
Electric Genetics PTY Ltd
Tel: +27 21 959 3645; Mobile: +2782 480 6097
www.egenetics.com



More information about the Bioperl-l mailing list