Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • BED Tools Tutorial

    Sorry I'm a total beginner in bioinformatics and really need some help.

    I'm on a mac and so far I've managed to install xcode and therefore was able to execute the installation process described here: http://code.google.com/p/bedtools/

    Now I've got the problem that I can't get anything started because I'm lacking the 'terminal' commands. So here I picked an example from the manual:

    ###########################

    5.4.1 Usage and option summary

    Usage: $ bamToBed [OPTIONS] -i <BAM>

    ###########################

    How would on use this command ?

    Thanx a lot in advance!

  • #2
    Do you have data in .bam format?

    Comment


    • #3
      Yes, but do I have to give the full path to the files or in which directory do I have to place them ?
      Which directory do I have to open in terminal, when I execute a command ?

      Comment


      • #4
        You can either give the full path to the files or work in the same directory as the files. You can call a command from any directory if it's installed somewhere in your PATH.

        Comment


        • #5
          I told you I'm an absolute beginner, but finally I've solved the problem.
          When I copied an example command like: $ bamToBed [OPTIONS] -i <BAM> I didn't remove the '$' so I had 2 '$' in the terminal which led to a 'command not known'. But now it works so I'll finally get back to you with some other basic usage questions

          Comment


          • #6
            Okay, as proposed here is my next question.
            When I use the coverageBed command, it gives me only the coverage for the specified regions in file b but not of the whole span. E.g. reference file b contains chr1: 1000-2000 , then I only get the mean coverage for this region. How can I get the per base coverage in this region ?

            Comment


            • #7
              Use the -d option for per base coverage. Available in version 2.9 and later.
              Aaron

              Comment


              • #8
                -d is what I used. The problem is, if in the reference b file are given e.g. 5 ranges like described before, I only get one value per each range. What I would like to have is for every base in the range. In this case: chr1: 1000-2000 I would expect 1000 lines of output.

                Comment


                • #9
                  That is odd. Could you post the command you used and the output?
                  Section 5.9.6 of the manual illustrates the way the output should loo.
                  Last edited by quinlana; 09-28-2010, 08:33 AM.

                  Comment

                  Latest Articles

                  Collapse

                  • seqadmin
                    The Impact of AI in Genomic Medicine
                    by seqadmin



                    Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
                    02-26-2024, 02:07 PM
                  • seqadmin
                    Multiomics Techniques Advancing Disease Research
                    by seqadmin


                    New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

                    A major leap in the field has
                    ...
                    02-08-2024, 06:33 AM

                  ad_right_rmr

                  Collapse

                  News

                  Collapse

                  Topics Statistics Last Post
                  Started by seqadmin, 02-28-2024, 06:12 AM
                  0 responses
                  22 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 02-23-2024, 04:11 PM
                  0 responses
                  70 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 02-21-2024, 08:52 AM
                  0 responses
                  78 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 02-20-2024, 08:57 AM
                  0 responses
                  67 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X