View Single Post
Old 12-19-2011, 03:43 AM   #1
alexbmp
Member
 
Location: Seoul, South Korea

Join Date: Oct 2011
Posts: 30
Question Setting @RG ID, PL, PU, LB, SM etc.

Hi all

When running bwa, I want the read group information to be implanted so that I can use the information in picard and GATK later.

I wrote down what I understood. Please point out if you think I'm wrong about something

-I see that RG ID is a read discriminator; even when BAM files are merged thanks to RG ID the reads can be discriminated.

-RG PL seems to be the sequencing machine platform (e.g. illumina).

-RG LB seems to be referred by picard MarkDuplicates so PCR duplicates in each sequencing library can be removed.


Now, I mention here that I read SAM-format spec.s.
( http://samtools.sourceforge.net/SAM1.pdf )
However, I'm confused about 2 more things.

1. Why does RG PU exist in the first place? Currently, the only reason I put in PU is to avoid picard and GATK errors later.

2. What's the difference between RG ID and RG SM? Is there any sutble difference between the two?

Thanks for your replies in advance!


Have a great day
alexbmp is offline   Reply With Quote