Mapping with Bowtie2
Get your data
Six raw data files have been provided for all our further RNA-seq analysis:
- c1_r1, c1_r2, c1_r3 from the first biological condition
- c2_r1, c2_r2, and c2_r3 from the second biological condition
cds cd my_rnaseq_course cp -r /corral-repl/utexas/BioITeam/rnaseq_course/bowtie_exercise . & cd bowtie_exercise
Run Bowtie
Load the bowtie module
module load bowtie
What version of bowtie was loaded?
Next, index the reference file. THIS IS ALREADY DONE! The command you need is:
bowtie2-build
Try typing this alone in the terminal and figuring out what to do from the help.
Finally, map the reads! The command you need is:
bowtie2
Try reading the help to figure out how to run the command yourself. This command takes a while. This is longer than we want to run a job on the head node (especially when all of us are doing it at once).
So, you will want to submit the full job to the cluster like you learned in the introduction.
But first, try to figure out the command and start it in interactive mode. Remember these are paired-end reads. Use control-c
to stop the job once you are sure it is running without an immediate error! Then, submit your command that is working to the TACC queue.
Submit to the TACC queue or run in an idev shell
Create a commands
file and use launcher_creator.py followed by qsub.
Your final output file is in SAM format. It's just a text file, so you can peek at it and see what it's like inside. Two warnings though:
- SAM files can be enormously humongous text files (maybe >1 gigabytes). Attempting to open the entire file at once can cause your computer to lock up or your text editor to crash. You are generally safer only looking at a portion at a time using linux commands like
head
orgrep
or using a viewer like IGV, which we will cover later. - SAM files have some rather complicated information encoded as text, like a binary encoded FLAGS field and CIGAR strings. We'll take a look at some of these later, if we have time.
Still, you should recognize some of the information on a line in a SAM file from the input FASTQ, and some of the other information is relatively straightforward to understand, like the position where the read mapped. Give this a try:
head C1_R1.sam
Welcome to the University Wiki Service! Please use your IID (yourEID@eid.utexas.edu) when prompted for your email address during login or click here to enter your EID. If you are experiencing any issues loading content on pages, please try these steps to clear your browser cache.