SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
GATK "MAPQ should be 0 for unmapped read" complaint for mapped read efoss Bioinformatics 12 10-18-2012 02:19 AM
Average Read Coverage for 454 paired end read data lisa1102 Core Facilities 8 10-18-2011 09:40 AM
Puzzling result from Illumina 150bp PE reads chkuo De novo discovery 8 06-06-2011 01:07 AM
Assembly of 150bp paired-end reads anli Illumina/Solexa 1 05-18-2011 04:01 AM

Reply
 
Thread Tools
Old 06-01-2011, 02:03 AM   #1
NGS_user
Junior Member
 
Location: Europe

Join Date: Nov 2010
Posts: 9
Default 150bp read?

Hi Folks,
I am in a situation where I need long reads but with as few sequencing errors as possible. I currently have 75bp reads but am looking to generate longer ones. Would you recommend 100bp or 150bp? Ideally the 150 bp reads would avoid assembling errors I am having but a high number of errors in the reads could also cause all sorts of bother! Any recommendations?
NGS_user is offline   Reply With Quote
Old 06-01-2011, 03:11 AM   #2
tonybolger
Senior Member
 
Location: berlin

Join Date: Feb 2010
Posts: 156
Default

Quote:
Originally Posted by NGS_user View Post
Hi Folks,
I am in a situation where I need long reads but with as few sequencing errors as possible. I currently have 75bp reads but am looking to generate longer ones. Would you recommend 100bp or 150bp? Ideally the 150 bp reads would avoid assembling errors I am having but a high number of errors in the reads could also cause all sorts of bother! Any recommendations?
Use 150bp if you have the choice.

It should be better out to 100bp, plus you should get decent quality (on some reads) out to 130-140 bases. Adaptive trimming is, as ever, a good idea to take out the poorer reads.
tonybolger is offline   Reply With Quote
Old 06-01-2011, 09:10 AM   #3
HESmith
Senior Member
 
Location: Bethesda MD

Join Date: Oct 2009
Posts: 509
Default

An alternative is the one used by the Broad Institute for de novo genome assembly. Construct libraries with 180bp inserts, performed paired-end 100bp sequencing, and use the 20bp overlap in the center to generate 180bp reads. In addition to increased length, there's the added advantage that the bases at the ends of the reads (which typically have the lowest q-scores) are sequenced twice. The critical parameter is a tight size distribution for your libraries so that the reads overlap.
HESmith is offline   Reply With Quote
Old 06-01-2011, 09:16 AM   #4
Andrew_Slatter
Member
 
Location: London

Join Date: Dec 2010
Posts: 12
Default

All depends on what quality you need at the end of the run. I've just been looking at the latest MySeq performance data and they report a mean qscore of >20 at 150 bases.
Andrew_Slatter is offline   Reply With Quote
Old 06-02-2011, 11:02 AM   #5
chkuo
Member
 
Location: Taiwan

Join Date: May 2010
Posts: 11
Default

We had a lane of 150bp PE reads for de novo assembly of a bacterial genome. Curiously, trimming the reads down to 100bp actually produced larger contigs. Not sure what's going on, I just posted this issue here:

http://seqanswers.com/forums/showthread.php?t=11784
chkuo is offline   Reply With Quote
Old 06-06-2011, 03:56 AM   #6
tonybolger
Senior Member
 
Location: berlin

Join Date: Feb 2010
Posts: 156
Default

Quote:
Originally Posted by HESmith View Post
An alternative is the one used by the Broad Institute for de novo genome assembly. Construct libraries with 180bp inserts, performed paired-end 100bp sequencing, and use the 20bp overlap in the center to generate 180bp reads. In addition to increased length, there's the added advantage that the bases at the ends of the reads (which typically have the lowest q-scores) are sequenced twice. The critical parameter is a tight size distribution for your libraries so that the reads overlap.
You can do the same trick with 150bp reads, using 250bp or so inserts. However, you have to make absolutely sure that the distribution of reads is tight, or you'll possibly end up collapsing a small repeat region.
tonybolger is offline   Reply With Quote
Old 06-08-2011, 08:32 AM   #7
kwaraska
Senior Member
 
Location: Boston,MA

Join Date: Nov 2008
Posts: 122
Default

it depends on what machine and chemistry version.

On the GA they can comnfortably go to 150 without too many errors.
On the current HiSeq chemistry(not the one released last week-v4) I wouldn't go beyond 100, maybe not even beyond 75
That has been our experience.
They claim the v4 chemitry will be accurate to 150 bases (each way) but I don't have first hand knowledge yet. We will not be switching until we use up our old stuff and we have 6-8 weeks worth.
kwaraska 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 06:53 PM.


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