bowtie2用法

 bowtie2的功能：短序列的比对

用法：bowtie2 [options]* -x <bt2-idx> {-1 <m1> -2 <m2> | -U <r>} [-S <sam>]

$bowtie2 -x ./index/GCF reads.fa -S ./bowtie2
#-x的写法：index文件夹下，bt2的文件名是GCF。
#-S的目录必须先建立，不会自动建立。

　　

-x <bt2-idx>：参考基因组的索引路径

{-1 <m1> -2 <m2> | -U <r>}：比对的reads路径。-1，-2是针对PE reads；-U是针对SE reads（unpaired）

[-S <sam>]：输出文件路径

 

常用options：

-q：输入文件为.fastq格式，默认选择。

-f：输入文件为.fa/.mfa文件

--end-to-end       entire read must align; no clipping (on) #比对是将整个read和参考序列进行比对. 该模式--ma的值为0，与--local 是非彼即此关系，默认选择。

 --local            local alignment; ends might be soft clipped (off) #该模式下对read进行局部比对, 即read两端的一些碱基不比对，使比对得分满足要求. 该模式下–ma默认为2. 与--end-to-end 是与非彼即此关系，默认缺省。

 

参考文章：

http://bio-infor.blog.163.com/blog/static/1072798152014397554959/