Unmapped reads from activated B cell and plasmablast datasets (healthy or mild COVID-19 (ref. 30 ), n = 6; Tdap, n = 7; Comirnaty, n = 15; Vaxzevria 57 , n = 8; MIS-C, n = 32; severe COVID-19 (ref. 30 ), n = 12; same cells sorted) tested against the EBV genome. EBV-specific UMI counts were compared (healthy or mild COVID-19, n = 12,716 cells, 0 EBV UMIs; Tdap, n = 32,677 cells, 0 EBV UMIs; Comirnaty, n = 11,600 cells, 0 EBV UMIs; Vaxzevria, n = 16,445 cells, 2 EBV UMIs; MIS-C, n = 16,849 cells, 59 EBV UMIs; severe COVID-19, n = 96,703 cells, 361 EBV UMIs).