Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273349,9990TG52.8% 27.2 / 50.5 36V528G (GTG→GGG) prpEpropionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (2/15);  new base G (19/0);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 2.44e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-07
Rejected as consensus: Frequency below/above cutoff threshold.

TGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAC  >  NZ_CP009273/349916‑350084
                                                                                   |                                                                                     
tGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGa                                                                                 >  2:866761/1‑90 (MQ=255)
tGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGa                                                                                 >  1:900007/1‑90 (MQ=255)
tGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGa                                                                                 >  1:516213/1‑90 (MQ=255)
 gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAt                                                                                <  1:484113/90‑1 (MQ=255)
 gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAAAt                                                                                >  2:503155/1‑90 (MQ=255)
          cGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGaaa                                                                       <  2:237679/90‑1 (MQ=255)
                     gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGGGCTGGAAGGGGAgggggg                                                            >  1:181015/1‑89 (MQ=255)
                     gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGGGCGGAAAGGGCAgggggg                                                            >  1:830616/1‑89 (MQ=255)
                        cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggcggt                                                         <  2:919275/90‑1 (MQ=255)
                                  gagagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTc                                               <  1:278100/90‑1 (MQ=255)
                                    gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGGTTGTCAt                                             >  2:278100/1‑90 (MQ=255)
                                      gTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAAGGGAGGGGGCGGGGGCGGTTGTCAttt                                           >  2:825571/1‑89 (MQ=255)
                                          ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGaa                                       <  2:287247/90‑1 (MQ=255)
                                            ccAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTTGGGGGAAAGAAGGGCTGGAAGGGCAGGGGGCGGTGGGGGTTGGCATTTCGCAag                                     >  2:303660/1‑90 (MQ=255)
                                            ccAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAag                                     <  1:833873/90‑1 (MQ=255)
                                            ccAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGGGCTTAAAGGGCAGGGGGCGGTGGGGTTTGTTATTGCGAAag                                     >  2:219112/1‑90 (MQ=255)
                                              aGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGAAGGGCTGAAAGGGCAGGGGGGGGGGGGGGTTGTCATTCCGGAaggg                                   >  2:981431/1‑88 (MQ=255)
                                              aGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGGGCTGGAAGGGGAGGGGGGGGGGGGGGTTGTCATTCCGAAagag                                   >  1:911912/1‑90 (MQ=255)
                                              aGCCATCCGGGCGTTGCCGAAGTGCGGGGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGGGGGGGCGGTTGTCATTTCGAAagaga                                  >  1:649642/1‑90 (MQ=255)
                                                    ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc                             <  2:396403/90‑1 (MQ=255)
                                                    ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc                             <  1:164690/90‑1 (MQ=255)
                                                    ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc                             <  2:756458/90‑1 (MQ=255)
                                                    ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAAATTCGCTGAAAGGGGAGGGGGCGGGGGCGTTTGTCATTCCCAAAGAGAGCGAc                             >  2:640839/1‑90 (MQ=255)
                                                            tGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGGGGGGGGGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGa                     >  2:164690/1‑90 (MQ=255)
                                                            tGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGa                     >  2:341107/1‑90 (MQ=255)
                                                            tGCCGAAGTGGCGGGGGGTGGGGGGGAAAATGCGCTGGAAGGGGAGGGGGGGGGGGGGGTTGTTATTCCGAAAGAGAGCGACAGTCTggg                     >  1:874685/1‑89 (MQ=255)
                                                                gAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGAt                 <  1:827520/90‑1 (MQ=255)
                                                                gAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGGAGGGGGGGGGGGGGTTTGTCATTCCGAAAGAGAGGGACAGTCTGGAAGAt                 >  1:228399/1‑90 (MQ=255)
                                                                gAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGGGGTTGTCATTCCGAAAAAGAGGGACAGTCTGGGAGAt                 >  1:173474/1‑90 (MQ=255)
                                                                  aGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCg               <  2:911912/90‑1 (MQ=255)
                                                                              tGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCa   <  2:33969/90‑1 (MQ=255)
                                                                               ggggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAc  <  1:113457/90‑1 (MQ=255)
                                                                               ggggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAc  <  2:874685/90‑1 (MQ=255)
                                                                               ggggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAc  >  2:591365/1‑90 (MQ=255)
                                                                                gggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTc                                               >  2:584337/1‑44 (MQ=255)
                                                                                gggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTc                                               <  1:584337/44‑1 (MQ=255)
                                                                                   |                                                                                     
TGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAC  >  NZ_CP009273/349916‑350084

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