454: Extremely general question

Protaeus · 03-29-2011, 07:37 AM

All,

I have 454 and PacBio SRA data (that are cleaned of adaptors) that I'd like to align to a reference. I've tried bwasw, but when I look at the alignments in samtools tview, they are rather subpar as compared to data that has been reported by others on these sequences. Any suggestions for general approaches to aligning long reads? bwasw should work, correct? Better tools for dealing with these data?

ulz_peter · 03-29-2011, 09:17 PM

I guess it depends on how many reads you've got. For a moderate number of seuqences BLAT does a good job as well.
I used bwasw for 454 data as well and it worked out nicely.

RCJK · 03-29-2011, 10:22 PM

How about Newbler for the 454 reads? Not sure how it would deal w/ PacBio data, but if you use fasta + qual I think it could handle it (depending on read length).

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
General Question: Is it possible to generate mRNA from genomic DNA (in vitro)?	SeqUser2	RNA Sequencing	3	01-07-2013 05:09 AM
General question about TopHat	mhayes	Bioinformatics	1	02-03-2012 02:00 PM
General Question regarding Tophat and Bioscope	piRWLan	Bioinformatics	0	04-25-2011 06:03 AM
general genetics question	hoffmann	General	3	10-14-2010 01:38 AM
Discussion the general pipeline to analyze 454, Illumina/Solexa data	edge	General	1	09-09-2009 06:29 PM

03-29-2011, 07:37 AM	#1
Protaeus Member Location: 86001 Join Date: Aug 2010 Posts: 21	454: Extremely general question All, I have 454 and PacBio SRA data (that are cleaned of adaptors) that I'd like to align to a reference. I've tried bwasw, but when I look at the alignments in samtools tview, they are rather subpar as compared to data that has been reported by others on these sequences. Any suggestions for general approaches to aligning long reads? bwasw should work, correct? Better tools for dealing with these data?

03-29-2011, 09:17 PM	#2
ulz_peter Senior Member Location: Graz, Austria Join Date: Feb 2010 Posts: 219	I guess it depends on how many reads you've got. For a moderate number of seuqences BLAT does a good job as well. I used bwasw for 454 data as well and it worked out nicely.

03-29-2011, 10:22 PM	#3
RCJK Senior Member Location: Australia Join Date: May 2009 Posts: 155	How about Newbler for the 454 reads? Not sure how it would deal w/ PacBio data, but if you use fasta + qual I think it could handle it (depending on read length).