Similar to #105, since rg_id_in_bam could be the same from different input BAMs, in bwa-mem-alignerhere using rg_id_in_bam to find corresponding read group in metadata JSON is no longer reliable.
The solution is similar to how it's done in payload-gen-dna-seq-qc at here