码迷,mamicode.com
首页 > 其他好文 > 详细

比对工具之 BWA 使用方法

时间:2016-06-20 20:33:51      阅读:668      评论:0      收藏:0      [点我收藏+]

标签:

BWA算法简介:

  • BWA-bactrack
  • BWA-SW
  • BWA-MEM

BWA安装:

# installing BWA
tar jxf ~/software/bwa-0.7.10.tar.bz2 -C /opt/biosoft/
cd /opt/biosoft/bwa-0.7.10/
make
echo PATH=$PATH:/opt/biosoft/bwa-0.7.10 >> ~/.bashrc
source ~/.bashrc

BWA使用步骤:

  1. 使用BWA构建参考基因组的index数据库
  2. BWA-MEM比对(BWA-SW 和 BWA-backtrack用的少)

具体命令:

# 创建运行目录并进入
mkdir -p ~/train/05.reads_aligment/bwa
cd 05.reads_aligment/bwa/

# 为原始数据创建软链接,方便数据处理
ln -s ~/data/00.incipient_data/data_for_gene_prediction_and_RNA-seq/genome.fasta
ln -s ~/data/00.incipient_data/data_for_variants_calling/F2-19.?.fastq ./

# 创建参考基因组的索引
bwa index genome.fasta -p genome

# 比对创建sam文件
bwa mem -t 4 genome F2-19.1.fastq F2-19.2.fastq > F2-19.mem.sam

 

 

参考资料:

BWA官方网站

NGS生物信息分析 V4.2

比对工具之 BWA 使用方法

标签:

原文地址:http://www.cnblogs.com/leezx/p/5601525.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!