samtools view -h recall.bam chr1:1116029-1116298 > get.sam # 获取指定位置reads
samtools view -bS get.sam | samtools sort - -o get.bam # sam--->bam,排序
samtools index get.sam # 建索引
# 深度统计
samtools depth -r chr1:1116029-1116298 -d 10000000 -a recall.bam > results.cov
因篇幅问题不能全部显示,请点此查看更多更全内容