bed 文件怎么生成（BED文件格式）

BED (Browser Extensible Data)格式文件就是通过规定行的内容来展示注释信息，今天小编就来说说关于bed 文件怎么生成?下面更多详细答案一起来看看吧!

bed 文件怎么生成

BED (Browser Extensible Data)格式文件就是通过规定行的内容来展示注释信息

BED文件每行至少包括chrom，chromStart，chromEnd三列；另外还可以添加额外的9列，这些列的顺序是固定的。

在自定义BED文件时，前面可以有注释行，以“browser”或“track”开头，可以设置一些参数便于浏览器更好展示BED文件信息。但是，下游的一些分析工具，例如bedToBigBed，是不接受有注释的BED文件的。

# BED文件必须的3列:

1. chrom - 染色体号; 例如，chr1，chrX。。。。。。。
2. chromStart - feature在染色体上起始位置. 从0开始算，染色体上第一个碱基位置标记为0。
3. chromEnd - feature在染色体上终止位置。染色体上前100个碱基片段的位置位置标记为：chromStart=0, chromEnd=100。 实际上，第100个碱基不属于当前片段中，当前片段的碱基应该是0-99。所以在BED文件中，起始位置从0开始，终止位置从1开始。

4.name - BED行名，在基因组浏览器左边显示；

5.score - 在基因组浏览器中显示的灰度设定，值介于0-1000；

6.strand - 正负链标记. Either "." (=no strand) or " " or "-".

7.thickStart - feature起始位置。绘制特征的起始位置（例如，基因显示中的起始密码子）。当没有这部分时，thickStart和thickEnd通常设置为chromStart位置。

8.thickEnd - feature编码终止位置

9.itemRgb - R,G,B (e.g. 255,0,0)值，当itemRgb 设置为 "On"，BED的行会显示颜色

10.blockCount - blocks (exons)数目

11.blockSizes - blocks (exons)大小列表，逗号分隔，对应于blockCount

12.blockStarts -blocks (exons)起始位置列表，逗号分隔，对应于blockCount.；这个起始位置是与chromStart的一个相对位置。

例如：

bio@ubuntu ~]$ less -S GRCh38.gene.bed chr3 124792319 124792562 ENSG00000276626 RF00100 - chr1 92700819 92700934 ENSG00000201317 RNU4-59P - chr14 100951856 100951933 ENSG00000200823 SNORD114-2 chr22 45200954 45201019 ENSG00000221598 MIR1249 - chr1 161699506 161699607 ENSG00000199595 RF00019