fix negative Percentage Length of Reference Sequence for rRNA gene va…

…riant model #237
arpcard · Sep 6, 2023 · 4ec630e · 4ec630e
1 parent 0e2565f
commit 4ec630e
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/app/ConvertRGIJsonToTSV.py b/app/ConvertRGIJsonToTSV.py
@@ -223,7 +223,7 @@ def run(self):
 									other_snps = "n/a"
 
 								if rgi_data[hsp][hit]["model_type_id"] in [40295]:
-									percentage_length_reference_sequence = format(((orf_end - orf_start) /\
+									percentage_length_reference_sequence = format((abs(orf_end - orf_start) /\
 										len(rgi_data[hsp][ordered[0]]["dna_sequence_from_broadstreet"]))*100, '.2f')
 								else:
 									percentage_length_reference_sequence = format((len(rgi_data[hsp][ordered[0]]["orf_prot_sequence"]) /\