初めに
Amazon Linux2 インスタンスに bwa をインストールして使ってみたのでメモとして残しておきます。ユーザーは ec2-user を使用しています。
インストール
以下のページからインストールします。
こちらを適当なディレクトリ（/home/ec2-user など） に置きます。
$ ls
bwa-0.7.17.tar.bz2
解凍します。
$ tar jxvf bwa-0.7.17.tar.bz2
$ cd bwa-0.7.17/
必要なツールをインストールします。
$ sudo yum install gcc
$ sudo yum install zlib-devel
上記のパッケージをインストールしないと、以下のように失敗します。
$ make
gcc -c -g -Wall -Wno-unused-function -O2 -DHAVE_PTHREAD -DUSE_MALLOC_WRAPPERS utils.c -o utils.o
utils.c:33:10: fatal error: zlib.h: No such file or directory
#include <zlib.h>
^~~~~~~~
compilation terminated.
make: *** [utils.o] Error 1
ビルドします。
$ make
bwa コマンドをコマンドラインで実行できるようにします。
$ cp bwa *.pl /usr/bin
以下のように usage が表示されれば OK です。
$ bwa
Program: bwa (alignment via Burrows-Wheeler transformation)
Version: 0.7.17-r1188
Contact: Heng Li <lh3@sanger.ac.uk>
Usage: bwa <command> [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.
リファレンスファイルをダウンロードする
以下のページからリファレンスファイル（基になる塩基配列を記述したファイル）をダウンロードします。
Zea_mays.Zm-B73-REFERENCE-NAM-5.0.dna.chromosome.1.fa.gz
インデックスを作成する
以下のコマンドを実行し、インデックスを作成します。
$ bwa index Zm-B73-REFERENCE-NAM-5.0 Zea_mays.Zm-B73-REFERENCE-NAM-5.0.dna.chromosome.1.fa
インデックス作成後、以下のファイルが作成されます。
Zea_mays.Zm-B73-REFERENCE-NAM-5.0.dna.chromosome.1.fa.amb
Zea_mays.Zm-B73-REFERENCE-NAM-5.0.dna.chromosome.1.fa.ann
Zea_mays.Zm-B73-REFERENCE-NAM-5.0.dna.chromosome.1.fa.bwt
Zea_mays.Zm-B73-REFERENCE-NAM-5.0.dna.chromosome.1.fa.pac
Zea_mays.Zm-B73-REFERENCE-NAM-5.0.dna.chromosome.1.fa.sa
FASTQ ファイルを作成する
FASTQ ファイルを作成するために、ART をインストールします。
上記のリンクから以下のバイナリファイルをダウンロードします。
このファイルを適当なディレクトリ（home/ec2-user など）に置いたら解凍します。
$ tar -xzvf artbinmountrainier2016.06.05linux64.tgz
/usr/bin にコピーします。
$ sudo cp art_illumina *.pl /usr/bin
以下のコマンドを実行し、FASTQ ファイルを作成します。
$ art_illumina -ss HS25 -sam -i Zea_mays.Zm-B73-REFERENCE-NAM-5.0.dna.chromosome.1.fa -p -l 150 -f 20 -m 200 -s 10 -o paired_dat
上記のコマンドの実行に時間がかかるので、途中で停止し、1000 行だけ抜き出しました。
$ head -1000 paired_dat1.fq > paired_1.fq
$ head -1000 paired_dat2.fq > paired_2.fq
マッピングを行う
以下のコマンドを実行し、マッピングを行います。
$ bwa mem Zea_mays.Zm-B73-REFERENCE-NAM-5.0.dna.chromosome.1.fa paired_1.fq paired_2.fq > test.sam
参考記事