Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • cliff
    Member
    • Oct 2009
    • 41

    #1

    Question about BED format - chromStart and End

    I am a bit confused about the chromStart and chromEnd positions in the BED format.

    According to UCSC:
    chromStart - The starting position of the feature in the chromosome or scaffold. The first base in a chromosome is numbered 0.
    chromEnd - The ending position of the feature in the chromosome or scaffold. The chromEnd base is not included in the display of the feature. For example, the first 100 bases of a chromosome are defined as chromStart=0, chromEnd=100, and span the bases numbered 0-99.

    Assuming I download a bed file for a gene from UCSC as below:

    chromStart: 300
    chromEnd: 500

    Now, I get a set of SNPs by mapping reads to hg18 and calling SNPs using whatever SNPcaller. I want to know how many SNPs were called within the above gene. Should I compare each SNP position with the gene range as

    300<=SNP_POSITION<=500

    or

    301<=SNP_POSITION<=500

    or

    301<=SNP_POSITION<=499

    ?

    Does anyone know which is correct?

    Thanks
  • ffinkernagel
    Senior Member
    • Oct 2009
    • 110

    #2
    300 <= x < 500

    so the first base is no. 300, the last base is no .499, and the range covers 200 bases.
    Just as the doc says. You may imagine the mark being between the bases.

    Comment

    • lh3
      Senior Member
      • Feb 2008
      • 686

      #3
      301st<=snp_position<=500th

      EDIT:

      BED is always 0-based. The first base in a sequence has coordinate 0 and therefore coordinate 300 denotes the 301st base. A more obvious example is

      0 1

      which denotes the first base.
      Last edited by lh3; 10-25-2010, 04:14 AM.

      Comment

      • Hena
        Member
        • Nov 2009
        • 19

        #4
        It depends on how the range is defined: 0 based or 1 based positions? If one based, then 301 <= SNP_POSITION <= 500 is the range you had in the bed file.
        Last edited by Hena; 10-25-2010, 03:43 AM.

        Comment

        • ndaniel
          Member
          • Feb 2009
          • 33

          #5
          300 <= x < 500

          Comment

          • ndaniel
            Member
            • Feb 2009
            • 33

            #6
            Originally posted by lh3 View Post
            301st<=snp_position<=500th

            EDIT:

            BED is always 0-based. The first base in a sequence has coordinate 0 and therefore coordinate 300 denotes the 301st base. A more obvious example is

            0 1

            which denotes the first base.
            In case that everything is zero-based one has for example for range 3 (instead of 300) to 5 (instead of 500):


            0123456789 <-- positions
            ---gg----- <-- gene in range [3,5)


            therefore the correct answer is 300<=x<500!

            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, 10:52 AM
            0 responses
            9 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 10-14-2026, 05:39 AM
            0 responses
            19 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 10-13-2026, 09:56 AM
            0 responses
            16 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 10-09-2026, 11:14 AM
            0 responses
            27 views
            0 reactions
            Last Post SEQadmin2  
            Working...