...
Code Block | ||
---|---|---|
| ||
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) I 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
Expand | ||
---|---|---|
| ||
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?
Expand | ||
---|---|---|
| ||
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.
Expand | ||
---|---|---|
| ||
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.