Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • How to correctly choose Read Groups ?

    Hello !

    I need to add read groups to my BAM files before merging and variant calling with GATK.
    Do you have any advice ?

    I will use AddOrReplaceReadGroups.jar (picard) and plan to fill in the following tags : ID, library (LB), sequencing platform=Illumina (PL), and sample (SM).
    Am I missing another important tag ?

    Also, if I understand well, the ID should be unique to a sample, lane and sequencing experiment ?

    Thanks,

    Muriel

  • #2
    I use a different read-group for each pair of (paired-end) fastq, even if they belong to the same lane/lib/sample etc...: in the end I can find the FASTQ from a given SAM-Record. My scripts look like this:

    Code:
    (...)
    (bwa...) -r "@RG	ID:${pair.generateId}	LB:${pair.sample.lib}	SM:${pair.sample.name}	PL:ILLUMINA	PU:${pair.lane}" \
    (...)

    Comment


    • #3
      OK thanks !

      But since I already used bwa to generate the BAM file for each sample which corresponds to the alignement of forward, reverse and single reads, how can I do this ?
      I don't think I can separate those three groups in a BAM file can I ?
      Because I juste wanted to use AddOrReplaceReadGroups and not to re-use bwa (which would be much longer...)

      Next time, I will know that I should have specified the RG before but what about those already mapped reads ?

      Thanks !

      Comment


      • #4
        if you already have one BAM per sample you can just add one distinct ID for each file.

        Comment


        • #5
          OK, I will re-map then :/

          So, according to your method, you define pair.generateId, pair.sample.lib etc.. at the beginning of the script so that after you can use this line right ?
          Sorry, I am not only new in NGS, I am also new in bash scripting

          Muriel

          Comment


          • #6
            define pair.generateId, pair.sample.lib etc.. at the beginning of the script so that after you can use this line right ?
            yes.

            I use https://github.com/lindenb/jvarkit/wiki/Illuminadir to generate a JSON description of my FASTQs and I then use https://github.com/lindenb/jsvelocity to generate my scripts.

            Comment


            • #7
              OK ! Thanks a lot for your help !!

              Comment


              • #8
                Hey I come back to you since I had a look on bwa.
                I seems that naming RG with bwa is during the process sampe and samse.
                But during sampe you give both forward and reverse fastq files at the same time so how can you give different read names to F and R ?

                Comment


                • #9
                  the option '-r' is set during bwa sampe: you can give only one group for both direction. You can later find the orientation using the sam flag.

                  Comment


                  • #10
                    Aaaah ok !
                    Thanks !

                    Comment

                    Latest Articles

                    Collapse

                    • seqadmin
                      Strategies for Sequencing Challenging Samples
                      by seqadmin


                      Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                      03-22-2024, 06:39 AM
                    • seqadmin
                      Techniques and Challenges in Conservation Genomics
                      by seqadmin



                      The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                      Avian Conservation
                      Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                      03-08-2024, 10:41 AM

                    ad_right_rmr

                    Collapse

                    News

                    Collapse

                    Topics Statistics Last Post
                    Started by seqadmin, 03-27-2024, 06:37 PM
                    0 responses
                    12 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 03-27-2024, 06:07 PM
                    0 responses
                    11 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 03-22-2024, 10:03 AM
                    0 responses
                    53 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 03-21-2024, 07:32 AM
                    0 responses
                    68 views
                    0 likes
                    Last Post seqadmin  
                    Working...
                    X