Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,563,4810TG29.0% 59.7 / 17.1 31V490G (GTG→GGG) glpDaerobic glycerol 3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/22);  new base G (8/1);  total (8/23)
Fisher's exact test for biased strand distribution p-value = 1.14e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGAT  >  NC_000913/3563415‑3563544
                                                                  |                                                               
ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggtg                                                                >  1:834898/1‑68 (MQ=255)
  cGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGa                                                              <  1:451783/68‑1 (MQ=255)
   gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggggggg                                                             >  1:777679/1‑66 (MQ=255)
   gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAg                                                             <  1:3845294/68‑1 (MQ=255)
   gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAg                                                             <  1:1720260/68‑1 (MQ=255)
   gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAg                                                             <  1:2957972/68‑1 (MQ=255)
   gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAg                                                             <  1:2071567/68‑1 (MQ=255)
   gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAg                                                             <  1:2357675/68‑1 (MQ=255)
    tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGt                                                            <  1:206588/68‑1 (MQ=255)
      gCACAAAACAAGGCATGTGGCTGAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtat                                                          <  1:257829/68‑1 (MQ=255)
      gCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtat                                                          <  1:3603566/68‑1 (MQ=255)
        acaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATAc                                                        <  1:3721907/68‑1 (MQ=255)
        acaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGCGGGTATAc                                                        >  1:3653408/1‑68 (MQ=255)
           aaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACgca                                                     <  1:17498/68‑1 (MQ=255)
             aCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAagcagc                                                   >  1:2157886/1‑68 (MQ=255)
                  ggaTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAgg                                              <  1:2558355/66‑1 (MQ=255)
                  gCATTTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAgg                                              <  1:3532862/68‑1 (MQ=255)
                  gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAgg                                              <  1:2859480/68‑1 (MQ=255)
                  gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAgg                                              <  1:2869813/68‑1 (MQ=255)
                  gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATAAGCAGACGAgg                                              >  1:1014946/1‑68 (MQ=255)
                  gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGGTACGGAGCcgaag                                              >  1:859005/1‑66 (MQ=255)
                  gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGGAgggg                                              >  1:84252/1‑68 (MQ=255)
                  gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGGAGAgg                                              >  1:3469130/1‑68 (MQ=255)
                  gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAAGCAGCAgggg                                              >  1:2277983/1‑68 (MQ=255)
                   cATCTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGt                                             <  1:2223326/68‑1 (MQ=255)
                       tgGCTAAATGCGGATCAAGAGTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATc                                         <  1:1494435/68‑1 (MQ=255)
                              aTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGcg                                  <  1:494194/68‑1 (MQ=255)
                               tGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAGGGTATCGCTTGcgg                                 >  1:2110792/1‑67 (MQ=255)
                               tGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGAATAGGGAGCAGAAGGGATCGCTggggg                                 >  1:3888679/1‑65 (MQ=255)
                                gCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAGGGTATTGCTGGggtc                                >  1:3940607/1‑68 (MQ=255)
                                    aTCAACAATCTCGTGTGAGTCAGTGGCTTGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTaa                            <  1:3681471/68‑1 (MQ=255)
                                         caaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTaa                       <  1:2097417/68‑1 (MQ=255)
                                             ctcGTGTGAGTCAGTGGCTGGGGGGGGATACGGATGAGCGGTGTTTCCTGGGGGCGTCAATTAACggg                   >  1:4152877/1‑66 (MQ=255)
                                               cGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGCGGGTATTGCTGGCGTCGGGAATTAACGTaaa                 >  1:2266662/1‑67 (MQ=255)
                                               cGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGACGGTTTTGCTGGCGTCGGTAATTAAAGTAAg                 >  1:280217/1‑68 (MQ=255)
                                                 tgtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGt               <  1:1111130/68‑1 (MQ=255)
                                                 tgtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGt               <  1:4006617/68‑1 (MQ=255)
                                                 tgtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGt               <  1:3534502/68‑1 (MQ=255)
                                                        taGTGTCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAgg        <  1:469815/67‑1 (MQ=255)
                                                             gCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGa   <  1:4140651/68‑1 (MQ=255)
                                                             gCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGa   <  1:1207930/68‑1 (MQ=255)
                                                              cTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGAt  <  1:3457317/68‑1 (MQ=255)
                                                                  |                                                               
GGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGAT  >  NC_000913/3563415‑3563544

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑