]> git.donarmstrong.com Git - samtools.git/commitdiff
Prevent division by zero
authorPetr Danecek <pd3@sanger.ac.uk>
Thu, 3 Feb 2011 08:13:48 +0000 (08:13 +0000)
committerPetr Danecek <pd3@sanger.ac.uk>
Thu, 3 Feb 2011 08:13:48 +0000 (08:13 +0000)
bcftools/vcfutils.pl

index 3c962eccb5f52b54628090a4d6f47025b39e1f54..12807c6d40d0af9b05fbdfd8d3065b68623ac0be 100755 (executable)
@@ -201,7 +201,7 @@ Note: This command discards indels. Output: QUAL #non-indel #SNPs #transitions #
          my @x;
          $x[0] = sprintf("%.4f", $c[1]-$c[2]? $c[2] / ($c[1] - $c[2]) : 100);
          $x[1] = sprintf("%.4f", $hsize? $c[3] / $hsize : 0);
-         $x[2] = sprintf("%.4f", $c[3] / $c[1]);
+         $x[2] = sprintf("%.4f", $c[1] ? $c[3] / $c[1] : 0);
          my $a = $c[1] - $lc[1];
          my $b = $c[2] - $lc[2];
          $x[3] = sprintf("%.4f", $a-$b? $b / ($a-$b) : 100);