SEQanswers (
-   Bioinformatics (
-   -   missing read groups for unaligned reads (

gfmgfm 01-19-2011 01:08 PM

missing read groups for unaligned reads
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?


aleidenroth 01-20-2011 05:10 AM

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 01-20-2011 06:49 AM

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

All times are GMT -8. The time now is 06:53 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.