Commit 54cdad06 authored by LANORE Vincent's avatar LANORE Vincent
Browse files

Fixed bug in nucleotide_fasta_gc

parent db6d9233
......@@ -25,7 +25,7 @@ let nucleotide_fasta_gc ?thirdpos:(tp=false) fa =
)
in
let total_len = Array.fold gc_counts ~init:0. ~f:(fun acc (n,_) -> acc +. n) in
Array.fold gc_counts ~init:0. ~f:(fun acc (len, k) -> acc +. k *. len /. total_len)
Array.fold gc_counts ~init:0. ~f:(fun acc (_, k) -> acc +. k /. total_len)
let main ~alignment () =
match Alignment.from_fasta alignment with
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment