SEQanswers

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 07:46 AM
Assembling De Novo 454 Transcriptome Contigs and Singletons with Illumina Short Reads Vickenstein Bioinformatics 7 03-05-2011 12:43 AM
error correction for RNA-seq reads fhb RNA Sequencing 6 11-18-2010 01:58 PM
PubMed: Correction of sequencing errors in a mixed set of reads. Newsbot! Literature Watch 0 04-10-2010 02:00 AM
Alignment of ABI solid reads and 454 reads baohua100 Bioinformatics 2 02-23-2009 04:58 PM

Reply
 
Thread Tools
Old 11-24-2010, 03:47 AM   #1
yvan.wenger
Member
 
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
yvan.wenger is offline   Reply With Quote
Old 11-24-2010, 04:08 AM   #2
colindaven
Senior Member
 
Location: Germany

Join Date: Oct 2008
Posts: 414
Default

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

http://seqanswers.com/forums/showthread.php?t=3635
colindaven is offline   Reply With Quote
Old 11-26-2010, 03:48 AM   #3
Brugger
Member
 
Location: Cambridge, UK

Join Date: Mar 2010
Posts: 21
Default

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, 04:17 AM   #4
jts
Member
 
Location: Cambridge

Join Date: Feb 2009
Posts: 22
Default

Have you seen this paper?

http://bioinformatics.oxfordjournals.../1284.abstract
jts is offline   Reply With Quote
Reply

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 12:35 PM.


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