使用SAMtools软件进行序列可视化展示

SAMtools序列可视化操作

1）准备 BAM 文件

# 使用参考序列构建索引，index为索引前缀

bowtie2-build reference.fa index

# 使用bowtie2进行比对

bowtie2 -x index -1 reads_1.fq -2reads_2.fq -p 32 -S example_pair.sam

# 使用samtools将sam文件转换为bam文件

samtools sort example_pair.sam >example_pair.bam

2）使用 SAMtools 对 BAM 文件进行排序和索引（这是可视化展示的前提）

# 排序 BAM 文件

samtools sort example_pair.bam -oexample_pair.sorted.bam

# 创建 BAM 索引

samtools index example_pair.sorted.bam

3）可视化

# 使用 Samtools tview查看 BAM 文件的比对情况

samtools tview example_pair.sorted.bamreference.fa

# 运行后会进入santools tview的交互界面：

#  按下“shift+?”即可显示帮助菜单栏，如下图所示：

① 按下 g ，则提示输入要到达基因组的某一个位点

② 使用H(左）J（上）K（下）L（右）移动显示界面。

③ Ctrl+H 向左移动1kb碱基； Ctrl+L 向右移动1kb碱基

④ 可以用颜色标注比对质量，碱基质量，核苷酸等。

     30～40的碱基质量或比对质量使用白色表示；

     20～30黄色；

     10～20绿色；

     0～10蓝色。

⑤ 使用点号’.'切换显示碱基和点号；

⑥ 使用r切换显示read name