Question: @RG ( read group ) in neandertals
0
gravatar for khatami.mahshid
4 weeks ago by
khatami.mahshid10 wrote:

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

snp rna-seq • 140 views
ADD COMMENTlink modified 4 weeks ago by WouterDeCoster41k • written 4 weeks ago by khatami.mahshid10
3
gravatar for WouterDeCoster
4 weeks ago by
Belgium
WouterDeCoster41k wrote:

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

ADD COMMENTlink modified 4 weeks ago by ATpoint24k • written 4 weeks ago by WouterDeCoster41k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2073 users visited in the last hour