使用bowtie2和samblaster一步到位的干净比对
bowtie2
以前都是和samtools组合，如下：
bowtie2 -x $index -U $id | samtools sort -@ 4 -o $sample.bam -
运行速度很慢，现在有高效工具啦，比如sambamba主要有filter，merge,slice和duplicate等七个功能来处理sam/bam文件，几乎可以替代 samtools啦，不过，这里要着重介绍的是samblaster
samblaster
主要参数：
-i --input 输入sam文件（必须包含header且按reads id排序）
-o --output 输出sam文件
-d --discordantFile 输出discordant read pairs
-s --splitterFile 输出split reads
-u --unmappedFile 输出unmapped/clipped reads
其他参数:
-a --acceptDupMarks 不去重
-e --excludeDups 去掉discordant, splitter, and/or unmapped等重复（具体定义详见samblaster主页）
-r --removeDups 去掉重复(-e --excludeDups类似)
--addMateTags 添加MC and MQ tags
-M 与bwa mem -M 类似
命令组合
bowtie2 -x $index -U $id | samblaster -e -d $sample.disc.sam -s $sample.split.sam | samtools view -Sb - > $sample.clean.bam
当然，也可以使用bwa啦，如下：
#自动输出discordant read pairs和split read alignments:
bwa mem <idxbase> samp.r1.fq samp.r2.fq | samblaster -e -d samp.disc.sam -s samp.split.sam | samtools view -Sb - > samp.out.bam
#从bam文件中提取 split reads和discordants read pairs
samtools view -h samp.bam | samblaster -a -e -d samp.disc.sam -s samp.split.sam -o /dev/null