Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,299,0170AG29.1% 75.1 / 39.6 55D15G (GAT→GGT) yhbOprotein/nucleic acid deglycase 2
Reads supporting (aligned to +/- strand):  ref base A (19/20);  new base G (11/5);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCCACGGCTCATGCACCAGGCTTAATAGACCATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAAAGGCAAAAAAGGAGAAGCCAGCGTGACCATCGATAAATCCATCGATGAAGT  >  NC_000913/3298880‑3299155
                                                                                                                                         |                                                                                                                                          
tCCACGGCTCATGCACCAGGCTTAATAGACCATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGAtt                                                                                                                                          <  1:985996/140‑1 (MQ=255)
  cACGGCTCATGCACCAGGCTTAATAGACCATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGGTTCa                                                                                                                                        >  2:886290/1‑140 (MQ=255)
       cTCATGCACCAGGCTTAATAGACCATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGGTTCAGAAtt                                                                                                                                   <  2:334804/140‑1 (MQ=255)
          aTGCACCAGGCTTAATAGACCATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGGTTCAGAATTTAc                                                                                                                                >  2:258761/1‑140 (MQ=255)
          aTGCACCAGGCTTAATAGACCATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTAc                                                                                                                                >  1:850092/1‑140 (MQ=255)
            gCACCAGGCTTAATAGACCATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGGTTCAGAATTTACtt                                                                                                                              >  1:776696/1‑140 (MQ=255)
            gCACCAGGCTTAATAGACCATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACtt                                                                                                                              >  1:172104/1‑140 (MQ=255)
            gCACCAGGCTTAATAGACCATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACtt                                                                                                                              <  1:508754/140‑1 (MQ=255)
             cACCAGGCTTAATAGACCATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTc                                                                                                                             <  1:631460/140‑1 (MQ=255)
                 aGGCTTAATAGACCATTAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCAccc                                                                                                                         >  1:977806/1‑140 (MQ=255)
                       aaTAGACCATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGGTTCAGAATTTACTTCACCCGCAGAc                                                                                                                   <  2:588598/140‑1 (MQ=255)
                       aaTAGACCATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGAc                                                                                                                   <  2:803246/140‑1 (MQ=255)
                       aaTAGACCATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGAc                                                                                                                   <  2:406203/140‑1 (MQ=255)
                            aCCATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGGTTCAGAATTTACTTCACCCTCAGACGAGtt                                                                                                              >  2:593001/1‑140 (MQ=255)
                            aCCATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGtt                                                                                                              >  1:897443/1‑140 (MQ=255)
                             ccATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGGTTCAGAATTTACTTCACCCGCAGACGAGTTc                                                                                                             >  2:907983/1‑140 (MQ=255)
                              cATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGGTTCAGAATTTACTTCACCCGCAGACGAGTTcc                                                                                                            >  1:125621/1‑140 (MQ=255)
                                 gAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTa                                                                                                         >  1:987952/1‑140 (MQ=255)
                                     tAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGc                                                                                                     <  1:779000/140‑1 (MQ=255)
                                       aTTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCg                                                                                                   >  1:572440/1‑140 (MQ=255)
                                       aTTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCg                                                                                                   >  2:16024/1‑140 (MQ=255)
                                          aTCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGac                                                                                                <  1:423318/140‑1 (MQ=255)
                                            ccccGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGacac                                                                                              >  1:50664/1‑140 (MQ=255)
                                             cccGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGGTTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACg                                                                                             <  2:221839/140‑1 (MQ=255)
                                              ccGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGa                                                                                            >  1:737615/1‑140 (MQ=255)
                                               cGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGGTTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGaa                                                                                           <  1:886290/140‑1 (MQ=255)
                                                  ttGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTg                                                                                        >  2:900456/1‑140 (MQ=255)
                                                   tgtgGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGa                                                                                       <  1:737539/140‑1 (MQ=255)
                                                     tgGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGAtt                                                                                     >  2:741045/1‑140 (MQ=255)
                                                     tgGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGAtt                                                                                     >  2:83963/1‑140 (MQ=255)
                                                       ggAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGGTTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTAc                                                                                   >  1:647632/1‑140 (MQ=255)
                                                       ggAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTAc                                                                                   <  1:678396/140‑1 (MQ=255)
                                                           aTTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCAtt                                                                               >  2:852688/1‑140 (MQ=255)
                                                                         aTGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGGTTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCggg                                                                 <  1:907983/140‑1 (MQ=255)
                                                                         aTGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCggg                                                                 <  2:572440/140‑1 (MQ=255)
                                                                          tGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGt                                                                <  1:869063/140‑1 (MQ=255)
                                                                              acaacGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTaaaa                                                            >  2:626991/1‑140 (MQ=255)
                                                                               caacGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAAc                                                           >  1:883904/1‑140 (MQ=255)
                                                                                 acGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGGTTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACgg                                                         >  1:285035/1‑140 (MQ=255)
                                                                                  cgcgGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGt                                                        >  1:245421/1‑140 (MQ=255)
                                                                                    cgGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGa                                                      >  1:855262/1‑140 (MQ=255)
                                                                                        ggaAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAAAgg                                                  <  2:850092/140‑1 (MQ=255)
                                                                                                 aGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGGTTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAAAGGCAAAAAAgg                                         >  2:592340/1‑140 (MQ=255)
                                                                                                 aGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAAAGGCAAAAAAgg                                         <  2:121678/140‑1 (MQ=255)
                                                                                                     aGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAAAGGCAAAAAAGGAGaa                                     <  2:987952/140‑1 (MQ=255)
                                                                                                       aaaaTTGCCGTTTTAATCACTGATGAATTTGAGGGTTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAAAGGCAAAAAAGGAGAAGc                                   >  2:573669/1‑140 (MQ=255)
                                                                                                       aaaaTTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAAAGGCAAAAAAGGAGAAGc                                   <  1:628484/140‑1 (MQ=255)
                                                                                                            tGCCGTTTTAATCACTGATGAATTAGAGGACTCAGAATTTAATTCACCCGCAGACGAATTCCGTAAAGAAGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAAAGGCAAAAAAGGAGAAGCCAGCg                              <  2:977806/140‑1 (MQ=255)
                                                                                                                gTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAAAGGCAAAAAAGGAGAAGCCAGCGTGAc                          <  1:900456/140‑1 (MQ=255)
                                                                                                                 ttttAATCACTGACGAATTGGAGGATTCAGAATTTACTTCGCCCGCAGACGAGTACCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAATGGCAAAAAAGGAGAAGCCAGCGGGAcc                         <  2:364633/140‑1 (MQ=255)
                                                                                                                      aTCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAAAGGCAAAAAAGGAGAAGCCAGCGTGACCATCGa                    <  2:172104/140‑1 (MQ=255)
                                                                                                                            gatgaATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAAAGGCAAAAAAGGAGAAGCCAGCGTGACCATCGATAAATc              <  1:741045/140‑1 (MQ=255)
                                                                                                                             atgaATTTGAGGGTTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAAAGGCAAAAAAGGAGAAGCCAGCGTGACCATCGATAAATcc             >  2:208809/1‑140 (MQ=255)
                                                                                                                             atgaATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAAAGGCAAAAAAGGAGAAGCCAGCGTGACCATCGATAAATcc             >  2:530682/1‑140 (MQ=255)
                                                                                                                                 aTTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAAAGGCAAAAAAGGAGAAGCCAGCGTGACCATCGATAAATCCATCg         <  1:83963/140‑1 (MQ=255)
                                                                                                                                  tttGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAAAGGCAAAAAAGGAGAAGCCAGCGTGACCATCGATAAATCCATCGa        >  2:539147/1‑140 (MQ=255)
                                                                                                                                        gATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAAAGGCAAAAAAGGAGAAGCCAGCGTGACCATCGATAAATCCATCGATGAAGt  <  1:530682/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                          
TCCACGGCTCATGCACCAGGCTTAATAGACCATGAGGTAATTATCCCCGATTGTGGGAAATTCGCCTCATCCAATGCAACAACGCGGAGGAAGCATGAGTAAGAAAATTGCCGTTTTAATCACTGATGAATTTGAGGATTCAGAATTTACTTCACCCGCAGACGAGTTCCGTAAAGCCGGACACGAAGTGATTACCATTGAAAAACAAGCGGGTAAAACGGTGAAAGGCAAAAAAGGAGAAGCCAGCGTGACCATCGATAAATCCATCGATGAAGT  >  NC_000913/3298880‑3299155

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑