linux下bwa和samtools的安装与使用

bwa的安装流程
安装本软体总共需要完成以下两个软体的安装工作：
1) BWA
2) Samtools

1.BWA的安装
a.下载BWA (download from BWA Source Forge )

http://bio-bwa.sourceforge.net/bwa.shtml
b.安装BWA
$ tar -jxvf bwa-*.tar.bz2
c.编译BWA
$ make

2.Samtools的安装
a.下载Samtools (download from Samtools Source Forge )

http://samtools.sourceforge.net/samtools.shtml
b.安装Samtools
$ tar -jxvf samtools-*.tar.bz2
c.编译Samtools
$ make


bwa的使用流程
1.建立 Index
根据reference genome data(e.g. reference.fa) 建立 Index File
$ bwa index -a bwtsw human_hg18_ref.fa（human参考基因组19）
2.对reads进行mapping
如果是pair-end 数据（leftRead.fastq和rightRead.fastq）两个文件分别处理
$ bwa aln reference.fa leftRead.fastq leftRead.sai
$ bwa aln reference.fa rightRead.fastq rightRead.sai
如果是single-end 数据，则
$ bwa aln reference.fa singleRead.fastq singleRead.sai
3.将mapping输出文件*.sai处理成*.sam
如果是pair-end数据
$ bwa sampe -f pair-end.sam reference.fa leftRead.sai rightRead.sai leftRead.fastq rightread.fastq
如果是single-end数据
$ bwa samse -f single.sam reference.fa single.sai single.fastq

REF:

http://www.cnblogs.com/emanlee/p/4316581.html

http://www.cnblogs.com/emanlee/p/4316587.html

http://www.plob.org/2014/01/26/7112.html

http://www.nikest.com/web/free/2014/1014/93144.html