Question: bwa samse: invalid option 'R'
0
gravatar for Max
3.2 years ago by
Max130
Max130 wrote:

When creating .bam files from single-end read fastq's, I ran into difficulties with mutect because the .bam files didn't contain read groups. I attempted to recreate the .bams by rerunning bwa samse with the -R <readgroup> option, but I get an error stating that -R is an invalid option with samse (in contrast, it works with sampe). Is there an option and argument that I should use with bwa sampe to generate read groups?

bwa • 1.1k views
ADD COMMENTlink modified 3.2 years ago by WouterDeCoster40k • written 3.2 years ago by Max130
0
gravatar for Devon Ryan
3.2 years ago by
Devon Ryan91k
Freiburg, Germany
Devon Ryan91k wrote:

Try -r in both samse and sampe...

ADD COMMENTlink written 3.2 years ago by Devon Ryan91k

Right, I think I was confusing the -R for read groups in bwa mem (-r is something else) with -r for reads in bwa sampe

ADD REPLYlink written 3.2 years ago by Max130
0
gravatar for WouterDeCoster
3.2 years ago by
Belgium
WouterDeCoster40k wrote:

Alternatively, you can use Picard AddOrReplaceReadGroups

java -jar AddOrReplaceReadGroups.jar I= infile.bam O= outfile.bam SORT_ORDER= coordinate RGID=foo RGLB=bar RGPL=illumina RGSM=Sample RGPU=NONE CREATE_INDEX=True
ADD COMMENTlink written 3.2 years ago by WouterDeCoster40k
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: 1615 users visited in the last hour