From 7079a9f306b9a4944cb7959c9d213edf2c1e5c1a Mon Sep 17 00:00:00 2001 From: Heng Li Date: Fri, 27 Mar 2020 16:55:28 -0400 Subject: [PATCH] working on mid-sized test example --- Correct.cpp | 3 --- anchor.cpp | 4 ++-- 2 files changed, 2 insertions(+), 5 deletions(-) diff --git a/Correct.cpp b/Correct.cpp index 4ca6da9..1733e94 100644 --- a/Correct.cpp +++ b/Correct.cpp @@ -7213,8 +7213,6 @@ void correct_overlap(overlap_region_alloc* overlap_list, All_reads* R_INF, haplotype_evdience_alloc* hap, Round2_alignment* second_round, int force_repeat, int is_consensus, int* fully_cov, int* abnormal) { - reverse_complement(g_read->seq, g_read->length); - clear_Correct_dumy(dumy, overlap_list); long long window_start, window_end; @@ -7261,7 +7259,6 @@ void correct_overlap(overlap_region_alloc* overlap_list, All_reads* R_INF, (*fully_cov) = check_if_fully_covered(overlap_list, R_INF, g_read, dumy, g, abnormal); - } diff --git a/anchor.cpp b/anchor.cpp index e4abb94..3b1cbca 100644 --- a/anchor.cpp +++ b/anchor.cpp @@ -108,8 +108,8 @@ void ha_get_new_candidates(ha_abuf_t *ab, int64_t rid, UC_Read *ucr, overlap_reg fprintf(stderr, "B\t%ld\t%ld\n", (long)overlap_list->length, (long)overlap_list->mapped_overlaps_length); for (int i = 0; i < (int)overlap_list->length; ++i) { overlap_region *r = &overlap_list->list[i]; - fprintf(stderr, "C\t%d\t%d\t%d\t%c\t%d\t%d\t%d\t%c\t%d\n", (int)r->x_id, (int)r->x_pos_s, (int)r->x_pos_e, "+-"[r->x_pos_strand], - (int)r->y_id, (int)r->y_pos_s, (int)r->y_pos_e, "+-"[r->y_pos_strand], (int)r->shared_seed); + fprintf(stderr, "C\t%d\t%d\t%d\t%c\t%d\t%d\t%d\t%c\t%d\t%d\n", (int)r->x_id, (int)r->x_pos_s, (int)r->x_pos_e, "+-"[r->x_pos_strand], + (int)r->y_id, (int)r->y_pos_s, (int)r->y_pos_e, "+-"[r->y_pos_strand], (int)r->shared_seed, r->is_match); } #endif }