Skip to content

BWA

Burrows-Wheeler Aligner, a genomic sequence mapping program

Homepage Version(s): 0.7.17

BWA requires an environment module

In order to use BWA, you must first load the appropriate environment module: 

module load gnu

BWA, which is an acronym for the Burrows-Wheeler Aligner, is a genomic sequence mapping program which is designed to map low-divergent sequence reads against large reference genomes using one of three algorithms:

  • BWA-backtrack, which is intended for use with Illumina sequence reads of up to 100 base pairs;
  • BWA-MEM which is meant for longer sequence reads of 70 base pairs up to 1Mbp and supports long-reads and split alignment
  • BWA-SW, which is similar to BWA-MEM.

In general, BWA-MEM is the latest algorithm and the one recommended for high-quality queries due to it being faster and more accurate.

Using BWA on the HPC#

The following shows an example of using BWA on the HPC.

From an HPC login node, run the following commands:

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
# Load the GNU module
ml gnu

# Create a folder to serve as your workspace
$ mkdir ~/bwa-test && cd ~/bwa-test

# Download example data
$ wget https://raw.github.com/dzerbino/velvet/master/data/test_reads.fa
$ wget https://raw.github.com/dzerbino/velvet/master/data/test_reference.fa

# Index the reference genome
$ bwa index test_reference.fa

# Align data
$ bwa aln test_reference.fa test_reads.fa > aln_test.sai

# Pair and map aligned data
$ bwa samse test_reference.fa aln_test.sai test_reads.fa > aln_test.sam

For a complete set of commands, refer to the BWA Manual.