ReadMGATK {Signac}

Read MGATK output

Description

Read output files from MGATK (https://github.com/caleblareau/mgatk).

Usage

ReadMGATK(dir, verbose = TRUE)

Arguments

dir

Path to directory containing MGATK output files

verbose

Display messages

Value

Returns a list containing a sparse matrix (counts) and two dataframes (depth and refallele).

The sparse matrix contains read counts for each base at each position and strand.

The depth dataframe contains the total depth for each cell. The refallele dataframe contains the reference genome allele at each position.

Examples

## Not run: 
data.dir <- system.file("extdata", "test_mgatk", package="Signac")
mgatk <- ReadMGATK(dir = data.dir)

## End(Not run)