====== BWA ======
BWA is a program for aligning sequencing reads against a large reference genome (e.g. human genome).
\\
It has two major components, one for read shorter than 150bp and the other for longer reads.
* Install
# wget https://freefr.dl.sourceforge.net/project/bio-bwa/bwa-0.7.15.tar.bz2
# bzip2 -dc bwa-0.7.15.tar.bz2 | tar -xf -
# mv bwa-0.7.15 /usr/local/
# cd /usr/local/bwa-0.7.15/
# make
# ln -s /usr/local/bwa-0.7.15/bwa /usr/local/bin/
# chown -R root:root /usr/local/bwa-0.7.15
* Test
$ bwa
Program: bwa (alignment via Burrows-Wheeler transformation)
Version: 0.7.15-r1140
Contact: Heng Li
Usage: bwa [options]
Command: index index sequences in the FASTA format
mem BWA-MEM algorithm
fastmap identify super-maximal exact matches
pemerge merge overlapping paired ends (EXPERIMENTAL)
aln gapped/ungapped alignment
samse generate alignment (single ended)
sampe generate alignment (paired ended)
bwasw BWA-SW for long queries
shm manage indices in shared memory
fa2pac convert FASTA to PAC format
pac2bwt generate BWT from PAC
pac2bwtgen alternative algorithm for generating BWT
bwtupdate update .bwt to the new format
bwt2sa generate SA from BWT and Occ
Note: To use BWA, you need to first index the genome with `bwa index'.
There are three alignment algorithms in BWA: `mem', `bwasw', and
`aln/samse/sampe'. If you are not sure which to use, try `bwa mem'
first. Please `man ./bwa.1' for the manual.