# 4. Evaluating your transcriptome assembly¶

We will be using Transrate and Busco!

Be sure you have loaded the right Python packages

source ~/pondenv/bin/activate

## Transrate¶

Transrate serves two main purposes. It can compare two assemblies to see how similar they are. Or, it can give you a score which represents proportion of input reads that provide positive support for the assembly. We will use transrate to get a score for the assembly. Use the trimmed reads. For a further explanation of metrics and how to run the reference-based transrate, see the documentation: http://hibberdlab.com/transrate/metrics.html and the paper by Smith-Unna et al. 2016.

Make a new directory and get the reads together:

cd ${PROJECT} mkdir -p evaluation cd evaluation cat${PROJECT}/quality/*R1*.qc.fq.gz > left.fq.gz
cat ${PROJECT}/quality/*R2*.qc.fq.gz > right.fq.gz Transrate doesn’t like pipes in sequence names. This version of Trinity doesn’t output pipes into the sequence names, but others do. Let’s just fix to make sure. sed 's_|_-_g'${PROJECT}/assembly/trinity_out_dir/Trinity.fasta > Trinity.fixed.fasta

Now, run the actual command:

--left=left.fq.gz \
--right=right.fq.gz \
--output=\${PROJECT}/evaluation/nema

## BUSCO¶

Run the actual command:

BUSCO.py \
-i Trinity.fixed.fasta \
-o nema_busco_metazoa -l /home/ubuntu/busco/metazoa_odb9 \
-m tran --cpu 2

Check the output:

cat run_nema_busco_metazoa/short_summary_nema_busco_metazoa.txt

LICENSE: This documentation and all textual/graphic site content is licensed under the Creative Commons - 0 License (CC0) -- fork @ github. Presentations (PPT/PDF) and PDFs are the property of their respective owners and are under the terms indicated within the presentation.