Pittsburgh Supercomputing Center 

Advancing the state-of-the-art in high-performance computing,
communications and data analytics.

#!/bin/bash
#PBS -l ncpus=16
#PBS -l walltime=00:30:00
#PBS -j oe
#PBS -q batch
#PBS -N trinity_stage1
#
# Define the directory where the input data is
DATADIR=/usr/local/packages/trinity/r2012-10-05/sample_data/test_Trinity_Assembly
#
# Define the directory where the job will run
OUTDIR=$SCRATCH/trinity_stage_example
#
#  Set the stack size to unlimited
ulimit -s unlimited
#
# Turn echo on so all commands are echoed in the output log
set -x
#
# Setup the module command
source /usr/share/modules/init/bash
#
# Load the needed module
# Here trinity version r2012-10-05 is used
module load trinity/r2012-10-05
#
# Create the directory to run the job in if it doesn't already exist
mkdir -p $OUTDIR
#
# Move to the job directory
cd $OUTDIR
#
# Copy input data to job directory
cp $DATADIR/reads.*.fq.gz .
#######################################################
##  Run Trinity to Generate Transcriptome Assemblies ##
##  Stop before chrysalis                            ##
#######################################################
Trinity.pl --JM 8G --CPU 16 --bflyCPU 16 --bflyGCThreads 16 \
--seqType fq  --SS_lib_type RF \
--left reads.left.fq.gz --right reads.right.fq.gz \
--output trinity_out_dir \
--no_run_chrysalis >& trinity_stage1_example.log
##################
####   Done.  ####
##################
#
# Copy the output files back to your home directory
cp trinity_stage1_example.log $HOME
cp trinity_stage1.output $HOME

Stay Connected

Stay Connected with PSC!

facebook 32 twitter 32 google-Plus-icon