SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Advice on RNA-Seq Depth KDS RNA Sequencing 3 02-19-2014 03:03 PM
Small RNA sequencing depth comparison to array or qPCR seqfan2011 RNA Sequencing 0 05-03-2012 08:08 AM
RNA-seq depth/plant Lucilia RNA Sequencing 3 10-27-2011 11:29 AM
RNA-Seq: Differential expression in RNA-seq: A matter of depth. Newsbot! Literature Watch 0 09-10-2011 02:12 AM
how to calculate ka/ks with RNA-seq data boyzoe Bioinformatics 0 11-03-2010 07:03 PM

Reply
 
Thread Tools
Old 11-02-2012, 09:06 AM   #1
wangli
Member
 
Location: Texas

Join Date: Apr 2012
Posts: 48
Default how to calculate the mean sequencing depth of RNA-seq

Dear users

I tried to calculate the mean sequencing depth of RNA-seq. Firstly, I utilized GATK to calculate that with the whole genome as reference. It turned out that some samples with around 15 million mapped reads only have around 2.5* mean sequencing depth. That looks to me weird. Some people suggested me to calculate it by hand, such as: (number of mapped reads * read length)/total length of the species' transcriptome

It looks reasonable. How do you calculate the sequencing depth of RNA-seq?
I will appreciate your comments very much!
wangli is offline   Reply With Quote
Old 02-18-2014, 05:38 PM   #2
woodydon
Member
 
Location: https://www.researchgate.net/profile/Woody_Lin

Join Date: Jan 2010
Posts: 52
Default

You should calculate the sequencing depth of a gene or an exon. It doesn't make sense to have a number representing the overall sequencing depth of RNA-Seq. DNA and RNA-Seq are different.
woodydon is offline   Reply With Quote
Old 02-19-2014, 02:56 PM   #3
puggie
Member
 
Location: Europe

Join Date: Nov 2011
Posts: 52
Default

You could use e.g. Cufflinks or DEseq to calculate the library depths, at least that would provide you a scaling factor for inter-library comparisons and normalization. I think the "sequencing depth" as per DNA-seq is not comparable in this respect at all.
puggie is offline   Reply With Quote
Old 02-19-2014, 04:44 PM   #4
Jeremy
Senior Member
 
Location: Pathum Thani, Thailand

Join Date: Nov 2009
Posts: 190
Default

The most informative coverage metric for RNA-seq is the number of genes with x number of reads, e.g. 80% of the known genes had at least 10 reads mapping to them, 95% of the known genes had at least one read map.
Jeremy 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:36 AM.


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