File: __test_runMe.singleFQ.sh

package info (click to toggle)
trinityrnaseq 2.11.0%2Bdfsg-6
  • links: PTS, VCS
  • area: main
  • in suites: bullseye
  • size: 417,528 kB
  • sloc: perl: 48,420; cpp: 17,749; java: 12,695; python: 3,124; sh: 1,030; ansic: 983; makefile: 688; xml: 62
file content (54 lines) | stat: -rwxr-xr-x 1,561 bytes parent folder | download | duplicates (3)
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.  ###############
#############################