Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
csfasta + qual to fastq moriah Bioinformatics 0 07-05-2012 04:11 AM
Convert csfasta to fastq yksikaksi Bioinformatics 2 10-30-2011 08:36 PM
Error Bowtie with fastq files (from Solid .csfasta & .qual) pepperoni Bioinformatics 0 10-12-2011 08:20 PM
fastq to csfasta and .qual samt SOLiD 15 10-29-2009 09:11 AM
using BWA to align SOLiD fastq files from 1000 Genomes tgenahmet Bioinformatics 1 10-15-2009 05:19 PM

Thread Tools
Old 10-23-2012, 08:32 AM   #1
Location: USA

Join Date: Apr 2010
Posts: 76
Default SOLiD csfasta to fastq using BFAST+BWA


I'm trying to convert my SOLiD reads to FASTQ using BFAST+BWA's utility. I have two sets of paired-end files
However, the utility only allows files labeled F3 and R3
Note: <in.title> is the string showed in the `# Title:' line of a
".csfasta" read file. Then <in.title>F3.csfasta is read sequence
file and <in.title>F3_QV.qual is the quality file. If
<in.title>R3.csfasta is present, this script assumes reads are
paired; otherwise reads will be regarded as single-end.

The read name will be <out.prefix>: panel_x_y/[12] with `1' for R3
tag and `2' for F3. Usually you may want to use short <out.prefix>
to save diskspace. Long <out.prefix> also causes troubles to maq.
The SOLiD documentation (see here), however, says that F5 is for paired end and R3 is for mate-pair. My reads are paired end and so they are labeled F5.

My question is how do you go about merging these two files? Is it as simple as renaming the F5-DNA files as R3? I also notice in the documentation that F5 and R3 are in opposite orientation. Is it then required to reverse-complement the F5 reads to make them R3? How is that done? Am I missing something?

If solid2fastq cannot accomplish this task, do you know of any other programs which could do this? Finally, I want to:
1) Trim adapters from each paired-end file individually using cutadapt (which takes a FASTQ file)
2) Map the filtered reads using BFAST (which also takes a FASTQ file)

flobpf is offline   Reply With Quote
Old 10-23-2012, 12:01 PM   #2
Location: SF Bay Area

Join Date: Feb 2012
Posts: 62

Cutadapt accepts csfasta and qual input. It will even output a fastq afterwards. should take it from there.

Last edited by jparsons; 10-23-2012 at 12:09 PM. Reason: add'l info
jparsons is offline   Reply With Quote

bfast, solid, solid paired end

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 11:23 AM.

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