This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Read groups information?

I am doing sequence analysis on the genomic data of E. Coli. My sequencing comes from Illumina sequencing. As I don't have the detailed information, I set the arguments of AddOrReplaceReadGroups method in Picard tools as follows:

RGID=SampleName RGLB=SampleName RGPU=illumina RGPL=illumina RGLB=SampleName RGSM=SampleName

Is it fine or should I look for the exact value for these fields?

sequencing sequence picard next-gen picard tools
• 1,847 views • link
written by Researcher • Researcher

There is no formal definition of what is a read group,

This was the first sentence of the tutorial that @Pierre had linked on the last question you had asked about this topic. Are you trying to use GATK for bacteria?

• 0 views • link
written by GenoMax GenoMax

Dear @genomax

yeah, I am using GATK, but, I skipped some steps like "Base quality score recalibration (BQSR)" for bacteria. So, now my question is for which steps of GATK, the read group information is needed?

• 0 views • link
written by Researcher • Researcher

Any answer from @genomax , @ATpoint or anyone else?

• 0 views • link
written by Researcher • Researcher

Any answer from @genomax , @ATpoint or anyone else?

• 0 views • link
written by Researcher • Researcher

I am not a regular GATK user but all GATK/Picard tools likely require RG.

That said snippy was a simpler tool recommended today for bacterial SNP calling.

• 0 views • link
written by GenoMax GenoMax

Do you plan to use any tools that rely on exact information for this like some of the Broad Institute/Picard/GATK tools family? if not then just leave it as is, I personally never care about the inconsistently-defined concept of RGs.

• 0 views • link
written by ATpoint ATpoint

0 answers

No answers yet.

Log in to answer this question.

More posts like this