Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,089,5710AG56.2% 0.5 / 14.4 16I83V (ATT→GTT) hisFimidazole glycerol phosphate synthase subunit HisF
Reads supporting (aligned to +/- strand):  ref base A (1/6);  new base G (9/0);  total (10/6)
Fisher's exact test for biased strand distribution p-value = 8.74e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCC  >  NZ_CP009273/2089484‑2089653
                                                                                       |                                                                                  
gCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGtt                                                                                  >  2:452167/1‑90 (MQ=255)
  ttCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGttta                                                                                >  2:386554/1‑88 (MQ=255)
  ttCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTaa                                                                                >  1:388044/1‑90 (MQ=255)
  ttCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTaa                                                                                >  2:79453/1‑90 (MQ=255)
       gCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGtctc                                                                           <  2:209110/90‑1 (MQ=255)
           tGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAATCTCTggg                                                                       >  1:221950/1‑89 (MQ=255)
                       aGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGaa                                                           <  1:164500/90‑1 (MQ=255)
                               gCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGtctc                                                                           >  2:181983/1‑66 (MQ=255)
                               gCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGtctc                                                                           <  1:181983/66‑1 (MQ=255)
                                               ggCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTAAGTCTCTGGAAGATGGGGCGCAAATTCTTTCCTTTTGCGCgggaaa                                   >  1:418185/1‑86 (MQ=255)
                                                         aTCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATc                         <  1:74488/90‑1 (MQ=255)
                                                                ttCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAAATTCCATCAACtctc                  >  1:299769/1‑90 (MQ=255)
                                                                ttCCGTTTTGTGTGGCGGGGGGGGTTTAGTCTCTGGGAGAAGCCCCGAAAAATTTTTCCTTTGGGGCGGGTAAAATTTCCATCAACtctc                  >  2:325842/1‑90 (MQ=255)
                                                                            tgGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAg      <  2:388044/90‑1 (MQ=255)
                                                                              gCGGGGGGGGTTAAGTCTCTGGAAGATGCCCCGAAAATTTTTTTCTTTGGCGCGGGTAAAATTTCCATCAACTCTCCTGCGCTGGCAGAc    >  2:133480/1‑90 (MQ=255)
                                                                                gggTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGGTGGCAGAccc  <  1:386554/90‑1 (MQ=255)
                                                                                       |                                                                                  
GCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCC  >  NZ_CP009273/2089484‑2089653

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: