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 character vector of barcode sequences to use for computing sequencing saturation. Defaults to use all barcodes. vector of proportions between 0 and 1 to downsample reads (defaults to seq(0, 1, by = 0.1)) if TRUE return plotting data instead of plotting distribution