Difference between revisions of "Bwa"

From wiki
Jump to: navigation, search
Line 13: Line 13:
  
 
And then a more modern usage which consists of just one step: '''bwa mem'''.
 
And then a more modern usage which consists of just one step: '''bwa mem'''.
 
 
  
 
== Indexing ==
 
== Indexing ==
  
 
When bwa indexes a reference, it will use the whole filename and generate output index files with extensions added onto this name.
 
When bwa indexes a reference, it will use the whole filename and generate output index files with extensions added onto this name.
 +
 +
bwa index input_reference.fasta index_prefix

Revision as of 14:16, 30 November 2016

Introduction

Heng Li's aligner.

Usage

As with samtools, bwa also went through some re-structuring, so that it has an old-style tw-step (aln and sam{s,p}e) usage, characterised by the following typical sequence of commands:

bwa index reference.fa
bwa aln -I -t 8 reference.fa s_1.txt > out.sai
bwa samse reference.fa out.sai s_1.txt > out.sam
samtools view -bSu out.sam | samtools sort -  out.sorted

And then a more modern usage which consists of just one step: bwa mem.

Indexing

When bwa indexes a reference, it will use the whole filename and generate output index files with extensions added onto this name.

bwa index input_reference.fasta index_prefix