SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Merging reads mapped on genome and CDS (SOLID data) jmandel Bioinformatics 1 01-04-2012 06:34 AM
Reads mapping of ends of chromosomes ETHANol Bioinformatics 0 11-03-2011 09:09 AM
Paired-end reads mapping to different chromosomes... rareaquaticbadger Bioinformatics 1 06-27-2011 04:24 AM
Finding exon-exon junction vincebrown General 10 11-16-2010 04:08 PM
tophat paired reads mapped to different chromosomes wenhuang Bioinformatics 3 06-23-2010 05:23 AM

Reply
 
Thread Tools
Old 09-17-2010, 05:41 AM   #1
sridharacharya
Member
 
Location: Institute, WV

Join Date: May 2010
Posts: 24
Arrow Exon-Junction mapping: re-assigning CDS-mapped reads to chromosomes

Hi,

I am working with RNA-SEQ data and want to align SOLiD reads using publicly available software. Tophat cannot handle SOLiD reads, so I am going with Bowtie. But the problem is the mapping of exon-junctions. For this, I downloaded the REF-SEQ cDNA sequences from UCSC.

So my strategy is to first use align the reads to the genome (chromosomes) and then align the remaining unmapped reads to REF-SEQ cDNA sequences. Subsequently, I want to re-assign the "REF-SEQ cDNA aligned reads" to the chromosomal sequences, so that I can use the alignments with program like Cufflinks, to find enrichment.

My question is regarding re-assigning the "REF-SEQ cDNA aligned reads" to the chromosomal sequences? Is there a tool available for this? Is there a program available apart from Cufflinks to find the enrichment of reads/compare two experiments?

Thankyou for your suggestions!
sridharacharya is offline   Reply With Quote
Old 10-21-2010, 04:07 PM   #2
glacerda
Member
 
Location: Brazil

Join Date: Aug 2008
Posts: 27
Default

There are many programas for differential expression analysis in RNA-seq data. My favorites are edgeR and DEseq. Both are R packages, but they are well documented and easy to run, even for a R begginer like me.

You should input the raw read counts to these programs, so you could sum the number of reads mapped to genome and mapped to reference transcripts, summarizing this information for each locus, without counting twice reads that map to both genome and transcript.
glacerda is offline   Reply With Quote
Reply

Tags
exon-junction mapping

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:42 AM.


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