diff --git a/strkit/call/call_locus.py b/strkit/call/call_locus.py index a7efb84..6a87d58 100644 --- a/strkit/call/call_locus.py +++ b/strkit/call/call_locus.py @@ -1109,6 +1109,7 @@ 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)) read_sc += round(read_offset_frac_from_starting_guess * read_sc) (read_cn, read_cn_score), n_read_cn_iters, new_offset_from_starting_count = get_repeat_count( start_count=max(read_sc, 0),