From e4a0313a7ad5aeced0fa45ded49bf2a8fb0a0dbe Mon Sep 17 00:00:00 2001
From: Coppini <guilherme.coppini@gmail.com>
Date: Tue, 7 Dec 2021 16:50:03 -0300
Subject: [PATCH] feat: prints information on how abundance_cutoff was
 calculated and what the minimum number of reads is

---
 NGSpeciesID | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/NGSpeciesID b/NGSpeciesID
index 2dcef92..d1b9c75 100755
--- a/NGSpeciesID
+++ b/NGSpeciesID
@@ -124,7 +124,10 @@ def main(args):
         print()
         work_dir = tempfile.mkdtemp()
         print("Temporary workdirektory for consensus and polishing:", work_dir)
-
+        print(
+            f"Forming draft consensus with abundance_cutoff >= {abundance_cutoff} "
+            f"({args.abundance_ratio * 100}% of {len(read_array)} reads)"
+        )
         centers = consensus.form_draft_consensus(clusters, representatives, sorted_reads_fastq_file, work_dir, abundance_cutoff, args)
 
         if args.primer_file or args.remove_universal_tails: