This post will show how to get the reference base of chr1 from 49999 to 500001 (Version: hg19). Please note: different tools has different coordinate (0 start or 1 start).

## 1, SAMtools

Index reference sequence in the FASTA format or extract subsequence from indexed reference sequence. If no region is specified, faidx will index the file and create .fai on the disk. If regions are specified, the subsequences will be retrieved and printed to stdout in the FASTA format.

## 5, UCSC UCSC Genome Browser

just input the chr and position in UCSC genome browser. http://genome.ucsc.edu/cgi-bin/hgTables

## 6，togows

http://togows.org/ TogoWS was originally developed to realize uniformed REST services over the public REST and SOAP services provided by NCBI, EBI, DDBJ, PDB, and KEGG. Because most of those SOAP services were discontinued in 2012, we then focus on providing intuitive REST APIs, inclusion of UCSC genome databases and support for the Semantic Web.

 1 2 3  \$wget http://togows.org/api/ucsc/hg19/chr1:49999-50001.fasta -O - -q >hg19:chr1:49999-50001 ATA 

And many other usful tools can make it. If you have a good idea, make a comment below.

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

#版权所有 转载请告知 版权归作者所有 如有侵权 一经发现 必将追究其法律责任

#Author: Jason

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