[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