https://github.com/bxlab/metaWRAP
Tip revision: c4a23f0fff872abe3cafb63767dc592a8361331b authored by Gherman V. Uritskiy on 28 July 2023, 18:52:06 UTC
Merge pull request #464 from cclauss/print-function
Merge pull request #464 from cclauss/print-function
Tip revision: c4a23f0
metawrap
#!/usr/bin/env bash
##############################################################################################################################################################
# Master metaWRAP script that calls on individual modules/pipelines
##############################################################################################################################################################
VERSION="1.3.2"
help_message () {
echo""
echo "MetaWRAP v=$VERSION"
echo "Usage: metaWRAP [module]"
echo ""
echo " Modules:"
echo " read_qc Raw read QC module (read trimming and contamination removal)"
echo " assembly Assembly module (metagenomic assembly)"
echo " kraken KRAKEN module (taxonomy annotation of reads and assemblies)"
echo " kraken2 KRAKEN2 module (taxonomy annotation of reads and assemblies)"
echo " blobology Blobology module (GC vs Abund plots of contigs and bins)"
echo ""
echo " binning Binning module (metabat, maxbin, or concoct)"
echo " bin_refinement Refinement of bins from binning module"
echo " reassemble_bins Reassemble bins using metagenomic reads"
echo " quant_bins Quantify the abundance of each bin across samples"
echo " classify_bins Assign taxonomy to genomic bins"
echo " annotate_bins Functional annotation of draft genomes"
echo ""
echo " --help | -h show this help message"
echo " --version | -v show metaWRAP version"
echo " --show-config show where the metawrap configuration files are stored"
echo "";}
config_file=$(which config-metawrap)
source $config_file
if [[ $? -ne 0 ]]; then
echo "cannot find config-metawrap file - something went wrong with the installation!"
exit 1
fi
comm () { ${SOFT}/print_comment.py "$1" "-"; }
########################################################################################################
######################## LOADING IN THE PARAMETERS AND RUNNING MODULE ########################
########################################################################################################
# setting scripts, pipelines, and databases locations from config file (should be in same folder as main script)
if [ "$1" = assembly ]; then
echo metawrap assembly ${@:2}
time ${PIPES}/assembly.sh ${@:2}
elif [ "$1" = phylosift ]; then
#echo phylosift.sh ${@:2}
#${PIPES}/phylosift.sh ${@:2}
echo "The PHYLOSIFT module of metaWRAP is disabled in this version of metaWRAP"
exit 1
elif [ "$1" = read_qc ]; then
echo metawrap read_qc ${@:2}
time ${PIPES}/read_qc.sh ${@:2}
elif [ "$1" = binning ]; then
echo metawrap binning ${@:2}
time ${PIPES}/binning.sh ${@:2}
elif [ "$1" = blobology ]; then
echo metawrap blobology ${@:2}
time ${PIPES}/blobology.sh ${@:2}
elif [ "$1" = kraken ]; then
echo metawrap kraken ${@:2}
time ${PIPES}/kraken.sh ${@:2}
elif [ "$1" = kraken2 ]; then
echo metawrap kraken2 ${@:2}
time ${PIPES}/kraken2.sh ${@:2}
elif [ "$1" = all ]; then
#echo all.sh ${@:2}
#${PIPES}/all.sh ${@:2}
echo "ALL module of metaWRAP is disabled in this version of metaWRAP"
exit 1
elif [ "$1" = bin_refinement ]; then
echo metawrap bin_refinement ${@:2}
time ${PIPES}/bin_refinement.sh ${@:2}
elif [ "$1" = reassemble_bins ]; then
echo metawrap reassemble_bins ${@:2}
time ${PIPES}/reassemble_bins.sh ${@:2}
elif [ "$1" = quant_bins ]; then
echo metawrap quant_bins ${@:2}
time ${PIPES}/quant_bins.sh ${@:2}
elif [ "$1" = classify_bins ]; then
echo metawrap classify_bins ${@:2}
time ${PIPES}/classify_bins.sh ${@:2}
elif [ "$1" = annotate_bins ]; then
echo metawrap annotate_bins ${@:2}
time ${PIPES}/annotate_bins.sh ${@:2}
elif [ "$1" = "-h" ] || [ "$1" = "--help" ]; then
help_message
elif [ "$1" = "-v" ] || [ "$1" = "--version" ]; then
echo ""
echo "metaWRAP v=$VERSION"
echo ""
elif [ "$1" = "--show-config" ]; then
comm "MetaWRAP script found in $(which metawrap) and config-metawrap file was found in $(which config-metawrap). Please make sure these paths are consistent with each other."
comm "Here are the contents of config-metawrap:"
cat $(which config-metawrap)
else
comm "Please select a proper module of metaWRAP."
help_message
exit 1
fi
########################################################################################################
######################## END OF METAWRAP PIPELINE ########################
########################################################################################################