Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273349,9990TG52.4% 34.1 / 62.1 42V528G (GTG→GGG) prpEpropionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (1/19);  new base G (22/0);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 4.48e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-07
Rejected as consensus: Frequency below/above cutoff threshold.

AACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAC  >  NZ_CP009273/349911‑350084
                                                                                         |                                                                                     
aaCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTg                                                                                      <  1:553373/90‑1 (MQ=255)
  cGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTTGGGGGGa                                                                                     >  2:44568/1‑90 (MQ=255)
    ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAg                                                                                  >  2:539307/1‑90 (MQ=255)
    ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGGGGGGGGGAAg                                                                                  >  2:487116/1‑90 (MQ=255)
        cGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATgc                                                                              <  1:48500/90‑1 (MQ=255)
          ggCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATgcgc                                                                            <  1:555237/90‑1 (MQ=255)
           gCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAAAGgcgcc                                                                           >  2:302544/1‑89 (MQ=255)
                    ggggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAGGCGGTGAAAagggg                                                                  >  1:385317/1‑88 (MQ=255)
                          gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGCGCTGAAAGGGCAgggggc                                                            >  1:329891/1‑90 (MQ=255)
                          gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGGAgggggg                                                            >  1:56975/1‑89 (MQ=255)
                          gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAggtggg                                                            >  1:93106/1‑89 (MQ=255)
                          gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGAAAGATGCGGTGAAAGGGCTgggggg                                                            >  1:399643/1‑89 (MQ=255)
                          gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGGGCTGAAAGGGCAgggggg                                                            >  1:598869/1‑89 (MQ=255)
                             cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggcggt                                                         <  1:487116/90‑1 (MQ=255)
                                       gagagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTc                                               <  1:434376/90‑1 (MQ=255)
                                         gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAAGGGAGGGGGCGGTGGCGGTTGTCAt                                             >  1:211759/1‑90 (MQ=255)
                                         gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGGAAGATGCGCTGGAAGGGGAGGGGGGGGGGGGGGTTGTTAt                                             >  2:452823/1‑90 (MQ=255)
                                                 ccAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAag                                     <  2:581102/90‑1 (MQ=255)
                                                 ccAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAag                                     <  1:368588/90‑1 (MQ=255)
                                                   aGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAagag                                   <  1:293467/90‑1 (MQ=255)
                                                   aGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAAGGCAGGGGGGGGGGGCGGTTGTCATTCCGAAagag                                   >  1:657867/1‑90 (MQ=255)
                                                   aGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGGAAGATGCGCTGAAAGGGCGGGGGGGGGGGGGGGTTTTCATTCCGAAagag                                   >  1:581102/1‑90 (MQ=255)
                                                   aGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAATGGGCTGGAAGGGGGGGGGGGGGGGGGGGTTGTCATTTCGAAaaag                                   >  2:571057/1‑90 (MQ=255)
                                                   aGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAAAGCGCTGGAAGGGGAGGGGGGGGGGGGGGTTGTTATTTCGGAaaag                                   >  2:368588/1‑90 (MQ=255)
                                                       aTCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGCGCTGAAAGGGGGGGGGGGGGGGGGGGTTTTCATTCCGAAagagaggg                               >  2:364130/1‑88 (MQ=255)
                                                       aTCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGGGGGGGCGTTTGTCATTCCGAAAGAGAGCg                               >  1:78263/1‑90 (MQ=255)
                                                       aTCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAAAAGCGCTGAAAGGGGAGGGGGCGGGGGGGGTTTTGCTTCCGGAAAAGAGCg                               >  2:263613/1‑90 (MQ=255)
                                                       aTCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAAAGCGCTGAAAAGGGGGGGGGGGGGGGGGGTTGTTATTTCGGAAGAGAGCg                               >  2:620587/1‑90 (MQ=255)
                                                        tCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGAAAGATGCGCTGGAAGGGGAGGGGGGGGGGGGGTTTGTCATTTCGAAAGAGAGGGa                              >  2:149544/1‑90 (MQ=255)
                                                         ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc                             <  2:370445/90‑1 (MQ=255)
                                                         ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc                             <  1:571057/90‑1 (MQ=255)
                                                         ccGGGCGGTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc                             <  2:399643/90‑1 (MQ=255)
                                                           gggCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAg                           <  1:641978/90‑1 (MQ=255)
                                                           gggCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAg                           <  1:44568/90‑1 (MQ=255)
                                                           gggCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAg                           <  2:598869/90‑1 (MQ=255)
                                                              cGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCt                        >  1:499577/1‑90 (MQ=255)
                                                                                   ttGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGgcgc    <  2:56975/90‑1 (MQ=255)
                                                                                    tGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCa   <  1:620587/90‑1 (MQ=255)
                                                                                    tGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCa   <  1:539307/90‑1 (MQ=255)
                                                                                     ggggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAc  <  1:149544/90‑1 (MQ=255)
                                                                                      gggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGAt                 <  2:315996/74‑1 (MQ=255)
                                                                                      gggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGAt                 >  1:315996/1‑74 (MQ=255)
                                                                                         |                                                                                     
AACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAC  >  NZ_CP009273/349911‑350084

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