[samtools] 文本查看语法,浏览SNP/INDEL位点

  santools可以作为文本查看工具,查看比对结果文件,下面做一简单介绍:

  1. 通过BWA比对获取sam比对文件,也可以将fastq文件转化为bam/sam文件;

  2. 转换sam文件为bam文件,samtools view -bS seq.sam > seq.bam

  3. 对bam文件进行排序,samtools sort seq.bam -o seq.sorted.bam

  4. 对bam文件进行index,samtools index seq.sorted.bam

  5. 查看比对结果文件,samtools tview seq.sorted.bam ref.fasta

注:

  所有的帮帮助说明都可以在官方网站查阅:

  http://www.htslib.org/doc/samtools.html

  https://github.com/samtools/hts-specs

posted on 2016-11-12 18:46  萧飞IDO  阅读(638)  评论(0编辑  收藏  举报

导航