Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,223,4100TG38.5% 32.8 / 13.0 26V409G (GTG→GGG) dldD‑lactate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (3/13);  new base G (10/0);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.08e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TACGAGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCAAAGCGTT  >  NC_000913/2223346‑2223476
                                                                |                                                                  
tACGAGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCtgggggg                                                                 >  1:1869821/1‑67 (MQ=255)
tACGAGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCtgggggg                                                                 >  1:3203132/1‑67 (MQ=255)
tACGAGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCtgggggg                                                                 >  1:1373479/1‑67 (MQ=255)
tACGAGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAAACGTGGCtgggggg                                                                 >  1:1790188/1‑67 (MQ=255)
tACGAGCATCATCTGCTGTTAAAAACGGCGGGCGAGGGCGTGGGCGACGCCCAATCGTGGCtgggggg                                                                 >  1:2658091/1‑67 (MQ=255)
  cGAGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGgtt                                                               >  1:401059/1‑68 (MQ=255)
       atcatcTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGGTTTTTTc                                                          >  1:1248444/1‑68 (MQ=255)
         catcTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCaa                                                        <  1:1350094/68‑1 (MQ=255)
           tcTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGGTTATTTTAAAc                                                      >  1:1694022/1‑68 (MQ=255)
                tgTTAAAAATGGCGGGCGGGGGTGTGGGGGAAGCCGAAAAGGGGGTGGGGGGTTATTTCAAACAggac                                                 >  1:1649003/1‑66 (MQ=255)
                     aaaaGGGCGGCGCATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAgg                                            <  1:2142439/68‑1 (MQ=255)
                        aTGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGGTTATTTTCAACAGGCCGAAggggg                                         >  1:810349/1‑65 (MQ=255)
                           gcgggcgATGGCTTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGAttt                                      <  1:3072998/68‑1 (MQ=255)
                           gcgggcgATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGAttt                                      <  1:2762837/68‑1 (MQ=255)
                                cgATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTg                                 <  1:965474/68‑1 (MQ=255)
                                   tGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCt                              <  1:1455909/68‑1 (MQ=255)
                                   tGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCt                              <  1:25365/68‑1 (MQ=255)
                                   tGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCt                              <  1:2118746/68‑1 (MQ=255)
                                    ggCGTGGGCGAAGCCAAATCGTGGCTGGGGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTg                             >  1:550630/1‑68 (MQ=255)
                                     gCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGt                            <  1:1259371/68‑1 (MQ=255)
                                      cGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTa                           >  1:218910/1‑68 (MQ=255)
                                      cGTGGGCGAAGCCAAATCGTGGCTGGGGGGTTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTa                           >  1:10995/1‑68 (MQ=255)
                                           gCGAAGCCAAATCGTGGCTGGTGGATTATTTCCAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCg                      >  1:1446888/1‑68 (MQ=255)
                                           gCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGGTTTCTTTGTCTGTACGCCg                      >  1:1174819/1‑68 (MQ=255)
                                           gCGAAGCCAAATCGTGGCTGGGGGATTATTTCAAACAGGCCGCAGGCGACTTCTTTGTCTGTACGCCg                      >  1:2309083/1‑68 (MQ=255)
                                           gCGAAGCCAAATCGTGGCTGGGGGATTATTTCAAACAGGCCGAAGGGGATTTCTTTGTCTGTACGCCg                      >  1:3445730/1‑68 (MQ=255)
                                           gCGAAGCCAAATCGTGGCTGGGGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCg                      >  1:2530152/1‑68 (MQ=255)
                                           gCGAAGCCAAATCGTGGCTGGGGGATTATTTCAAAAAGGCCGAAGGCGATTTCTTTGTCTTTACGCCg                      >  1:1107233/1‑68 (MQ=255)
                                           gCGAAGCCAAATCGTGGCTGGAGGATTATTTCAAACAGGCCGGAGGCGATTTCTTTTTCTGTATGCCg                      >  1:1512214/1‑68 (MQ=255)
                                                 ccAAATCGTGGCTGGTGGGTTTTTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCggaggct                >  1:1342980/1‑66 (MQ=255)
                                                 ccAAATCGTGGCTGGTGGGTTATTTTAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGaa                >  1:1007523/1‑68 (MQ=255)
                                                 ccAAATCGTGGCTGGGGGATTATTTCAAAAAGGCCGAAGGCGATTTCTTTGTCTGTACGCCggaggga                >  1:2432592/1‑66 (MQ=255)
                                                  cAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAg               <  1:770747/68‑1 (MQ=255)
                                                  cAAATCGTGGCTGGTGGAGTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAg               <  1:3270361/68‑1 (MQ=255)
                                                     aTCGTGGCTGGGGGATTATTTCAAACAGGCCGAAGGGGGTTTCTTTGTCTGTACGCCGGGGGAAGgca            >  1:1084353/1‑68 (MQ=255)
                                                      tCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGgcag           <  1:2429006/68‑1 (MQ=255)
                                                       cGTGGCTGGGGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGGAGgcagc          >  1:2358232/1‑68 (MQ=255)
                                                          ggCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCaaa       <  1:740035/68‑1 (MQ=255)
                                                               gtggATTATTTCAACCAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCAAAGCGtt  <  1:2025566/68‑1 (MQ=255)
                                                               gtggATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCAAAGCGtt  <  1:1016484/68‑1 (MQ=255)
                                                               gtggATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCAAAGCGtt  <  1:3133024/68‑1 (MQ=255)
                                                               gtggATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCAAAGCGtt  <  1:227776/68‑1 (MQ=255)
                                                               gtggATTATTGCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCAAAGCGtt  <  1:2956022/68‑1 (MQ=255)
                                                                |                                                                  
TACGAGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCAAAGCGTT  >  NC_000913/2223346‑2223476

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