Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273349,9990TG57.1% 21.8 / 46.0 35V528G (GTG→GGG) prpEpropionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (3/12);  new base G (20/0);  total (23/12)
Fisher's exact test for biased strand distribution p-value = 5.45e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.55e-06
Rejected as consensus: Frequency below/above cutoff threshold.

TTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACT  >  NZ_CP009273/349915‑350085
                                                                                    |                                                                                       
ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAg                                                                                    >  2:93283/1‑90 (MQ=255)
ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAg                                                                                    >  2:579476/1‑90 (MQ=255)
  gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAt                                                                                  <  1:444013/90‑1 (MQ=255)
  gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAt                                                                                  >  2:117345/1‑90 (MQ=255)
    cGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATgc                                                                                <  1:494035/90‑1 (MQ=255)
       gCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCt                                                                             <  1:231386/90‑1 (MQ=255)
                ggggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAAAGGGCTGGAAAgggg                                                                    >  2:708604/1‑89 (MQ=255)
                     cgcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAggggg                                                              >  2:52945/1‑90 (MQ=255)
                      gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAATGGGCTGAAAGGGGAgggggg                                                             >  1:88263/1‑89 (MQ=255)
                      gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAAAGGGCTGGAAGGGCggggggg                                                             >  1:405019/1‑89 (MQ=255)
                       cgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAAAGCGCTGGAAAGGGAggggggg                                                            >  1:582705/1‑88 (MQ=255)
                            gaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAATGCGCTGAAAAGGGAgggggggggggg                                                       >  1:175114/1‑89 (MQ=255)
                                     gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGGGGGGGCGGTTGTCAt                                              >  2:148001/1‑90 (MQ=255)
                                     gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAAATGCGCTGAAAGGGCAGGGGGGGGGGGGGGTTGTCAt                                              >  2:485176/1‑90 (MQ=255)
                                     gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGAAAGATGCGCTGAAAAGGCAGGGGGCGGGGGCGGTTGGCAt                                              >  2:317722/1‑90 (MQ=255)
                                     gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAATGCGCCGGAAGGGGGGGGGGGGGGGGCGGTTGTTAt                                              >  1:603696/1‑90 (MQ=255)
                                     gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAAATGCGCTGAAAGGGGAGGGGGCGGGGGGGGTTGTCAt                                              >  1:35167/1‑90 (MQ=255)
                                       gTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGGGCTGAAAAGGGAGGGGGCGGGGGGGGTTGTCATTc                                            >  2:565937/1‑90 (MQ=255)
                                          tctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAATGCGCTGAAAAGGGCGGGGGGGGGGGCGGTTTGCATTTCGa                                         >  2:319356/1‑90 (MQ=255)
                                               aGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGGAGGTGGGGGGGGCGTTTGTCATTTCGAAagag                                    >  2:91912/1‑90 (MQ=255)
                                                   aTCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGGGTTTGTCATTCCGAAAGAGAGCg                                >  2:10804/1‑90 (MQ=255)
                                                   aTCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAATGGGCTGAAAGGGGAGGGGGGGGGGGGGGTTTTCCTTTCGAAaaagaagg                                >  1:69330/1‑87 (MQ=255)
                                                    tCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGTGGAAAAAGCGCGGAAAAGGCAGGGGGGGGGGGCGGTTTTGTTTCCGGAagagagggg                               >  1:636465/1‑87 (MQ=255)
                                                     ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc                              <  2:289791/90‑1 (MQ=255)
                                                     ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc                              <  1:91912/90‑1 (MQ=255)
                                                     ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc                              <  2:568341/90‑1 (MQ=255)
                                                       gggCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAg                            <  1:479955/90‑1 (MQ=255)
                                                       gggCGTTGCCGAAGTGGCGGTGGTTGGGATGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAg                            <  1:285601/90‑1 (MQ=255)
                                                          cGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCt                         >  2:26504/1‑90 (MQ=255)
                                                            ttGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTgg                       <  2:223364/90‑1 (MQ=255)
                                                             tGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAAGGGAGGGGGGGGGGGGGTTTGTCATTCCGAAAGAGAGCGACAGTCTggg                      >  2:528705/1‑89 (MQ=255)
                                                               ccGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAg                    <  2:34626/90‑1 (MQ=255)
                                                                   aGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCg                <  2:476750/90‑1 (MQ=255)
                                                                                ggggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAc   <  1:10804/90‑1 (MQ=255)
                                                                                 gggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACt  >  2:66258/1‑90 (MQ=255)
                                                                                    |                                                                                       
TTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACT  >  NZ_CP009273/349915‑350085

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