SEQanswers

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

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

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

Hi,

I'm trying to convert my SOLiD reads to FASTQ using BFAST+BWA's solid2fastq.pl utility. I have two sets of paired-end files
Quote:
*_F3.csfasta
*_F3_QV.qual
*_F5-DNA.csfasta
*_F5-DNA.QV.qual
However, the solid2fastq.pl utility only allows files labeled F3 and R3
Quote:
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)

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

Join Date: Feb 2012
Posts: 62
Default

Cutadapt accepts csfasta and qual input. It will even output a fastq afterwards.
http://seqanswers.com/forums/showthread.php?t=7100 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
Reply

Tags
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