This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Platform Unit For Sam Read Group

I'm wondering exactly what the "Platform unit" means in the read group header in the SAM format. This is what I found in the specification:

Platform unit (e.g. flowcell-barcode.lane for Illumina or slide for SOLiD). Unique identi er.

Since I have Illumina data, does this mean that I should use something like this: "FC706VJ.1" - assuming that the flowcell barcode is "FC706VJ" and the lane is 1. (The example names are from the FASTQ enty at wikipedia)? And If this is the case, does any one know if it is possible to extract the flowcell name from the report.xml generated Casava? I have a hunch that this is the same as the RunFolder attribute - but I might be wrong.

• 6,181 views • link
written by Johan Johan

1 answer

Ok. I will answer my own question with what I have figured out this far. The flowcell-barcode is a unique identifier (as described in the documentation quoted above) of the flowcell. In the report.xml file it should be included as an attribute called "FlowCellId" in the MetaData tag. However, this does not always seem to be included. This is then coupled with the lane(s) that the sample was run on to form the platform unit.

• 185 views • link
written by Johan Johan

Log in to answer this question.

More posts like this
  • Bcl2fastq conversion
    written by abiuma • abiuma
    Hi, I am trying to run bcl2fastq bash file in linux and I keep getting this error. How could I fix this. "Failed to parse …
  • Changing @RG header information in *.bam files using picard
    written by kirannbishwa01 kirannbishwa01
    I am working with multiple bam files which is aligned to my reference (using Bowtie2). I am planning to move forward to variant calling using …
  • Demultiplexing fastq files. Undocumented Illumina feature.
    written by BioApps BioApps
    I am building a full blown fastq/sff editor: https://www.biostars.org/p/98851 While building the demultiplexing module I ran into a FastQ with strange Illumina encoding. First line …
  • Bedops Bam2Bed Error With Gatk Recalibrated Bam-File
    written by Irsan Irsan
    <p>Whey I try to convert one of my <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a>-files with bam2bed from the bedops suite (version 2.3.0) I get an error saying that my …
  • Advice On Adding Readgroups
    written by Jamie • Jamie
    <p>Hi there,</p> <p>Ok so I have some Illumina whole genome bams. They have been aligned by illumina using Casava but we wanted to re-run the …
  • Tophat: Proper-Paired Reads
    written by ai • ai
    <p>Hi Everyone,</p> <p>I tried Tophat v 2.0.6 and GSNAP v 2013-02-05 to align Illumina paired-end reads. However, the numbers of proper-paired reads were very different. …
  • What Are The (Non-Standard) Bd And Bi Tags In Samformat?
    written by Johan Johan
    <p>I'm running the <a href='https://github.com/broadgsa/gatk/blob/master/public/scala/qscript/org/broadinstitute/sting/queue/qscripts/DataProcessingPipeline.scala'>dataprocessing pipeline qscript</a> provided by the GATK people to do dataprocessing according to their <a href='http://gatkforums.broadinstitute.org/discussion/1186/best-practice-variant-detection-with-the-gatk-v4-for-release-2-0'>best practices</a>. Looking at the files …
  • [Casava] Error: Barcode Cttgta For Lane 1 Has Length 6: Expected Barcode Lenth Is 7
    written by Pierre Lindenbaum Pierre Lindenbaum
    <p>That's the first time I'm using <strong>Casava</strong> to generate my <strong>fastq</strong> files.</p> <p>My command line was:</p> <pre><code>/path/to/CASAVA/bin/configureBclToFastq.pl \ --input-dir ../120414_SNL149_0002_AC0KR2ACXX/Data/Intensities/BaseCalls \ --output-dir 20120720.out \ --sample-sheet …
  • How To Set Rg Header For Tablerecalibration?
    written by PeterPan • PeterPan
    <p>hi, everyone~ I am using GATK recently. And I also use Picard's AddOrReplaceReadGroups to add RG header. Also I checked <a href='http://picard.sourceforge.net/command-line-overview.shtml#AddOrReplaceReadGroups'>http://picard.sourceforge.net/command-line-overview.shtml#AddOrReplaceReadGroups</a>, I still don't …
  • Is Casava Multithreaded?
    written by Misha • Misha
    <p>I built a copy of CASAVA from source. </p> <p>I am using CASAVA to perform a BCL to FASTQ conversion on a single lane (my …