處理基因組數(shù)據(jù),很多時(shí)候我們會(huì)覺得直接看序列文件不夠直觀,如果繪圖的話,把n多G把數(shù)據(jù)用畫圖出來不僅費(fèi)勁,就算操作也不方便。因此我們可以用UCSC開發(fā)出的genome browser,可以直接把數(shù)據(jù)信息寫成track,連上genome browser 上查看,它還支持安裝到本地服務(wù)器上(genome browser in box ,簡(jiǎn)稱GBIB),genome browser 支持的格式有bedGraph, GTF, PSL, BED, bigBed, WIG, bigGenePred, bigMaf, bigChain, bigPsl, bigWig, BAM, CRAM, VCF, MAF, BED detail, Personal Genome SNP, broadPeak, narrowPeak, and microarray (BED15),GFF和GTF文件必須tab分隔。 廢話少說,直接入門。本文主要講SAM,BAM,WIG,bigWig,VCF,BED文件上傳及使用。
一、格式的前期處理
1.1 WIG 和 bigWig
WIG 文件格式,有兩種可選的格式,variableStep和fixedStep。variableStep用于區(qū)間變化的,fixedStep用于區(qū)間固定的。
variableStep WIG文件以variableStep 開頭,chrom染色體,可選參數(shù)span(默認(rèn)span=1),指定每一行的位置區(qū)間,比如2,區(qū)間就是chromStart~chromStart+2。chromStart染色體位置,dataValue染色體位置上的值。
1 variableStep chrom=chrN2 [span=windowSize]3 chromStartA dataValueA4 chromStartB dataValueB5 ... etc ... ... etc ...
延伸閱讀
- ssh框架 2016-09-30
- 阿里移動(dòng)安全 [無線安全]玩轉(zhuǎn)無線電——不安全的藍(lán)牙鎖 2017-07-26
- 消息隊(duì)列NetMQ 原理分析4-Socket、Session、Option和Pipe 2024-03-26
- Selective Search for Object Recognition 論文筆記【圖片目標(biāo)分割】 2017-07-26
- 詞向量-LRWE模型-更好地識(shí)別反義詞同義詞 2017-07-26
- 從棧不平衡問題 理解 calling convention 2017-07-26
- php imagemagick 處理 圖片剪切、壓縮、合并、插入文本、背景色透明 2017-07-26
- Swift實(shí)現(xiàn)JSON轉(zhuǎn)Model - HandyJSON使用講解 2017-07-26
- 阿里移動(dòng)安全 Android端惡意鎖屏勒索應(yīng)用分析 2017-07-26
- 集合結(jié)合數(shù)據(jù)結(jié)構(gòu)來看看(二) 2017-07-26