diff --git a/src/walt/mapping.cpp b/src/walt/mapping.cpp index 3df4a74..1f168b1 100644 --- a/src/walt/mapping.cpp +++ b/src/walt/mapping.cpp @@ -236,7 +236,6 @@ void SingleEndMapping(const string& org_read, const Genome& genome, continue; IndexRegion(read_seed, genome, hash_table, seed_len, region); - if (region.second - region.first + 1 > b) { continue; } diff --git a/src/walt/paired.cpp b/src/walt/paired.cpp index 5da0f08..b32f5e9 100644 --- a/src/walt/paired.cpp +++ b/src/walt/paired.cpp @@ -121,11 +121,10 @@ void PairEndMapping(const string& org_read, const Genome& genome, continue; IndexRegion(read_seed, genome, hash_table, seed_len, region); -#ifdef SMALLBUCKETS - if (region.second - region.first + 1 > 5000) { + if (region.second - region.first + 1 > b) { continue; } -#endif + for (uint32_t j = region.first; j <= region.second; ++j) { uint32_t genome_pos = hash_table.index[j]; uint32_t chr_id = getChromID(genome.start_index, genome_pos);