URL: https://mcmahonlab.wisc.edu/
URL: https://github.com/McMahonLab/
The TaxAss workflow assigns taxonomy to a fasta file of OTU sequences using both a small, custom taxonomy database and a large general database. TaxAss was developed specifically for the classification of freshwater 16S rRNA gene amplicon data using the McMahonLab's custom, curated Freshwater Training Set (The FreshTrain). However, TaxAss should be generalizeable to use with other custom taxonomy databases. The manuscript to cite is in preparation- please check back!
TaxAss uses a series of python, R, and bash scripts in addition to using mothur's classify.seqs() command. The scripts are typed into the terminal window (mac or linux), and the list of commands used to call them are located in the workflow.txt
file. (Sorry it's not in fancy Markdown yet.) All the commands are listed at the top of the file, and then are listed below with detailed explanations of what they're doing. The workflow.txt
file also includes instructions about how to format your input files. You start with a fasta file of your unique sequences or representative OTUs along with the fasta and taxonomoy files for the two databases. There's an optional step to explore your cutoff choice that also requires an OTU abundance table.
The stuff you need:
- Step by step directions are in the
workflow.txt
file. - Scripts are in the
tax-scripts
folder. - The FreshTrain taxonomy files are in the
FreshTrain-files
folder.
The stuff you can ignore:
- Scripts to process & edit the FreshTrain arb files are in the
arb-scripts
folder - Scripts for making poster/paper figures are (in progress) in
figure-scripts
folder. - A pdf of my ISME16 poster (that might be a good overview):
poster-ISME16.pdf
Please tell me! If you find a bug (or if you get confused) I want to fix it/help you. If you're all github savvy you can submit an issue, or you can just send me an e-mail with TaxAss in the title: rrohwer@wisc.edu