Skip to content

Commit

Permalink
fix bug AFreq calculation on recomb analysis
Browse files Browse the repository at this point in the history
  • Loading branch information
Chienchi Lo committed Jun 9, 2022
1 parent 411f086 commit a20cf6d
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion ramifi/ramifi.py
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down

0 comments on commit a20cf6d

Please sign in to comment.