名称 bwa – Burrows-Wheeler Alignment Tool 内容摘要描述命令行与选项SAM 比对格式短序列比对注意事项 比对精确性 估计插入大小分布 内存需求 速度Bwa-0.6中的改变其他作者引用与授权历史 摘要 b w a i n d e x r e f . f ab w a m e m r e f . f a r e a d s . f q > a l n - s e . s a mb w a m e m r e…
一.bwa比对软件的使用 1.对参考基因组构建索引 bwa index -a bwtsw hg19.fa # -a 参数:is[默认] or bwtsw,即bwa构建索引的两种算法,两种算法都是基于BWT的(BWT search while the CIGAR string by Smith-Waterman alignment.).-a bwtsw对于短的参考序列是不工作的,必须要大于等于10Mb:-a is 不适用于大的参考序列,必须要小于等于2G: output:hg19.fa.am…
一 建立索引 比对之前,需要对fasta文件构建FM-index索引:bwa index -a bwtsw hg19.fasta 生成 hg19.fasta.amb.hg19.fasta.ann.hg19.fasta.bwt.hg19.fasta.pac.hg19.fasta.sa四个文件. 起可参数 -a 有两种构建index算法: bwtsw 大的基因组数据,必须大于10MB,比如人的全基因组. is 默认的算法,速度较快,需要较大的内存,不能构建大于2GB的数据库 -P str 输出数…
bwa的使用需要两中输入文件: Reference genome data(fasta格式 .fa, .fasta, .fna) Short reads data (fastaq格式 .fastaq, .fq)step 1: 建立 Index根据reference genome data(e.g. reference.fa) 建立 Index File bwa index -a bwtsw reference.fabwa index 指令更多的用法及 options,通过以下的…
好吧,这是本周(2016.10.21-28)的学习任务之一:安装bowtie2并学习其使用方法&参数设置 所以,啃文档咯,官方文档Version 2.2.9 http://bowtie-bio.sourceforge.net/bowtie2/manual.shtml 以下是我的整理.我不生产文档,我只是文档的搬运工么么哒- Bowtie2适合将长度50-1000bp的reads比对到长的参考序列上.Bowtie 2 indexes the genome with an FM Index (bas…