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

Try -r in both samse and sampe...

ADD COMMENTlink written 4.8 years ago by Devon Ryan98k

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.8 years ago by Max140
0
gravatar for WouterDeCoster
4.8 years ago by
Belgium
WouterDeCoster45k 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.8 years ago by WouterDeCoster45k
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: 1646 users visited in the last hour
_