<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=utf-8">
</head>
<body bgcolor="#FFFFFF" text="#000066">
Hello Will,<br>
<br>
I'm interested to print AC and AN count for each population in ExAC
VEP plugin output. I've modified the following lines to store the
data:<br>
<br>
<meta http-equiv="content-type" content="text/html; charset=utf-8">
<meta http-equiv="content-type" content="text/html; charset=utf-8">
<pre style="background-color:#2b2b2b;color:#a9b7c6;font-family:'Source Code Pro';font-size:13.5pt;"> <span style="color:#cc7832;">foreach my </span>$a(@vcf_alleles) {
<span style="color:#cc7832;">my </span>$ac = <span style="color:#ab51ba;">shift </span>@ac;
$total_ac += $ac;
$data->{$a}->{<span style="color:#6a8759;">'ExAC_'</span>.$afh.<span style="color:#6a8759;">'_AC'</span>} = $ac;
$data->{$a}->{<span style="color:#6a8759;">'ExAC_'</span>.$afh.<span style="color:#6a8759;">'_AN'</span>} = $an;
$data->{$a}->{<span style="color:#6a8759;">'ExAC_'</span>.$afh} = sprintf(<span style="color:#6a8759;">"%.3g"</span>, $ac / $an);
}
<span style="color:#808080;font-style:italic;"># use total to get ref allele freq
</span><span style="color:#808080;font-style:italic;"> </span>$data->{$ref_allele}->{<span style="color:#6a8759;">'ExAC_'</span>.$afh.<span style="color:#6a8759;">'_AC'</span>} = $total_ac;
$data->{$ref_allele}->{<span style="color:#6a8759;">'ExAC_'</span>.$afh.<span style="color:#6a8759;">'_AN'</span>} = $an;
$data->{$ref_allele}->{<span style="color:#6a8759;">'ExAC_'</span>.$afh} = sprintf(<span style="color:#6a8759;">"%.3g"</span>, <span style="color:#6897bb;">1 </span>- ($total_ac / $an));
}</pre>
However I'm getting a bit messy when it comes to header function. It
loops over ExAC header, but I don't know how to modify
get_header_info() to output the desired columns.<br>
<br>
<br>
Regards,<br>
Guillermo.<br>
</body>
</html>