[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