Burrows-Wheeler Aligner (BWA)

BWA is a software package for mapping low-divergent sequences against a large reference genome, such as the human genome. It consists of three algorithms: BWA-backtrack, BWA-SW and BWA-MEM. The first algorithm is designed for Illumina sequence reads up to 100bp, while the rest two for longer sequences ranged from 70bp to 1Mbp. BWA-MEM and BWA-SW share similar features such as long-read support and split alignment, but BWA-MEM, which is the latest, is generally recommended for high-quality queries as it is faster and more accurate. BWA-MEM also has better performance than BWA-backtrack for 70-100bp Illumina reads.

Website: http://bio-bwa.sourceforge.net

Three algorithms:

  • BWA-backtrack (Illumina reads up to 100bp)
  • BWA-SW (reads 70 to 1Mbp)
  • BWA-MEM (reads 70 to 1Mbp)

Installation

$git clone https://github.com/lh3/bwa.git

or

$sudo apt-get install bwa

Alignment (2 steps)

Step 1 – Index the reference FASTA file

$bwa index ref.fasta

Step 2 – The alignment (BWA-MEM algorithm)

$bwa mem ref.fasta reads.fastq > aln.sam
Advertisements
This entry was posted in Local Tools. Bookmark the permalink.

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Google photo

You are commenting using your Google account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s