SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
BWA mapping for SOLiD single-end reads PedroA Bioinformatics 7 12-13-2012 11:07 AM
Mapping SOLiD colorspace paired end reads NestorNotabilis SOLiD 10 12-12-2012 06:14 PM
SOLiD paired-end mapping yaira04 SOLiD 6 12-12-2012 04:41 PM
BWA paired end mapping quality pparg Bioinformatics 9 11-14-2011 06:51 PM
BWA align SOLiD PE data gives poor mapping & 0.5% properly paired alig Bioinformatics 3 07-08-2011 08:44 AM

Reply
 
Thread Tools
Old 09-16-2013, 03:12 AM   #1
Kath
Member
 
Location: Oxford

Join Date: Nov 2010
Posts: 24
Default SOLiD paired end data and BWA - cleaning/mapping

Hello everyone,

I notice there has been the odd post about this general topic before but I am still confused!

I have SOLiD single and paired end data which I would like to take through the following steps:
1) Quality assessment - Fastqc
2) Data cleaning - Quality trimming/removing adapters if necessary
3) Mapping using BWA

At the moment I am converting the csfasta and qual files in to double-encoded (fakespace) fastq files for BWA. I reverse the F5 reads so BWA treats them correctly. However, can the same fastq files be used with Fastqc? Also, can CutAdapt and fastq_quality_trimmer handle this format? If anyone else out there still uses SOLiD data and can give me some details of their workflows I would be very grateful.

Best wishes,
Kath
Kath is offline   Reply With Quote
Old 09-18-2013, 01:42 PM   #2
BrianJames
Member
 
Location: San Diego

Join Date: May 2010
Posts: 13
Default

I used to work with SOLiD data, and I have to ask why you are going to all of this trouble. In my experience, the best way to map SOLiD data is to stay in colorspace, and use Bfast, or mapreads (part of lifescope/bioscope).


Brian
BrianJames is offline   Reply With Quote
Old 09-19-2013, 01:25 AM   #3
Kath
Member
 
Location: Oxford

Join Date: Nov 2010
Posts: 24
Default

Thanks for your reply. I have Lifescope mapped BAM files but I'm trying to dedup using Picard and am getting some very high duplication rates. I've been told that Picard seems to overestimate duplication rates in Lifescope mapped data and that BWA provides better compatibility with Picard. I did look into using Bfast but previous benchmarking studies suggested BWA would be better.

Anyway, it looks like there is a script - solid-trimmer.py - that will do what I want.

https://github.com/brentp/bio-playgr...lid-trimmer.py
Kath is offline   Reply With Quote
Old 09-19-2013, 08:53 AM   #4
BrianJames
Member
 
Location: San Diego

Join Date: May 2010
Posts: 13
Default

OK. Are you working with DNA-seq, or RNA-seq? If it is DNA seq, I would be worried about not working in colorspace, if it is RNA-seq, removing duplicates from SOLiD short read data may not be a good idea. Just my experience.

Brian
BrianJames is offline   Reply With Quote
Old 09-20-2013, 04:24 AM   #5
colindaven
Senior Member
 
Location: Germany

Join Date: Oct 2008
Posts: 415
Default

We tested all SOLiD mappers very thoroughly, and NovoalignCS (commercial but not too expensive) won hands-down. I would really recommend it.

BWA dropped color space support a long time ago and did not perform very well in terms of # mapped reads. Bfast wasn't much fun either.
colindaven is offline   Reply With Quote
Reply

Tags
bwa, cutadapt, fastqc, fastq_quality_filter, solid 5500

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 03:25 PM.


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