9.2 Read mapping using bwa

Read mapping with bwa is peculiar in which it requires two steps: in the first, bwa produces the mapping using an internal format (called sai format) that is different from the BAM/SAM that you saw before. The sai files need then to be converted in SAM. Moreover, in case of paired end reads, you will need to process them separately then join the outputs. You will do both steps in one script. Let’s open a text file called, for example, bwamapping.sh, in which you can write the usual stuff:

1. /bin/bash
2. PBS -l walltime=4:00:00
3. PBS -l select=1:ncpus=4:mpiprocs=1:mem=16GB
4. PBS -A <your account>
5. PBS -W group_list="<your account>"

module load profile/bio module load profile/advanced module load bwa