Hi,

I am using DSK to count the k-mers that are present in multiple fasta files. I would like to get the list of all solid k-mers that were found in all files.

The command line interface has a "-solid-kmers-out" argument which, I assume, should output this list in a file. However, even when I set a value, it does not seem to output a file. Am I using it correctly?

Also, is there a programmatic way to go from the integer representation of k-mers to their string representation? 

Thanks for your help,
Alex
Bonjour Alexandre,

To output the solid kmers in string representation, you can either use the dsk2ascii tool provided with DSK, or parse the hdf5 file using h5dump. See https://github.com/GATB/dsk#results-visualization

To integrate DSK results in a C++ program, it is currently recommended to use the [GATB-core][1] library. Have a look at [the doxygen page][2] and also the [example code][3] in the GATB-core repository that demonstrate some custom processing:

* [kmer10.cpp][4] is the most up to date version of that snippet.
* [kmer16.cpp][5] shows how to get separate kmer counts (one abundance per file) over multiple files (*multi-bank kmer counting*)
* kmer11-15 are a bit more advanced

Thanks for pointing this out by the way, the "-solid-kmers-out" argument was actually made for GATB-core graph construction (to output the graph and the solid kmers in different files), and is not used in DSK. I've added a note for the future versions of DSK.


  [1]: https://github.com/GATB/gatb-core/
  [2]: http://gatb-core.gforge.inria.fr/doc/api/snippets_kmer.html#snippets_kmer_snippet10
  [3]: https://github.com/GATB/gatb-core/tree/master/gatb-core/examples/kmer
  [4]: https://github.com/GATB/gatb-core/blob/master/gatb-core/examples/kmer/kmer10.cpp
  [5]: https://github.com/GATB/gatb-core/blob/master/gatb-core/examples/kmer/kmer16.cpp