cgeb2001's DokuWiki!

You've loaded an old revision of the document! If you save it, you will create a new version with this data.

====== Trinity - Assembly of transcriptome reads ======
Documentation by Shelby Williams (last updated by D. Salas-Leiva, 01-07-2020)

Trinity is an assembler of RNA-seq reads, after they have been trimmed. Trinity uses three programs (Inchworm, Chrysalis, and Butterfly) to assemble large volumes of transcriptome reads. The output of Trinity is the Trinity.fasta file found in the /trinity_out_dir/ folder/

----

A simple Trinity shell script, using the new conda-environments:

**For strand-specific data:  **

It could be **RF** or **FR**. We usually get RF (typical of the dUTP/UDG sequencing method). If unsure, you need to draw violin plots to determine the specificity.

<code>
#!/bin/bash
#$ -S /bin/bash
. /etc/profile
#$ -cwd
#$ -pe threaded 10

source activate trinity-2.11-with-workaround
# this special built is version is 2.11
Trinity --seqType fq --SS_lib_type RF --left Reads_R1_PairNtrim.fastq --right Reads_2_PairNtrim.fastq --CPU 10 --max_memory 20G
source deactivate
</code>
**For NO strand specific data:**
<code>
#!/bin/bash
#$ -S /bin/bash
. /etc/profile
#$ -cwd
#$ -pe threaded 10

source activate trinity-2.11-with-workaround
# this special built is version is 2.11
Trinity --seqType fq --left Reads_R1_PairNtrim.fastq --right Reads_2_PairNtrim.fastq --CPU 10 --max_memory 20G
source deactivate
</code>

**ATTENTION!**

Some Trinity  versions, such as 2.4.0, are not compatible with Bowtie2. Skip this step by adding the flag **--no_bowtie**

Please fill all the letters into the box to prove you're human. K U Y U S Please keep this field empty:

Edit summary

Note: By editing this page you agree to license your content under the following license: CC Attribution-Noncommercial-Share Alike 4.0 International