-#!/usr/bin/env perl -w
+#!/usr/bin/env perl
# Copyright (C) 2007-2009 Martin A. Hansen.
# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
+use warnings;
use strict;
use Maasha::Biopieces;
use Maasha::Plot;
close $fh;
-Maasha::Filesys::dir_remove( $tmp_dir );
-
Maasha::Biopieces::close_stream( $in );
Maasha::Biopieces::close_stream( $out );
foreach $match ( @{ $matches } )
{
- if ( $match->{ "DIR" } =~ /^f/ )
+ if ( ($match->{ "DIR" } and $match->{ "DIR" } =~ /^f/) or ($match->{ "STRAND" } and $match->{ "STRAND" } eq '+') )
{
print $fh_forward join( "\t", $match->{ "Q_BEG" } + 1, $match->{ "S_BEG" } + 1 ), "\n";
print $fh_forward join( "\t", $match->{ "Q_END" } + 1, $match->{ "S_END" } + 1 ), "\n";
close $fh_forward;
close $fh_backward;
- $cmd = "gnuplot";
+ $cmd = "gnuplot -persist";
$pid = open2( $fh_out, $fh_in, $cmd );