Kraken 2016

From HPC users
Revision as of 14:13, 18 June 2021 by Schwietzer (talk | contribs) (→‎Loading / Using Kraken2)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigationJump to search


Kraken 2 is the newest version of Kraken, a taxonomic classification system using exact k-mer matches to achieve high accuracy and fast classification speeds. This classifier matches each k-mer within a query sequence to the lowest common ancestor (LCA) of all genomes containing the given k-mer. The k-mer assignments inform the classification algorithm. (see: Kraken 1's Webpage for more details). 1

Installed version(s)

The following versions are installed and currently available...

... on environment hpc-env/8.3:

  • name/version-GCC-8.3.0
  • Kraken2/2.1.2-gompi-2019b
  • Kraken2/2.1.1-gompi-2019b

Loading / Using Kraken2

To load the desired version of the module, use the module load command, e.g.

module load hpc-env/8.3
module load Kraken2

Always remember: this command is case sensitive!

To find out on how to use Kraken2you can just type in kraken2 without any arguments to print out a help text to get you started:

$ kraken2
Need to specify input filenames!
Usage: kraken2 [options] <filename(s)>

  --db NAME               Name for Kraken 2 DB
                          (default: none)
  --threads NUM           Number of threads (default: 1)
  --quick                 Quick operation (use first hit or hits)
  --unclassified-out FILENAME
                          Print unclassified sequences to filename
  --classified-out FILENAME
                          Print classified sequences to filename
  --output FILENAME       Print output to filename (default: stdout); "-" will
                          suppress normal output
  --confidence FLOAT      Confidence score threshold (default: 0.0); must be
                          in [0, 1].
  --minimum-base-quality NUM
                          Minimum base quality used in classification (def: 0,
                          only effective with FASTQ input).
  --report FILENAME       Print a report with aggregrate counts/clade to file
  --use-mpa-style         With --report, format report output like Kraken 1's
  --report-zero-counts    With --report, report counts for ALL taxa, even if
                          counts are zero
  --report-minimizer-data With --report, report minimizer and distinct minimizer
                          count information in addition to normal Kraken report
  --memory-mapping        Avoids loading database into RAM
  --paired                The filenames provided have paired-end reads
  --use-names             Print scientific names instead of just taxids
  --gzip-compressed       Input files are compressed with gzip
  --bzip2-compressed      Input files are compressed with bzip2
  --minimum-hit-groups NUM
                          Minimum number of hit groups (overlapping k-mers
                          sharing the same minimizer) needed to make a call
                          (default: 2)
  --help                  Print this message
  --version               Print version information

Kraken2 also comes with a hand full of scripts which can be fount at the bin directory:

$ ls $EBROOTKRAKEN2/bin                  estimate_capacity  lookup_accession_numbers          kraken2    kraken2-build         kraken2-inspect
build_db                   classify                 dump_table         


The full documentation can be found here and at the projects' github page.