Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,101,5140GT23.2% 143.5 / 21.4 56intergenic (+84/‑176)sodA/kdgTsuperoxide dismutase (Mn)/2‑dehydro‑3‑deoxy‑D‑gluconate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (9/34);  new base T (13/0);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 2.63e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAG  >  NC_000913/4101379‑4101645
                                                                                                                                       |                                                                                                                                     
gAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGAt                                                                                                                                   <  2:812687/140‑1 (MQ=255)
gAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGAt                                                                                                                                   >  2:412796/1‑140 (MQ=255)
     tggtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAAtt                                                                                                                              >  2:711393/1‑140 (MQ=255)
         gAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCaa                                                                                                                          <  2:828663/140‑1 (MQ=255)
            cTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCTGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTaaaaa                                                                                                                       >  2:1208423/1‑139 (MQ=255)
               ggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGTTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGGTTAATTTCACaaaaaa                                                                                                                    >  1:855108/1‑140 (MQ=255)
               ggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTCAAAATaaa                                                                                                                    >  1:1719264/1‑140 (MQ=255)
                gACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGTTTTTTTTTGATCAATTTCTAAATaaaa                                                                                                                   >  2:102782/1‑140 (MQ=255)
                    aaGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGGTTAATTTCAAAATaaaaaaaa                                                                                                               >  2:671801/1‑140 (MQ=255)
                    aaGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTGATCAATTTTAAAATAAAAACaa                                                                                                               >  2:1690375/1‑140 (MQ=255)
                          cGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTTATCAAATTCAAAAAAAAAACAATAATcc                                                                                                         >  1:1157559/1‑140 (MQ=255)
                             cACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGat                                                                                                      <  1:919122/140‑1 (MQ=255)
                             cACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGat                                                                                                      <  1:929409/140‑1 (MQ=255)
                             cACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGat                                                                                                      <  2:1719264/140‑1 (MQ=255)
                               cGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGAATCGataa                                                                                                    >  2:400539/1‑140 (MQ=255)
                               cGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAAAAAAAAAAAAGATCCGataa                                                                                                    >  2:1589743/1‑140 (MQ=255)
                               cGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGataa                                                                                                    >  1:505720/1‑140 (MQ=255)
                                       ggcgAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGAGGGTTTTTTTTGATCAATTTTAAAAAAAAAAAAAAGATCCGataaaaaaaaaa                                                                                            >  2:1687933/1‑140 (MQ=255)
                                             aaaaTAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGt                                                                                      <  1:1690375/140‑1 (MQ=255)
                                             aaaaTAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGt                                                                                      <  1:1589743/140‑1 (MQ=255)
                                                ataatCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTc                                                                                   <  1:1414500/140‑1 (MQ=255)
                                                  aatCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCaa                                                                                 <  2:556131/140‑1 (MQ=255)
                                                      aTTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGa                                                                             <  2:1289464/140‑1 (MQ=255)
                                                        ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtg                                                                           <  1:616658/140‑1 (MQ=255)
                                                         tGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCaa                                                                                 >  1:983415/1‑133 (MQ=255)
                                                         tGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCaa                                                                                 <  2:983415/133‑1 (MQ=255)
                                                          gccgccTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGGTCAATTTTAAAATAAAAAAAATGATCCGATAAAAATAAAAAAGCGGTTCAATTTAtgtg                                                                         >  2:296018/1‑140 (MQ=255)
                                                                tgctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTg                                                                  <  2:1646078/140‑1 (MQ=255)
                                                                  ctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAAt                                                                                <  2:1661016/125‑1 (MQ=255)
                                                                  ctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAAt                                                                                >  1:1661016/1‑125 (MQ=255)
                                                                  ctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTTAAATAAAAAAAAAGATCCGATAAAAAAAAAAAAACGTTTTAATTTATGTGGttttttt                                                                >  2:1130578/1‑137 (MQ=255)
                                                                    gcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCa                                                             <  1:1565311/140‑1 (MQ=255)
                                                                               tataGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGAt                                                  <  2:180450/140‑1 (MQ=255)
                                                                               tataGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGAt                                                  >  1:180450/1‑140 (MQ=255)
                                                                               tataGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAAATATCCGATAAAAAAAAAAAAGCGGTTTAATTGATGTGGTTTTTATCACTTttatttat                                                  >  2:1514113/1‑140 (MQ=255)
                                                                                 taGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTa                                                <  2:505720/140‑1 (MQ=255)
                                                                                    gctgCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaat                                             <  1:1208423/137‑1 (MQ=255)
                                                                                      cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAACAATGATCCGATAAAAATAAAAAAGCGGTTCAATTGATGTGGGTTTGATCACTTTTTTTGATTAaatga                                           >  1:272903/1‑140 (MQ=255)
                                                                                      cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTAGTGTTTTTTTGATGAATTTTAAAATAAAAACAAATATTCGCTAAAAATAAAACAGGGTTTCAAATTGTGGGGGTTTGGTCGCTTTTATTATTTAaaaga                                           >  2:1510555/1‑140 (MQ=255)
                                                                                         atatCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTTATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGATTCTATTTATGTGGGTTTTTTCACTTTTATTGATTAaaagaaag                                        >  2:212886/1‑138 (MQ=255)
                                                                                         atatCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAATTTCAAAATAAAAAAAATGATCCGATAAAAATAAAACAGCGGTTCAATTGATGTGGGTTTGATCACTTTTTTTTATTAaatgaatg                                        >  1:913194/1‑140 (MQ=255)
                                                                                          tatCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATAAATTTTAAAATAAAAACAAAGAATCGATAAAAATAAAAAAACGCTTTAATTGTGAGAGTTTTTATACCTTTTTATTGTTTAaataatat                                        >  2:851576/1‑135 (MQ=255)
                                                                                            tCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCt                                     <  1:105487/140‑1 (MQ=255)
                                                                                            tCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCt                                     <  1:1223034/140‑1 (MQ=255)
                                                                                              aGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAATTTTAAAAAAAAAAAAAAGATCCGATAAAAAAAAAAAAGCGGTTCAATTGATGTGGTTTTTATCACTTTTTTTGATTAAATGAATGTCTAt                                   >  2:1674680/1‑140 (MQ=255)
                                                                                                     aaaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGtt                            <  2:272903/140‑1 (MQ=255)
                                                                                                     aaaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGtt                            <  2:1293257/140‑1 (MQ=255)
                                                                                                     aaaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGtt                            <  1:1547948/140‑1 (MQ=255)
                                                                                                      aaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGttt                           <  2:855108/140‑1 (MQ=255)
                                                                                                        aTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTcc                         <  1:671801/140‑1 (MQ=255)
                                                                                                          gAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTCAAAATAAAAAAAATGATCCGATAAAAATAAAACAACGTTTCAATTGATGTGGGTTTTATCACTTTTTTTTATTAAATGAAAGTCTTTCTTCGTTTTCAt                       >  2:1547948/1‑140 (MQ=255)
                                                                                                          gAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCAt                       <  2:1157559/140‑1 (MQ=255)
                                                                                                               aTCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAaca                  <  1:1514113/140‑1 (MQ=255)
                                                                                                                 cGCCAACGGCGGTGGGTTTTTTGTGATCAATTTTAAAATAAAAACAATGATCCGATAAAAATAAAAAAGCGTTTCAATTGATGTGGGTTTGATCACTTTTATTTATTAAATGAATGTCTATCTTTGTTTCCATCAACACt                >  1:1745372/1‑140 (MQ=255)
                                                                                                                  gCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTg               <  1:412796/140‑1 (MQ=255)
                                                                                                                     aaCGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATg            <  1:678413/140‑1 (MQ=255)
                                                                                                                     aaCGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATg            <  2:58834/140‑1 (MQ=255)
                                                                                                                     aaCGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATg            <  1:400539/140‑1 (MQ=255)
                                                                                                                        ggcggTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTc         <  2:942259/140‑1 (MQ=255)
                                                                                                                        ggcggTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTc         <  1:296018/140‑1 (MQ=255)
                                                                                                                          cggTGGTTTTTTTGTGATCAATTTCAAAATAAAAAAAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCa       >  2:313834/1‑140 (MQ=255)
                                                                                                                               gtTTTTTTGTGATCAAATTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAg  <  1:1086980/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
GAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAG  >  NC_000913/4101379‑4101645

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