染色體共線性可視化


摘要

基因組組裝質量的評估手段之一就是用現有的參考基因組和組裝的基因組進行共線性分析,級別有全基因組的共線性分析和染色體級別的共線性分析,全基因組級別的共線性可以用circlize進行可視化,現在進行染色體級別的共線性分析可視化。

准備工作

  • NUCMER進行染色體級別的比對
numcer --mum -p chr1 -t 20 chr1_g1.fasta chr1_g2.fasta
  • 准備格式
show-coords -c -l chr1.delta > chr1.delta.coords

可視化

python LINKVIEW.py -t 2 -s --svg2png_dpi 1200 -k k_1.txt --no_dash --svg_width 2500 --chro_axis --chro_axis_density 2 -o chr1 chr1.delta.coords

三個基因組染色體之間的共線性可視化

  • 比對
numcer --mum -p chr1_2 -t 20 chr1_g1.fasta chr1_g2.fasta
numcer --mum -p chr1_3 -t 20 chr1_g1.fasta chr1_g3.fasta
  • 准備格式
show-coords -c -l chr1_2.delta > chr1_2.delta.coords
show-coords -c -l chr1_3.delta > chr1_3.delta.coords
cat chr1_2.delta.coords chr1_3.delta.coords > chr1.coords
  • 可視化
python LINKVIEW.py -t 2 -s --svg2png_dpi 1200 -k k_1.txt --no_dash --svg_width 2500 --chro_axis --chro_axis_density 2 -o chr1 chr1.coords

出圖

image


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM