Skip variants with population AN=0 instead of dividing by zero

abhisheknrl · Aug 18, 2016 · 4e3e994 · 4e3e994
1 parent f2c5f1c
commit 4e3e994
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/ExAC.pm b/ExAC.pm
@@ -213,6 +213,7 @@ sub run {
       my $ref_allele = shift @vcf_alleles;
 
       # iterate over required headers
+      HEADER:
       foreach my $h(@{$self->{headers} || []}) {
         my $total_ac = 0;
 
@@ -242,6 +243,7 @@ sub run {
             foreach my $a(@vcf_alleles) {
               my $ac = shift @ac;
               $an = shift @an if @an;
+              next HEADER unless $an;
 
               $total_ac += $ac;
               if ($self->{display_ac}){