Compute exome coverage from a region file

This function will take a bed file, and return the sum of the lengths of unique regions

compute_exome_coverage(targets_bed_file, out_file = NULL)

Arguments

targets_bed_file

Path to a bed file with exome target regions

out_file

A file name to which the number of covered bases will be written, instead of returning the value

Value

An integer value of the sum of the length of the covered regions

Author

Mayank Tandon, Ashish Jain

Examples

library(MAFDash)
bedFile <- system.file("extdata", "test.bed", package = "MAFDash")
coverage<-compute_exome_coverage(bedFile)