5转载 bwa的使用方法

Posted 风中之铃

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了5转载 bwa的使用方法相关的知识,希望对你有一定的参考价值。

http://www.plob.org/?p=25

bwa的使用需要两中输入文件:

  1. Reference genome data(fasta格式 .fa, .fasta, .fna)
  2. Short reads data (fastaq格式 .fastaq, .fq)

step 1: 建立 Index
根据reference genome data(e.g. reference.fa) 建立 Index File

bwa index -a bwtsw reference.fa

bwa index 指令更多的用法及 options,通过以下的命令来查看

bwa index

step 2: 寻找 SA coordinates
如果是pair-end 数据(leftRead.fastq和rightRead.fastq)两个文件分别处理

bwa aln reference.fa leftRead.fastq > leftRead.sai
bwa aln reference.fa rightRead.fastq > rightRead.sai
bwa aln reference.fa singleRead.fastq > singleRead.sai

如果希望多线程运行,在其中加入 -t这个参数,另外-f这个参数可以指定结果输出文件,如:

bwa aln -c -t 3 -f leftreads.sai reference.fa leftreads.fastq

step 3:转换SA coordinates输出为sam
如果是pair-end数据

bwa sampe -f pair-end.sam reference.fa leftRead.sai rightRead.sai leftRead.fastq rightread.fastq

如果是single reads数据

bwa samse -f single.sam reference.fa single.sai single.fastq

值此Reads的mapping工作已经完成,关于bwa更详细的用法以及输出结果SAM文件的格式说明,可以参考官方文档.

以上是关于5转载 bwa的使用方法的主要内容,如果未能解决你的问题,请参考以下文章

10.12 bwa使用 安装文件路径与使用 sh权限

GATK--使用转载

二代下机文件在Integrative Genomics Viewer中实现reads可视化的流程

比对软件BWA使用

比对工具之 BWA 使用方法

精心收集的 48 个 JavaScript 代码片段,仅需 30 秒就可理解!(转载)