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:

add plot depth and pabio read mapping of bam file, as out3 in example/work.sh , inspired by https://github.com/nspies/genomeview

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