Lighter

Described in:

Song, L., Florea, L. and Langmead, B., Lighter: Fast and Memory-efficient Sequencing Error Correction without Counting. Genome Biol. 2014 Nov 15;15(11):509.

Lighter includes source code from the bloom C++ Bloom filter library. bloom is distributed under the Mozilla Public License (MPL) v1.1.

What is Lighter?

Lighter is a kmer-based error correction method for whole genome sequencing data. Lighter uses sampling (rather than counting) to obtain a set of kmers that are likely from the genome. Using this information, Lighter can correct the reads containing sequence errors.

Install

Clone the GitHub repo, e.g. with git clone https://github.com/mourisl/Lighter.git
Run make in the repo directory

Lighter is small and portable, with pthreads and zlib being the only library dependency.

Usage

Usage: ./lighter [OPTIONS]
OPTIONS:
Required parameters:
    -r seq_file: seq_file is the path to the sequence file. Can use multiple -r to specifiy multiple sequence files
                     The file can be fasta and fastq, and can be gzip'ed with extension *.gz. 
                     When the input file is *.gz, the corresponding output file will also be gzip'ed.
    -k kmer_length genome_size alpha
Other parameters:
    -od: output_file_directory (default: ./)
    -t: number of threads to use (default: 1)
    -noQual: ignore the quality socre (default: false)
    -stable: sequentialize the sampling stage, output the same result with different runs (default: false)
    -trim: allow trimming (default: false)
    -discard: discard unfixable reads. Will LOSE paired-end mataching when discarding (default: false)
    -maxcor: the maximum number of correction for within a kmer_length window (default: 4)

NOTICE: genome_size does not need to be accurate, but it should be at least as large as the size of the sequenced genome.
        alpha is the sampling rate and decided by the user. A rule of thumb: alpha=(7/C), where C is the average coverage of the data set.

Example

Suppose the data set is from E.Coli whose genome size is about 4.7M. If the data set has 3.3M reads in total with read length 100bp, then the average coverage is about 70x and we can set the alpha to be 7/70=0.1.

Single-end data set:

./lighter -r read.fq -k 17 5000000 0.1 -t 10

Paired-end data set:

./lighter -r left.fq -r right.fq -k 17 5000000 0.1 -t 10

Terms of use

This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received (LICENSE.txt) a copy of the GNU General Public License along with this program; if not, you can obtain one from http://www.gnu.org/licenses/gpl.txt or by writing to the Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA

Support

Create a GitHub issue.