Question: @RG ( read group ) in neandertals
0
gravatar for khatami.mahshid
12 months ago by
khatami.mahshid20 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 • 283 views
ADD COMMENTlink modified 12 months ago by WouterDeCoster44k • written 12 months ago by khatami.mahshid20
3
gravatar for WouterDeCoster
12 months ago by
Belgium
WouterDeCoster44k 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 12 months ago by ATpoint39k • written 12 months ago by WouterDeCoster44k
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: 1182 users visited in the last hour