SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
ENCODE RNA-Seq cdias RNA Sequencing 1 09-25-2012 11:00 PM
ENCODE RNA-seq Quality Dario1984 RNA Sequencing 2 08-15-2012 01:13 PM
Question: Information about the Caltech RNA-Seq data on ENCODE and its corresponding tmy1018 Bioinformatics 1 08-15-2012 01:04 PM
is is possible to visualize RNA-seq data in UCSC genome browser(ENCODE)? HSV-1 RNA Sequencing 5 07-21-2012 07:24 PM

Reply
 
Thread Tools
Old 12-12-2012, 06:58 AM   #1
Karenj
Member
 
Location: Aarhus

Join Date: Dec 2012
Posts: 13
Default How to work with ENCODE RNA-seq fastq files?

I have downloaded the fastq files for paired end RNA-sequencing for a single cell line from the ENCODE project, as I want to map the reads my self in the search for chimeric read pairs.
My problem is that for each sample, the read1 and read 2 fastq.tgz folders contain 5 separate fastq-files. Should these files be merged before or midway in my mapping? If yes - how?
Karenj is offline   Reply With Quote
Old 12-12-2012, 07:04 AM   #2
alexdobin
Senior Member
 
Location: NY

Join Date: Feb 2009
Posts: 161
Default

Hi @Karenj,

which particular files did you download and where did you download them?
My guess is you are seeing the files separate by Illumina lanes.
The lanes are supposed to be merged on the UCSC server for most of the samples. If not, you can merge them before or while mapping. It's important that the order of merging is exactly the same for both reads 1 and 2.

Alex
alexdobin is offline   Reply With Quote
Old 12-12-2012, 07:10 AM   #3
Karenj
Member
 
Location: Aarhus

Join Date: Dec 2012
Posts: 13
Default

Hi Alex,

The files are RNA-seq for K562 from the Caltech data set, 2x75 paired end. I found few minutes ago, that for the CSHL RNA-seq data set the fastq files are merged and the data set is newer, so perhaps it is easier to work with this.
Karenj is offline   Reply With Quote
Old 12-12-2012, 07:17 AM   #4
alexdobin
Senior Member
 
Location: NY

Join Date: Feb 2009
Posts: 161
Default

I would also recommend working with CSHL dataset
It's indeed newer and is more uniform in terms of quality and sequencing depth.
Note, that CSHL's dataset is stranded dUTP protocol, while Caltech's is mostly standard unstranded.
alexdobin is offline   Reply With Quote
Old 12-12-2012, 08:18 AM   #5
Karenj
Member
 
Location: Aarhus

Join Date: Dec 2012
Posts: 13
Default

Thanks for your answer. I will try to see if the data from CSHL behaves better :-)
Karenj 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 04:32 PM.


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