Mapping SOliD reads to a Newbler 454 alignment to correct errors

Bukowski · 03-09-2010, 03:20 AM

I have a 454 assembly from Newbler, which I know has frameshifts/unexpected stop codons in the middle of ORF's.

I would like to use a SOLiD run (which gives me 150x + coverage over the Newbler contigs) in order to correct for potential errors.

I've been using MAQ successfully to map the SOLiD reads to the Newbler contigs. However I can't see the correct - if there is one - method to extract a sequence that represents the 'overloading' of the SOLiD data. I'm using the 454 contigs as a scaffold, but I would like to extract the 'best possible' sequence over the scaffold from only the SOLiD reads.

Is there any way to do this in MAQ or in another package?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Using solexa to correct 454 homopolymer errors	coldturkey	Bioinformatics	21	01-26-2016 02:39 AM
454 reads correct with illumina	biocomfun	454 Pyrosequencing	6	02-12-2012 04:00 AM
How to find out if mapping is correct/not	mitochy	Bioinformatics	3	01-17-2012 02:09 AM
Newbler mapping: 50% chimeric reads?	yvan.wenger	RNA Sequencing	1	05-10-2011 07:05 AM
Alignment of ABI solid reads and 454 reads	baohua100	Bioinformatics	2	02-23-2009 05:58 PM

03-09-2010, 03:20 AM	#1
Bukowski Senior Member Location: UK Join Date: Jan 2010 Posts: 390	Mapping SOliD reads to a Newbler 454 alignment to correct errors I have a 454 assembly from Newbler, which I know has frameshifts/unexpected stop codons in the middle of ORF's. I would like to use a SOLiD run (which gives me 150x + coverage over the Newbler contigs) in order to correct for potential errors. I've been using MAQ successfully to map the SOLiD reads to the Newbler contigs. However I can't see the correct - if there is one - method to extract a sequence that represents the 'overloading' of the SOLiD data. I'm using the 454 contigs as a scaffold, but I would like to extract the 'best possible' sequence over the scaffold from only the SOLiD reads. Is there any way to do this in MAQ or in another package?