我们一般会熟悉sam/bam格式文件，就是把测序reads比对到参考基因组后的文件！bam或者bed格式的文件主要是为了追踪我们的reads到底比对到了参加基因组的什么区域，而UCSC规定的这几个文件格式(wig、bigWig和bedgraph)用处不一样，仅仅是为了追踪参考基因组的各个区域的覆盖度，测序深度！而且这些定义好的文件，可以无缝连接到UCSC的Genome Browser工具里面进行可视化！
这个网站提供了这几种数据格式的构造及转换脚本：http://barcwiki.wi.mit.edu/wiki/SOPs/coordinates
具体文件格式定义请直接看UCSC的官网，下面是我基于自己的理解来翻译的，没什么特殊的，建议大家看原文，然后自己翻译一个，跟我比较！
- Wiggle Track Format (WIG)：http://genome.ucsc.edu/goldenPath/help/wiggle.html
- bigWig Track Format ：http://genome.ucsc.edu/goldenPath/help/bigWig.html
- BedGraph Track Format ：http://genome.ucsc.edu/goldenPath/help/bedgraph.html
这3种文件格式都是UCSC规定的，所以它提供了系列工具进行互相转换，可以直接下载可执行版本程序：http://hgdownload.cse.ucsc.edu/admin/exe/
常见的工具如下：
- bigWigToBedGraph — this program converts a bigWig file to ASCII bedGraph format.
- bigWigToWig — this program converts a bigWig file to wig format.
- bigWigSummary — this program extracts summary information from a bigWig file.
- bigWigAverageOverBed — this program computes the average score of a bigWig over each bed, which may have introns.
- bigWigInfo — this program prints out information about a bigWig file.
其实对我们的bam文件，用samtools软件也可以很容易得到基因组区域的覆盖度和测序深度，比如：
samtools depth -r chr12:126073855-126073965 Ip.sorted.bam chr12 126073855 5 chr12 126073856 15 chr12 126073857 31 chr12 126073858 40 chr12 126073859 44 chr12 126073860 52
~~~~~~~~~其余省略输出~~~~~~~~~
这其实就是wig文件的雏形，但是wig文件会更复杂一点！
首先它不需要第一列了，因为全部是重复字段，只需要在每个染色体的第一行定义好染色体即可。
首先需要设置这个wig文件在UCSC的Genome Browser工具里面显示的属性：
track type=wiggle_0 name=track_labeldescription=center_labelvisibility=display_modecolor=r,g,baltColor=r,g,bpriority=priorityautoScale=on|offalwaysZero=on|offgridDefault=on|offmaxHeightPixels=max:default:mingraphType=bar|pointsviewLimits=lower:upperyLineMark=real-valueyLineOnOff=on|offwindowingFunction=mean+whiskers|maximum|mean|minimumsmoothingWindow=off|2-16
type=wiggle_0 这个是默认的， 而且到目前为止，必须是这样的！其余的都是可选参数，自己读官网说明
这些参数一般不用管，除非你很熟悉了UCSC的Genome Browser工具
然后需要设置每条染色体的属性，几个比较重要的参数是：
fixedStepchrom=chrNstart=positionstep=stepInterval[span=windowSize]
下面是wig的一个具体例子：
track type=print wiggle_0 name=hek description=hek variableStep chrom=chr1 span=10 10008 7 10018 14 10028 27 10038 37 10048 45 10058 43 10068 37 10078 26
~~~~~~~~~其余省略输出~~~~~~~~~
UCSC也提供了一个wig文件：http://genome.ucsc.edu/goldenPath/help/examples/wiggleExample.txt
可以看到我设置的参数很少很少，而且我是直接对sort后的bam文件用脚本变成wig文件的。
那么bigwig格式文件就没什么好讲的了，它就是wig格式文件的二进制压缩版本，这样更加节省空间。
我们只需要用UCSC提供的工具把自己的wig文件转换一下即可，步骤如下：
Save this wiggle file to your machine (this satisfies steps 1 and 2 above).
Save this text file to your machine. It contains the chrom.sizes for the human (hg19) assembly (this satisfies step 4 above).
Download the wigToBigWig utility (see step 3).
Run the utility to create the bigWig output file (see step 5):
wigToBigWig wigVarStepExample.gz hg19.chrom.sizes myBigWig.bw
最后我们讲一下BedGraph格式文件，它是BED文件的扩展，是4列的BED格式，但是需要添加UCSC的Genome Browser工具里面显示的属性，但是一般就定义有限的几个属性即可。
track type=bedGraph name=track_labeldescription=center_label visibility=display_modecolor=r,g,baltColor=r,g,b priority=priorityautoScale=on|offalwaysZero=on|off gridDefault=on|offmaxHeightPixels=max:default:min graphType=bar|pointsviewLimits=lower:upper yLineMark=real-valueyLineOnOff=on|off windowingFunction=maximum|mean|minimumsmoothingWindow=off|2-16
有一点需要注意：
- These coordinates are zero-based, half-open.
- Chromosome positions are specified as 0-relative.
- The first chromosome position is 0.
- The last position in a chromosome of length N would be N – 1.
- Only positions specified have data.
- Positions not specified do not have data and will not be graphed.
- All positions specified in the input data must be in numerical order.
我这里有一个MACS对CHIP-seq数据call peaks附带的BedGraph文件，也可以用工具直接从bam格式文件得到：
track type=bedGraph name=”hek_treat_all” description=”Extended tag pileup from MACS version 1.4.2 20120305″ chr1 9997 9999 1 chr1 9999 10000 2 chr1 10000 10001 4 chr1 10001 10003 5 chr1 10003 10007 6 chr1 10007 10010 7 chr1 10010 10012 8 chr1 10012 10015 9 chr1 10015 10016 10 chr1 10016 10017 11 chr1 10017 10018 12
原文来自：http://www.bio-info-trainee.com/1815.html