added tests for create_seq_index and get_seq

[biopieces.git] / bp_bin / plot_matches

diff --git a/bp_bin/plot_matches b/bp_bin/plot_matches
index 3fc25ad2f871c11799007229d859c3e9673ecaa1..b5d5ed11c70405f8bbd622655e228e9212775e99 100755 (executable)
--- a/bp_bin/plot_matches
+++ b/bp_bin/plot_matches
@@ -113,7 +113,7 @@ sub dotplot_matches
 
     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";
@@ -136,7 +136,7 @@ sub dotplot_matches
     close $fh_forward;
     close $fh_backward;
 
-    $cmd  = "gnuplot";
+    $cmd  = "gnuplot -persist";
 
     $pid = open2( $fh_out, $fh_in, $cmd );