Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273349,9930TG52.3% 47.6 / 56.0 44V526G (GTT→GGT) prpEpropionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (7/14);  new base G (23/0);  total (30/14)
Fisher's exact test for biased strand distribution p-value = 1.01e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-06
Rejected as consensus: Frequency below/above cutoff threshold.

GTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATG  >  NZ_CP009273/349905‑350076
                                                                                          |                                                                                   
gTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGtt                                                                                    <  1:384416/90‑1 (MQ=255)
  aaTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTgg                                                                                  >  2:116135/1‑90 (MQ=255)
      aaCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTgggggg                                                                              >  2:88914/1‑88 (MQ=255)
      aaCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTgggggg                                                                              >  1:417207/1‑88 (MQ=255)
        cGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTgggggg                                                                              >  1:150485/1‑88 (MQ=255)
        cGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTggggggga                                                                            >  1:170986/1‑86 (MQ=255)
          ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTTGGGGGaaaa                                                                          >  1:234344/1‑89 (MQ=255)
          ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAg                                                                          >  2:79262/1‑90 (MQ=255)
            gCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCGGAAGTGGCGGTGGGTGGGGGGGAAGaa                                                                        >  1:423385/1‑89 (MQ=255)
            gCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAt                                                                        <  1:358085/90‑1 (MQ=255)
            gCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAt                                                                        >  1:65983/1‑90 (MQ=255)
            gCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAt                                                                        >  1:85580/1‑90 (MQ=255)
                    tCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAATGGGCTgga                                                                >  1:300808/1‑88 (MQ=255)
                                gcgcgAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGCGCTGGAAGGGGAgggggg                                                    >  2:18816/1‑89 (MQ=255)
                                gcgcgAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGCGCTGAAAGGGGAgggggg                                                    >  2:126343/1‑89 (MQ=255)
                                    gagaTTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggcggtg                                                <  1:116135/90‑1 (MQ=255)
                                         tGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGtt                                           <  1:128876/90‑1 (MQ=255)
                                          gA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGGAGGTGGCGGGGGCGgtt                                          >  1:391153/1‑90 (MQ=255)
                                           a‑‑aGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggt                                                       >  1:174659/1‑76 (MQ=255)
                                           a‑‑aGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggt                                                       <  2:174659/76‑1 (MQ=255)
                                                 gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAATGCGCTGAAAGGGGAGGGGGGGGGGGGGGTTGTCAt                                     >  1:394461/1‑90 (MQ=255)
                                                 gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGGAAGGGGAGGGGGGGGGGGGGGTTGTCAt                                     >  2:125873/1‑90 (MQ=255)
                                                 gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTTGGGGGGAAGGTGGGCCGGAAGGGGAGGGGGGGGGGGGGGTTgggct                                     >  2:201283/1‑86 (MQ=255)
                                                 gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGCGCTGAAAGGGGAGGGGGGGGGGGGGGTTGTCAt                                     >  2:96808/1‑90 (MQ=255)
                                                 gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAAGCGCTGGAAGGGGAGGGGGGGGGGGGGGTTGTCAt                                     >  2:427104/1‑90 (MQ=255)
                                                 gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAAAAGCGCTGAAAGGGCGGGGGGGGGGGGGGGTTGtgct                                     >  2:446470/1‑87 (MQ=255)
                                                 gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGAGGGTGGGGGGGAAGATGGGCTGAAAAGGGAGGGGGGGGGGGGGGTTTTCAt                                     >  2:316228/1‑90 (MQ=255)
                                                   gTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCGGGGGGGGGGGGCGGTTTTCATTc                                   >  1:350831/1‑90 (MQ=255)
                                                   gTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGGAAAATGCGCTGAAAGGGGAGGGGGGGGGGGGGGTTGTCATTc                                   >  2:28561/1‑90 (MQ=255)
                                                       ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGaa                               <  1:18331/90‑1 (MQ=255)
                                                       ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGaa                               <  1:99641/90‑1 (MQ=255)
                                                       ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGaa                               <  1:308672/90‑1 (MQ=255)
                                                           aGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGCGGGGGCGGTTGTCATTCCGAAagag                           >  1:166983/1‑90 (MQ=255)
                                                           aGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAAGGGCCGGAAGGGGAGGGGGGGGGGGCGGTTGTTGTTTCGAAagag                           >  2:285598/1‑90 (MQ=255)
                                                               aTCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAAATGCGCTGCAAGGGGAGGGGGCGGGGGGGGTTGTCATTTCGAAAGAGAGCg                       >  1:462129/1‑90 (MQ=255)
                                                                tCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGTGAAAGATGCGCTGAAAAGGGAGGGGGGGGGGGGGGTTGTCCTTTCGCAAAAAAGCGa                      >  1:342777/1‑90 (MQ=255)
                                                                 ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc                     <  1:208745/90‑1 (MQ=255)
                                                                 ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc                     <  1:245441/90‑1 (MQ=255)
                                                                 ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc                     <  2:242035/90‑1 (MQ=255)
                                                                 ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc                     <  1:271700/90‑1 (MQ=255)
                                                                 ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc                     <  2:394461/90‑1 (MQ=255)
                                                                   gggCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAg                   <  2:279210/90‑1 (MQ=255)
                                                                           ccGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAg           <  1:316228/90‑1 (MQ=255)
                                                                                    cGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCgtggtg  >  1:264671/1‑90 (MQ=255)
                                                                                    cGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGGGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGAtg  >  2:248511/1‑90 (MQ=255)
                                                                                          |                                                                                   
GTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATG  >  NZ_CP009273/349905‑350076

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