Genome Coverage Calculator

Estimate the average sequencing depth for your targeted experiment based on a BED file’s target region size.

Average Sequencing Coverage

0X

Coverage vs. Read Count

This chart illustrates how the average coverage (X-depth) increases as the total number of sequencing reads increases, keeping read length and target size constant.

---

Coverage Projection Table

Projected average coverage at different read counts.

Total Reads	Projected Coverage

What is Genome Coverage?

In genomics, genome coverage (or sequencing depth) refers to the average number of times a specific nucleotide in a genome or target region is sequenced. For instance, a coverage of 50X means that, on average, each base in the region of interest was read 50 times. When you ‘calculate genome coverage using a BED file’, you are specifically calculating this metric for the targeted regions defined in that file, which is common in exome sequencing or with targeted gene panels.

This is a critical quality control metric in Next-Generation Sequencing (NGS). Higher coverage increases the statistical confidence in identifying genetic variants like SNPs and indels, reducing the likelihood that a detected variant is a sequencing error. The required coverage depends on the application; for example, detecting rare somatic mutations in cancer tissue requires much higher coverage than calling common germline variants. This genome coverage calculator is designed to help you plan your sequencing experiments to achieve the desired depth.

The Formula to Calculate Genome Coverage

The calculation is based on the Lander-Waterman equation, which provides a straightforward way to estimate the average coverage across a defined area. The formula is:

Coverage (C) = (Number of Reads (N) × Average Read Length (L)) / Total Size of Target Region (T)

This formula helps you understand the relationship between the three main parameters of a sequencing experiment. To achieve higher coverage, you can either increase the number of reads (more sequencing) or use longer reads, assuming the target size remains constant.

Formula Variables

Variable	Meaning	Unit	Typical Range
N	Total number of sequencing reads	(Unitless count)	Millions to billions
L	Average length of a single read	base pairs (bp)	50 – 300 bp (for short-read)
T	Total size of the genomic regions	bp, kb, Mb	1 kb – 3,000 Mb

Practical Examples

Understanding how to calculate genome coverage is easier with realistic scenarios.

Example 1: Human Exome Sequencing

• Inputs:
  • Number of Reads (N): 80,000,000
  • Average Read Length (L): 150 bp
  • Total Size of Target Region (T): 45 Mb
• Calculation:
  • Total Sequenced Bases = 80,000,000 × 150 = 12,000,000,000 bases
  • Target Size = 45,000,000 bases
  • Coverage = 12,000,000,000 / 45,000,000 = 266.7X

Example 2: Targeted Gene Panel

• Inputs:
  • Number of Reads (N): 5,000,000
  • Average Read Length (L): 100 bp
  • Total Size of Target Region (T): 500 kb
• Calculation:
  • Total Sequenced Bases = 5,000,000 × 100 = 500,000,000 bases
  • Target Size = 500,000 bases
  • Coverage = 500,000,000 / 500,000 = 1000X

How to Use This Genome Coverage Calculator

1. Enter Total Reads: Input the total number of reads your sequencing run produced. This is often provided in the run summary from the sequencing facility.
2. Enter Read Length: Provide the average length of your reads in base pairs (e.g., 100 for 100bp reads).
3. Enter Target Region Size: Input the cumulative size of all target regions from your BED file. You can select the most convenient unit (bp, kb, or Mb). To get this value, you might need a simple script. For example, using a command-line tool like `awk` on your BED file: awk -F'\t' '{{sum += $3 - $2}} END {{print sum}}' yourfile.bed.
4. Review Results: The calculator will instantly provide the average sequencing coverage (X-depth). It also shows intermediate values and projects coverage at different sequencing depths, helping you plan future experiments. For more analysis options, consider our BED file format guide.

Key Factors That Affect Genome Coverage

• Sequencing Depth (Total Reads): The most direct factor. More reads lead to higher coverage.
• Read Length: Longer reads contribute more bases per read, increasing coverage for the same number of reads.
• Target Region Size: A larger target region requires more sequencing to achieve the same level of coverage.
• Library Quality: Poor library preparation can lead to a high number of PCR duplicates, which are often removed during analysis, effectively reducing the number of useful reads and thus lowering coverage. Our guide on PCR duplicate removal provides more context.
• Enrichment Efficiency: In targeted sequencing, the efficiency of capturing the desired regions is crucial. Low efficiency means more reads will map “off-target,” wasting sequencing data and reducing on-target coverage.
• Sequencing Uniformity: Reads are rarely distributed perfectly evenly. GC-rich or repetitive regions are often harder to sequence, leading to lower coverage in those areas and “valleys” in your coverage plot. Good read quality control (QC) is essential.

Frequently Asked Questions (FAQ)

1. What is considered “good” sequencing coverage?

It’s highly dependent on the application. For germline variant calling in diploid genomes, 30X-50X is a common standard. For detecting low-frequency somatic mutations in tumors, coverage of 500X or even higher is often required. For a deeper dive, see our comparison of somatic vs. germline variant calling.

2. How do I calculate the ‘Total Size of Target Region’ from my BED file?

You need to sum the lengths of all intervals. Each line in a standard BED file has a start (column 2) and end (column 3) coordinate. The length of that interval is `end – start`. You can use a simple script or command-line tool to sum these lengths for all lines in the file. A common `awk` command is: awk -F'\t' '{{sum += $3 - $2}} END {{print sum}}' your_panel.bed

3. Does this calculator account for off-target reads?

No, this calculator provides a theoretical average coverage assuming all reads map to the target region. In reality, some percentage of reads will be off-target. Your actual on-target coverage will be lower, depending on the capture efficiency of your experiment.

4. Why is my actual coverage from `samtools depth` different from the calculated value?

This calculator estimates the *average* coverage. Actual coverage varies position by position due to sequencing biases (e.g., GC content) and random sampling. Tools like `samtools depth` report the per-base depth, and the average of those values should be close to the estimate from this genome coverage calculator, but will be affected by real-world factors like mapping quality and duplicate removal.

5. What is the difference between sequencing depth and coverage?

The terms are often used interchangeably. “Depth” and “coverage” both refer to the number of times a base is sequenced (e.g., 50X depth or 50X coverage). Sometimes, “coverage” can also refer to the percentage of the target region sequenced at least once (e.g., 99% coverage breadth at 20X depth), but “depth” almost always refers to the “X” value.

6. Can I use this for Whole-Genome Sequencing (WGS)?

Yes. For WGS, the “target region” is the entire genome. Simply enter the haploid genome size of your organism in the “Total Size of Target Region” field. For humans, this is approximately 3,100 Mb. Check a genome size database for other species.

7. Does read pairing (paired-end sequencing) affect this calculation?

Not directly in this formula. You should use the total number of reads (e.g., if you have 50 million pairs, that’s 100 million reads in total) and the length of a single read (e.g., 150bp for a 2x150bp run). Do not double the read length.

8. How do I choose the right unit for Target Region Size?

Choose the unit that is most convenient for you. The calculator handles the conversion automatically. Genetic panels are often measured in kilobases (kb), exomes in megabases (Mb), and whole genomes in megabases (Mb) or gigabases (Gb).

Related Tools and Internal Resources

Expand your bioinformatics toolkit with these related resources:

• Genome Size Database: Find the haploid genome size for thousands of species to use in WGS calculations.
• BED File Format Guide: A deep dive into the BED format and how to manipulate it.
• Introduction to NGS Analysis: A beginner’s guide to the entire Next-Generation Sequencing workflow.
• Read Quality Control (QC) Metrics: Understand the key metrics for assessing the quality of your sequencing data.
• Somatic vs. Germline Variant Calling: Learn about the different requirements and challenges for each type of analysis.
• PCR Duplicate Removal Tool: An example tool demonstrating how duplicate reads can affect your analysis.