Reads per cell are downsampled and the sequencing saturation is computed. Sequencing saturation is defined as 1 - UMIs / total reads.

plot_saturation(molecules, bcs_to_use = NULL, proportions = seq(0, 1,
  by = 0.1), return_data = FALSE)

Arguments

molecules

path molecule.tsv.gz flatfile produced by scrunchy pipeline or molecule_df data.frame

bcs_to_use

character vector of barcode sequences to use for computing sequencing saturation. Defaults to use all barcodes.

proportions

vector of proportions between 0 and 1 to downsample reads (defaults to seq(0, 1, by = 0.1))

return_data

if TRUE return plotting data instead of plotting distribution