void querySeqc(BloomFilter & myFilter, const string & seq, size_t & fHit) { if (seq.size() < opt::kmerLen) return; for (size_t i = 0; i < seq.size() - opt::kmerLen + 1; i++) { string kmer = seq.substr(i, opt::kmerLen); getCanon(kmer); if(myFilter.containsCit(kmer.c_str())) { #ifdef _OPENMP #pragma omp atomic #endif ++fHit; } } }