Skip to content

Latest commit

 

History

History
77 lines (68 loc) · 4.05 KB

README.md

File metadata and controls

77 lines (68 loc) · 4.05 KB

ClustersPloter

Usage:

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

Dependence

samtools >=1.7 # add this to linux PATH 
perl >=v5.10.1
sort

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
  • illunimate pair-end or mate-pair reads or pacbio/nanopore long reads mapping with varation(indel, sv) in bam, support snpindel in vcf file

Update:

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. test is example/out7.sh

2018-12-20:

  • support indel of cigar or indel in vcf now, but not support snp yet. test is example/out11.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

ToDo:

  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
    3. embed svg in svg to make heatmap or scatter or line or histogram more easy to product, like example/out.svg*svg or <foreignObject>
 

sh example/out11.sh # output out11.svg and convert out11.svg to out11.png

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: [email protected]