From a20cf6d34745e9aeb201006baa57329b03c6464b Mon Sep 17 00:00:00 2001 From: Chienchi Lo Date: Thu, 9 Jun 2022 09:23:09 -0600 Subject: [PATCH] fix bug AFreq calculation on recomb analysis --- ramifi/ramifi.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/ramifi/ramifi.py b/ramifi/ramifi.py index fdbe0f3..1ea9c1c 100755 --- a/ramifi/ramifi.py +++ b/ramifi/ramifi.py @@ -151,7 +151,7 @@ def parse_vcf(argvs,nt_to_variants): depth=m[1] m = re.search(r'DP4=(\d+),(\d+),(\d+),(\d+)', line) (ref_foward,ref_reverse,alt_forward,alt_reverse) = m.groups() - AFreq = float((int(alt_forward) + int(alt_reverse)/int(depth))) + AFreq = float((int(alt_forward) + int(alt_reverse))/int(depth)) ref_bases = content[3] alt_bases = content[4] mutations_count += 1