SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
RNA-Seq: SAVoR: a server for sequencing annotation and visualization of RNA structure Newsbot! Literature Watch 0 04-12-2012 03:00 AM
RNA-Seq: Detection, annotation and visualization of alternative splicing from RNA-Seq Newsbot! Literature Watch 0 01-10-2012 04:00 AM

Reply
 
Thread Tools
Old 03-09-2013, 03:30 AM   #1
tamari
Member
 
Location: Israel

Join Date: Oct 2012
Posts: 12
Red face RNA-seq data annotation

Hi,

I have an alingment of PE reads (100bp) from RNA-seq (BAM file) and I would like to know what is the best way to estimate the percantege of the genomic features in the library that were aligned, i.e exon, intron, 5UTR, 3UTR etc... as well as types of RNAs like coding, non-coding ( miRNA, lincRNA ,tRNA/rRNA...)

I know that there is something in bedtools package (perhaps intersect or coverage) but I am not sure how to use it in the best way so I get the most accurate picture because I expect overlaps between different features.

I tried to look for something simillar in R (bioconductor ) but I can't find anything.

I would really appreciate your help,

Tamari
tamari is offline   Reply With Quote
Old 03-09-2013, 03:51 AM   #2
Bukowski
Senior Member
 
Location: Aberdeen, Scotland

Join Date: Jan 2010
Posts: 388
Default

http://picard.sourceforge.net/comman...tRnaSeqMetrics
Bukowski is offline   Reply With Quote
Old 03-09-2013, 06:45 AM   #3
pbluescript
Senior Member
 
Location: Boston

Join Date: Nov 2009
Posts: 224
Default

This software can do many of the things you want:
https://code.google.com/p/rseqc/

If you want to get really specific about the types of RNA that are represented in your sample, then using something like bedtools intersect with custom GTF files would be useful.
pbluescript is offline   Reply With Quote
Old 03-10-2013, 11:29 PM   #4
tamari
Member
 
Location: Israel

Join Date: Oct 2012
Posts: 12
Default

Thanks for your help!
tamari is offline   Reply With Quote
Reply

Tags
annotation, annotation mouse, bam, rna seq

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 09:49 AM.


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