Rather than a list, we could consider returning a SummarizedExperiment object directly from digestFastqs. These could then be combined across samples with the combineCols() function. We would need to provide a sample name to digestFastqs, to use as the column name in the assays and to name the parameters in the metadata (since these will be combined without name checking by combineCols()). We can have multiple assays (reads, maxNReads, UMIs if applicable). This would also mean that we could use the existing QC plot functions directly on each sample as well.