搜索
您的当前位置:首页正文

从bam中截取指定位置的 reads 或深度

来源:欧得旅游网

 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

因篇幅问题不能全部显示,请点此查看更多更全内容

Top