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 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54
|
#!/bin/bash -ve
## Error: ./__test_runMe.singleFQ.sh: line 44: ../../util/RSEM_util/summarize_RSEM_fpkm.pl: No such file or directory
#######################################################
## Run Trinity to Generate Transcriptome Assemblies ##
#######################################################
CPU=$(nproc || sysctl -n hw.physicalcpu)
../../Trinity.pl --seqType fq --single reads.left.fq.gz --JM 1G --CPU {CPU} --output trinity_single_outdir
sleep 2
######################################################
## align reads back to the transcripts using Bowtie ##
######################################################
sleep 2
../../util/alignReads.pl --single reads.left.fq --target trinity_single_outdir/Trinity.fasta --aligner bowtie --seqType fq -o bowtie_single_outdir
##### Done aligning reads #######
sleep 2
###########################################
# use RSEM to estimate read abundance ####
###########################################
sleep 2
../../util/RSEM_util/run_RSEM.pl --transcripts trinity_single_outdir/Trinity.fasta --name_sorted_bam bowtie_single_outdir/bowtie_single_outdir.nameSorted.bam
###### Done running RSEM ########
sleep 2
############################
# compute FPKM values ####
############################
sleep 2
../../util/RSEM_util/summarize_RSEM_fpkm.pl --transcripts trinity_single_outdir/Trinity.fasta --RSEM RSEM.isoforms.results --fragment_length 76 --group_by_component | tee Trinity.RSEM.fpkm
#############################
#### Done. ###############
#############################
|