This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Ucsc Dbsnp131 Table Format

Hi, I want to create a bed file that contains leftflank-rightflank positions of a subset of dbsnp131 table I downloaded from UCSC. The problem is when I look at the data there are two fileds chromStart and ChromEnd but in most cases that I saw these were sequencial numbers like Chr1 134567 134568 .So what do those chrStart and chromEnd positions exactly mean and how can I get a leftflank-right flank position data for the snps that I am interested in? Thanks

mysql> select * from snp131 limit 10;

bin | chrom | chromStart | chromEnd | name| score | strand | refNCBI | refUCSC | observed | molType

| 585 | chr1 | 10433 | 10433 | rs56289060 | 0 | + | - | - | -/C | genomic | insertion | unknown | 0 | 0 | near-gene-5 | between | 1 |

| 585 | chr1 | 10491 | 10492 | rs55998931 | 0 | + | C | C | C/T | genomic | single | unknown | 0 | 0 | near-gene-5 | exact | 1 |

| 585 | chr1 | 10518 | 10519 | rs62636508 | 0 | + | G | G | C/G | genomic | single | unknown | 0 | 0 | near-gene-5 | exact | 1 |

dbsnp bed ucsc
• 3,788 views • link
updated by Stephane Plaisance Stephane Plaisance • written by Biomed Biomed

2 answers

It is an half-open interval.

chromStart: is the starting position ( 0-based) of the SNP on the chromosome

chromEnd: is the position where there is NO MORE snp.

so chromEnd-chromStart = length of the variation. e.g. if length=0, the SNP is a deletion

• 0 views • link
updated by Ram Ram • written by Pierre Lindenbaum Pierre Lindenbaum

SO for a Single nucleotide variation chromstart is actually where the variation occures and the leftflank would be chromstart-1 while the rightflank would be chromEnd, right?

• 0 views • link
written by Biomed Biomed

and

BedTools has commands to retreive (boundary) sequences from the fasta ref genome and a BED file

the current SNP track first need be processed to add few bases each side using yet another bed tool. Then you get the fasta sequences with fastaFromBed

have a look to the version v2.9.0 manual about the command details slopBed (get bordering bases) fastaFromBed (extract resulting sequences)

I love BedTools !

Cheers, Stephane, BITS-VIB

• 1 views • link
written by Stephane Plaisance Stephane Plaisance

Log in to answer this question.

More posts like this
  • News: A Json-Based Annotation Format
    written by Pierre Lindenbaum Pierre Lindenbaum
    FYI, I wrote a java [FeatureCodec][1] for JSON with a the [tribble][2] library. See: https://github.com/lindenb/jsontribble The library contains the tools to **sort**, **index** and **query** …
  • Snp Position Problem
    written by Biomed Biomed
    <p>I need help regarding snp positions on dbsnp130. I downloaded the data files for dbsnp130 from UCSC download area <a href="http://hgdownload.cse.ucsc.edu/goldenPath/hg18/database/snp130.txt.gz">dbsnp 130 txt</a> and dumped …
  • Multiple Entries With The Same Rs Number In Dbsnp131
    written by Biomed Biomed
    <p>when I run the query </p> <pre><code>mysql&gt; select * from snp131 where name='rs67465876'; +------+-------+------------+-----------+------------+-------+--------+---------+---------+----------+---------+--------+-------------------------------------+-------+---------+----------------+---------+--------+ | bin | chrom | chromStart | chromEnd | name | …