SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
bwa sampe -c option cband Bioinformatics 1 08-07-2012 03:50 AM
-A option in bwa chsfredrick Bioinformatics 0 07-31-2012 01:37 PM
BWA -n option seq_lover Bioinformatics 2 05-16-2012 05:47 AM
velveth_de running time and mem Ori Bioinformatics 10 11-03-2011 04:39 PM
BWA samse -n option kellyv Bioinformatics 0 01-19-2010 11:15 AM

Reply
 
Thread Tools
Old 05-02-2013, 07:13 AM   #1
dakin
Junior Member
 
Location: US

Join Date: Dec 2009
Posts: 3
Default bwa mem -R option

Trying to use bwa mem with the -R option to add read group headers to the alignment output in one step, but it keeps coming back with the error:

[E::bwa_set_rg] no ID at the read group line

Below was what I tried for paired-end reads and then for all single reads:

bwa mem -aM -R @RG\tID:Seq01p\tSM:Seq01\tPL:ILLUMINA\tPI:330 Refgen Seq01pair1.fastq Seq01pair2.fastq > Seq01pairs.sam

bwa mem -aM -R @RG\tID:Seq01a\tSM:Seq01\tPL:ILLUMINA\tPI:330 Refgen Seq01all.fastq > Seq01all.sam

If I take out the entire -R option, bwa mem appears to run just fine. I suppose I can always add the read group header on later with Picard or samtools, but it would be nice to add it now....

Wondering if there is something wrong with my read group ID line? Has anyone else successfully run bwa mem with the -R option and what did you use?
dakin is offline   Reply With Quote
Old 05-02-2013, 07:32 AM   #2
Richard Finney
Senior Member
 
Location: bethesda

Join Date: Feb 2009
Posts: 700
Default

I"m not sure but try quotes for your mem argument ...

bwa mem -aM -R "@RG\tID:Seq01p\tSM:Seq01\tPL:ILLUMINA\tPI:330" Refgen Seq01pair1.fastq Seq01pair2.fastq > Seq01pairs.sam
Richard Finney is offline   Reply With Quote
Old 05-02-2013, 07:42 AM   #3
dakin
Junior Member
 
Location: US

Join Date: Dec 2009
Posts: 3
Default

Thanks for the quick (and correct!) response - I tried it with the quotes, and it worked!
dakin is offline   Reply With Quote
Reply

Tags
-r option, bwa, mem, options, read group header

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 01:31 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO