计算测序深度和覆盖度

# 概念：

深度（depth）与覆盖度（coverage）

假想实验

对长100bp的目标区域进行捕获测序：采用单端测序，每个read长5bp；总共得到了200个reads；把所有的reads比对到目标区域后，100bp的目标区域中有98bp的位置至少有1个read覆盖到，换言之，剩余的2bp没有1个read覆盖。

深度（depth）

200 x 5 / 100 = 10

我们说这此测序的深度为10X。

覆盖度（coverage）

98 / 100 × 100% = 98%

我们说这次测序的覆盖度为98%

ref： 深度（depth）与覆盖度（coverage） | Public Library of Bioinformatics

# 使用genomecov：

ref： genomecov — bedtools 2.27.0 documentation

samtools view -b xxx.bam | genomeCoverageBed -ibam stdin -g hg18.genome -bga > cov.bedgraph

hg19:hg19.chrom.sizes

mm10:mm10.chrom.sizes

或者samtools查看bam：

samtools view -H xxx.bam  

-bg Reporting genome coverage in BEDGRAPH format

-bga Reporting genome coverage for all positions in BEDGRAPH format

另一个ref：测序数据基本信息统计 | reads,coverage,depth -

# 使用samtools的depth：

ref：[samtools]depth命令简介 - CSDN博客

 samtools depth xxx.bam -a > xxx.bam.depth 

结果： 一共得到3列以指标分隔符分隔的数据，第一列为染色体名称，第二列为位点，第三列为覆盖深度。