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)
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
if TRUE return plotting data instead of plotting distribution