SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
IGV .genome files: where is the sequence? pmiguel Bioinformatics 7 12-22-2011 07:03 PM
Pattern location, Please help me happy2000 Bioinformatics 2 05-16-2011 02:00 AM
Ensembl human genome naming pattern rboettcher General 3 03-03-2011 06:32 AM
qseq files versus sequence.txt files drio Illumina/Solexa 3 11-09-2009 10:02 AM
Can anyone make sense of the quality scores in the qseq.txt files? TylerBackman Bioinformatics 2 04-29-2009 10:23 AM

Reply
 
Thread Tools
Old 08-26-2010, 09:27 AM   #1
brentp
Member
 
Location: salt lake city, UT

Join Date: Apr 2010
Posts: 72
Default qseq files pattern of reads with no sequence

i dont believe this is a problem, just something i'm curious about:

normally we get our data from a sequencing center that returns a fastq file. this time, we have a bunch of qseq files like: s_3_1_0001_qseq.txt
in all of them, the first and last few hundred reads in the file have no sequence (just .'s for the length of the read). as one moves down from the top of the file (or up from the bottom), it starts that sequence appears on either end of the read, so it will look like:
Code:
G.TCTG..............................................................ACCCGGA.
(and the quality is all B's). going to the center of the file, the reads are all sequence.
i'm curious as to why this occurs at the ends of the file. is there some correspondence between position in the file and location on the cell.

thanks.
brentp is offline   Reply With Quote
Old 08-26-2010, 09:44 AM   #2
kmcarr
Senior Member
 
Location: USA, Midwest

Join Date: May 2008
Posts: 1,178
Default

Your guess is correct. Reads are sorted by their x-coordinate (but apparently there's no sorting by y-coordinate). This means that reads at the start and end of the qseq are from the extreme edges of the tile which have much lower image quality. What you are observing in these qseq files is perfectly normal and expected.

It is also likely that the fastq file they typically send you only contains filter passed reads whereas the qseq files contain all reads. Reads like the one shown would not end up in the filter passing fastq file.

Last edited by kmcarr; 08-26-2010 at 09:47 AM.
kmcarr is offline   Reply With Quote
Old 08-27-2010, 05:08 AM   #3
simonandrews
Simon Andrews
 
Location: Babraham Inst, Cambridge, UK

Join Date: May 2009
Posts: 871
Default

Incidentally, you can request that your sequencing service provide fastq files still. The latest Illumina pipeline doesn't generate them by default, but there's an option to turn them back on.
simonandrews is offline   Reply With Quote
Reply

Tags
flow cell, ngs, qseq

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 08:05 PM.


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