Skip to content
This repository has been archived by the owner on Mar 15, 2024. It is now read-only.
/ GSSplayground Public archive

Lightweight single-html-file-based Genome Segments playground for Visualize genome features cluster(gene arrow map or other features), add synteny among genome fragments or add crosslink among features, add short(PE/MP)/long reads(pacbio or nanopore) mapping or snpindel in vcf(not support complex sv yet), support all CIGAR of sam alignment, dire…

Notifications You must be signed in to change notification settings

orangeSi/GSSplayground

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

36 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

ClustersPloter

Usage:

%  git clone http://github.com/orangeSi/ClustersPloter.git
%  cd  ClustersPloter/example/
%  cat work.sh # output svg and html
#  or go to https://github.com/orangeSi/ClustersPloter/wiki or https://clustersploter.readthedocs.io/en/latest/# 



gene cluster image

Dependence

#method1: install by docker

#method2: manually install, tested only in Centos yet
* perl >=v5.10.1
optional:
  samtools >=1.7 # if you give sam or bam file, add this to linux PATH 
  Imager::Font  # perl package for cacualte text height/width by font size 

Features:

  • plot gene clusters of many samples, one track means one sample, one track contain more than one fragments. one fragment contain gene cluster. you can defined every gene or feature(rotation,color,label,order depth,font size) in clusters. And add crossing link for any pair of genes.
  • every track mean one sample , one sample can has more than one fragments. you can defind the feature color/lable font size/label color/label rotaion in feature.color.label.conf
  • you can draw crosslink or sysnteny among features of different tracks or same track
  • illunimate pair-end or mate-pair reads or pacbio/nanopore long reads mapping with varation(indel, sv) in bam, support snpindel in vcf file, not support complex sv yet

Update:

  • support gff/bed in main.conf(instead of track.list)
  • support GC plot
  • support directly modify almost features attribution in html with mouse-click (color, opacity, postion, font_size)
  • support minimap2/blast/Mummer and commom alignment format for synteny
  • support both reverse fragment and reverse alignment hit in synteny

2018-12-28:

  • add html format output, so you can pan or zoom the svg in html

2018-12-21:

  • add synteny between genomes, only support paf(like minimap2 ouput) format yet. example is out8 in example/work.sh

2018-12-20:

  • support indel of cigar or indel in vcf now, but not support snp yet. test is out3 in example/work.sh

2018-12-10:

2018-11-12:

  • rewrite feature.crossing.link.conf so that defined color and opactity or order or anchor position of every pair of links.
  • add tracks_reorder to plot tracks by new order, otherwise you must adjust --list file

2018-11-07:

  • feature_shaple now support circle_point/rect/arrow, add feature_shift_y, as example/out10.svg. combine with feature_height_ratio or feature_shift_y, try to plot line or scatter or heatmap or histogram is possible but hard for much data points, I will try do this by embed svg in svg, the embedded svg is producted by python or R tool(such as https://github.com/nspies/genomeview or https://github.com/svviz/svviz), and the linux convert not support embed svg, but cairosvg support it (should use the absoluted embedded svg path)!

2018-10-29:

  • redesign the the main.conf and feature.crossing.link.conf to make more freely to reset every feature, reset every link color and opacity one by one, and reset order depth of every feature or crosslink or track, reset feature height one by one as in out2.svg
  • add rect for feature shape, now have arrow and rect, you can use arrow and rect at the same time
  • remove the legend border line

Tips:

$samtools calmd or SamFixCigar.jar of jvarkit will convert M to X/= in cigar of sam file, if do this, will see snp of reads in read_mapping

ToDo:

     add gff/bed to reads_mapping      add reads alignment start/end for reads_mapping, only have target alignment start/end yet      support Hi-C data/reads mapping      support julia/python api to call clustersploter   2018-12-28:
    1. add gtf fomrat, not only gff format
    2. add blat_spl and mummer4 fomrat, not only blast_m8 or paf format
  2018-11-09:
    1. sort by feature,so same feature of different tracks can align centre
    2. add ratio scale in the bottom to recognise the approximate length of every feature

cat example/work.sh

gene cluster image


gene cluster image


gene cluster image


gene cluster image


gene cluster image


gene cluster image


gene cluster image


gene cluster image


gene cluster image


gene cluster image

gene cluster image

gene cluster image


contact:
    QQ: 1522051171
    mail: ilikeorangeapple@gmail.com

License

License: MIT

About

Lightweight single-html-file-based Genome Segments playground for Visualize genome features cluster(gene arrow map or other features), add synteny among genome fragments or add crosslink among features, add short(PE/MP)/long reads(pacbio or nanopore) mapping or snpindel in vcf(not support complex sv yet), support all CIGAR of sam alignment, dire…

Topics

Resources

Stars

Watchers

Forks

Packages

No packages published

Languages