SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Using bedgraph to transcript assemble Feng Tian Bioinformatics 1 09-05-2014 12:05 AM
Coverage of a transcript and accuracy of DE bruce01 General 0 01-08-2014 09:03 AM
RNA-seq : 5 to 3 transcript coverage lnzzz Bioinformatics 7 01-07-2014 10:25 AM
coverage calculation - Bedgraph files AnushaC Bioinformatics 1 11-20-2013 08:09 AM
Transcript coverage estimates? tboothby General 2 12-14-2011 03:42 AM

Reply
 
Thread Tools
Old 01-22-2015, 10:38 PM   #1
manauwer
Junior Member
 
Location: new delhi, india

Join Date: Dec 2014
Posts: 3
Default get coverage of a particular transcript from bedgraph

i have a bedgraph..

1 0 10540 0
1 10540 10564 4
1 10564 10569 5
1 10569 10590 7
1 10590 10614 3
1 10614 10619 2
1 10619 12045 0

from reads mapped to the genome. now i want to extract a particular transcript
1 925741 944581
from the bedgraph, which would give me only those coordinates and coverage of the transcript provided.

are there any tools to do this job.

i can write perl script to do this, but i have many transcript to find out, please suggest me with a tool, if available.
manauwer is offline   Reply With Quote
Old 01-23-2015, 12:38 AM   #2
sarvidsson
Senior Member
 
Location: Berlin, Germany

Join Date: Jan 2015
Posts: 137
Default

You can extract the overlapping coverage data with bedtools intersect, e.g.

Code:
bedtools intersect -wb -a coverage.bdg -b transcripts.bed
(-wb is helpful to see which transcript in the "b" file overlapped the features in the "a" file - you could add transcript identifiers to the "b" for easier parsing)

But if you are actually interested in the actual number of reads overlapping with transcript features, you're probably better of using the original BAM with a GFF3/GTF file and HTSeq-count (http://www-huber.embl.de/users/ander...doc/count.html) - it has some further interesting features for RNA-Seq (if that's what you're into).
sarvidsson is offline   Reply With Quote
Reply

Tags
bedgraph, bioninformatics tool, tools

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 08:45 AM.


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