mlst icon indicating copy to clipboard operation
mlst copied to clipboard

Add option to report all matched schemes

Open hmontenegro opened this issue 4 years ago • 2 comments

It would be nice to have an option to report all matched schemes above the score, instead of reporting just the best one.

Case in point: for an isolate genome which is not truly an isolate, the output is:

contigs.fa	abaumannii_2	15	Pas_cpn60(6)	Pas_fusA(6)	Pas_gltA(8)	Pas_pyrG(2)	Pas_recA(3)	Pas_rplB(5)	Pas_rpoB(4)

However, --debug shows:

$VAR1 = [
          'abaumannii_2',
          15,
          '6/6/8/2/3/5/4',
          100
        ];
$VAR2 = [
          'bcereus',
          1280,
          '33/8/13/19/8/17/201',
          100
        ];
$VAR3 = [
          '-',
          '-',
          '-/-/-/-/-/-/-',
          0
        ];

Showing multiple matching schemes would allow quick screening of contaminated isolates. The default for this option would be to report the best match, so no changes to current behaviour.

hmontenegro avatar May 19 '20 16:05 hmontenegro