SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Planning an RNA-Seq Experiment gavin.oliver RNA Sequencing 20 02-12-2014 12:08 AM
rna-seq experiment not strand-specific or coverage symmetry laupl Introductions 2 05-09-2012 07:19 PM
regarding depth and coverage of rna seq jyoshna.jo RNA Sequencing 0 11-18-2011 09:11 PM
RNA-Seq: ReCount: A multi-experiment resource of analysis-ready RNA-seq gene count da Newsbot! Literature Watch 0 11-18-2011 02:20 AM
Depth of coverage variation in RNA Seq jay RNA Sequencing 19 08-11-2010 11:26 PM

Reply
 
Thread Tools
Old 03-17-2011, 06:46 PM   #1
PFS
Member
 
Location: USA

Join Date: Mar 2010
Posts: 55
Default coverage depth, bedtools and RNA-Seq experiment

I want to estimate the % of a transcriptome covered by an RNA-Seq experiment.

I tried to use BEDTools: one file consisted of exons only coordinates (from UCSC), the other file consisted of TopHat alignments.

The results show that 70% of the transcriptome has coverage depth 0.

I was expecting a much smaller percentage, given that I started with 100,000,000 reads x 2 (PE) and the total length of the transcriptome is 17,470,227.

Am I making a big mistake here?

ps. when I say "transcriptome", I mean the collection of all exons.
PFS is offline   Reply With Quote
Old 03-18-2011, 12:04 AM   #2
natstreet
Member
 
Location: Sweden

Join Date: Nov 2009
Posts: 83
Default

Can you show us the command you used to generate that value?
natstreet is offline   Reply With Quote
Old 03-18-2011, 07:27 AM   #3
PFS
Member
 
Location: USA

Join Date: Mar 2010
Posts: 55
Default

Quote:
Originally Posted by natstreet View Post
Can you show us the command you used to generate that value?
coverageBed -abam my_accepted_hits.bam -b rn4_UCSC_kg_exons.bed -hist -s -split
PFS 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 07:08 AM.


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