SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Bedtools Intersection files bnfoguy Bioinformatics 1 06-17-2014 12:31 AM
Bedtools and VCF files morangal Bioinformatics 4 07-11-2012 10:55 PM
Tophat for multiple files. oxydeepu Bioinformatics 1 04-02-2012 10:52 AM
Genome files for BEDTools GenomeCoverage command qnc Bioinformatics 4 10-25-2011 04:30 AM
Multiple index files in Bowtie? spb2003 Bioinformatics 3 10-13-2009 04:22 PM

Reply
 
Thread Tools
Old 07-23-2015, 01:02 PM   #1
asapp
Junior Member
 
Location: Cambridge, MA

Join Date: Jun 2015
Posts: 8
Default Intersecting multiple bed files

Hello,

I have several RNA-Seq transcript data sets assembled and in BED format. I have masked out previously identified transcripts from the data, leaving me with "unidentified" transcripts. A visual of this data in the UCSC Genome Browser is attached.

I need to computationally gather the locations where all of the resulting transcripts intersect. In looking though the bedtools intersect function, it seems that there is no option to specify regions where ALL the files intersect, only those where at least 2 of the supplied files intersect.

Have any of you encountered this issue before? What are your suggested workarounds? Would I need to continually compare 2 files and then keep comparing that output to the next file and so on? Or is there an option I am unaware of/another software function that can handle this?

Thank you all very much!
-Anna
Attached Images
File Type: png seqanswers.png (142.6 KB, 5 views)
asapp is offline   Reply With Quote
Old 07-23-2015, 06:57 PM   #2
blancha
Senior Member
 
Location: Montreal

Join Date: May 2013
Posts: 367
Default

bedops --intersect first.bed second.bed third.bed > result.bed

https://bedops.readthedocs.org/en/la...ct-i-intersect
blancha is offline   Reply With Quote
Old 07-24-2015, 05:18 AM   #3
asapp
Junior Member
 
Location: Cambridge, MA

Join Date: Jun 2015
Posts: 8
Default

Perfect, thanks!
asapp is offline   Reply With Quote
Reply

Tags
bed, bedtools, bioinformatics, intersect, 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 12:12 AM.


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