chore(call): rm print debugging + add iters to read log

davidlougheed · Dec 19, 2024 · e02f86d · e02f86d
1 parent 9dbe371
commit e02f86d
Showing 1 changed file with 3 additions and 7 deletions.
diff --git a/strkit/call/call_locus.py b/strkit/call/call_locus.py
@@ -1110,17 +1110,14 @@ def get_read_length_partition_mean(p_idx: int) -> float:
         # Set initial integer copy number guess based on aligned TR size, plus the previous read offset (how much the
         # last guess was wrong by, as a delta.)
         read_sc = round(tr_len / motif_size)
-        print(round(read_offset_frac_from_starting_guess * read_sc))
         if (read_sc_offset := round(read_offset_frac_from_starting_guess * read_sc)) < -1 * read_sc:
             # If our new guess is negative, it was probably a vastly different read, so we should ignore the offset
             # fraction and just use a new starting guess and start again with our offset.
-            read_offset_frac_from_starting_guess = 0
+            read_offset_frac_from_starting_guess = 0.0
         else:
             # Otherwise, use the offset.
             read_sc += read_sc_offset
 
-        print(read_sc)
-
         (read_cn, read_cn_score), n_read_cn_iters, new_offset_from_starting_count = get_repeat_count(
             start_count=read_sc,  # should always be >= 0 with above logic
             tr_seq=tr_read_seq_wc,
@@ -1146,7 +1143,7 @@ def get_read_length_partition_mean(p_idx: int) -> float:
         read_adj_score: float = match_score if tr_len == 0 else read_cn_score / tr_len_w_flank
 
         logger_.debug(
-            "%s - %s | start=%d, of=%.4f | rct=%fs, cn=%d, s=%d[%f]",
+            "%s - %s | start=%d, of=%.4f | rct=%fs, cn=%d, s=%d[%f], i=%d",
             locus_log_str,
             rn,
             read_sc,
@@ -1155,6 +1152,7 @@ def get_read_length_partition_mean(p_idx: int) -> float:
             read_cn,
             read_cn_score,
             read_adj_score,
+            n_read_cn_iters,
         )
 
         if rc_time >= really_bad_read_alignment_time:
@@ -1207,8 +1205,6 @@ def get_read_length_partition_mean(p_idx: int) -> float:
                         tr_read_seq_wc[:20],
                     )
 
-                    print(tr_read_seq_wc)
-
                     return {
                         **locus_result,
                         "peaks": None,