Parses and tidies DRAGEN coverage outputs (per-contig mean coverage, coverage
metrics, fine histograms, and coverage-report BEDs). Some files fan out into
multiple output tables: tidy_metricsmain splits the coverage metrics into
metricsmain/metricsbins/metricscumu, and tidy_reportbedmain splits the
coverage-report BED into reportbedmain/reportbedcumu. Every other table has
tidy_name == parser, so its name is unchanged. The per-method docs below cover
each table's specifics.
The coverage region (wgs / tmb /qc-coverage-region-{<region>}) and
phenotype (normal / tumor) are folded from the filename into the output
prefix by DragenTool's refine_files() hook, so one schema table serves
every region/phenotype variant.
Super classes
nemo::Tool -> tidydragen::DragenTool -> DragenCov
Public fields
flat_tidy_names(
logical(1))TRUE: fanned-out sub-tables are named<tool>_<tidy_name>directly (parser token dropped), so each sub-table'snameis its final output table. See the class description for the fan-out map.
Methods
Method new()
Create a new DragenCov object.
Usage
DragenCov$new(path = NULL, files_tbl = NULL)Arguments
path(
character(1))
Output directory of tool. Iffiles_tblis supplied, this is ignored.files_tbl(
tibble(n))
Tibble of files fromnemo::list_files_dir().
Method tidy_metricsmain()
Tidy *_coverage_metrics.csv into metricsmain (wide summary),
metricsbins (bucketed depth histogram), and metricscumu (cumulative).
Examples
cls <- DragenCov; tool <- "dragencov"
indir <- system.file("extdata", tool, package = "tidydragen")
odir <- tempdir()
obj <- cls$new(indir)
obj$run(output_dir = odir, format = "parquet", input_id = "run1")
(lf <- list.files(odir, pattern = "dragencov_.*parquet", full.names = FALSE))
#> [1] "sampleA_exon_dragencov_metricsbins.parquet"
#> [2] "sampleA_exon_dragencov_metricscumu.parquet"
#> [3] "sampleA_exon_dragencov_metricsmain.parquet"
#> [4] "sampleA_target_bed_dragencov_metricsbins.parquet"
#> [5] "sampleA_target_bed_dragencov_metricscumu.parquet"
#> [6] "sampleA_target_bed_dragencov_metricsmain.parquet"
#> [7] "sampleA_umccr_dragencov_readreportbed.parquet"
#> [8] "sampleA_umccr_dragencov_reportbedcumu.parquet"
#> [9] "sampleA_umccr_dragencov_reportbedmain.parquet"
#> [10] "sampleA_wgs_dragencov_contigmean.parquet"
#> [11] "sampleA_wgs_dragencov_finehist.parquet"
#> [12] "sampleA_wgs_dragencov_metricsbins.parquet"
#> [13] "sampleA_wgs_dragencov_metricscumu.parquet"
#> [14] "sampleA_wgs_dragencov_metricsmain.parquet"
#> [15] "sampleA_wgs_normal_dragencov_contigmean.parquet"
#> [16] "sampleA_wgs_tumor_dragencov_contigmean.parquet"
