File: README.md

package info (click to toggle)
python-pymummer 0.9.0-1~bpo8%2B1
  • links: PTS, VCS
  • area: main
  • in suites: jessie-backports
  • size: 280 kB
  • sloc: python: 1,008; sh: 55; makefile: 8
file content (81 lines) | stat: -rw-r--r-- 2,106 bytes parent folder | download | duplicates (4)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
pymummer
========

Python3 wrapper for running MUMmer and parsing the output. 

[![Build Status](https://travis-ci.org/sanger-pathogens/pymummer.svg?branch=master)](https://travis-ci.org/sanger-pathogens/pymummer)

Installation
------------

###Homebrew/LinuxBrew###
```
brew tap homebrew/python
brew install pymummer
```

##Pip
###Pre-requisites###

The MUMmer package must be installed.
Instructions to install MUMmer can be found [here](http://mummer.sourceforge.net/manual/#installation)
		
###Installation###

Install with

    pip3 install pymummer


Usage (for developers)
----------------------

Example showing how pymummer can be used to run nucmer on a fasta file and parse the output file to produce a set of alignment objects:

	from pymummer import coords_file, alignment, nucmer
	...
	runner = nucmer.Runner(reference_file, query_file, results_file) 
	runner.run()
	file_reader = coords_file.reader(results_file)
	alignments = [coord for coord in file_reader if not coord.is_self_hit()] #Remove self hits
	...

###pymummer nucmer class###

Wraps the `nucmer`, `delta-filter`, `show-coords` and `show-snps` commands. 

Arguments:

__ref__			reference file  
__query__			query file  
__outfile__		output file  
__min\_id__		min\_id for delta-filter command (default None)  
__min\_length__	min\_length for delta-filter command (default None)  
__breaklen__		breaklen for nucmer command (nucmer's default is 200)   
__coords\_header__	print header in show-coords output (default True)  
__maxmatch__		maxmatch for nucmer (default False)  
__show\_snps__		run show-snps (default False)  
__snps\_header__ 	print header in show-snps output (default True)  

###pymummer promer class###

[TODO]

###pymummer coords_file class###

Parses the nucmer output and populate an alignment object for each hit in the output

  
###pymummer alignment class###

Check attributes of a hit, swap the reference and query, check if it's a self hit and so on

Contact
-------

Authors: Martin Hunt, Nishadi De Silva

Affiliation: Wellcome Trust Sanger Institute, Hinxton, UK

Email: path-help@sanger.ac.uk