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/

05-11 18:12