Difference between revisions of "Kraken"
Moskalenko (talk | contribs) |
|||
| (3 intermediate revisions by 2 users not shown) | |||
| Line 1: | Line 1: | ||
| − | [[Category:Software]][[Category: | + | [[Category:Software]][[Category:Biology]][[Category:NGS]] |
{|<!--CONFIGURATION: REQUIRED--> | {|<!--CONFIGURATION: REQUIRED--> | ||
|{{#vardefine:app|kraken}} | |{{#vardefine:app|kraken}} | ||
| Line 22: | Line 22: | ||
In its fastest mode of operation, for a simulated metagenome of 100 bp reads, Kraken processed over 4 million reads per minute on a single core, over 900 times faster than Megablast and over 11 times faster than the abundance estimation program MetaPhlAn. Kraken's accuracy is comparable with Megablast, with slightly lower sensitivity and very high precision. | In its fastest mode of operation, for a simulated metagenome of 100 bp reads, Kraken processed over 4 million reads per minute on a single core, over 900 times faster than Megablast and over 11 times faster than the abundance estimation program MetaPhlAn. Kraken's accuracy is comparable with Megablast, with slightly lower sensitivity and very high precision. | ||
<!--Modules--> | <!--Modules--> | ||
| − | == | + | ==Environment Modules== |
| − | + | Run <code>module spider {{#var:app}}</code> to find out what environment modules are available for this application. | |
| − | |||
==System Variables== | ==System Variables== | ||
* HPC_{{uc:{{#var:app}}}}_DIR - installation directory | * HPC_{{uc:{{#var:app}}}}_DIR - installation directory | ||
| Line 43: | Line 42: | ||
$ module load kraken/2 | $ module load kraken/2 | ||
$ kraken2-inspect -db fungi | $ kraken2-inspect -db fungi | ||
| + | |||
| + | If you need a custom database [https://github.com/DerrickWood/kraken2/wiki/Manual#custom-databases built] let us know and provide the fasta file with taxonomic IDs. We'd be happy to build and host a commonly used database. | ||
| + | |||
| + | If you get the ''Loading database information...classify: Error reading in hash table'' error allocate more memory to the job. Keep incrementing the memory request until you see a status message in the job log like | ||
| + | Loading database information... done. | ||
| + | Processed 2488462 sequences (226450042 bp) ... | ||
|}} | |}} | ||
| Line 70: | Line 75: | ||
<!--Turn the Table of Contents and Edit paragraph links ON/OFF--> | <!--Turn the Table of Contents and Edit paragraph links ON/OFF--> | ||
__NOTOC____NOEDITSECTION__ | __NOTOC____NOEDITSECTION__ | ||
| − | |||
| − | |||
Latest revision as of 19:28, 12 August 2022
Description
Kraken is a system for assigning taxonomic labels to short DNA sequences, usually obtained through metagenomic studies. Previous attempts by other bioinformatics software to accomplish this task have often used sequence alignment or machine learning techniques that were quite slow, leading to the development of less sensitive but much faster abundance estimation programs. Kraken aims to achieve high sensitivity and high speed by utilizing exact alignments of k-mers and a novel classification algorithm.
In its fastest mode of operation, for a simulated metagenome of 100 bp reads, Kraken processed over 4 million reads per minute on a single core, over 900 times faster than Megablast and over 11 times faster than the abundance estimation program MetaPhlAn. Kraken's accuracy is comparable with Megablast, with slightly lower sensitivity and very high precision.
Environment Modules
Run module spider kraken to find out what environment modules are available for this application.
System Variables
- HPC_KRAKEN_DIR - installation directory
- HPC_KRAKEN_BIN - executable directory
- KRAKEN2_DB_PATH - directory with 'db' builds.
Additional Information
After loading the kraken/2 module you can check the available databases with
$ ls $KRAKEN2_DB_PATH
Use any of the standard or custom databases we host with the '-db DBNAME' argument. E.g.
$ module load kraken/2 $ kraken2-inspect -db fungi
If you need a custom database built let us know and provide the fasta file with taxonomic IDs. We'd be happy to build and host a commonly used database.
If you get the Loading database information...classify: Error reading in hash table error allocate more memory to the job. Keep incrementing the memory request until you see a status message in the job log like
Loading database information... done. Processed 2488462 sequences (226450042 bp) ...