bwa samse: invalid option 'R'
2
0
Entering edit mode
8.0 years ago
Max ▴ 150

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 • 2.1k views
ADD COMMENT
0
Entering edit mode
8.0 years ago

Try -r in both samse and sampe...

ADD COMMENT
0
Entering edit mode

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 REPLY
0
Entering edit mode
8.0 years ago

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 COMMENT

Login before adding your answer.

Traffic: 2380 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6