Skip to content

Bam header RG conflict #12

@arnavaz

Description

@arnavaz

Describe the bug

The header in the bam files that are created by bwa aligner have two instances of @rg written in the bam header and GATK tools have issues with that. Example:
@rg ID:1 SM:SE17-1161_N_cfDNA-010-R1.fastq.gz PL:Illumina
@pg ID:bwa PN:bwa VN:0.7.15-r1140 CL:/mnt/work1/software/bwa/0.7.15/bwa mem -M -t4 -R @rg ID:1 SM:SE17-1161_N_cfDNA-010-R1.fastq.gz

To Reproduce

Steps to reproduce the behavior:

  1. Go to '...'
  2. Use commands '....'
  3. Use inputs available at '...'
  4. See error

Expected behavior

A clear and concise description of what you expected to happen.

Screenshots

If applicable, add screenshots to help explain your problem.

System (please complete the following information):

  • OS: [e.g. iOS, Linux, etc]
  • Software versions

Additional context

Add any other context about the problem here.

Relevant inputs

Please make sure to provide relevant input files

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions