Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,101,5140GT28.6% 84.4 / 10.5 35intergenic (+84/‑176)sodA/kdgTsuperoxide dismutase, Mn/2‑keto‑3‑deoxy‑D‑gluconate transporter
Reads supporting (aligned to +/- strand):  ref base G (2/23);  new base T (10/0);  total (12/23)
Fisher's exact test for biased strand distribution p-value = 3.60e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCCGGCTTCCCGATTATGGGCCTGGATGTGTGGGAACATGCTTACTACCTGAAATTCCAGAACCGCCGTCCGGACTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAAT‑GAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTA‑CGCATC‑AGCCCTTAAAAATATGCCG‑ACAGGTGATGGAAATGCAGATAAAACGCTCGATTG  >  NC_000913/4101287‑4101711
                                                                                                                                                                                                                                   |                                                                                                                                                                                                         
tCCGGCTTCCAGATTATGGACCTGGATGTGGTGGAACATGCTTACTACCTGAAATTCCAGAACCGCCGTCCGGACTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTc                                                                                                                                                                                                <  2:417680/239‑1 (MQ=255)
  cGGCTTCCCGATTATGGGCCTGGATGGGTGGGAACATGCTTACTACCTGAAATTCCAGAACCGCCGTCCTGACTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCaa                                                                                                                                                                                              <  2:447130/239‑1 (MQ=255)
                         aTGTGTGGGAACATGCTTACTACCTGAAATTCCAGAACCGCCGTCCGGACTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGataa                                                                                                                                                                        <  2:672967/238‑1 (MQ=255)
                                  aaCATGCTTACTACCTGAAATTCCAGAACCGCCGTCCGGACTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAAc                                                                                                                                                               <  1:637483/238‑1 (MQ=255)
                                                   aaaTTCCAGAACCGCCGTCCGGACTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTGATTAAATTCTAAATAAAAAAAAAAATCCGATAAAAATAAAACAAAGTTTCAACTGAtgtt                                                                                                                                             >  2:580691/1‑238 (MQ=255)
                                                                   gTCCGGACTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTa                                                                                                                             <  1:580691/239‑1 (MQ=255)
                                                                         aCTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTa                                                                                                                             <  1:237863/233‑1 (MQ=255)
                                                                           tACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGataa                                                                                                                                                                        <  1:808/188‑1 (MQ=255)
                                                                                   agagTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTTGTCAATTTTAAAATAAAAACAAAAATCCGATAAAAAAAAAAAACAGCGTTAATTTAAGGGGGTTTTATTACATTTATTGATTAaatac                                                                                                                 >  2:237863/1‑233 (MQ=255)
                                                                                      gTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAAT‑GAATGTCt                                                                                                          <  1:86409/238‑1 (MQ=255)
                                                                                         cTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGTTTAATTTCAAAAAAAAAACAAAAGTCCGGTAAAAAAAAAAAAACGATTTAATTTATATGGTTTTGGTCACTTTTATTTATTAaat‑gaatggtgatat                                                                                                      >  2:267520/1‑232 (MQ=255)
                                                                                              aCGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAAT‑GAATGTCTATCTTCGtt                                                                                                 <  1:267520/239‑1 (MQ=255)
                                                                                               cGTGGTGAACTGGGACGAAGCGGCGGCACGTTTTGCGGCGAAAAAATAATCATTCGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAAT‑GAATGTCTATCTTCGtt                                                                                                 <  2:274216/238‑1 (MQ=255)
                                                                                                      aaCTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGAtt                                                                                                                       <  1:199457/210‑1 (MQ=255)
                                                                                                      aaCTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAAATTTAAAAAAAAAAAAAAGAAACCGTAAAAAAAAAAAAACGCTTTAAATTATTTGGGTTTGTTTACTTTTTTTTATTTAATAGAATGT‑TATCTTCTGTTCCTTCaa                                                                                         >  2:566235/1‑239 (MQ=255)
                                                                                                                aaGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGGTCAATTTTCAAATAAAAAAAAAGATCCCGTAAAAATAAAAAAAGGGTTTAATTTATTTGGGTTTTTTCACCTTTTTTTTTTATAT‑GAAAGt                                                                                                            >  2:199457/1‑210 (MQ=255)
                                                                                                                     gcGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAAT‑GAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCAtt                                                                          <  2:143880/239‑1 (MQ=255)
                                                                                                                      cGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCaa                                                                                                                                                                                              >  2:242503/1‑123 (MQ=255)
                                                                                                                      cGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCaa                                                                                                                                                                                              <  1:242503/123‑1 (MQ=255)
                                                                                                                       ggCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAAT‑GAATGTCTATCTTCGTTTCCATCAACACTGATGCt                                                                               <  2:465475/232‑1 (MQ=255)
                                                                                                                           cGTTTGGCAGCGAAAAAATAATCATTTGCATCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAAT‑GAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGaa                                                                    <  2:511297/239‑1 (MQ=255)
                                                                                                                              tttGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAAT‑GAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAAtt                                                                  <  1:566235/238‑1 (MQ=255)
                                                                                                                                       aaaaaaTAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAAT‑GAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTA‑CGCATC‑Ag                                                       <  1:79211/238‑1 (MQ=255)
                                                                                                                                          atataatCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAAT‑GAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTA‑CGCATC‑AGCCCt                                                   <  1:56004/237‑1 (MQ=255)
                                                                                                                                                    ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTGGTCAATTTTCAAATAAAAAAAAAGAACCGATAAAAAAAAAACAACGGTTCCATTTATGAGGGTTTTATTCCTTTTTTTTATTAAAT‑TAAAGTCTTACTTCGTTTCCTTTAAAACTCATGCTCCATTTATGAATTA‑ACGATCAAGCCCTTACAAata                                           >  2:690047/1‑238 (MQ=255)
                                                                                                                                                               tgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGAtt                                                                                                                       <  1:629589/153‑1 (MQ=255)
                                                                                                                                                                      aGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaat‑ga                                                                                                                <  1:664762/152‑1 (MQ=255)
                                                                                                                                                                         cGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTCAAAATAAAAAAAAAAATCCGATAAAAAAAAAAAAACGTTTCTATTGATGTGGGTTTTTTCACTTTTTTTTTTTAAAT‑TAAAGt                                                                                                            >  2:629589/1‑153 (MQ=255)
                                                                                                                                                                                          aGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATTAAAATTAAAAAAAAAAAAAAGATCCGATAAAAATAAAAAAGCGGTTTAATTGAAGTGGGTTTGATCACATTTATTGAATAAAC‑GAATGTCTATCTTCGTTTTCATTAACACTTATGGTGCAATTAGGAATTA‑AGCATA‑AGGCCCTAAAAAAAAGACGAACAGGGAGTGGAAAGGCAG‑TAAAAAGGTCGAttt  >  2:69073/1‑238 (MQ=255)
                                                                                                                                                                                          aGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAAATTTCAAATAAAAAAAAAAATCCGAAAAAAATAAAAAAGCGGTTTTATTTATGTGGTTTTTATTCCTTTTATTGATTAAAA‑GAAAATCTATTTTTCTTTTCATCAAAACCTTAGCTCCATTGTGGAAGTATTGCTCA‑ACCCCTAAAAAATATGCTG‑ACACGGAATGGAAATCAAGATAAAAACCTCGCtt   >  2:267046/1‑239 (MQ=255)
                                                                                                                                                                                            cTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAAT‑GAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTA‑CGCATC‑AGCCCTTAAAAATATGCCG‑a                                    <  1:439684/203‑1 (MQ=255)
                                                                                                                                                                                                         ccATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAAT‑GAATGTCTATCTTCg                                                                                                   <  1:502542/130‑1 (MQ=255)
                                                                                                                                                                                                             cGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGAt                                                                                                                                     <  2:179335/93‑1 (MQ=255)
                                                                                                                                                                                                             cGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGAt                                                                                                                                     >  1:179335/1‑93 (MQ=255)
                                                                                                                                                                                                                   cggcggTGGTTTTTTTTTGATCAATTTCAAAATAAAAACAAAGATCCGATAAAAATAAAAAAACGTTTTAATTGATGTGGGTTTTATCACTTTTTTTGATTAAAT‑GAATGTCTATCTTTGTTTTCATCaa                                                                                         >  2:502542/1‑130 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                         
TCCGGCTTCCCGATTATGGGCCTGGATGTGTGGGAACATGCTTACTACCTGAAATTCCAGAACCGCCGTCCGGACTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAAT‑GAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTA‑CGCATC‑AGCCCTTAAAAATATGCCG‑ACAGGTGATGGAAATGCAGATAAAACGCTCGATTG  >  NC_000913/4101287‑4101711

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