dotPlotly 比較兩個(gè)不同組裝的效果

可視化dotPlotly工具蜕青,是利用R語言寫的可視化Minimap2/mummer比對輸出的結(jié)果文件覆履。

安裝



# 1 激活環(huán)境后輸入R乖仇,進(jìn)入R環(huán)境

# 2 指定安裝包的來源(這些步驟都是用的百度到的
options(repos=structure(c(CRAN="https://mirrors.tuna.tsinghua.edu.cn/CRAN/")))
# 3 R安裝模塊
install.packages("optparse", "plotly","ggplot2")

4  在Linux服務(wù)器終端從github上下載dotPlotly的兩個(gè)R腳本
git clone https://github.com/piyixing/dotPlotly.git

5 根據(jù)參數(shù)進(jìn)行運(yùn)行腳本,先使用例子的paf文件也行
/public/home/fengting/task/5.12ragtag/7.22/out/test/dotPlotly/pafCoordsDotPlotly.R -i test.paf -o 1 -s -t -m 500 -q 500000 -k 7 -l

minimap2 生成paf文件命令,-ax 生成sam格式,-cx生成paf格式鹿寨。

minimap2 -cx asm5 HHZ.fasta rice_MH63.fa > test.paf

參數(shù)說明:

Usage: pafCoordsDotPlotly.R -i alignments.coords -o out [options]
Options:
-i INPUT, --input=INPUT # 輸入文件,minimap paf文件
coords file from mummer program 'show.coords' [default NULL]
-o OUTPUT, --output=OUTPUT #輸出文件前綴
output filename prefix [default out]
-v, --verbose
Print out all parameter settings [default]
-q MIN-QUERY-LENGTH, --min-query-length=MIN-QUERY-LENGTH # 最短的查詢序列長度薪夕,用于過濾
filter queries with total alignments less than cutoff X bp [default 4e+05]
-m MIN-ALIGNMENT-LENGTH, --min-alignment-length=MIN-ALIGNMENT-LENGTH #最短的比對長度
filter alignments less than cutoff X bp [default 10000]
-p PLOT-SIZE, --plot-size=PLOT-SIZE # 繪圖圖片大小
plot size X by X inches [default 15]
-l, --show-horizontal-lines # 顯示水平線脚草,默認(rèn)關(guān)上,建議打開
turn on horizontal lines on plot for separating scaffolds [default FALSE]
-k NUMBER-REF-CHROMOSOMES, --number-ref-chromosomes=NUMBER-REF-CHROMOSOMES # 選擇繪圖時(shí)染色體數(shù)目寥殖,默認(rèn)繪制所有序列
number of sorted reference chromosomes to keep [default all chromosmes]
-s, --identity # -s 不同比對一致性顯示不同的顏色玩讳,默認(rèn)關(guān)上,建議打開
turn on color alignments by % identity [default FALSE]
-t, --identity-on-target # 僅對目標(biāo)對比時(shí)計(jì)算一致性嚼贡,可以理解為比對很好時(shí)才會計(jì)算一致性熏纯,如果加上-s 參數(shù),那就是只有比對很好的序列會顯示顏色粤策,其他都是黑白
turn on calculation of % identity for on-target alignments only [default FALSE]
-x, --interactive-plot-off #是否打開互動展示樟澜,默認(rèn)開,建議關(guān)上,打開報(bào)錯(cuò)
turn off production of interactive plotly [default TRUE]
-r REFERENCE-IDS, --reference-ids=REFERENCE-IDS #要保留的引用id列表(以逗號分隔)
comma-separated list of reference IDs to keep [default NULL]
-h, --help
Show this help message and exit

結(jié)果展示

經(jīng)過安裝dotPlotly秩贰,就得到了共線性圖霹俺,會自動排序,調(diào)整正負(fù)方向毒费,和mummer軟件出來的圖有所差別丙唧,更美觀,同時(shí)也會損失正負(fù)鏈信息觅玻。

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末想际,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子溪厘,更是在濱河造成了極大的恐慌胡本,老刑警劉巖,帶你破解...
    沈念sama閱讀 217,084評論 6 503
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件畸悬,死亡現(xiàn)場離奇詭異侧甫,居然都是意外死亡,警方通過查閱死者的電腦和手機(jī)蹋宦,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,623評論 3 392
  • 文/潘曉璐 我一進(jìn)店門披粟,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人妆档,你說我怎么就攤上這事僻爽。” “怎么了贾惦?”我有些...
    開封第一講書人閱讀 163,450評論 0 353
  • 文/不壞的土叔 我叫張陵,是天一觀的道長敦捧。 經(jīng)常有香客問我须板,道長,這世上最難降的妖魔是什么兢卵? 我笑而不...
    開封第一講書人閱讀 58,322評論 1 293
  • 正文 為了忘掉前任习瑰,我火速辦了婚禮,結(jié)果婚禮上秽荤,老公的妹妹穿的比我還像新娘甜奄。我一直安慰自己,他們只是感情好窃款,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,370評論 6 390
  • 文/花漫 我一把揭開白布课兄。 她就那樣靜靜地躺著,像睡著了一般晨继。 火紅的嫁衣襯著肌膚如雪烟阐。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 51,274評論 1 300
  • 那天,我揣著相機(jī)與錄音蜒茄,去河邊找鬼唉擂。 笑死,一個(gè)胖子當(dāng)著我的面吹牛檀葛,可吹牛的內(nèi)容都是我干的玩祟。 我是一名探鬼主播,決...
    沈念sama閱讀 40,126評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼屿聋,長吁一口氣:“原來是場噩夢啊……” “哼空扎!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起胜臊,我...
    開封第一講書人閱讀 38,980評論 0 275
  • 序言:老撾萬榮一對情侶失蹤勺卢,失蹤者是張志新(化名)和其女友劉穎,沒想到半個(gè)月后象对,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體黑忱,經(jīng)...
    沈念sama閱讀 45,414評論 1 313
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,599評論 3 334
  • 正文 我和宋清朗相戀三年勒魔,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了甫煞。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點(diǎn)故事閱讀 39,773評論 1 348
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡冠绢,死狀恐怖抚吠,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情弟胀,我是刑警寧澤楷力,帶...
    沈念sama閱讀 35,470評論 5 344
  • 正文 年R本政府宣布,位于F島的核電站孵户,受9級特大地震影響萧朝,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜夏哭,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,080評論 3 327
  • 文/蒙蒙 一检柬、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧竖配,春花似錦何址、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,713評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至龄减,卻和暖如春项钮,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 32,852評論 1 269
  • 我被黑心中介騙來泰國打工烁巫, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留署隘,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 47,865評論 2 370
  • 正文 我出身青樓亚隙,卻偏偏與公主長得像磁餐,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個(gè)殘疾皇子阿弃,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,689評論 2 354

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