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

TTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACC  >  NZ_CP009273/2089486‑2089652
                                                                                     |                                                                                  
ttCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTaa                                                                                >  1:461280/1‑90 (MQ=255)
        aTGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAGTCTCTgg                                                                        >  1:714222/1‑90 (MQ=255)
         tGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGa                                                                       >  1:736990/1‑90 (MQ=255)
                                gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGtc                                                                             >  1:486006/1‑61 (MQ=255)
                                gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGtc                                                                             <  2:486006/61‑1 (MQ=255)
                                                     tGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCa                          <  2:173468/90‑1 (MQ=255)
                                                     tGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCa                          <  2:461280/90‑1 (MQ=255)
                                                      gATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAAGCTTTGGGGGAAGGCCCGGAAAATTTTTTCTTTTGGGCGGGTAAAAATTTCAt                         >  2:699095/1‑90 (MQ=255)
                                                       aTCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGATGCCGCGAAAAATCTTTTCTTTTGCGCGGGTAAAAATTCCATc                        >  1:542879/1‑90 (MQ=255)
                                                          gACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAAc                     <  2:624035/90‑1 (MQ=255)
                                                           aCATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTTGCGCGGATAAAATTTCCATCAAct                    >  1:347657/1‑90 (MQ=255)
                                                              ttCCGTTTTGTGTGGCGGGGGGGGTTAAATCTCTGGAAGATGCCGCGAAAATTTTTTTCTTTTGGGGGGGTAAAAATTTCATCAACtctc                 >  2:594313/1‑90 (MQ=255)
                                                              ttCCGTTTTGTGTGGCGGGGGGGGTTAAATCTCTGGAAGATGCCGCGAAAATTTTTTTCTTTTGGGCGGGTAAAATTTTCATCAACtctc                 >  2:709882/1‑90 (MQ=255)
                                                               tCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTcc                <  1:594313/90‑1 (MQ=255)
                                                                  gTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAAAAGCCCCGAAAATTTTTTTCTTTTGGGGGGGTAAAAATTTCATCAACCCTTCtct             >  1:658280/1‑88 (MQ=255)
                                                                    tttGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTgcgc           <  1:828404/90‑1 (MQ=255)
                                                                            gCGGGGGGGGTTAAGTCTCTGGGAGATGCCGCGAAAATTTTTTTCTTTGGGGCGGGTAAAAATTTCATCAACTCTCCTGCGCTGGGAGAc   >  2:8090/1‑90 (MQ=255)
                                                                             cGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGAcc  <  2:736990/90‑1 (MQ=255)
                                                                                     |                                                                                  
TTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACC  >  NZ_CP009273/2089486‑2089652

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