phac-nml / staramr

Scans genome contigs against the ResFinder, PlasmidFinder, and PointFinder databases.
Apache License 2.0
111 stars 25 forks source link

A Full Set of pbp5 Ampicillin Mutations are on the Threshold of the Default Percent Identity #208

Open emarinier opened 3 months ago

emarinier commented 3 months ago

When identifying pbp5 PointFinder mutations using the included test data, the full set of 20 pbp5 mutations results in a BLAST percent identity of 98.233%. The default StarAMR BLAST percent identity threshold is 98%. After manually adding 2 codon mutations to the test file, the percent identity drops to 97.987% and no mutations (individual point or single complex) are reported.

Considering the mutations we're looking for in this region necessarily amount to a lower percent identity, I wonder if there is something we should do by default to accommodate the lower percent identity.