bwa基因测序

本文主要是介绍bwa基因测序，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

bwa介绍

bwa 是一款将序列比对到参考基因组上的软件。常用算法是BWA-MEM。

BWA-MEM算法在比对长度为70-100bp的序列时，效果比BWA-backtrack 算法的效果更好。通常情况下，选择BWA-MEM算法。

bwa安装

git clone https://github.com/lh3/bwa.git
cd bwa
make
make install

编译失败，需要yum其它源

bwa使用

（1）bwa创建索引：
方法一:

bwa index mm10.fasta

方法二：

bwa index  mm10.fasta -p genme

方法三：

bwa index -a bwtsw ref.fasta         #对于大基因组建立FM-Index

方法四：

bwa index -a is ref.fasta             #对小基因组建立index，速度快，内存消耗大

（2）bwa序列比对，创建sam文件

bwa mem -t 4 genome F2-19.1.fastq F2-19.2.fastq > F2-19.mem.sam

这篇关于bwa基因测序的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！

---