Error in output when a gene contains a parentheses

[nvlachos]

I am reaching out in response to the recent communications with Alison Laufer Halpin ( and Jill Hagey) about trying to use srst2 in our pipeline. We have run into a pair of errors that we have mitigated locally, but would like to pass along to you to include, if desired. These could be just fixes for our specific environment but hopefully they may help others as well. Both occur in srst2.py in the 0.2.0 branch

1. Output error if a '(' is in the description of the gene in the database header (ex:>980__aph(3')__aph(3')-Id_NG_047445.1__06702__aminoglycoside__NCBI [NCBI]aph(3')-Id:1:NG_047445.1:aminoglycoside_O-phosphotransferase_APH(3')-Id|WP_010891085.1|499193545:aminoglycoside;NG_047445.1;aminoglycoside;NCBI)
   Source of the error is

   srst2/scripts/srst2.py

   Line 1502 in 73f885f

   header_string = os.popen(" ".join(["grep",allele,fasta]))

   
   Our fix was to change the line to (note: there are backslashes where the big X's are since the markup took them away)
   command = "grep X""+allele+"X" "+fasta
   header_string = os.popen(command)

2. Rounding Error resulting in the script to completely failing with RuntimeError: floating point number truncated to an integer
   However, this has been resolved since release 0.2.0, but not part of any official release/tag (RuntimeWarning: floating point number trucated to an integer #69). This prevents us from properly documenting tool versions for any validation purposes.

We also would like to bump issue#141 to see if anyone was able to look into it yet.

Thank you so much for your help! Nick