Merge pull request sanger-pathogens#204 from martinghunt/mlst_hetmin_bug

Mlst hetmin bug

ariba/mlst_reporter.py

@@ -47,7 +47,7 @@ def _call_gene(self, gene):
                     depths = [int(x) for x in d['smtls_nts_depth'].split(',')]
                     depths.sort()
                     het_pc = round(100.0 * depths[-1] / sum(depths), 2)
-                    if results['hetmin'] == '.' or results['hetmin'] < het_pc:
+                    if results['hetmin'] == '.' or results['hetmin'] > het_pc:
                         results['hetmin'] = het_pc
         if len(het_data):
             results['hets'] = '.'.join(het_data)

ariba/tests/data/mlst_reporter.het_snps.out.details.tsv

@@ -1,3 +1,3 @@
 gene	allele	cov	pc	ctgs	depth	hetmin	hets
-gene1	1*	100.0	100.0	1	42.2	75.0	30,10.25,10,5
+gene1	1*	100.0	100.0	1	42.2	62.5	30,10.25,10,5
 gene2	2	100.0	100.0	1	40.2	.	.

setup.py

@@ -55,7 +55,7 @@
 setup(
     ext_modules=[minimap_mod, fermilite_mod, vcfcall_mod],
     name='ariba',
-    version='2.10.2',
+    version='2.10.3',
     description='ARIBA: Antibiotic Resistance Identification By Assembly',
     packages = find_packages(),
     package_data={'ariba': ['test_run_data/*']},