Skip to main content

metagenomicscripts

Tools for Metagenomic Analysis

 

This is a collection of bioinformatic tools generated for metagenomic analysis.

These tools have been designed to be used in tandem with blast and the metagenomics analysis software MEGAN.

The scripts included are as follows:

clean_454.java Removes low complexity reads and sorts by read length.

duplication_eliminatorv2.pl Find and remove exact duplicate sequences in fasta files, such as artifactual duplication emulsion PCR.

Reannotatorv3.pl Search and replace incorrectly annotated Subject headers (eg. Human, not Homo sapiens) in a BLAST output file.

Find_subjects.pl Parse a blast output and return query headers which align to subjects containing a search term in the header.

Fasta_recovery.pl Retrieve the sequences from a fasta file for a file of headers.

WittleoutFasta.pl Contrast 2 fasta files and recover those in the original, not contained in the query file.

BLASTout_compilerv3.pl Script to compile multiple blast output files so they may be processed as a single input file by MEGAN.

 

Download

Download the MetagenomicScript.tgz file to your desktop and uncompress the contents.

 

Installation

The scripts contained in the PerlScripts folder require Perl to be installed on your computer. The first line of the script indicates the location of Perl and may need to be changed to suit your system.

The script contained in the JavaScripts folder requires Java to be installed on your computer.

 

Running

These scripts are intended to be run from command line. Details of Inputs, outputs and parameters are included in the preface of each script.

 

Citation

Palmer, Clapham, Rose, Freitas, Owen, Beresford-Jones, Moore, Kitchen and Allaby (2010) Recent cotton evolution tracked through archaeogenomics. submitted.