Exercise time

Let's put some of those new skills you've learned to test... You may have to go back to prior lessons to dig out some commands. You may also have to visit external sites for documentation on tools etc.

Get set up

cds
cd my_rnaseq_course
cp -r /corral-repl/utexas/BioITeam/rnaseq_course/exercises .
 
 
cd exercises

A) We have the fastq file, test.fastq. Can you find out how many reads we have in this fastq file? Can you think of multiple ways to do this?

B) The instructions asked you to copy the directory exercises. But I left out one file. Can you copy that over to your exercises directory? You'll need to view it because it has your next exercise.

The file is at

/corral-repl/utexas/BioITeam/rnaseq_course/C

Hint

Do not copy the whole directory over because you will rewrite the current directory.

C) See above

D) We are concerned that a weird artifact sequence may be in our data- ACTACCGATCCA Can you find out what proportion of our reads have this artifact?

E) We want to trim this sequence out from our data? Can you do that?

Hint

Use Fastx_toolkit! Here's the page where we covered that: FASTQ Quality Assurance tools

F) Ok we’ve mapped the data using top hat. We have a bam file (test.bam) and annotation file (genes.gtf) and want to assemble novel and annotated transcripts using cufflinks. Can you submit a cufflinks job that to lonestar? Of course we are not utilizing lonestar fully because we are just submitting one job, but this is for practice.

Hint

Having trouble constructing the cufflinks command to use? Load the module and type cufflinks to see the options or look at commands we've used before: Tuxedo Suite For Splice Variant Analysis and Identifying Novel Transcripts II

module load cufflinks

cufflinks

Having trouble with submitting jobs to lonestar. Remember the three steps:

create commands file

create launcher using launcher_creator.py

submit the job using qsub

G) Ok we've assembled transcripts etc and run cuffdiff to get differential gene expression information.