     timer.stop();
     const float unsorted_time = timer.seconds();
 
+    timer.start();
+
     // and now repeat the same operation with the sorted q-grams
-    timer.start();
-
-    // and search the genome q-grams in the index
     thrust::transform(
         d_sorted_qgrams.begin(),
         d_sorted_qgrams.begin() + n_queries,

    log_info(stderr, "  querying q-group index... started\n");
}

template <typename qgram_index_type, typename genome_string>
void test_qgram_index_query(
          qgram_index_type& qgram_index,

    log_info(stderr, "    throughput : %.2f B q-grams/s\n", (1.0e-9f * float( stats.queries )) / stats.filter_time);
}

// main test entry point
//
int qgram_test(int argc, char* argv[])
{
    uint32 n_qgrams      = 10000000;
