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