genomediff-python parses files in the GenomeDiff format generated by the breseq variant caller for haploid microbial organisms.
python setup.py install
Only Python 3.x is tested.
GenomeDiff files are read using GenomeDiff.read(file). The GenomeDiff object contains a metadata dict with
the metadata, as well as mutations, evidence and validation lists—each containing records of that type.
Records can be accessed through this list or by id. GenomeDiff is iterable and iterating will return all records of all types.
>>> from genomediff import *
>>> document = GenomeDiff.read(open('MyDiff.gd', 'r', encoding='utf-8'))
>>> document.metadata
{'GENOME_DIFF': '1.0', 'AUTHOR': ''}
>>> document.mutations[0]
Record('SNP', 1, [191], new_seq='A', seq_id='NC_000913', snp_type='intergenic', position=12346)
>>> document.mutations[0].parent_ids
[191]
>>> document[191]
Record('RA', 191, None, tot_cov='46/42', new_base='A', insert_position=0, ref_base='G', seq_id='NC_000913', quality=252.9, position=12345)
>>> document.mutations[0].parents
[Record('RA', 191, None, tot_cov='46/42', new_base='A', insert_position=0, ref_base='G', seq_id='NC_000913', quality=252.9, position=12345)]
>>> document.write(open('NewDiff.gd', 'w', encoding='utf-8'))