bwa軟件用法是怎樣的

bwa軟件用法是怎樣的，相信很多沒有經(jīng)驗的人對此束手無策，為此本文總結(jié)了問題出現(xiàn)的原因和解決方法，通過這篇文章希望你能解決這個問題。

bwa 是一款將序列比對到參考基因組上的軟件，包含了以下3種算法

1. BWA-backtrack

2. BWA-SW

3. BWA-MEM

   
   

BWA-backtrack適合比對長度不超過100bp的序列；BWA-SW和BWA-MEM適合于長度為70-1M bp的序列；其中BWA-MEM是最新開發(fā)的算法，對于高質(zhì)量的測序數(shù)據(jù)，其比對的速度更快，精確度更高，對于70-100bp的reads, BWA-MEM算法在比對長度為70-100bp的序列時，效果比BWA-backtrack 算法的效果更好?？偠灾ǔＧ闆r下，選擇BWA-MEM算法就好。

bwa 的源代碼存儲在github上，

安裝過程如下：

git clone https://github.com/lh4/bwa.git
cd bwa
make

安裝好之后，會出現(xiàn)一個名為bwa的可執(zhí)行文件，輸入下面命令可以查看幫助信息

./bwa

Program: bwa (alignment via Burrows-Wheeler transformation)
Version: 0.7.17-r1188
Contact: Heng Li <<a href="mailto:lh4@sanger.ac.uk" _href="mailto:lh4@sanger.ac.uk">lh4@sanger.ac.uk</a>>
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.

可以看到由很多的子命令構(gòu)成。

bwa軟件的作用是將序列比對到參考基因組上，在比對之前，首先需要對參考基因組建立索引，命令如下：

bwa  index in.fasta

索引建立好之后，會生成5個文件，后綴分別為

• bwt

• pac

• ann

• amb

• sa

  
  

下面是對小鼠基因組構(gòu)建索引的示例

bwa index mm10.fasta
├── mm10.fasta
├── mm10.fasta.amb
├── mm10.fasta.ann
├── mm10.fasta.bwt
├── mm10.fasta.pac
└── mm10.fastq.sa

建立好參考基因組之后，就可以進(jìn)行比對了。不同的比對算法，命令不同。

1. BWA-backtrack 算法

對應(yīng)的子命令為aln/samse/sample
單端數(shù)據(jù)用法如下：

bwa aln ref.fa reads.fq > aln_sa.sai
bwa samse ref.fa aln_sa.sai reads.fq > aln-se.sam

雙端數(shù)據(jù)用法如下：

bwa aln ref.fa read1.fq > aln1_sa.sai
bwa aln ref.fa read2.fq > aln2_sa.sai
bwa sampe ref.fa aln_sa1.sai aln_sa2.sai read1.fq read2.fq > aln-pe.sam

2. BWA-SW 算法

對應(yīng)的子命令為bwasw, 基本用法如下

bwa bwasw ref.fa reads.fq > aln-se.sam
bwa bwasw ref.fa read1.fq read2.fq > aln-pe.sam

3. BWA-MEM` 算法

對應(yīng)的子命令為mem, 基本用法如下

bwa mem ref.fa reads.fq > aln-se.sam
bwa mem ref.fa read1.fq read2.fq > aln-pe.sam

對于超長讀長的reads,比如PacBio和Nanopore測序儀產(chǎn)生的reads, 用法如下

bwa mem -x pacbio ref.fa reads.fq > aln.sam
bwa mem -x ont2d ref.fa reads.fq > aln.sam

上述的代碼中, ref.fa指的是參考基因組索引的名字，對于上面提供的小鼠的例子而言，參考基因組索引的名字為mm10.fasta, 注意是不包含后綴的。

默認(rèn)用法都非常簡單，有時還需要對參數(shù)進(jìn)行調(diào)整，以mem子命令為例，常用的參數(shù)包括以下幾個:

-t指定線程數(shù)，默認(rèn)為1，增加線程數(shù)，會減少運行時間；-p 忽略第二個輸入序列，默認(rèn)情況下，輸入一個序列文件認(rèn)為是單端測序，輸入兩個序列文件則是雙端測序，加上這個參數(shù)后，會忽略第二個輸入序列文件，把第一個文件當(dāng)做單端測序的數(shù)據(jù)進(jìn)行比對;-Y參數(shù)對數(shù)據(jù)進(jìn)行soft clipping, 當(dāng)錯配或者gap數(shù)過多比對不上時，會對序列進(jìn)行切除，這里的切除并只是在比對時去掉這部分序列，最終輸出結(jié)果中序列還是存在的，所以稱為soft clipping。

看完上述內(nèi)容，你們掌握bwa軟件用法是怎樣的的方法了嗎？如果還想學(xué)到更多技能或想了解更多相關(guān)內(nèi)容，歡迎關(guān)注億速云行業(yè)資訊頻道，感謝各位的閱讀！