SEQanswers

Go Back   SEQanswers > General



Similar Threads
Thread Thread Starter Forum Replies Last Post
Converting FPKM from Cufflinks to raw counts for DESeq jebe Bioinformatics 34 02-05-2014 08:19 AM
Raw read counts for RNAseq biofreak Introductions 13 01-16-2013 05:28 AM
Extracting 'counts' from v values crh Bioinformatics 0 12-31-2011 09:29 AM
Raw read counts for RNAseq biofreak RNA Sequencing 2 06-15-2011 05:56 AM
How to convert cufflinks output to raw counts jebe RNA Sequencing 0 01-26-2010 11:29 AM

Reply
 
Thread Tools
Old 04-05-2010, 09:03 AM   #1
gen2prot
Member
 
Location: Hyderabad, India

Join Date: Apr 2010
Posts: 64
Default How can one get raw read counts from RPKM values

Hi All,

I am new to the transcriptome sequencing field and I am learning to use Bowtie, Tophat and Cufflinks. I find that Cufflinks spits out RPKM and FPKM values which though very informative does not help my purpose since I am interested in the number of raw reads per gene. That is to say without normalization. Is there a way of figuring this out?

Thanks
Abhijit
gen2prot is offline   Reply With Quote
Old 04-05-2010, 12:35 PM   #2
mgogol
Senior Member
 
Location: Kansas City

Join Date: Mar 2008
Posts: 193
Default

One way I've done it is using bedtools coverageBed with the bam file generated from the sam file from the tophat alignment and a bed file describing the genes.

I am curious to hear other people's techniques, because I'm sure there are more efficient ways to go about this.
mgogol is offline   Reply With Quote
Old 08-25-2010, 07:19 AM   #3
lexa
Member
 
Location: MPI

Join Date: Jun 2010
Posts: 17
Default

hi,

I`m also using bedtools but does anybody has a better solution?
lexa is offline   Reply With Quote
Old 06-22-2011, 08:04 PM   #4
edge
Senior Member
 
Location: China

Join Date: Sep 2009
Posts: 199
Default

Hi,

Do you mind to share the command that you type in bedtools in order to extract the read count per transcripts?
I'm facing the same problem as well.
Thanks.
edge is offline   Reply With Quote
Old 06-22-2011, 08:04 PM   #5
edge
Senior Member
 
Location: China

Join Date: Sep 2009
Posts: 199
Default

Hi,

Do you mind to share the command that you type in bedtools in order to extract the read count per transcripts?
I'm facing the same problem as well.
Thanks.
edge is offline   Reply With Quote
Old 06-23-2011, 07:07 AM   #6
gen2prot
Member
 
Location: Hyderabad, India

Join Date: Apr 2010
Posts: 64
Default Getting raw read counts

I ended up using Simon Ander's program for computing read counts per gene. However you will have to provide the annotation of the genome in GTF format in order to view that. I find that this is quite accurate. Cheers...
gen2prot is offline   Reply With Quote
Old 06-24-2011, 11:08 AM   #7
biofreak
Member
 
Location: USA

Join Date: Jun 2011
Posts: 44
Default

I tried Htseq-count software.
python -m HTSeq.scripts.count accepted_hits.sam hg19RefGene.gtf

it works. But I want to do this in R. Maybe functions in the GenomicRanges packages can help? e.g. countOverlaps function?
biofreak 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 01:44 AM.


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