UCSC常见格式:http://genome.ucsc.edu/FAQ/FAQformat.html UCSC工具:http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/
1 Wig
示例:
fixedStep chrom=chr3 start=400601 step=100
11
22
33
2 bigWig
http://genome.ucsc.edu/goldenPath/help/bigWig.html wig生成的带索引的二进制文件,eg. bigWigExample.bw eg. hg19.chrom.sizes
生成命令:
wigToBigWig input.wig chrom.sizes myBigWig.bw
查看命令:
bigWigToBedGraph
bigWigToWig
bigWigSummary
bigWigAverageOverBed
bigWigInfo
3 bedGraph
http://genome.ucsc.edu/goldenPath/help/bedgraph.html 提取bigWig文件
bigWigToBedGraph byBigWig.bw out.bedGraph
示例:
browser position chr19:49302001-49304701
browser hide all
browser pack refGene encodeRegions
browser full altGraph
# 300 base wide bar graph, autoScale is on by default == graphing
# limits will dynamically change to always show full range of data
# in viewing window, priority = 20 positions this as the second graph
# Note, zero-relative, half-open coordinate system in use for bedGraph format
track type=bedGraph name="BedGraph Format" description="BedGraph format" visibility=full color=200,100,0 altColor=0,100,200 priority=20
chr19 49302000 49302300 -1.0
chr19 49302300 49302600 -0.75
chr19 49302600 49302900 -0.50
chr19 49302900 49303200 -0.25
chr19 49303200 49303500 0.0
chr19 49303500 49303800 0.25
chr19 49303800 49304100 0.50
chr19 49304100 49304400 0.75
chr19 49304400 49304700 1.00