用法如下:

bedtools coverage -a diff_peak.bed -b clean.bam -counts >.peak.count.xls # 计算区间的平均测序深度(常用)

samtools bedcov gene.bed sample.bam # 计算bed区间所有位点测序深度的加和

samtools depth test.bed sample.bam # 计算单个位点或给定bed文件区间内所有单位点的深度

bedtools coverage diff_peak.bed -b clean.bam > samp.depth.txt
samp.depth.txt 统计结果的输出文件,格式为7列:染色体、区间起始位点、区间结束位点、该区间内的reads数、该区间内的碱基数、区间大小、该区间的平均覆盖度

Logo

有“AI”的1024 = 2048,欢迎大家加入2048 AI社区

更多推荐