Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,101,5140GT25.9% 131.0 / 21.0 54intergenic (+84/‑176)sodA/kdgTsuperoxide dismutase (Mn)/2‑dehydro‑3‑deoxy‑D‑gluconate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (7/33);  new base T (14/0);  total (21/33)
Fisher's exact test for biased strand distribution p-value = 3.58e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAG  >  NC_000913/4101382‑4101645
                                                                                                                                    |                                                                                                                                   
cGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCaa                                                                                                                              >  2:905533/1‑140 (MQ=255)
      gAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCaa                                                                                                                        <  1:757435/140‑1 (MQ=255)
      gAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCaa                                                                                                                        <  2:330942/140‑1 (MQ=255)
           gggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATaa                                                                                                                   <  2:693936/140‑1 (MQ=255)
            ggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATaaa                                                                                                                  <  2:991855/140‑1 (MQ=255)
              aCGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAAt                                                                                                                     >  1:430110/1‑135 (MQ=255)
              aCGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAAt                                                                                                                     <  2:430110/135‑1 (MQ=255)
                 aaGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTCAAAATAAAAACaa                                                                                                             >  1:744631/1‑140 (MQ=255)
                      gcGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTTGTTAATTTTAaaaaaaaaaaaaaaatc                                                                                                        >  2:296851/1‑135 (MQ=255)
                         gCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGa                                                                                                     >  1:476208/1‑140 (MQ=255)
                          cACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGat                                                                                                    <  2:389531/140‑1 (MQ=255)
                               tttGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGataaaaa                                                                                               <  1:282983/140‑1 (MQ=255)
                                    ggcgAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGataaaaataaaa                                                                                          <  1:623073/140‑1 (MQ=255)
                                        aaaaaaTAATCA‑TTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCCATTTCAAAAAAAAAAAAAATAACCGgtaaaaataaaaaaacg                                                                                     >  2:285852/1‑135 (MQ=255)
                                          aaaaTAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGt                                                                                    <  2:476208/140‑1 (MQ=255)
                                          aaaaTAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGt                                                                                    <  1:470022/140‑1 (MQ=255)
                                                atCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAAt                                                                              <  2:1169449/140‑1 (MQ=255)
                                                atCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAAt                                                                              <  1:658419/140‑1 (MQ=255)
                                                atCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAAt                                                                              >  1:709135/1‑140 (MQ=255)
                                                atCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAAt                                                                              <  2:998805/140‑1 (MQ=255)
                                                     ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTTAAAATAAAAAAAATGATCCGATAAAAATAAAACAGCGGTTCAATTGAtg                                                                         >  1:769605/1‑140 (MQ=255)
                                                     ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtg                                                                         <  2:561449/140‑1 (MQ=255)
                                                     ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtg                                                                         <  2:312703/140‑1 (MQ=255)
                                                      tGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATCCGATAAAAATAAAACAGCGTTTTAATTGAtgt                                                                        >  2:507096/1‑140 (MQ=255)
                                                      tGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAATTTTAAAATAAAAAAAAAGATCCGATAAAAATAAAAAAACGTTTCAATTTAtgt                                                                        >  2:989471/1‑140 (MQ=255)
                                                             agctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTg                                                                 <  1:592976/139‑1 (MQ=255)
                                                                   aaTGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACt                                                           <  1:296851/140‑1 (MQ=255)
                                                                       aGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTa                                                       <  2:182773/140‑1 (MQ=255)
                                                                       aGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAATGATCCGATAAAAATAAAAAAGCGTTTTAATTGATGTGGGTTTGATCACttttt                                                       >  2:990812/1‑139 (MQ=255)
                                                                       aGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAAATTCAAAATAAAAAAAAAAATCCGATAAAAAAAAAAAAACGTTTCAATTTATGTGGTTTTTATCACttttt                                                       >  2:975742/1‑139 (MQ=255)
                                                                        ggCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTAt                                                      <  2:709135/140‑1 (MQ=255)
                                                                          cGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTCAAAATAAAAAAAATGATCCGATAAAAATAAAACAGCGTTTTAATTGATGTGGTTTTTATCACTTTTAttt                                                    >  2:38442/1‑139 (MQ=255)
                                                                          cGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAAAAAAAAAAAAAATCCGATAAAAAAAAAAAAACGTTTTAATTTATGTGGGTTTTATCACTttttttt                                                    >  2:592976/1‑139 (MQ=255)
                                                                          cGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTCAAAATAAAAACAAAGATCCGATAAAAAAAAAAAAGCGGTTCAATTTATGTGGTTTTGATCACTTTTAttt                                                    >  1:8133/1‑139 (MQ=255)
                                                                          cGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAAATTCAAAAAAAAAAAAATGATCCGATAAAAATAAAAAAGCGGTTCAATTGATGTGGGTTTGATTACTttttttt                                                    >  2:1038847/1‑139 (MQ=255)
                                                                           gTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGa                                                   <  1:719719/140‑1 (MQ=255)
                                                                              taGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATCCGATAAAAATAAAACAACGTTTTAATTTATGTGGGTTTGATCACTTTTTTttatta                                                >  2:363339/1‑140 (MQ=255)
                                                                                   cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatga                                           <  1:465080/140‑1 (MQ=255)
                                                                                   cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAATTTTAAAATAAAAAAAAAGATCCGATAAAAAAAAAAAAACGTTTTCATTGATGTGGTTTTTATTACTTTTTTTTATTAaatga                                           >  2:691169/1‑140 (MQ=255)
                                                                                    gCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGAGGTTTTTTTGTGATCAAATTCAAAAAAAAAAAAATGATTCCATAAAAAAAAAAAAACACTTTCATTGATTTGGTTTTGATTACTTTTATTTATTTaatgaa                                          >  2:982931/1‑140 (MQ=255)
                                                                                     cATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatga                                           >  2:444889/1‑138 (MQ=255)
                                                                                     cATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatga                                           <  1:444889/138‑1 (MQ=255)
                                                                                     cATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGGTCAATTTCAAAATAAAAAAAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatgaat                                         >  1:1027286/1‑140 (MQ=255)
                                                                                           aGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGGTTTTATCACTTTTATTTATTAAATGAATGTCTAt                                   >  2:446206/1‑140 (MQ=255)
                                                                                               tAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATCCGATAAAAAAAAAAAAGCGTTTTAATTGATGTGGGTTTTATCACTTTTTTTTATTAAATGAATGTCTATCttt                               >  2:607666/1‑139 (MQ=255)
                                                                                                aaaaaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCg                              <  1:989471/140‑1 (MQ=255)
                                                                                                      tGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCa                        <  1:905533/140‑1 (MQ=255)
                                                                                                       gAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCAt                       <  1:579289/140‑1 (MQ=255)
                                                                                                          ccATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCaa                    <  1:923628/140‑1 (MQ=255)
                                                                                                          ccATCGCCAACGGCGGTGGGTTTTTTTTGATTAATTTTAAAAAAAAAACAAAGATCCGATAAAAATAAAAAAGCGGTTTAATTTATGTGGTTTTTATTACTTTTTTTTATTAAATTAATTTCTATCTTTGTTTTCCTCaa                    >  2:148864/1‑140 (MQ=255)
                                                                                                            aTCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAaca                  <  2:8133/140‑1 (MQ=255)
                                                                                                             tCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAacac                 <  1:363339/140‑1 (MQ=255)
                                                                                                             tCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAacac                 <  1:607666/140‑1 (MQ=255)
                                                                                                                  aaCGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATg            <  1:1038847/140‑1 (MQ=255)
                                                                                                                    cggcggTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCt          <  1:691169/140‑1 (MQ=255)
                                                                                                                    cggcggTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCt          <  1:975742/140‑1 (MQ=255)
                                                                                                                     ggcggTGGTTTTTTTTTGATCAATTTCAAAATAAAAAAAATGATCCCATAAAAATAAAACAAGGTTTCAAATGATGTGGGTTTTATCACTTTTTTTTATTAAATGAATTTCTATCTTCGTTTCCCTCAAAACTGATGCTc         >  2:394353/1‑140 (MQ=255)
                                                                                                                            gtTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAg  <  2:769605/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                   
CGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAG  >  NC_000913/4101382‑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: ‑