X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=misc%2Fplot-bamcheck;h=d5857f15798d25995913f4f517048a88c4d64ef0;hb=859feefb171ac2bd02e20e73a3054bd469c69913;hp=adaebe6e92d218533925754758e7c46ba4e852b3;hpb=8928103fafd05b9b555f9ae1710031666725867a;p=samtools.git

diff --git a/misc/plot-bamcheck b/misc/plot-bamcheck
index adaebe6..d5857f1 100755
--- a/misc/plot-bamcheck
+++ b/misc/plot-bamcheck
@@ -803,14 +803,17 @@ sub plot_indel_cycles
         set style line 1 linetype 1  linecolor rgb "red"
         set style line 2 linetype 2  linecolor rgb "black"
         set style line 3 linetype 3  linecolor rgb "green"
+        set style line 4 linetype 4  linecolor rgb "blue"
         set style increment user
         set ylabel "Indel count" 
         set xlabel "Read Cycle"
         set title "$$args{title}"
-        plot '-' w l ti 'Insertions', '' w l ti 'Deletions'
+        plot '-' w l ti 'Insertions (fwd)', '' w l ti 'Insertions (rev)', '' w l ti 'Deletions (fwd)', '' w l ti 'Deletions (rev)'
     ];
     for my $len (@{$$opts{dat}{IC}}) { print $fh "$$len[0]\t$$len[1]\n"; } print $fh "end\n";
     for my $len (@{$$opts{dat}{IC}}) { print $fh "$$len[0]\t$$len[2]\n"; } print $fh "end\n";
+    for my $len (@{$$opts{dat}{IC}}) { print $fh "$$len[0]\t$$len[3]\n"; } print $fh "end\n";
+    for my $len (@{$$opts{dat}{IC}}) { print $fh "$$len[0]\t$$len[4]\n"; } print $fh "end\n";
     close($fh);
     plot($$args{gp});
 }