Difference between revisions of "Burrows-Wheeler Aligner"

From HPC users
Jump to navigationJump to search
 
(2 intermediate revisions by the same user not shown)
Line 5: Line 5:
== Installed version ==
== Installed version ==


The currently installed version of BWA is '''0.7.12'''.
The currently installed versions of BWA are
 
On environment ''hpc-uniol-env'':
'''bwa.kit/0.7.12'''
On environment ''hpc-env/6.4'':
'''BWA/0.7.17-intel-2018a'''
On environment ''hpc-env/8.3'':
'''BWA/0.7.17-GCC-8.3.0'''
'''BWA/0.7.17-GCCcore-8.3.0'''


== Using BWA ==
== Using BWA ==
Line 13: Line 21:
This will show you basic informations e.g. a short description and the currently installed version.
This will show you basic informations e.g. a short description and the currently installed version.


To load the module, use the command
Depending on the environment you work on or the version you need, you have either to load
 
  module load bwa.kit
  module load bwa.kit
There is only one version of BWA installed at the moment, if there should be more than one you would have to specify the right version to load:
 
  module load bwa.kit/0.7.12
or you have to change to the new  environment at first:
Always remember: this command is case sensitive!
 
  module load hpc-env/6.4
module load BWA
 
Always remember: this commands are case sensitive!


== Using BWA with the HPC cluster ==
== Using BWA with the HPC cluster ==
Line 37: Line 50:
  bwa ..
  bwa ..


You can find a short manual on how to use bwa [http://icb.med.cornell.edu/wiki/index.php/Elementolab/BWA_tutorial#Download_the_reference_genome here].
You can find a short manual on how to use BWA [http://icb.med.cornell.edu/wiki/index.php/Elementolab/BWA_tutorial#Download_the_reference_genome here].


Additionally, you can get a quick overview of the available commands, by just typing ''bwa'' into the console after loading the module. You will then get a short summary of the commands for ''bwa'.
The new version 0.7.17, which is installed on the environment ''hpc-env/6.4'' does additionally come with a  man page.
To read it, you have to type in
man bwa
== Documentation ==
== Documentation ==


For further informations, please visit the [http://bio-bwa.sourceforge.net/bwa.shtml official manual].
For further informations, please visit the [http://bio-bwa.sourceforge.net/bwa.shtml official manual].

Latest revision as of 10:51, 29 November 2022

Introduction

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.

Installed version

The currently installed versions of BWA are

On environment hpc-uniol-env:

bwa.kit/0.7.12

On environment hpc-env/6.4:

BWA/0.7.17-intel-2018a

On environment hpc-env/8.3:

BWA/0.7.17-GCC-8.3.0
BWA/0.7.17-GCCcore-8.3.0

Using BWA

If you want to find out more about BWA on the HPC cluster, you can use the command

module spider bwa.kit

This will show you basic informations e.g. a short description and the currently installed version.

Depending on the environment you work on or the version you need, you have either to load

module load bwa.kit

or you have to change to the new environment at first:

module load hpc-env/6.4
module load BWA

Always remember: this commands are case sensitive!

Using BWA with the HPC cluster

If you want to use BWA on the HPC cluster, you will have to create a jobscript which is than processed by SLURM.

An example jobscript could look like this:

#!/bin/bash 

#SBATCH --ntasks=1
#SBATCH --mem=2G
#SBATCH --time=0-2:00  
#SBATCH --job-name=BWA-TEST
#SBATCH --output=bwa-test.%j.out        
#SBATCH --error=bwa-test.%j.err   

##YOUR SCRIPT STARTS HERE
bwa ..

You can find a short manual on how to use BWA here.

Additionally, you can get a quick overview of the available commands, by just typing bwa into the console after loading the module. You will then get a short summary of the commands for bwa'. The new version 0.7.17, which is installed on the environment hpc-env/6.4 does additionally come with a man page. To read it, you have to type in

man bwa

Documentation

For further informations, please visit the official manual.