feat: prints information on how abundance_cutoff was calculated and w…

…hat the minimum number of reads is
ksahlin · Dec 7, 2021 · e4a0313 · e4a0313
1 parent 604445d
commit e4a0313
Showing 1 changed file with 4 additions and 1 deletion.
diff --git 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: