[Bioperl-l] issue writing to pipe with bio sequoia
Stephane Plaisance | VIB |
stephane.plaisance at vib.be
Tue May 31 14:06:48 UTC 2016
Dear,
I try to bgzip my fasta outpmut fbut have an issue
my $seq_out;
if ( defined($zipit) ) {
my $bgzip = `which bgzip`;
chomp($bgzip);
die "No bgzip command available\n" unless ( $bgzip );
$seq_out = Bio::SeqIO -> newFh( -format => 'Fasta', -file => " | $bgzip -c > $outfile\.gz");
} else {
$seq_out = Bio::SeqIO -> newFh( -format => 'Fasta', -file => ">$outfile" );
}
the code fails although bgzip is in my path and the $bgzip variable sets it right
Any help is very welcome
Stephane
More information about the Bioperl-l
mailing list