missing read groups for unaligned reads

[gfmgfm]

Hi everyone,

I am using bwa to align Illumina reads. I used the -r option in order to specify the read group.
In the output sam file I am getting read groups only for the aligned reads. Reads that are not mapped (with '*' in the chromosome and position) don't have a read group.
In the downstream analysis I use GATK, which gives errors because of this.

Any suggestions how to overcome this?

Thanks

[aleidenroth]

I have the same problem, I also use bwa -r for alignment.

When I re-run TableRecalibrate with the added parameter

--default_read_group name_of_my_read_group

it seems to work and ignores the unmapped reads. However, I'm not sure this will cause any trouble in downstream analysis?

[lh3]

you should join the bwa mailing list. already fixed a while ago.

https://github.com/lh3/bwa