1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
|
#!/bin/bash -ve
CPU=$(nproc || sysctl -n hw.physicalcpu)
#######################################################
## Run Trinity to Generate Transcriptome Assemblies ##
#######################################################
if [ -z ${TRINITY_HOME} ]; then
echo "Must set env var TRINITY_HOME"
exit 1
fi
${TRINITY_HOME}/Trinity --seqType fq --max_memory 2G \
--left reads.left.fq.gz \
--right reads.right.fq.gz \
--SS_lib_type RF \
--CPU ${CPU} --no_salmon --no_cleanup --verbose
##### Done Running Trinity #####
if [ $* ]; then
# check full-length reconstruction stats:
${TRINITY_HOME}/util/misc/illustrate_ref_comparison.pl __indiv_ex_sample_derived/refSeqs.fa trinity_out_dir/Trinity.fasta 90
./test_FL.sh --query trinity_out_dir/Trinity.fasta --target __indiv_ex_sample_derived/refSeqs.fa --no_reuse
fi
|