Created
March 7, 2011 23:19
-
-
Save arq5x/859487 to your computer and use it in GitHub Desktop.
Make a BED file from the raw UCSC file
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
export GENOME=hg19 | |
export SNPBUILD=131 | |
curl -s http://hgdownload.cse.ucsc.edu/goldenPath/$GENOME/database/snp$SNPBUILD.txt.gz | \ | |
zcat | \ | |
cut -f 2,3,4,5,6,7,10,16 > dbsnp.$SNPBUILD.$GENOME.bed | |
head dbsnp.$SNPBUILD.$GENOME.bed | |
chr1 10433 10433 rs56289060 0 + -/C near-gene-5 | |
chr1 10491 10492 rs55998931 0 + C/T near-gene-5 | |
chr1 10518 10519 rs62636508 0 + C/G near-gene-5 | |
chr1 10582 10583 rs58108140 0 + A/G near-gene-5 | |
chr1 10827 10828 rs10218492 0 + A/G near-gene-5 | |
chr1 10903 10904 rs10218493 0 + A/G near-gene-5 | |
chr1 10926 10927 rs10218527 0 + A/G near-gene-5 | |
chr1 10937 10938 rs28853987 0 + A/G near-gene-5 | |
chr1 11001 11002 rs79537094 0 + A/C near-gene-5 | |
chr1 11013 11014 rs28484712 0 + A/G near-gene-5 |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment