1.JupiterPlot下載地址：https://github.com/JustinChu/JupiterPlot
一種基于Circos的工具，與參考基因組相比，可以可視化基因組組裝的一致性。依賴軟件circos，minimap2 ，samtools。
2.JupiterPlot安裝
直接github 下載，下載方式很多，下載解壓即可。
接著需要將minimap2加入當(dāng)前環(huán)境：
export PATH=/share/nas2/genome/biosoft/minimap2/current/:$PATH
3.參考命令如下：
分析前，建議提取兩個物種染色體序列進行比較畫圖。
export PATH=/share/nas2/genome/biosoft/minimap2/current/:$PATH
export PATH=/share/nas2/genome/biosoft/circos/circos-0.69/bin/:$PATH
reference=Chr.ref.fa
scaffolds=Chr.query.fa
/share/nas1/user/01.software/JupiterPlot/JupiterPlot-master/jupiter name=test ref=$reference fa=$scaffolds labels=both
主要使用的參數(shù)說明：
name 輸出文件前綴
ref 參考基因組，寫的時候在前面
fa 自己的物種基因組，寫的時候再后面
labels=both 表示顯示ref 和fa兩者的染色體id
4.使用中可能存在的問題：
如果使用過程中發(fā)現(xiàn)id不全，需要在日志文件nohup.out中一步步debug，手動調(diào)整整理好。例如需要加上-g 10000，ref 的染色體id才會抓全。
nohup.out中:
perl /share/nas1/user/01.software/JupiterPlot/JupiterPlot-master/bin/generateConf.pl -n 200 -m -1 -r /share/nas1/user/01.software/JupiterPlot/JupiterPlot-master/config/rawConf.conf -p test -s test_scaffolds.fa -b test-agp.bed -a test.agp -k test_reference.karyotype -g 100000
Optional commands:
sam= #Specify an existing alignment of scafftigs to if they already exist (naming convention that fatoagp.pl produces must be consistent)
######General Parameters
t=4 #number of threads to use for minimap2
>######Karyotype options
m=100000 #only use genomic reference chromosomes larger than this value
ng=75 #use largest scaffolds that are equal to 75% of the genome
maxScaff=-1 #Instead of ng filter by this number of scaffolds
i=0 #increment for colouring chromosomes (HSV colour shift by setting 0-360), when set to >360 it generates random colours
g=1 #minimum gap size in reference to render
gScaff=100000 #minimum gap size in scaffolds to render
labels = ref #Shows reference chromosome name "ref", scaffolds "scaf" or "both".
######Link options
maxGap=100000 #maximum alignment gap allowed to consider a region contiguous
minBundleSize=50000 #minimum size of a contiguous region to render
MAPQ=50 #maximum mapping quality allowed when filtering
linkAlpha=5 #alpha of links 1 = 17% , 2 = 33%, 3 = 50%, 4 = 67% and 5 = 83%.
結(jié)果如下：
image.png
5.結(jié)果調(diào)整
1.同時顯示兩者的染色體id需要注釋掉配置文件
#注釋掉test.conf中的label_format = eval( var(chr) =~ /scaf(\d+)$/ ? "": var(label) )
#radius = 0.80r #染色體id太多，調(diào)小半徑
