Skip to content

kristyhoran/maskrc-svg

 
 

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

16 Commits
 
 
 
 
 
 
 
 

Repository files navigation

maskrc-svg

Masks recombinant regions in an alignment based on ClonalFrameML or Gubbins output
Option to draw SVG of recombinant regions

Authors

Jason Kwong (@kwongjc) :: kwongj
Torsten Seemann (@torstenseemann) :: tseemann

Dependencies

Usage

$ maskrc-svg.py -h

usage: 
  maskrc-svg.py --aln FASTA --out OUTFILE [--gubbins] <PREFIX>

Mask recombination from ClonalFrameML/Gubbins output and draw SVG of recombinant regions

positional arguments:
  PREFIX               prefix used for CFML/Gubbins input files (required)

optional arguments:
  -h, --help           show this help message and exit
  --gubbins            parse as Gubbins instead of ClonalFrameML
  --aln FASTA          multiFASTA alignment used as input for CFML (required)
  --out OUTFILE        output file for masked alignment (default="maskrc.aln")
  --symbol CHAR        symbol to use for masking (default="?")
  --regions FILE       output recombinant regions to file
  --svg FILE           draw SVG output of recombinant regions and save as specified file
  --svgsingle          draw a single SVG for each isolate in input, with isolate label as filename
  --svgsize WIDExHIGH  specify width and height of SVG in pixels (default="800x600")
  --svgorder FILE      specify file containing list of taxa (1 per line) in desired order
  --svgcolour COLOUR   specify colour of recombination regions in HEX format (default=black)
  --consensus          add consensus row of recombination hotspots
  --version            show program's version number and exit

Requires:

  • Output from ClonalFrameML or Gubbins
  • MultiFASTA genome alignment used as input

Options:

  • Specify output file using --out OUTFILE
  • Specify symbol to use in the alignment for masking --symbol ?
  • Save tab-separated file of recombinant region coordinates --regions FILE
  • Draw SVG of recombinant regions --svg FILE
  • Specify size of SVG in pixels eg. 800x600 --svgsize WIDTHxHEIGHT
  • Specify desired order of taxa in SVG --svgorder FILE
  • Specify colour to show extant recombination (ancestral recombination is shown in grey) --svgcolour COLOUR
  • Add consensus row of recombination hotspots to SVG --consensus

Bugs

Please submit via the GitHub issues page: https://github.com/kwongj/maskrc-svg/issues

Software Licence

GPLv3: https://github.com/kwongj/maskrc-svg/blob/master/LICENSE

Other links

About

Masks recombination as detected by ClonalFrameML or Gubbins and draws an SVG

Resources

License

Stars

Watchers

Forks

Packages

No packages published

Languages

  • Python 100.0%