SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Reply
 
Thread Tools
Old 12-08-2011, 09:41 AM   #1
cnicolet
Member
 
Location: Los Angeles

Join Date: Dec 2008
Posts: 35
Default Casava1.8.2

Hi:
How are people feeling about CASAVA 1.8.2? We are just starting to use it and so far are underwhelmed. The new formats are cumbersome and the way it delivers multiple files is complicating everything we do downstream. We also benchmarked the new ELAND aligner and it performed exactly the same as the old aligner (which in our hands only delivers about 80% of the mapped reads as bwa), despite claims of superiority. Maybe we're not using it correctly? Any tips?
cnicolet is offline   Reply With Quote
Old 12-08-2011, 11:49 PM   #2
simonandrews
Simon Andrews
 
Location: Babraham Inst, Cambridge, UK

Join Date: May 2009
Posts: 870
Default

There's supposed to be an option in 1.8.2 to generate only a single fastq.gz file which might make things a bit easier.

I'm afraid we gave up trying to work with Eland since 1.8 due to the extra complexity of the run commands. We were perfectly happy with the performance of the old Eland, but have now moved to bowtie for everything (which is also a pain as you can't run paired end gzipped data through it without uncompressing first).
simonandrews is offline   Reply With Quote
Old 12-09-2011, 03:21 PM   #3
cnicolet
Member
 
Location: Los Angeles

Join Date: Dec 2008
Posts: 35
Default

Thanks Simon! We did find some info on that using the --fastq-cluster-count parameter. But it's a little confusing, the manual says "Specify 0 to ensure creation of a single FASTQ file" then later it says "If you need to generate one unique fastq gzipped file for use in a third-party tool, you can set the --fastq-cluster-count option to -1" But I agree with you about Eland, I think our days of using that are over too.
cnicolet is offline   Reply With Quote
Old 12-09-2011, 05:32 PM   #4
BAMseek
Senior Member
 
Location: St. Louis, MO, USA

Join Date: Apr 2011
Posts: 124
Default

Quote:
Originally Posted by simonandrews View Post
... but have now moved to bowtie for everything (which is also a pain as you can't run paired end gzipped data through it without uncompressing first).
Hi Simon,

If the pain is about unzipping the files into temporary storage first, I think there might be some ways to uncompress and feed into Bowtie on the fly.

One way that I think works would be with named pipes like this:
Code:
rm -f pipe1
rm -f pipe2
mkfifo pipe1
mkfifo pipe2
gunzip -c reads/e_coli_1000_1.fq.gz > pipe1 & \
gunzip -c reads/e_coli_1000_2.fq.gz > pipe2 & \
./bowtie e_coli -1 pipe1 -2 pipe2
You could also try uncompressing and merging the paired files on-the-fly into the special tab-delimited format accepted by Bowtie, where the columns are

Quote:
read_name[tab]sequence_1[tab]quality_1[tab]sequence_2[tab]quality_2
and pipe that into Bowtie.

best,
Justin
BAMseek is offline   Reply With Quote
Old 12-12-2011, 06:10 AM   #5
simonandrews
Simon Andrews
 
Location: Babraham Inst, Cambridge, UK

Join Date: May 2009
Posts: 870
Default

Thanks for the suggestions Justin. I hadn't realised there was a paired end format already - I might just make up a wrapper for that.

I'm sorely tempted to just patch in gzip support to bowtie, but don't want to have the hassle of maintaining a fork in perpetuity.
simonandrews is offline   Reply With Quote
Old 01-20-2012, 07:25 AM   #6
earonesty
Member
 
Location: United States of America

Join Date: Mar 2011
Posts: 52
Default

Bowtie 0.12.7 patch for gzip support:

http://code.google.com/p/ea-utils/so...tie-gzip.patch
earonesty is offline   Reply With Quote
Reply

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 12:27 AM.


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