Linux | Mac OSX |
---|---|
DSK is a k-mer counter for reads or genomes.
It takes as input a set of sequences in FASTA or FASTQ format (see "Input" section). DSK outputs a set of solid kmers, i.e. kmers which occur more than a minimal amount of times in the input (-abundance-min parameter). It also outputs the number of times these kmers occur. See "Results visualization" section to learn how to use its output.
CMake 2.6+; see http://www.cmake.org/cmake/resources/software.html
c++ compiler; compilation was tested with gcc and g++ version>=4.5 (Linux) and clang version>=4.1 (Mac OSX).
# get a local copy of source code
git clone --recursive https://github.com/GATB/dsk.git
# compile the code an run a simple test on your computer
cd dsk
sh INSTALL
If you retrieved the software from a git clone, do the following AFTER compiling it:
cd ../scripts # we suppose you are in the build directory
./simple_test.sh
If you retrieved the software from a binary distribution:
cd test
./simple_test.sh
- type
./dsk
for usage instructions.
File input format can be fasta, fastq, either gzipped or not.
To pass several files as input, separate file names by a comma (","), for example:
./dsk -file A1.fa,A2.fa,A3.fa -kmer-size 31
Alternatively, input can be a list of files (one file per line):
ls -1 *.fastq > list_reads
./dsk -file list_reads
DSK uses a binary format for its output.
By default, the format is HDF5, so you can use HDF5 tools for getting ascii output from the HDF5 output (such tools are provided with DSK distribution).
In particular, the HDF5 output of DSK holds two sets of data:
- the couples of (kmer,abundance)
- the histogram of abundances
You can get content of a dataset (kmers in partition 0, and whole histogram) with:
h5dump -y -d dsk/solid/0 output.h5
h5dump -y -d histogram/histogram output.h5
To see the results as a list of "[kmer] [count]\n", use the dsk2ascii
binary with the output file from dsk
dsk2ascii -file output.h5 -out output.txt
To plot kmer coverage distribution,
h5dump -y -d histogram/histogram output.h5 | grep "^\ *[0-9]" | tr -d " " | paste - - | gnuplot -p -e 'plot "-" with lines'
DSK converts all kmers to their canonical representation with respect to reverse-complementation.
In other words, a kmer and its reverse complement are considered to be the same object. For example, with k=3 and assuming the kmer AAA and its reverse complement TTT are both present the input dataset, DSK will consider that one of them is the canonical kmer, e.g. AAA. If AAA is present 2 times and TTT 3 times, then DSK outputs that the count of AAA is 5 and won't return the count of TTT at all.
Note: a canonical kmer is not necessarily the lexicographically smallest one! DSK uses a different ordering for faster performance. Specifically, DSK considers tha A<C<T<G and returns the lexicographically smaller kmer using this alphabet order. So, in the example above, AAA is indeed the canonical kmer. For the GTA/TAC pair, the lexicographically smallest is GTA however the canonical kmer is TAC (as DSK considers that T<G).
DSK supports arbitrary large k-mer lengths. Just compile from the source, to support k-mer lengths up to, say, 160, type this in the build folder:
rm -Rf CMake* && cmake -Dk4=160 .. && make
https://gatb.inria.fr/software/dsk/
To contact a developer: https://gatb.inria.fr/contact/