SEQanswers

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



Similar Threads
Thread Thread Starter Forum Replies Last Post
MiSeq gDNA reads still fail "Kmer content" and "per base seq content" after trimming" ysnapus Illumina/Solexa 4 11-12-2014 08:25 AM
AVA-difference betweeen "computation" reads and "global align" reads CCBIO 454 Pyrosequencing 0 01-17-2014 04:13 AM
DEXSeq error in estimateDispersions: match.arg(start.method, c("log(y)", "mean")) fpadilla Bioinformatics 14 07-03-2013 03:11 PM
The position file formats ".clocs" and "_pos.txt"? Ist there any difference? elgor Illumina/Solexa 0 06-27-2011 08:55 AM
SEQanswers second "publication": "How to map billions of short reads onto genomes" ECO Literature Watch 0 06-30-2009 12:49 AM

Reply
 
Thread Tools
Old 11-14-2014, 01:39 AM   #1
ymc
Senior Member
 
Location: Hong Kong

Join Date: Mar 2010
Posts: 498
Default What does it mean by "150/350 PE reads" from MiSeq?

http://www.ncbi.nlm.nih.gov/sra/SRX573008[accn]

I downloaded this SRA data. I noticed that the R1 reads are 150bp but the R2 reads are 350bp. How could that be? Aren't they supposed to be the same length???
ymc is offline   Reply With Quote
Old 11-14-2014, 01:58 AM   #2
nucacidhunter
Jafar Jabbari
 
Location: Melbourne

Join Date: Jan 2013
Posts: 1,238
Default

Unequal length on MiSeq is permitted.
nucacidhunter is offline   Reply With Quote
Old 11-14-2014, 09:57 AM   #3
ymc
Senior Member
 
Location: Hong Kong

Join Date: Mar 2010
Posts: 498
Default

Quote:
Originally Posted by nucacidhunter View Post
Unequal length on MiSeq is permitted.
So this data was generated by a 2x250 run?

Why do people want to do unequal length?

Thanks in advance!
ymc is offline   Reply With Quote
Old 11-14-2014, 10:43 AM   #4
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

All else equal... longer reads are more useful. 350bp reads can span repeats that 250bp reads can't, for example.

However, all else is NOT equal, and I seriously doubt the quality at the end of those 350bp reads is very good.
Brian Bushnell is offline   Reply With Quote
Old 11-14-2014, 10:50 AM   #5
ymc
Senior Member
 
Location: Hong Kong

Join Date: Mar 2010
Posts: 498
Default

Quote:
Originally Posted by Brian Bushnell View Post
All else equal... longer reads are more useful. 350bp reads can span repeats that 250bp reads can't, for example.

However, all else is NOT equal, and I seriously doubt the quality at the end of those 350bp reads is very good.
Thanks a lot for your reply. Are there any limitations to adjusting the lengths? Can you do 1/499 or even 0/500???
ymc is offline   Reply With Quote
Old 11-14-2014, 11:23 AM   #6
HESmith
Senior Member
 
Location: Bethesda MD

Join Date: Oct 2009
Posts: 509
Default

Well, you need 4-5 cycles in the first read to call clusters, and at least 25 cycles for base-calling. Plus, as Brian noted, the quality drops precipitously with the longer reads. But I think that the instrument could be programmed for 25/475.
HESmith is offline   Reply With Quote
Old 11-14-2014, 11:38 AM   #7
ymc
Senior Member
 
Location: Hong Kong

Join Date: Mar 2010
Posts: 498
Default

Quote:
Originally Posted by HESmith View Post
Well, you need 4-5 cycles in the first read to call clusters, and at least 25 cycles for base-calling. Plus, as Brian noted, the quality drops precipitously with the longer reads. But I think that the instrument could be programmed for 25/475.
Thanks for your reply. It is good to know that this is possible anyway. 25/575 for the 2x300 kit can come in handy sometimes.
ymc is offline   Reply With Quote
Old 11-14-2014, 12:58 PM   #8
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,087
Default

Or you could do a 600 cycle SE run.
GenoMax is online now   Reply With Quote
Old 11-14-2014, 05:15 PM   #9
ymc
Senior Member
 
Location: Hong Kong

Join Date: Mar 2010
Posts: 498
Default

Quote:
Originally Posted by GenoMax View Post
Or you could do a 600 cycle SE run.
So you can turn PE into SE in V3 chemistry but not V2?
ymc is offline   Reply With Quote
Old 11-14-2014, 05:59 PM   #10
nucacidhunter
Jafar Jabbari
 
Location: Melbourne

Join Date: Jan 2013
Posts: 1,238
Default

It can be done with v2 too and possibly on HiSeq. In MiSeq after cycle 200 the quality drops and is in lowest value at 300 cycle. Increasing cycle beyond that results in low quality so it may not be useful depending on intended use of data.
nucacidhunter is offline   Reply With Quote
Old 11-25-2014, 03:10 AM   #11
matth431
Member
 
Location: UK

Join Date: Jan 2012
Posts: 36
Default

Quote:
Originally Posted by ymc View Post
So this data was generated by a 2x250 run?

Why do people want to do unequal length?

Thanks in advance!
One reason we've done this is where there's a high proportion of short fragments in the library - we've had a user make Nextera libraries from very low amounts of viral DNA which tend to over-tagment, so lots of ~200bp inserts. If the bulk of the reads sequence past the end of the molecules, the machine loses focus and the run will abort prior to the index reads.

So we often do 100/400 reads to ensure we get index data for all clusters, then hopefully get as much of a reverse read as possible for the small percentage of larger fragments in the library - if the run aborts on read 2, we can at least recover the data, but if it aborts on read 1, the data is useless as it can't be demultiplexed.
matth431 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 10:31 AM.


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