Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • find overlaps/common in multiple bed file

    Does anyway have experience to find commonalies (or overlaps) between bed files. I have 4 bed files and want to find what intervals are common to all.

    And I don't have any rigid criteria for overlap, any intersection will do.

    Appreciate any answers.

  • #2
    Have a look at bedtools, specifically the intersectBed utility.

    Comment


    • #3
      http://pypi.python.org/pypi/pybedtools is a python extension of the bedtools that was mentioned in the post #2 by dpryan.

      Comment


      • #4
        Originally posted by epi View Post
        Does anyway have experience to find commonalies (or overlaps) between bed files. I have 4 bed files and want to find what intervals are common to all.

        And I don't have any rigid criteria for overlap, any intersection will do.

        Appreciate any answers.
        I mentioned the wrong function I meant to mention Galaxy with the Intersect function in the operate on genomic intervals functions
        here is example


        Last edited by husamia; 01-23-2012, 07:03 AM. Reason: corrected intersect instead of merege for overlapping bed files

        Comment


        • #5
          Originally posted by dpryan View Post
          Have a look at bedtools, specifically the intersectBed utility.
          Thanks for response everyone !

          I could implement intersectBed as pairwise. And to find common among 4, I can do multiple pairwise, but it seems there is some chance for false negatives.

          I dint get chance to look at the python script so far, I wonder if that adresses this issue.

          Comment


          • #6
            I mentioned the wrong function, I modified my response above to reflect the correct reference function. Sorry for confusion

            Comment


            • #7
              Originally posted by epi View Post
              Thanks for response everyone !

              I could implement intersectBed as pairwise. And to find common among 4, I can do multiple pairwise, but it seems there is some chance for false negatives.

              I dint get chance to look at the python script so far, I wonder if that adresses this issue.
              The python interface isn't very different from direct command line usage and I would suspect produces the same results. I also don't see how you'd get a false negative, provided you actually want regions existing in all of the bed files.

              Comment


              • #8
                Originally posted by dpryan View Post
                The python interface isn't very different from direct command line usage and I would suspect produces the same results. I also don't see how you'd get a false negative, provided you actually want regions existing in all of the bed files.
                Thanks again for responding. I realize i did not state my objective well enough.
                This is chip-seq analysis for which I have bed files (peaks). There could be a situation when peak 1 intersect peak 3 at 5' and peak 2 intersect peak 3 at 3`. but peak 1 and peak 2 do not intersect.
                Read Peak1, Peak2 and Peak3 coming form Samples 1,2 and 3 please.
                intersectBed will not reveal these peaks on my first paiwise comparison (peak1 and peak2) so it will be gone. Logically, it seems they come from same region so I was wondering if there is a tool that could capture those.

                I must mention that this is not a real example, just i theoretical possibility which crossed my mind. May b I am just too obsessed over it

                Comment


                • #9
                  You want multiIntersectBed...

                  Comment


                  • #10
                    Originally posted by epi View Post
                    Read Peak1, Peak2 and Peak3 coming form Samples 1,2 and 3 please.
                    intersectBed will not reveal these peaks on my first paiwise comparison (peak1 and peak2) so it will be gone. Logically, it seems they come from same region so I was wondering if there is a tool that could capture those.
                    You can include both A and B regions in the output from intersectBed (e.g. with -wo), merge them with mergeBed and then intersect the merged output to the next sample (to "grow" the overlapping regions).

                    Comment


                    • #11
                      Originally posted by mgogol View Post
                      bullseye !!

                      Looks like can the job, will try it out ...
                      appreciate your response

                      Comment


                      • #12
                        If you prefer a much more scalable solution that can do this simple intersection (and any other set-like operation) on any number of bedfiles at once, check out BEDOPS.

                        Comment

                        Latest Articles

                        Collapse

                        • seqadmin
                          Strategies for Sequencing Challenging Samples
                          by seqadmin


                          Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                          03-22-2024, 06:39 AM
                        • seqadmin
                          Techniques and Challenges in Conservation Genomics
                          by seqadmin



                          The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                          Avian Conservation
                          Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                          03-08-2024, 10:41 AM

                        ad_right_rmr

                        Collapse

                        News

                        Collapse

                        Topics Statistics Last Post
                        Started by seqadmin, 03-27-2024, 06:37 PM
                        0 responses
                        12 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 03-27-2024, 06:07 PM
                        0 responses
                        11 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 03-22-2024, 10:03 AM
                        0 responses
                        53 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 03-21-2024, 07:32 AM
                        0 responses
                        68 views
                        0 likes
                        Last Post seqadmin  
                        Working...
                        X