Hello!
I'm using GATK to process some BAM files using their Best Practices. However, I forgot to add platform information under the @RG tag in the BAM header when I did the alignment (using BWA).
So, I've done dedup marking, indel realignment and the next step is base score recalibration. GATK won't do this since the platform information is missing from the BAM file. My first thought was to just use AddOrReplaceReadGroups from Picard, but as I understand it, this removes all other read groups in the BAM file and replace them with the one specified when running the program.
My question is: Can I just add platform information to the existing BAM file without removing the read groups in it (and replacing them with a single one)?
To clarify further, my header now looks something like this:
...
@RG ID:ABC0021 SM:SAMPLE123
@RG ID:ABC0026 SM:SAMPLE123
@RG ID:ABC0013 SM:SAMPLE123
...
And I want it to look something like this:
...
@RG ID:ABC0021 SM:SAMPLE123 PL:illumina PU:1 LB:Lib1
@RG ID:ABC0026 SM:SAMPLE123 PL:illumina PU:1 LB:Lib1
@RG ID:ABC0013 SM:SAMPLE123 PL:illumina PU:1 LB:Lib1
...
Also with the platform info added to the reads.
Is this possible to achieve using Picard?
it is much like
How can i edit some rows in .bam header file?
(closing)