有用的UCSC tools推薦

先上網(wǎng)站鏈接奖年,這些都不需要安裝,直接下載下來就可以用,比較便捷
http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/

1. bedGeneParts - Given a bed, spit out promoter, first exon, or all introns.

bedGeneParts part in.bed out.bed
Where part is either 'exons' or 'firstExon' or 'introns' or 'promoter' or 'firstCodingSplice' or 'secondCodingSplice'
options:
   -proStart=NN - start of promoter relative to txStart, default -100
   -proEnd=NN - end of promoter relative to txStart, default 50

2. bedGraphToBigWig v 4 - Convert a bedGraph file to bigWig format.

bedGraphToBigWig in.bedGraph chrom.sizes out.bw
where in.bedGraph is a four column file in the format:
      <chrom> <start> <end> <value>
and chrom.sizes is a two-column file/URL: <chromosome name> <size in bases>
and out.bw is the output indexed big wig file.
If the assembly <db> is hosted by UCSC, chrom.sizes can be a URL like
  http://hgdownload.soe.ucsc.edu/goldenPath/<db>/bigZips/<db>.chrom.sizes
or you may use the script fetchChromSizes to download the chrom.sizes file.
If not hosted by UCSC, a chrom.sizes file can be generated by running
twoBitInfo on the assembly .2bit file.
The input bedGraph file must be sorted, use the unix sort command:
  sort -k1,1 -k2,2n unsorted.bedGraph > sorted.bedGraph
options:
   -blockSize=N - Number of items to bundle in r-tree.  Default 256
   -itemsPerSlot=N - Number of data points bundled at lowest level. Default 1024
   -unc - If set, do not use compression.

3. bigWigAverageOverBed v2 - Compute average score of big wig over each bed, which may have introns.

bigWigAverageOverBed in.bw in.bed out.tab
The output columns are:
   name - name field from bed, which should be unique
   size - size of bed (sum of exon sizes
   covered - # bases within exons covered by bigWig
   sum - sum of values over all bases covered
   mean0 - average over bases with non-covered bases counting as zeroes
   mean - average over just covered bases
Options:
   -stats=stats.ra - Output a collection of overall statistics to stat.ra file
   -bedOut=out.bed - Make output bed that is echo of input bed but with mean column appended
   -sampleAroundCenter=N - Take sample at region N bases wide centered around bed item, rather
                     than the usual sample in the bed item.
   -minMax - include two additional columns containing the min and max observed in the area.

4. bigWigCorrelate - Correlate bigWig files, optionally only on target regions.

bigWigCorrelate a.bigWig b.bigWig
or
   bigWigCorrelate listOfFiles
options:
   -restrict=restrict.bigBed - restrict correlation to parts covered by this file
   -threshold=N.N - clip values to this threshold
   -rootNames - if set just report the root (minus directory and suffix) of file
                names when using listOfFiles
   -ignoreMissing - if set do not correlate where either side is missing data
                Normally missing data is treated as zeros

5. bigWigMerge v2 - Merge together multiple bigWigs into a single output bedGraph.

You'll have to run bedGraphToBigWig to make the output bigWig.
The signal values are just added together to merge them

bigWigMerge in1.bw in2.bw .. inN.bw out.bedGraph
options:
   -threshold=0.N - don't output values at or below this threshold. Default is 0.0
   -adjust=0.N - add adjustment to each value
   -clip=NNN.N - values higher than this are clipped to this value
   -inList - input file are lists of file names of bigWigs
   -max - merged value is maximum from input files rather than sum

6. fetchChromSizes

fetchChromSizes <db> > <db>.chrom.sizes
    used to fetch chrom.sizes information from UCSC for the given <db>
<db> - name of UCSC database, e.g.: hg38, hg18, mm9, etc ...

This script expects to find one of the following commands:
    wget, mysql, or ftp in order to fetch information from UCSC.
Route the output to the file <db>.chrom.sizes as indicated above.
This data is available at the URL:
  http://hgdownload.soe.ucsc.edu/goldenPath/<db>/bigZips/<db>.chrom.sizes

Example:   fetchChromSizes hg38 > hg38.chrom.sizes
?著作權歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市虑稼,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌势木,老刑警劉巖蛛倦,帶你破解...
    沈念sama閱讀 206,839評論 6 482
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場離奇詭異啦桌,居然都是意外死亡溯壶,警方通過查閱死者的電腦和手機,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,543評論 2 382
  • 文/潘曉璐 我一進店門甫男,熙熙樓的掌柜王于貴愁眉苦臉地迎上來且改,“玉大人,你說我怎么就攤上這事板驳∮瞩耍” “怎么了?”我有些...
    開封第一講書人閱讀 153,116評論 0 344
  • 文/不壞的土叔 我叫張陵若治,是天一觀的道長慨蓝。 經(jīng)常有香客問我感混,道長,這世上最難降的妖魔是什么礼烈? 我笑而不...
    開封第一講書人閱讀 55,371評論 1 279
  • 正文 為了忘掉前任弧满,我火速辦了婚禮,結果婚禮上此熬,老公的妹妹穿的比我還像新娘谱秽。我一直安慰自己,他們只是感情好摹迷,可當我...
    茶點故事閱讀 64,384評論 5 374
  • 文/花漫 我一把揭開白布疟赊。 她就那樣靜靜地躺著,像睡著了一般峡碉。 火紅的嫁衣襯著肌膚如雪近哟。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 49,111評論 1 285
  • 那天鲫寄,我揣著相機與錄音吉执,去河邊找鬼。 笑死地来,一個胖子當著我的面吹牛戳玫,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播未斑,決...
    沈念sama閱讀 38,416評論 3 400
  • 文/蒼蘭香墨 我猛地睜開眼咕宿,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了蜡秽?” 一聲冷哼從身側(cè)響起府阀,我...
    開封第一講書人閱讀 37,053評論 0 259
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎芽突,沒想到半個月后试浙,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 43,558評論 1 300
  • 正文 獨居荒郊野嶺守林人離奇死亡寞蚌,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 36,007評論 2 325
  • 正文 我和宋清朗相戀三年田巴,在試婚紗的時候發(fā)現(xiàn)自己被綠了。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片挟秤。...
    茶點故事閱讀 38,117評論 1 334
  • 序言:一個原本活蹦亂跳的男人離奇死亡壹哺,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出煞聪,到底是詐尸還是另有隱情斗躏,我是刑警寧澤逝慧,帶...
    沈念sama閱讀 33,756評論 4 324
  • 正文 年R本政府宣布昔脯,位于F島的核電站啄糙,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏云稚。R本人自食惡果不足惜隧饼,卻給世界環(huán)境...
    茶點故事閱讀 39,324評論 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望静陈。 院中可真熱鬧燕雁,春花似錦、人聲如沸鲸拥。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,315評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽刑赶。三九已至捏浊,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間撞叨,已是汗流浹背金踪。 一陣腳步聲響...
    開封第一講書人閱讀 31,539評論 1 262
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留牵敷,地道東北人胡岔。 一個月前我還...
    沈念sama閱讀 45,578評論 2 355
  • 正文 我出身青樓,卻偏偏與公主長得像枷餐,于是被迫代替她去往敵國和親靶瘸。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 42,877評論 2 345

推薦閱讀更多精彩內(nèi)容