1）samtools简介
--------------------------------------------------------------------------
背景：前面我们讲过sam/bam格式，sam文件虽然是可读的文本文件形式，但是通常是非常大，因此一般会对其压缩来节省磁盘空间，且对于很多软件来说，相比于对sam文件，对bam文件进行处理更加有效。SAMtools 是一款优秀的用以解析、处理sam/bam格式文件的一种软件包工具。其详细的文档可以在其官网里面找到。它主要包含以下功能：
1、格式转换(convert from other alignment formats)
2、排序 ，合并(sort and merge alignments）
3、去除PCR重复( remove PCR duplicate)
4、call SNP和indels( call SNPs and short indel variants)
5、数据提取(data extraction)
6、可视化（show alignments in a text-based viewer)
2）用法和命令
----------------------------------------------------------------------
像许多linux命令一样，SAMtool命令遵循流模型，允许将多个命令组合到数据处理管道中。标准流（stdin，stdout和stderr）， 发送到stdout的数据默认打印到屏幕上，但很容易使用普通的Unix重定向器（>和>>）重定向到另一个文件，或通过管道（|）重定向到另一个命令。
SAMtools 命令的基本用法是： samtools <command> [options]
2.1) view
view命令用来过滤SAM或BAM格式的数据。输入通常是指定为参数的sam或bam文件，但可以是来自任何其他命令的sam或bam数据。 可选功能包括查看原始文件内容、BAM和SAM格式之间进行转换、将数据子集提取到新文件中， 提取读取的顺序得以保留。
-- Viewing
flags explain BAM flags
tview text alignment viewer
view SAM<->BAM<->CRAM conversion
depad convert padded BAM to unpadded BAM
示例：
2.1.1) 将bam文件转化为sam文件,h 参数表示header
samtools view -h RNA-seq.bam >RNA-seq.sam
2.1.2) 将sam文件转化为bam文件
samtools view -bS RNA-seq.sam > sample.bam
2.1.3) 提取某个区间的比对结果（需要索引文件）
samtools view -hb RNA-seq.bam 'chr17:40465680-40465850' >out.bam
2.14）tview命令启动基于ascii的交互式查看器，可用于可视化读取与参考基因组的指定区域。
Usage: samtools tview <aln.bam> [ref.fasta]
2.15）提取比对上的区域
samtools view -h -F 4 -b RNA-seq.bam > aln_only_mapped.bam
2.2）index
index命令创建一个新的索引文件，生成的* .sam.sai或* .bam.bai文件，允许快速查找（已排序）SAM或BAM中的数据。
-- Indexing
dict create a sequence dictionary file
faidx index/extract FASTA
fqidx index/extract FASTQ
index index alignment
示例：
2.2.1）samtools faidx longreads.fa #对fasta建立索引，会生成.fai文件
2.2.2）samtools faidx longreads.fa r1 >r1.fas #建立索引后就可以提取指定的序列，这里提r1
2.2.3）samtools fqidx longreads.fq #对fastq建立索引
2.2.4）samtools index RNA-seq.bam #对比对文件建立索引
2.3)格式转换
-- File operations
collate shuffle and group alignments by name
cat concatenate BAMs
merge merge sorted alignments
mpileup multi-way pileup
sort sort alignment file
split splits a file by read group
quickcheck quickly check if SAM/BAM/CRAM file appears intact
fastq converts a BAM to a FASTQ
fasta converts a BAM to a FASTA
示例
2.3.1）samtools sort -n RNA-seq.bam >RNA-seq_sorted_out.bam #将bam文件进行排序
2.3.2）samtools fastq RNA-seq.bam >RNA-seq.fq #将bam 文件转化为fastq
2.3.3）samtools fasta RNA-seq.bam >RNA-seq.fa #将bam 文件转化为fasta
2.4）编辑
-- Editing
calmd recalculate MD/NM tags and '=' bases
fixmate fix mate information
reheader replace BAM header
targetcut cut fosmid regions (for fosmid pool only)
addreplacerg adds or replaces RG tags
markdup mark duplicates
示例：
2.4.1）samtools markdup RNA-seq.bam RNA-seq. markdup.bam #进行duplicates标记
2.6）统计
-- Statistics
bedcov read depth per BED region
depth compute the depth
flagstat simple stats
idxstats BAM index stats
phase phase heterozygotes
stats generate stats (former bamcheck)
示例：
2.6.1）samtools depth RNA-seq.bam >depth.out #统计depth信息
2.6.2）samtools flagstat RNA-seq.bam >flagstat.out # 统计大致信息
2.6.1）samtools idxstats RNA-seq.bam >out # 统计index信息
3）参考资源
----------------------------------------------------------------------
1）https://en.wikipedia.org/wiki/SAMtools
2）The Sequence Alignment/Map format and SAMtools
3）https://davetang.org/wiki/tiki-index.php?page=SAMTools
4）https://github.com/davetang/learning_bam_file