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

                    • kinpict.com
                      Reply to 16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
                      by kinpict.com
                      The part about this topic feels surprisingly grounded, and the details on the page is the bit that made me slow down. I’d keep <a href="https://kinpict.com/">ai family photo generator</a> bookmarked for later.
                      Today, 06:37 PM
                    • SEQadmin2
                      16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
                      by SEQadmin2





                      NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
                      10-30-2026, 06:46 AM

                    ad_right_rmr

                    Collapse

                    News

                    Collapse

                    Topics Statistics Last Post
                    Started by SEQadmin2, 11-12-2026, 07:59 AM
                    0 responses
                    26 views
                    0 reactions
                    Last Post SEQadmin2  
                    Started by SEQadmin2, 11-05-2026, 09:25 AM
                    0 responses
                    35 views
                    0 reactions
                    Last Post SEQadmin2  
                    Started by SEQadmin2, 11-03-2026, 11:54 AM
                    2 responses
                    54 views
                    0 reactions
                    Last Post kinpict.com  
                    Started by SEQadmin2, 10-28-2026, 09:57 AM
                    0 responses
                    175 views
                    0 reactions
                    Last Post SEQadmin2  
                    Working...