jvarkit

WesCnvSvg

Last commit

SVG visualization of bam DEPTH for multiple regions

Usage

This program is now part of the main jvarkit tool. See jvarkit for compiling.

Usage: java -jar dist/jvarkit.jar wescnvsvg  [options] Files

Usage: wescnvsvg [options] Files
  Options:
    -css, --css
      custom svg css stylesheet
    -x, --extend
      Extending interval. The following syntaxes are supported: 1000; 1kb; 
      1,000; 30%(shrink); 150% (extend); 0.5 (shrink); 1.5 (extend)
      Default: 0
    -height, --height
      Sample Track height
      Default: 100
    -h, --help
      print help and exit
    --helpFormat
      What kind of help. One of [usage,markdown,xml].
    -u, --url, --hyperlink
      creates a hyperlink an area is 'clicked'. creates a hyperlink when 
      'click' in an area. The URL must contains __CHROM__, __START__ and 
      __END__ that will be replaced by their values. Predefined values are 
      'hg19','hg38','igv'. IGV : 
      "http://localhost:60151/goto?locus=__CHROM__%3A__START__-__END__" , 
      UCSC: "http://genome.ucsc.edu/cgi-bin/hgTracks?org=Human&db=hg19&position=__CHROM__%3A__START__-__END__"
      Default: <empty string>
    -B, --bed, -b, --capture, -r, -rgn, --region, --interval
      A source of intervals. The following suffixes are recognized: vcf, 
      vcf.gz bed, bed.gz, gtf, gff, gff.gz, gtf.gz.Otherwise it could be an 
      empty string (no interval) or a list of plain interval separated by '[ 
      \t\n;,]' 
      Default: (unspecified)
    -Q, --mapq
      Min mapping quality
      Default: 1
    -o, --output
      Output file. Optional . Default: stdout
  * -R, --ref, --reference
      Indexed fasta Reference file. This file must be indexed with samtools 
      faidx and with picard/gatk CreateSequenceDictionary or samtools dict
    -smooth, --smooth
      how to smooth data
      Default: AVERAGE
      Possible Values: [AVERAGE, MEDIAN]
    --title
      document title
      Default: WesCnvSvg
    --vcf
      plot VCF data
    --version
      print version and exit
    -w, --width
      Page width
      Default: 1000
    -D
      other parameters. '-Dkey=value'. Undocumented.
      Syntax: -Dkey=value
      Default: {}

Keywords

Creation Date

20180726

Source code

https://github.com/lindenb/jvarkit/tree/master/src/main/java/com/github/lindenb/jvarkit/tools/bam2svg/WesCnvSvg.java

Unit Tests

https://github.com/lindenb/jvarkit/tree/master/src/test/java/com/github/lindenb/jvarkit/tools/bam2svg/WesCnvSvgTest.java

Contribute

License

The project is licensed under the MIT license.

Citing

Should you cite wescnvsvg ? https://github.com/mr-c/shouldacite/blob/master/should-I-cite-this-software.md

The current reference is:

http://dx.doi.org/10.6084/m9.figshare.1425030

Lindenbaum, Pierre (2015): JVarkit: java-based utilities for Bioinformatics. figshare. http://dx.doi.org/10.6084/m9.figshare.1425030

Input

input is a set of bam file or a file with the ‘*.list’ suffix containing the path to the bam files.

Example

$ find dir -name "*.bam"  > bam.list
$ java -jar dist/wescnvsvg.jar -R ref.fasta -B cnv.bed bam.list > cnv.svg 

Screenshots

https://twitter.com/yokofakun/status/1022503372669300738 :

ScreenShot

https://twitter.com/yokofakun/status/1022805656905150464

ScreenShot

https://twitter.com/yokofakun/status/1023953180927963137

ScreenShot

https://twitter.com/yokofakun/status/1024315948847849472

ScreenShot

https://twitter.com/yokofakun/status/1024315948847849472

ScreenShot

https://twitter.com/yokofakun/status/1025330308193779712

ScreenShot

https://twitter.com/yokofakun/status/1040592885786263554

ScreenShot

https://twitter.com/yokofakun/status/1040577235856580608

ScreenShot

https://twitter.com/yokofakun/status/1057625407913111557

ScreenShot

https://twitter.com/yokofakun/status/1180046139502059521

ScreenShot