diff --git a/CommandLines.cpp b/CommandLines.cpp index 3568aa7..85a37e6 100644 --- a/CommandLines.cpp +++ b/CommandLines.cpp @@ -101,7 +101,7 @@ void init_opt(hifiasm_opt_t* asm_opt) asm_opt->thread_num = 1; asm_opt->k_mer_length = 51; asm_opt->mz_win = 51; - asm_opt->sample_dist = 200; + asm_opt->mz_sample_dist = 500; asm_opt->bf_shift = 37; asm_opt->max_kmer_cnt = 2000; asm_opt->high_factor = 5.0; diff --git a/CommandLines.h b/CommandLines.h index 082e141..bd61e93 100644 --- a/CommandLines.h +++ b/CommandLines.h @@ -3,7 +3,7 @@ #include -#define HA_VERSION "0.13-r330-dirty" +#define HA_VERSION "0.13-r332-dirty" #define VERBOSE 0 @@ -34,7 +34,7 @@ typedef struct { int thread_num; int k_mer_length; int mz_win; - int sample_dist; + int mz_sample_dist; int bf_shift; int max_kmer_cnt; double high_factor; // coverage cutoff set to high_factor*hom_cov diff --git a/anchor.cpp b/anchor.cpp index 9f0e83d..fc7fa56 100644 --- a/anchor.cpp +++ b/anchor.cpp @@ -73,7 +73,7 @@ void ha_get_new_candidates(ha_abuf_t *ab, int64_t rid, UC_Read *ucr, overlap_reg rlen = Get_READ_LENGTH(R_INF, rid); // read length // get the list of anchors - ha_sketch(ucr->seq, ucr->length, asm_opt.mz_win, asm_opt.k_mer_length, 0, !(asm_opt.flag & HA_F_NO_HPC), &ab->mz, ha_flt_tab, asm_opt.sample_dist, k_flag, dbg_ct); + ha_sketch(ucr->seq, ucr->length, asm_opt.mz_win, asm_opt.k_mer_length, 0, !(asm_opt.flag & HA_F_NO_HPC), &ab->mz, ha_flt_tab, asm_opt.mz_sample_dist, k_flag, dbg_ct); // minimizer of queried read if (ab->mz.m > ab->old_mz_m) { ab->old_mz_m = ab->mz.m;