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

GCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGAT  >  NC_000913/3563416‑3563544
                                                                 |                                                                
gCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggtgg                                                                <  1:2854153/68‑1 (MQ=255)
gCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggtgg                                                                <  1:2259013/68‑1 (MQ=255)
 cGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGa                                                               <  1:3275149/68‑1 (MQ=255)
  gTCGCANAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggggggg                                                              >  1:2029796/1‑66 (MQ=255)
  gTCGCACAAGACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAg                                                              <  1:2961642/68‑1 (MQ=255)
  gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggggggg                                                              >  1:870940/1‑66 (MQ=255)
  gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggggggg                                                              >  1:2947206/1‑66 (MQ=255)
  gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAg                                                              <  1:996744/68‑1 (MQ=255)
  gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAg                                                              <  1:2335659/68‑1 (MQ=255)
  gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAg                                                              <  1:2052370/68‑1 (MQ=255)
     gCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtat                                                           <  1:1185670/68‑1 (MQ=255)
          aaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCGGTGGCTGGTGGAGTATACgca                                                     <  1:2477479/68‑1 (MQ=255)
          aaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgca                                                     <  1:2549853/68‑1 (MQ=255)
            gacaGGGATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcagc                                                   <  1:508567/65‑1 (MQ=255)
            gacaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcagc                                                   <  1:1689504/65‑1 (MQ=255)
            aCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTAAGTCAGTGGCTGGGGGAGTATACgccac                                                   >  1:2960392/1‑65 (MQ=255)
                 gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAgg                                              <  1:1280864/68‑1 (MQ=255)
                 gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAgg                                              <  1:3178202/68‑1 (MQ=255)
                 gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAAGCAGCAGAgg                                              >  1:431989/1‑68 (MQ=255)
                 gCAGGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAgg                                              <  1:2773606/68‑1 (MQ=255)
                  cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGt                                             <  1:2563754/68‑1 (MQ=255)
                  cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAAGCGGCAGAggg                                             >  1:3058669/1‑67 (MQ=255)
                  cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAAGCCGCAGAggg                                             >  1:2590025/1‑67 (MQ=255)
                      tgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATc                                         <  1:687383/68‑1 (MQ=255)
                        gCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCCGAGGGGATCGc                                       >  1:1762201/1‑68 (MQ=255)
                            aaTGCGGATCAACAATCTCGCGTGAGTCAGTGGCTGGGGGGGTAAAACCAGCACCAGGTATAGCTGGc                                   >  1:494792/1‑68 (MQ=255)
                              tGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAGGGTATTGCTggggg                                 >  1:812488/1‑65 (MQ=255)
                               gCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAAGGCGGCGAAGGTTTTCCCTGGcgt                                 >  1:3134344/1‑68 (MQ=255)
                               gCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCCGCAGAGGGTATTGCTGGcggg                                >  1:2243861/1‑66 (MQ=255)
                               gCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAAGCAGCAGAGGGTATTGCCGGggtg                                >  1:1786786/1‑67 (MQ=255)
                                  gATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAGGGTATCGCTGgcggcgga                             >  1:147546/1‑66 (MQ=255)
                                        caaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTaa                       <  1:2427182/68‑1 (MQ=255)
                                        caaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTaa                       <  1:2336989/68‑1 (MQ=255)
                                         aaTCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAGGGTATCGCTGGCGGCGGTAATTaaa                      >  1:1319059/1‑67 (MQ=255)
                                            ctcGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAGGGTATCGCTGGGGGCGTAAAATAAAGTa                   >  1:1727269/1‑68 (MQ=255)
                                             tcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTaa                  <  1:2603260/68‑1 (MQ=255)
                                             tcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTaa                  <  1:3232363/68‑1 (MQ=255)
                                              cGTGTGAGTCAGTGGCTGGGGGGGGATAAGCAGCAGAGGGTATCGCTTGGGGCGGAAATTaaagaaag                 >  1:3305001/1‑68 (MQ=255)
                                              cGTGTGAGTCAGTGGCTGGGGGGGGATAAGCAGCAGAGGGTATCGCTGGGGGCGTAAATTAAAGTAAg                 >  1:1471683/1‑68 (MQ=255)
                                               gtgtGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAgg                <  1:2594010/68‑1 (MQ=255)
                                                tgtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGt               <  1:1081104/68‑1 (MQ=255)
                                                tgtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGAGTCGTAAATTAACGTAAGGt               <  1:1965496/68‑1 (MQ=255)
                                                       cAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAGTTAACGTAAGGTGATCAgg        <  1:54610/68‑1 (MQ=255)
                                                       cAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAgg        <  1:596891/68‑1 (MQ=255)
                                                            gCTGGTGGAGTATGCGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGa   <  1:157422/68‑1 (MQ=255)
                                                            gCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGa   <  1:3161723/68‑1 (MQ=255)
                                                            gCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGa   <  1:815666/68‑1 (MQ=255)
                                                             cTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTGAGAt  <  1:1728453/68‑1 (MQ=255)
                                                             cTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGAt  <  1:409638/68‑1 (MQ=255)
                                                                 |                                                                
GCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGAT  >  NC_000913/3563416‑3563544

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