/fastx_toolkit

FASTA/FASTQ pre-processing programs

Primary LanguageCOtherNOASSERTION

FASTX-Toolkit
=============


Short Summary
===============

The FASTX-Toolkit is a collection of command line tools for Short-Reads 
FASTA/FASTQ files preprocessing.



More Details
==============

Next-Generation sequencing machines usually produce FASTA or FASTQ files, 
containing multiple short-reads sequences (possibly with quality information).

The main processing of such FASTA/FASTQ files is mapping (aka aligning)
the sequences to reference genomes or other databases using specialized
programs. 

Example of such mapping programs are:
Blat (http://www.kentinformatics.com/index.asp), 
SHRiMP (http://compbio.cs.toronto.edu/shrimp),
LastZ (http://www.bx.psu.edu/miller_lab),
MAQ (http://maq.sourceforge.net/)
And many many others.

However, 
It is sometimes more productive to preprocess the FASTA/FASTQ files before 
mapping the sequences to the genome - manipulating the sequences to 
produce better mapping results.

The FASTX-Toolkit tools perform some of these preprocessing tasks.



Available Tools
===============

FASTQ-to-FASTA - Converts a FASTQ file to a FASTA file..

FASTQ-Statistics - scans a FASTQ file, and produces some statistics about the 
	quality and the sequences in the file.
	
FASTQ-Quality-BoxPlot, and
FASTQ-Nucleotides-Distribution - Generates charts based on the statistics 
	generated by FASTQ-Statistics. These charts can be used to quickly
	see the quality of the sequenced library.
	
FASTQ-Quality-Converter - Converts from ASCII to numeric quality scores.

FASTQ-Quality-Filter - removes low-quality sequences from FASTQ files.

FASTX-Artifacts-Filter - removes some sequencing artifacts from FASTA/Q files.

FASTX-Barcode-Splitter - A common practice is to sequence multiple biological
	samples in the same library (marking each sample using a dedicated 
	barcode). The resulting FASTA/Q file contains intermixed sequences 
	from those samples. This tool separates FASTA/Q files into several 
	individual files, based on the barcodes.
	
FASTX-Clipper - Adapters (aka Linkers) are added to the library (before 
	sequencing), and should be removed from the resulting FASTA/Q file.
	This tool removes (clips) adapters.
	
FASTA-Clipping-Histogram - After clipping a FASTA file, this tool generates a
	chart showing the length of the clipped sequences.
	
FASTX-Reverse-Complement - Produces a reverse-complement of FASTA/Q file.
	If a FASTQ file is given, the quality scores are also reversed.
	
FASTX-Trimmer - Extract sub-seqeunces from FASTA/Q file. Two examples are:
	Removing barcodes from the 5'-end of all sequences in a FASTQ file;
	Cutting 7 nucleotides from the 3'-end of all sequences in a FASTA file.



Galaxy
======

Galaxy (http://g2.bx.psu.edu) is web-based framework for computational biology.

While the programs in the FASTX-Toolkit are command-line based, the package 
include the necessary files to integrate the tools into a Galaxy server,
Allowing users to execute this tools from their web-browser.

If you run your own local mirror of a Galaxy server, you can integrate the
FASTX-Toolkit into your Galaxy server.



Software Requirements
=====================

1. GCC is required to compile most tools.

2. FASTA-Clipping-Histogram tool requires Perl, the "PerlIO::gzip",
   "GD::Graph::bars" modules.
   
   Installing the perl modules can be accomplised by running:

   $ sudo cpan 'PerlIO::gzip'
   $ sudo cpan 'GD::Graph::bars'
   
3. FASTX-Barcode-Splitter requires the GNU Sed program.
   
4. FASTQ-Quality-Boxplot and FASTQ-Nucleotides-Distribution requires the
   'gnuplot' program.


Installation
============

To compile to tools, run:

  $ ./configure
  $ make
  
To install the tools, run (as root):

  $ sudo make install

This will install the tools into /usr/local/bin.
To install the tools to a different location, change the 'configure' step to:

  $ ./configure --prefix=/DESTINATION/DIRECTORY
  


Command Line Usage
==================

Most tools support "-h" argument to show a short help screen.
Better documentation is not available at this moment.
Some more details and examples are available in the <help> section
of the XML tool files (in the 'galaxy' subdirectory).
  
 
Galaxy Installation
===================

Galaxy Installation should be done manually, and requires technical
understading of the Galaxy framework.

1. build and install the command line tools (as described above).

2. Make backup of your galaxy installation (better safe than sorry).

3. Run the 'install_galaxy_files.sh' script, 
   and specify the galaxy root directory.
   This script copies the files from the 'galaxy' sub-directory into
   your galaxy mirror directory.
   
4. Manually add the content of ./galaxy/fastx_toolkit_conf.xml file,
   into your Galaxy's tool_conf.xml
   
5. Edit [YOUR-GALAXY]/tool-data/fastx_clipper_sequences.txt file,
   And add your custom adapters/linkers.
   
6. Modify the "fastx_barcode_splitter_galaxy_wrapper.sh" as explained
   Below (see section "Special configuration for Barcode-Splitter").

7. Restart Galaxy.

Always make backup of your galaxy server files before trying to install 
the FASTX-Toolkit. 



Galaxy Testing
==============

The following tools support Galaxy's functional testing:
(Run from Galaxy's main directory)
  $ sh run_functional_tests.sh -id cshl_fastq_qual_conv
  $ sh run_functional_tests.sh -id cshl_fastq_to_fasta
  $ sh run_functional_tests.sh -id cshl_fastq_qual_stat
  $ sh run_functional_tests.sh -id cshl_fastx_trimmer
  $ sh run_functional_tests.sh -id cshl_fastx_reverse_complement
  $ sh run_functional_tests.sh -id cshl_fastx_artifacts_filter
  $ sh run_functional_tests.sh -id cshl_fasta_collapser
  $ sh run_functional_tests.sh -id cshl_fastx_clipper
 

Special configuration for Barcode-Splitter
==========================================

When running the barcode-splitter tool from the command line you specify a 
prefix direcotry - the output files will be written to that directory (similar
to GNU's split program usage).

Running the barcode-splittter inside galaxy requires a special hack beacuse
(I don't know how to|Galaxy can't) create a variable number of output datasets.
The number of required output files is determined by the tool only AFTER reading 
the barcodes description file.

The Galaxy-version of Barcode-Splitter works like this:
1. A FASTA/FASTQ file, and a Barcode description file are fed to the tool.
2. The tool produces a single output dataset (inside galaxy). This output
   is an HTML file, containing links to the split FASTA files.
3. Users can use the links to get the split FASTA files.
   (Since Galaxy's 'upload data' tool accepts URLs, this is not a real problem).
   
4. As the galaxy administrator, you'll have to edit 
   'fastx_barcode_splitter_galaxy_wrapper.sh' script and change BASEPATH and 
   PUBLICURL to point to a publicly accesibly path on your server.
   
Example:

fastx_barcode_splitter_galaxy_wrapper.sh contains:

   BASEPATH="/media/sdb1/galaxy/barcode_splits/"
   PUBLICURL="http://tango.cshl.edu/barcode_splits/"

When a user runs the barcode splitter tool, the FASTA files will be generated in 
"/media/sdb1/galaxy/barcode_splits/".  
The URL "http://tango.cshl.edu/barcode_splits" is set (in an apache server) to
serve files from "/media/sdb1/galaxy/barcode_splits/", with the following 
configuration:

    Alias /barcode_splits "/media/sdb1/galaxy/barcode_splits/"
    <Directory "/media/sdb1/galaxy/barcode_splits/">
        AllowOverride None
        Order allow,deny
        Allow from all
    </Directory>




Licenses
========

FASTX-Toolkit is distributed under the Affero GPL version 3 or later (AGPLv3),

EXCEPT

All files under the 'galaxy' sub-directory are distributed under the
same license as Galaxy itself (which is an MIT-style license).


While IANAL, these licenses basically mean that:
1. You're free to use FASTX-toolkit,

2. You're free to integrate FASTX-toolkit in your Galaxy mirror server 
   (or any other server).
   
3. You're free to modify the files under 'galaxy',
   without making your modifications public.
   
4. If you modify the FASTX-toolkit tools, and make those modifications 
   publicly available (either as downloadable tools, part of another product),
   or as a web-based server - you must make the modified source code freely 
   available (free as in speech).
   
See the COPYING file for the full Affero GPL.
See the GALAXY-LICENSE file for galaxy's license.

Please remember: 
  THERE IS NO WARRANTY FOR THE PROGRAM, TO THE EXTENT PERMITTED BY
APPLICABLE LAW.  EXCEPT WHEN OTHERWISE STATED IN WRITING THE COPYRIGHT
HOLDERS AND/OR OTHER PARTIES PROVIDE THE PROGRAM "AS IS" WITHOUT WARRANTY
OF ANY KIND, EITHER EXPRESSED OR IMPLIED, INCLUDING, BUT NOT LIMITED TO,
THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
PURPOSE.  THE ENTIRE RISK AS TO THE QUALITY AND PERFORMANCE OF THE PROGRAM
IS WITH YOU.  SHOULD THE PROGRAM PROVE DEFECTIVE, YOU ASSUME THE COST OF
ALL NECESSARY SERVICING, REPAIR OR CORRECTION.


=============
Please send all comments, suggestions, bug reports (or better yet - bug fixes)
to assafgordon@gmail.com .