SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Help about setting up TopHat lewewoo Illumina/Solexa 2 11-14-2011 06:52 PM
bowtie parameter setting nunu_ping RNA Sequencing 0 09-28-2010 02:59 PM
Setting up and efficiency ColNYC General 0 09-14-2010 09:26 AM
Mismatch setting in BWA genelab Bioinformatics 4 08-04-2010 05:02 AM
tophat parameter setting hollandorange Bioinformatics 0 05-14-2010 04:12 PM

Reply
 
Thread Tools
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
Old 12-19-2011, 04:06 AM   #2
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,543
Default

SM is sample, and one sample might be sequenced twice, for instance with 454 and Illumina, which would be two read groups with different platform. At least, that is my understanding.
maubp is offline   Reply With Quote
Old 12-20-2011, 11:20 PM   #3
alexbmp
Member
 
Location: Seoul, South Korea

Join Date: Oct 2011
Posts: 30
Default

Thank you for your reply. So although I don't need the information right now, it can be important in some other pipeline I might take, right?

Quote:
Originally Posted by maubp View Post
SM is sample, and one sample might be sequenced twice, for instance with 454 and Illumina, which would be two read groups with different platform. At least, that is my understanding.
alexbmp 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:57 PM.


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