Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,492,5780G.62.5% 29.1 / 82.9 48intergenic (‑53/+9)idnT/idnOL‑idonate/5‑ketogluconate/gluconate transporter/5‑keto‑D‑gluconate 5‑reductase
Reads supporting (aligned to +/- strand):  ref base G (7/11);  new base . (16/14);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGG  >  NC_000913/4492439‑4492714
                                                                                                                                           |                                                                                                                                        
aaCGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAg                                                                                                                                          >  1:834451/1‑140 (MQ=255)
   gAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAgggg                                                                                                                                       <  1:1457781/140‑1 (MQ=255)
    aGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAggggg                                                                                                                                      >  1:764399/1‑140 (MQ=255)
     gggCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAgggggc                                                                                                                                     >  2:464480/1‑139 (MQ=255)
      ggCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAggg                                                                                                                                        >  1:1087480/1‑136 (MQ=255)
      ggCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAggg                                                                                                                                        <  2:1087480/136‑1 (MQ=255)
          aTAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTa                                                                                                                               >  2:989362/1‑140 (MQ=255)
            aaaGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTaaa                                                                                                                             <  1:989362/140‑1 (MQ=255)
             aaGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTaaa                                                                                                                             <  1:543779/140‑1 (MQ=255)
             aaGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTaaa                                                                                                                             <  2:1360367/140‑1 (MQ=255)
             aaGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAAc                                                                                                                            <  2:455602/140‑1 (MQ=255)
                           aaaGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATg                                                                                                               <  1:1129768/140‑1 (MQ=255)
                                  atcatcAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATc                                                                                                       >  2:443234/1‑140 (MQ=255)
                                  atcatcAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATc                                                                                                       >  1:1056921/1‑140 (MQ=255)
                                  atcatcAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATc                                                                                                       >  1:617300/1‑140 (MQ=255)
                                        aGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCACCATCaacaaa                                                                                                 <  1:1259374/140‑1 (MQ=255)
                                               gaagCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAg                                                                                           >  1:867454/1‑140 (MQ=255)
                                                  gcagcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTgg                                                                                       <  2:617300/140‑1 (MQ=255)
                                                    agcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCaacaaa                                                                                                 >  1:1167253/1‑129 (MQ=255)
                                                    agcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCaacaaa                                                                                                 <  1:1254164/129‑1 (MQ=255)
                                                    agcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCaacaaa                                                                                                 >  2:1254164/1‑129 (MQ=255)
                                                    agcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCaacaaa                                                                                                 <  2:1167253/129‑1 (MQ=255)
                                                    agcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGcc                                                                                     >  1:704414/1‑140 (MQ=255)
                                                        cgACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGttt                                                                                 >  2:697110/1‑140 (MQ=255)
                                                              cTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGaa                                                                            <  2:490330/140‑1 (MQ=255)
                                                                 ccGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTgg                                                                                       >  1:798451/1‑125 (MQ=255)
                                                                 ccGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTgg                                                                                       <  2:798451/125‑1 (MQ=255)
                                                                    cAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCaga                                                                      <  2:867454/140‑1 (MQ=255)
                                                                              tAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTgaa                                                            <  1:930713/140‑1 (MQ=255)
                                                                               aaTGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTgaaga                                                          >  1:534080/1‑140 (MQ=255)
                                                                                aTGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTgaagaa                                                         <  1:1001971/140‑1 (MQ=255)
                                                                                     aTAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGaa                                                    <  2:704414/140‑1 (MQ=255)
                                                                                     aTAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGaa                                                    <  2:534080/140‑1 (MQ=255)
                                                                                      tAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGaa                                                    <  1:955408/140‑1 (MQ=255)
                                                                                      tAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAac                                                   >  1:727698/1‑140 (MQ=255)
                                                                                             tATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCgcag                                            >  1:789352/1‑140 (MQ=255)
                                                                                             tATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCgcag                                            <  1:464480/140‑1 (MQ=255)
                                                                                              aTCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCgcag                                            <  1:1062656/140‑1 (MQ=255)
                                                                                                   ttattaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCa                                       >  2:1182539/1‑140 (MQ=255)
                                                                                                   ttattaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCa                                       >  1:556854/1‑140 (MQ=255)
                                                                                                   ttattaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCa                                       >  1:1244774/1‑140 (MQ=255)
                                                                                                   ttattaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCaa                                      >  1:1021667/1‑140 (MQ=255)
                                                                                                     attaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTTAAGAAAGGAACACCGCAGCACCAATc                                    >  2:1404284/1‑140 (MQ=255)
                                                                                                     attaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCg                                               >  2:362341/1‑129 (MQ=255)
                                                                                                     attaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCg                                               <  1:362341/129‑1 (MQ=255)
                                                                                                        aGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGt                                 <  2:789352/140‑1 (MQ=255)
                                                                                                             cAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCt                             >  1:263959/1‑140 (MQ=255)
                                                                                                              aGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTg                            <  2:263959/140‑1 (MQ=255)
                                                                                                               gAAAAGCAGTGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCg                          >  2:667958/1‑140 (MQ=255)
                                                                                                               gAAAAGCAGTGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCg                          <  1:667958/140‑1 (MQ=255)
                                                                                                                aaaaGCAGCGCTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCgg                         >  1:677342/1‑140 (MQ=255)
                                                                                                                       gcgcTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTccc                  <  1:443234/140‑1 (MQ=255)
                                                                                                                         gcTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCa                >  2:643398/1‑140 (MQ=255)
                                                                                                                         gcTTTCGCGCTGCACGTAGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCa                <  2:727698/140‑1 (MQ=255)
                                                                                                                                        agaGGGGGGTGGTTAACCAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAACCAGTTCCTGCGGATCTCCCCAGCGTGCGGCggggg  <  2:417431/138‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                        
AACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGG  >  NC_000913/4492439‑4492714

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