Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
short reads from amplicon for 454 using Titanium chemistry pseudorabies 454 Pyrosequencing 9 08-15-2011 08:46 AM
Assembling De Novo 454 Transcriptome Contigs and Singletons with Illumina Short Reads Vickenstein Bioinformatics 7 03-05-2011 01:43 AM
error correction for RNA-seq reads fhb RNA Sequencing 6 11-18-2010 02:58 PM
PubMed: Correction of sequencing errors in a mixed set of reads. Newsbot! Literature Watch 0 04-10-2010 03:00 AM
Alignment of ABI solid reads and 454 reads baohua100 Bioinformatics 2 02-23-2009 05:58 PM

Thread Tools
Old 11-24-2010, 04:47 AM   #1
Location: Switzerland

Join Date: Aug 2009
Posts: 30
Default 454 Reads correction with Short reads ?

Hello everybody,

I can see that my 1Mio 454 reads are full of insertions/deletions in homopolymer regions, and of course this disrupts my ORFs in later stages (annotation step). Fortunately, I have ca 60mio Illumina reads that I can use to correct these errors. Both dataset are cDNA.

Has anyone a preferred method to correct short sequencing error indels? If tractable, I would prefer correcting my reads rather than correcting my draft assembly obtained from those reads. I was going to use ssaha (or segemehl, qpalma, other?) to map the illumina reads on my 454 reads, and then extract a consensus by parsing the pileup file generated with samtools

I have seen a few published, more advanced, methods available but I*do not now if one of them performs better (e.g. iCorn). Do you?

Thanks a lot!

yvan.wenger is offline   Reply With Quote
Old 11-24-2010, 05:08 AM   #2
Senior Member
Location: Germany

Join Date: Oct 2008
Posts: 415

I don't know if any programs will correct your reads!

See here for another effort .. I have got it to work but am still interpreting results
colindaven is offline   Reply With Quote
Old 11-26-2010, 04:48 AM   #3
Location: Cambridge, UK

Join Date: Mar 2010
Posts: 21

I did this once. My approach was to map the reads to the contigs (BWA -e5), realigned the indels using GATK and polished the contigs based on the the indel report. I did this multiple times as some corrections made it possible to correct other neighbouring regions.

I dont have the polishing script any more, but I recall it as fairly simple. It should be possible to polish the reads instead, with some perl regex magic.
Brugger is offline   Reply With Quote
Old 11-26-2010, 05:17 AM   #4
Location: Cambridge

Join Date: Feb 2009
Posts: 22

Have you seen this paper?

jts is offline   Reply With Quote

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 09:00 AM.

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