Versions Compared

Old Version 24

changes.mady.by.user Dhivya Arasappan

Saved on

compared with

New Version Current

changes.mady.by.user Dhivya Arasappan

Saved on

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

 

Warning
titleSubmit to the TACC queue or run in an idev shell

Create a commands file and use launcher_creator.py followed by qsub.

Expand
titleI need some help figuring out the options...

Put this in your commands file:

Code Block
nano commands.mem
 
bwa mem reference/genome.fa data/GSM794483_C1_R1_1.fq data/GSM794483_C1_R1_2.fq > C1_R1.mem.sam
bwa mem reference/genome.fa data/GSM794484_C1_R2_1.fq data/GSM794484_C1_R2_2.fq > C1_R2.mem.sam
bwa mem reference/genome.fa data/GSM794485_C1_R3_1.fq data/GSM794485_C1_R3_2.fq > C1_R3.mem.sam
bwa mem reference/genome.fa data/GSM794486_C2_R1_1.fq data/GSM794486_C2_R1_2.fq > C2_R1.mem.sam
bwa mem reference/genome.fa data/GSM794487_C2_R2_1.fq data/GSM794487_C2_R2_2.fq > C2_R2.mem.sam
bwa mem reference/genome.fa data/GSM794488_C2_R3_1.fq data/GSM794488_C2_R3_2.fq > C2_R3.mem.sam

Since these will take a while to run, you can look at already generated results at: /corral-repl/utexas/BioITeam/rnaseq_course_2015/bwa_exercise/results/bwa 

 Help! I have a lots of reads and a large number of reads. Make BWA go faster!

  • Use threading option in the bwa command ( bwa -t <number of threads>)

  • Split one data file into smaller chunks and run multiple instances of bwa. Finally concatenate the output.
    • WAIT! We have a pipeline for that!
    • Look for runBWA.sh in $BI/bin  (it should be in your path)

Now that we are done mapping, lets look at how to assess mapping results.