Question: Finding intergenic transcripts from a bedgraph file or bam file
0
gravatar for harris
6.3 years ago by
harris0
United States
harris0 wrote:

Hi, I have  a series of bedgraph files made from .bam files using bedtools-2.17.0/genomeCoverageBed -ibam my.sorted.bam -bg > my.sorted.bg command. I can load into UCSC browser and works fine. Now my question is, I did paired end strand specific libraries and analysis.  From my .bg files I want to find the intergenic transcripts Ie ncRNAs, possibe eRNAs( a complete list with positions that I can then work with).  I have another file with the chromsome and start and stop of all the transcriptions units, about 23,000 (from a cufflink file). How do I use bedtools or any tool to find the transcripts that are NOT in the transcription units? Thanks for any help.  I know this is basic, but I have to start somewhere. Thanks for any help. steve

rna-seq • 1.8k views
ADD COMMENTlink modified 6.3 years ago by Sukhdeep Singh10k • written 6.3 years ago by harris0
0
gravatar for Sukhdeep Singh
6.3 years ago by
Netherlands
Sukhdeep Singh10k wrote:

For the second question, use intersectBed with -v parameter on the bedGraph file and file transcript location.

-v Only report those entries in A that have no overlap in B. Restricted by -f and -r.

From my .bg files I want to find the intergenic transcripts Ie ncRNAs, possibe eRNAs( a complete list with positions that I can then work with)

You have to fetch files for the locations you want to work with (intergenic, exon, intron etc) and again use intersectBed

ADD COMMENTlink modified 8 months ago by RamRS30k • written 6.3 years ago by Sukhdeep Singh10k

Thanks so much Sukhdeep, will give your thoughts a try. steve

ADD REPLYlink modified 6.3 years ago by Sukhdeep Singh10k • written 6.3 years ago by harris0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 593 users visited in the last hour