Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Changing header of BAM

    Hello,

    I apologize if my question is redundant....

    I created bam files from fastq files by ELAND_standalone and tried to feed them to GATK. Then I got the following error message.
    ##### ERROR MESSAGE: SAM/BAM file reanalysis.bam is malformed: SAM file doesn't have any read groups defined in the header. The GATK no longer supports SAM files without read groups
    ##### ERROR -

    How can I change the header? Samtools?

    Thank you very much.

    Hiroki

  • #2
    you could easily change the header using samtools, but that won't solve the problem as the readgroup is necessary in each alignment section as well.
    You could have a look in the AddOrReplaceReadGroup Utility of Picard:


    I don't know if that changes the header as well but if not have a look here:

    Comment


    • #3
      Thanks ulz_peter!

      I was told to use Picard from another guy.
      Now I'm trying it.

      Thank you so much anyways!

      Comment


      • #4
        I've run into similar issues using GATK for bam files, and Picard. The option that always has helped me when I knew hte BAM file was "ok,"

        java [args to jvm] GATK.jar [GAKT opts] VALIDATION_STRINGENCY=SILENT

        Comment


        • #5
          Hi, I have the same problem as Hiroki.

          I made a sam file using command:
          Code:
          bwa samse ref.fa my.sai my.fastq > my.sam
          then, I have creaded .bam file:
          Code:
          samtools view -bS my.sam > my.bam
          next, I have sorted my.bam file:
          Code:
          samtools sort my.bam my_sorted
          Using GATK, I need to put a .bam file as an input.
          According to the instructions from http://www.broadinstitute.org/gsa/wi...s_for_the_GATK I have used Picard's ReordereSam, to reorder reads:
          Code:
          java -jar ReordereSam.jar I=/path/my_sorted.bam O=/path/my_reordered.bam R=/path/ref.fa
          Trying i.e.
          Code:
          java -jar GenomeAnalysisTK.jar -T DepthOfCoverage -R /paht/ref.fa -I /path/aln_reordered.bam
          an error is returned:
          HTML Code:
          ##### ERROR MESSAGE: SAM/BAM file /path/my_reordered.bam is malformed: SAM file doesn't have any read groups defined in the header.  The GATK no longer supports SAM files without read groups
          I have tried Picards AddOrReplaceReadGroups.jar and ReplaceSamHeader.jar but it didn't help me.

          Any suggestions how to bulid proper .bam file for GATK?
          thx

          Comment


          • #6
            hi the damian,
            when creating your sam file, you have to add read group header.
            this should work:

            bwa samse -r @RG\tID:IDa\tSM:SM\tPL:Illumina ref.fa my.sai my.fastq > my.sam
            bests

            colin

            Comment


            • #7
              Yes, I've already realised how to do this. Another option is to use Picard's AddOrReplaceReadGropus.jar.
              Code:
              java -jar AddOrReplaceReadGroups I=my.bam O=myGr.bam LB=whatever PL=illumina PU=whatever SM=whatever
              Of course prior using AddOrReplaceReadGroups your .sam file needs to pass Picard's validation (ValidateSamFile.jar)

              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, Yesterday, 06:37 PM
              0 responses
              11 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, Yesterday, 06:07 PM
              0 responses
              10 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-22-2024, 10:03 AM
              0 responses
              51 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