Tools to calculate coverage of a bam file window by window · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

Tools to calculate coverage of a bam file window by window

Hi there,

I am looking for a tool that scans an input bam file window by window (by window I mean consecutive counting bins of constant size, e.g. 10kb) and outputs the total or average coverage in percent per window. What would be the best way to do this? I've been looking at gatk CollectReadCounts, samtools coverage and DeepTools bamcoverage.

Thanks!

samtools deeptools bam gatk

• 1,993 views • link

updated 3.1 years ago by Ram

• written 3.7 years ago by marion.prieler •

1 answer

mosdepth (https://github.com/brentp/mosdepth) is a fast tool to do this; and you can specify -b 10000 to run on 10kb windows, instead of (like GATK or FeatureCounts etc) needig to create and pass in a binned .bed file.

• 0 views • link

written 3.7 years ago by LChart

thanks so much, I'll have a look!

• 0 views • link

written 3.7 years ago by marion.prieler •

Log in to answer this question.

More posts like this

ChIP-seq visualization: Is it valid to do a coverage normalization in addition to applying a spike-…
written 2.9 years ago by kalavattam •

Let's say a ChIP-seq sample is scaled in the [following way](https://www.biostars.org/p/9572653/#9572655): 1. Calculate the percent spike-in reads in the immunoprecipitate (IP) 2. Calculate the percent …
Spike-in normalization via deepTools multiBamSummary --scalingFactors: What's an appropriate --binS…
written 3.0 years ago by kalavattam •

To normalize ChIP-seq data with a separate-species spike-in, I'm following the advice described in this [GitHub issue](https://github.com/deeptools/deepTools/issues/1073): 1. Align to spiked-in species. 2. Compute scaling …
Disable automatic merging of consecutive windows in DeepTools bamcoverage
written 3.6 years ago by marion.prieler •

Hi all, I'm using deeptools bamcoverage to output the coverage of my bam files window by window. For my analysis the windows have to be …
Get average base quality of bases in a read
written 3.8 years ago by pablosolar.r •

Hi all! Given a BAM, I want to plot the insert size of a read pair vs average base quality of bases in a read …
BWA alignment generating errors
written 5.8 years ago by MAPK

I am trying to align interleaved fastq files (multiple runs) for SRA data I have downloaded. These fastq files are splitted by Read groups. I …
Difference in .bed target region base count and mapped base count for different tools
written 6.5 years ago by gavehan •

Hello! I was initially doing performance comparison for tools to obtain sample wide **average depth of coverage**. I compared several tools listed below *Qualimap, GATK …
Different bigwig file sizes between BedgraphToBigwig and bamCoverage
written 8.0 years ago by m98

I am trying to convert BAM files to Bigwig format. I was originally using a combination of genomeCoverageBed and BedgraphToBigWig to go from BAM to …
Where, When and How to use the "breadth of coverage" values reported per bin in BEDtools coverageBE…
written 8.3 years ago by Zee_S

Hello Everyone, I use BEDtools extensively for my data analyses and I have a confusion about the BEDtools coverage tool. here is a summary: Aim …
Detecting expression in intergenic regions
written 8.7 years ago by A. Domingues

**Background** We observed some RNA-seq coverage in regions outside annotated genes - let's called them intergenic. This expression appears to be more pronounced, or unique, …
Calculate percentage of bases covered less than nX in targeted sequencing experiment
written 11.6 years ago by Ram

<p>Hi,</p> <p>I'm looking to calculate some coverage statistics on the alignment in my targeted sequencing experiment. I have a BAM file and a BED file …