Category
Format Conversion
Usage
twoBitToFa input.2bit output.fa
Manual
This tool is part of UCSC Genome Browser's Utility Tools.
Options
- -seq=name: restrict this to just one sequence
- -start=X: start at given position in sequence (zero-based)
- -end=X: end at given position in sequence (non-inclusive)
- -seqList=file: file containing list of the desired sequence names in the format seqSpec[:start-end], e.g. chr1 or chr1:0-189 where coordinates are half-open zero-based, i.e. [start,end)
- -noMask: convert sequence to all upper case
- -bpt=index.bpt: use bpt index instead of built in one
- -bed=input.bed: grab sequences specified by input.bed. Will exclude introns
- -bedPos: with -bed, to use chrom:start-end as the fasta ID in output.fa
- -udcDir=/dir/to/cache: place to put cache for remote bigBed/bigWigs
Sequence and range may also be specified as part of the input file name using the syntax:
/path/input.2bit:name or /path/input.2bit:name or /path/input.2bit:name:start-end
Example
Assume you have a sequence file in 2bit format, and you want to extract sequence for chr1 from that file:
twoBitToFa -seq=chr1 -udcDir=. http://your-website.edu/~user/genome.2bit stdout > genome.chr1.fa
Share your experience or ask a question