From 9da737246359c7be8431f6125cb495e09a5d4c0b Mon Sep 17 00:00:00 2001
From: William McLaren <wm2@ebi.ac.uk>
Date: Fri, 29 Jul 2016 11:19:16 +0100
Subject: [PATCH] fix for ANs with commas

---
 ExAC.pm | 8 ++++++--
 1 file changed, 6 insertions(+), 2 deletions(-)

diff --git a/ExAC.pm b/ExAC.pm
index f194a9fd..0f681958 100644
--- a/ExAC.pm
+++ b/ExAC.pm
@@ -234,11 +234,14 @@ sub run {
           if(/$anh\=([0-9\,]+)/) {
             
             # grab AN
-            my $an = $1;            
-            next unless $an;
+            my @an = split /\,/, $1;
+            next unless @an;
+            my $an;
 
             foreach my $a(@vcf_alleles) {
               my $ac = shift @ac;
+              $an = shift @an if @an;
+
               $total_ac += $ac;
               if ($self->{display_ac}){
                 $data->{$a}->{'ExAC_'.$ach} = $ac;
@@ -246,6 +249,7 @@ sub run {
               if ($self->{display_an}){
                 $data->{$a}->{'ExAC_'.$anh} = $an;
               }
+
               $data->{$a}->{'ExAC_'.$afh} = sprintf("%.3g", $ac / $an);
             }