Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • fanwei
    Member
    • Mar 2013
    • 21

    #1

    Error with BEDTools

    Dear all,
    I want your help!
    I want to use bedtools to calculate the overlap between two bed format files.It recorded chromosome,startposition and endposition informations.
    When i run bedtools, there was no error indicated and no results displayed. But i'm sure overlaps exist!
    Maybe bedtools restricted the file size?
    I feel confused. Can you kindly help me?
  • Jetse
    Member
    • Nov 2012
    • 38

    #2
    What is the command you used?
    How do the files look like?

    Comment

    • fanwei
      Member
      • Mar 2013
      • 21

      #3
      Originally posted by Jetse View Post
      What is the command you used?
      How do the files look like?
      I used bedtools intersect -a A.bed -b.B.bed
      The first few lines are as follows,
      Click image for larger version

Name:	A.png
Views:	1
Size:	10.7 KB
ID:	304292

      Click image for larger version

Name:	B.png
Views:	1
Size:	13.5 KB
ID:	304293

      Comment

      • Jetse
        Member
        • Nov 2012
        • 38

        #4
        What is your output when you use the command:
        Code:
        bedtools intersect -a A.bed -b B.bed -wao | less
        Do you see a lot of output, all with a zero at the end? (this means they don't have overlaps)

        Comment

        • fanwei
          Member
          • Mar 2013
          • 21

          #5
          Originally posted by Jetse View Post
          What is your output when you use the command:
          Code:
          bedtools intersect -a A.bed -b B.bed -wao | less
          Do you see a lot of output, all with a zero at the end? (this means they don't have overlaps)
          There was nothing even zero.
          But it works when i write some data in new files to test that.
          Is it necessary to sort my data?

          Comment

          • Jetse
            Member
            • Nov 2012
            • 38

            #6
            It could be your computer ran out of memory, then sorting would solve your problem.

            From the manual:
            Note
            If you are trying to intersect very large files and are having trouble with excessive memory usage, please presort your data by chromosome and then by start position (e.g., sort k1,1 -k2,2n in.bed > in.sorted.bed for BED files) and then use the -sorted option. This invokes a memory-efficient algorithm designed for large files.

            Comment

            • fanwei
              Member
              • Mar 2013
              • 21

              #7
              Originally posted by Jetse View Post
              It could be your computer ran out of memory, then sorting would solve your problem.

              From the manual:
              I have tried it,but it doesn't work.
              Any other suggestions?

              Comment

              • fanwei
                Member
                • Mar 2013
                • 21

                #8
                Originally posted by Jetse View Post
                It could be your computer ran out of memory, then sorting would solve your problem.

                From the manual:
                It has been resolved.Thank you!

                Comment

                • Jetse
                  Member
                  • Nov 2012
                  • 38

                  #9
                  How did you solve it?
                  So other people who find this thread with the same problem can try your solution too...

                  Comment

                  • fanwei
                    Member
                    • Mar 2013
                    • 21

                    #10
                    Originally posted by Jetse View Post
                    How did you solve it?
                    So other people who find this thread with the same problem can try your solution too...

                    It is due to file format error caused by newline character.
                    And your suggestions are also really important!

                    Comment

                    Latest Articles

                    Collapse

                    • SEQadmin2
                      Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
                      by SEQadmin2



                      Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

                      Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
                      ...
                      10-14-2026, 09:55 AM
                    • SEQadmin2
                      New Genomics Technologies Take Aim at Long-Standing Limits
                      by SEQadmin2


                      Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

                      We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
                      ...
                      09-28-2026, 10:25 AM

                    ad_right_rmr

                    Collapse

                    News

                    Collapse

                    Topics Statistics Last Post
                    Started by SEQadmin2, Yesterday, 11:13 AM
                    0 responses
                    19 views
                    0 reactions
                    Last Post SEQadmin2  
                    Started by SEQadmin2, Yesterday, 10:48 AM
                    0 responses
                    18 views
                    0 reactions
                    Last Post SEQadmin2  
                    Started by SEQadmin2, 10-15-2026, 10:52 AM
                    0 responses
                    36 views
                    0 reactions
                    Last Post SEQadmin2  
                    Started by SEQadmin2, 10-14-2026, 05:39 AM
                    0 responses
                    42 views
                    0 reactions
                    Last Post SEQadmin2  
                    Working...