SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   missing read groups for unaligned reads (http://seqanswers.com/forums/showthread.php?t=8973)

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?

Thanks

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.

https://github.com/lh3/bwa


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.