Question: bwa samse: invalid option 'R'
0
gravatar for Max
4.1 years ago by
Max140
Max140 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.3k views
ADD COMMENTlink modified 4.1 years ago by WouterDeCoster44k • written 4.1 years ago by Max140
0
gravatar for Devon Ryan
4.1 years ago by
Devon Ryan95k
Freiburg, Germany
Devon Ryan95k wrote:

Try -r in both samse and sampe...

ADD COMMENTlink written 4.1 years ago by Devon Ryan95k

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 4.1 years ago by Max140
0
gravatar for WouterDeCoster
4.1 years ago by
Belgium
WouterDeCoster44k 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 4.1 years 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: 1113 users visited in the last hour