When run in multi-sample mode, the pipeline assigns sample-id = patient-id. I think plot_vaf is looking for the actual sample id within the VCF and fails. Either skip plot_vaf for multiiple or update to plot VAF for each sample (always only mutect2 for multiple samples).