-
Notifications
You must be signed in to change notification settings - Fork 8
Home
This documentation is still in an early stage. If you have any problems/questions/comments please don’t hesitate to contact [email protected] or [email protected]!
NextGenMap is a flexible and fast read mapping program that is more than twice as fast as BWA while achieving a mapping sensitivity similar to Stampy. NextGenMap uses a memory efficient index structure (hash table) to store the positions of all 13-mers present in the reference genome. This index enables a quick identification of potential mapping regions for every read. Unlike other methods, NextGenMap dynamically determines for each read individually how many of the potential mapping regions have to be evaluated by a pairwise sequence alignment. Moreover, NextGenMap uses fast SIMD instructions (SSE) to accelerate the alignment calculations on the CPU. If available NextGenMap calculates the alignments on the GPU (using OpenCL/CUDA) resulting in a runtime reduction of another 20 – 50 %.
Current release: 0.4.10 (tar, zip)
Older releases: see here
Fritz J. Sedlazeck, Philipp Rescheneder, and Arndt von Haeseler NextGenMap: fast and accurate read mapping in highly polymorphic genomes. Bioinformatics (2013) 29 (21): 2790-2791 first published online August 23, 2013 doi:10.1093/bioinformatics/btt468
see here for recent releases or here for older releases
| CPU: | SSE enabled 64 bit dual-core (quad-core or more recommended) |
| RAM: | 8 GB for human (4 GB for genomes < 500 MBp) |
| GPU (optional): | CUDA (Nvidia) or ATI Stream Technology (ATI) enabled |
| OS: | Linux (OpenSUSE with gcc 4.3.4 or Ubuntu 12.04 with gcc 4.6 recommended) |
| Software: | cmake |
see Installation
see here for the documentation and FAQ for frequently asked questions.
Single-end mapping using four CPU threads for the candidate search:
ngm -q reads.fastq -r reference-sequence.fa -o results.sam -t 4
Single-end mapping using four CPU threads for the candidate search and the GPU for alignment computation:
ngm -q reads.fastq -r reference-sequence.fa -o results.sam -t 4 -g
Paired-end mapping using four CPU threads for the candidate search:
For preparing your paired end data see here Only required for version <=0.4.4
ngm -r reference-sequence.fa -1 forward_reads.fq -2 reverse_reads.fq -o results.sam -t 4 -p
Paired-end mapping using four CPU threads for the candidate search and the GPU for alignment computation:
For preparing your paired end data see here Only required for version <=0.4.4
ngm -r reference-sequence.fa -1 forward_reads.fq -2 reverse_reads.fq -o results.sam -t 4 -p -g
Index the genome only:
ngm -r reference-sequence.fa
- NextGenMap falsely interprets negative values on the command line as parameters. e.g. “—score-mismatch -1” will cause an “Argument missing for parameter” error.
Workarround: create a config file (e.g. config.txt) add “score_mismatch = -1” and call NextGenMap with -c config.txt
- The length of a line in a input FASTA file must not exceed 4096 bp.
Solved:
- Error message “error while loading shared libraries: libOpenCL.so.1: cannot open shared object file: No such file or directory” when running NextGenMap
This can be fixed by changing the permissions of the opencl/lib folder: chmod a+rx ngm-0.4.6/opencl/lib/
-
NextGenMap doesn’t work with AMD-APP-SDK v2.7. Please use v2.8 instead.AMD APP SDK installation is now done automatically when building NextGenMap
-
The error message “FATAL: Module fglrx not found.” can be ignored. It will be removed in one of the next AMD APP SDK releases (http://devgurus.amd.com/message/1286600)NextGenMap now uses APP SDK 2.8.1
This documentation is still in an early stage. If you have any problems/questions/comments please don’t hesitate to contact [email protected] or [email protected]!