SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
picard add read groups HGENETIC Bioinformatics 0 01-25-2012 04:19 AM
best-practices for assigning read-groups for snp-caller brentp Bioinformatics 2 12-20-2011 08:34 AM
basic question about read groups efoss Bioinformatics 2 10-19-2011 05:32 PM
unaligned reads of tophat ae_ucla RNA Sequencing 1 04-07-2011 11:06 AM
Reorganize unaligned read output lew13 Bioinformatics 2 03-14-2010 07:56 PM

Reply
 
Thread Tools
Old 01-19-2011, 01:08 PM   #1
gfmgfm
Member
 
Location: il

Join Date: Jun 2010
Posts: 64
Default 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
gfmgfm is offline   Reply With Quote
Old 01-20-2011, 05:10 AM   #2
aleidenroth
Junior Member
 
Location: Nottingham, UK

Join Date: Nov 2009
Posts: 8
Default

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?
aleidenroth is offline   Reply With Quote
Old 01-20-2011, 06:49 AM   #3
lh3
Senior Member
 
Location: Boston

Join Date: Feb 2008
Posts: 693
Default

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

https://github.com/lh3/bwa
lh3 is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:59 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO