@RG ( read group ) in neandertals

I have some bam files from Neandertals I downloaded from UCSC. I want to first merge them and then convert to VCF file. but when I search for @RG in bam files to see if they are from the same library, it shows they don't have @RG. I use this command to find @RG in bam files:

samtools view -H SLFeld1.hg19.bam | grep "^@RG" > RG_SLFeld

why these bam files don't have @RG?

thanks a lot

rna-seq snp

• 1,482 views • link

updated 6.9 years ago by WouterDeCoster

• written 6.9 years ago by khatami.mahshid •

1 answer

Because it is not required. You can set or change a read group using samtools addreplacerg, or with Picard and probably others too.

• 0 views • link

updated 6.9 years ago by ATpoint

• written 6.9 years ago by WouterDeCoster

More posts like this

Samtools merge bam issue with header tags
written 3.2 years ago by sam •

Hi, I have a few bam files after mapping using Bowtie2. I want to extract the multi-mapped reads from these bam files. For doing so …
How to split multi run bam files and add read groups
written 5.8 years ago by MAPK

I have a few hundred WGS samples that came without read groups. I need to run GATK on these samples, but before running GATK, I …
Merge bam files
written 6.9 years ago by khatami.mahshid •

I already ask a question about this, but I didn't get my answer complete. https://www.biostars.org/p/399693/ I have 6 bam file that one of them is …
Neandertals VCF file
written 6.9 years ago by khatami.mahshid •

Hi all, I have Bam files of Neandertals from UCSC database. I tried to convert these bam files to VCF files with `samtools mpileup` and …
Merge bam files and get VCF file
written 6.9 years ago by khatami.mahshid •

I have 2 sorted Bam files and I want to merge them and get a VCF file from them using samtools pileup. first i merge …
Chimeric read insert length in SAM $9
written 7.2 years ago by Yijun Tian •

Hello everyone, I am recently doing a script to dissect chimeric alignment from BAM file, here is my code: #!/bin/bash samtools=`which samtools` bn=$(basename $1 .bam) …
Replace @RG group by the one of another BAM
written 7.4 years ago by da.de •

Hi, I want to replace the @RG info from one bam file (output.bam) by the @RG info from another (input.bam) creating a new.bam I tried …
BWA (mem, sampe, or aln) used in bam file
written 8.2 years ago by rborges •

I'm trying to figure out which version/type of bwa was used in three aligned bam files (sampe, mem or aln). Is there a way of …
Sam file header and truncated file
written 9.8 years ago by TriS

I posted this on the samtools groups with no luck, so hopefully here i get better luck I'm changing part of the header in a …
samtools merge respecting original read groups
written 10.2 years ago by abascalfederico

Hi, I'm struggling to merge many bam files respecting the original read groups with samtools. Some of the input files have the same read group, …