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

TGGATGTGTGGGAACATGCTTACTACCTGAAATTCCAGAACCGCCGTCCGGACTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCG  >  NC_000913/4101309‑4101738
                                                                                                                                                                                                              |                                                                                                                                                                                                                                 
tGGATGTGTGGGAACATGCTTACTACCTGAAATTCCAGAACCGCCGTCCGGACTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGC‑GGGGTTTTTTTGTTTCAATTTCACATTaaaaa                                                                                                                                                                                                             >  2:542870/1‑228 (MQ=255)
        tgGGAACATGCTTACTACCTGAAATTCCAGAACCGCCGTCCGGACTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGGTTTTTTTTG‑TACAATTTATAAATaaaaaaatagtac                                                                                                                                                                                                     >  2:206788/1‑220 (MQ=255)
             aCATGCTTACTACCTGAAATTCCAGAACCGCCGTCCGGACTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGA‑GGGGTTTTTTTG‑TTCAATTTTCTACTAAAAACAATgtaacgcta                                                                                                                                                                                                >  2:443705/1‑219 (MQ=255)
              cATGCTTAACACCTGAAATTCCAGACCCGCCGGCCGGACTACATTAAAGAGTTCTTGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGataa                                                                                                                                                                                               <  1:510097/227‑1 (MQ=255)
                tGCTTACTACCTGAAATTCCAGAACCGCCGTCCGGACTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGGGTTTTTTTG‑TTCAAATTTTAAAAAAAAACAAAGAtacctgtaa                                                                                                                                                                                              >  2:577959/1‑218 (MQ=255)
                                                  gACTACATTAAAGGGTTCTGGAACGTTGTAAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCa                                                                                                                                                          <  2:356660/228‑1 (MQ=255)
                                                                 ttCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGGTTTTTTTTG‑TTCAATTTTCAAATTAAAAAAAAAAACCGCTAAAAAAAAAAAAAAGATTTCATTGATTTGGGTTTTTTCACAtttttttttttaa                                                                                                                                           >  2:351895/1‑221 (MQ=255)
                                                                 ttCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGGTTTTTTGTG‑ATCAATTTCAAAATa                                                                                                                                                                                                                 >  1:468571/1‑158 (MQ=255)
                                                                 ttCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGC‑GGTTTTTTTTGG‑ATCAATTTTTAAATAAAAAACATGATTCGCTAAAAATCAAAAAGCGTTTCAATTTATGTAGGTTTTATCACTTTTa                                                                                                                                                    >  1:89611/1‑219 (MQ=255)
                                                                  tCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATaa                                                                                                                                                                                                                <  2:467803/158‑1 (MQ=255)
                                                                  tCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGat                                                                                                                                                                                                 <  1:292964/173‑1 (MQ=255)
                                                                           tggtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatgaa                                                                                                                                       >  1:431381/1‑222 (MQ=255)
                                                                            ggtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatgaat                                                                                                                                      <  2:430679/222‑1 (MQ=255)
                                                                                              agcGGCACGTTTTGCGGCGAAAAAATAATCATTTGACGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAacac                                                                                                              <  1:543746/228‑1 (MQ=255)
                                                                                                      gTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTc                                                                                                                                                                                                                       <  1:485112/115‑1 (MQ=255)
                                                                                                             ggcgAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCAAATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGGTTTTTTGTG‑ATCAATTTTAAAATAAAAAAAAAGATCCGACAAAAAAAAAAAAGCGGTTCAAATGATGTGGGTTTGa                                                                                                                                                             >  2:309820/1‑166 (MQ=255)
                                                                                                              gcgAAAAAATAATCATTTGCAGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGATTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGAt                                                                                                                                                            <  1:310311/166‑1 (MQ=255)
                                                                                                                  aaaaaTAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTg                                                                                                 <  1:578888/221‑1 (MQ=255)
                                                                                                                        aatCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCAt                                                                                    <  1:506139/228‑1 (MQ=255)
                                                                                                                                        ctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGGTTTTTTTTG‑GTCAATTTCAAAATAAAAAAAAAGATCCGATAAAAAAAAAAAAACGTTTTAAATGATTTGGGTTTGATTACTTTTTTTTATTAAATGAATTGCTATTTTCGTTTtcatca                                                                                                                  >  2:60442/1‑182 (MQ=255)
                                                                                                                                         tgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCaa                                                                                                                 <  1:60532/182‑1 (MQ=255)
                                                                                                                                             aTGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGa                                                               <  1:218812/228‑1 (MQ=255)
                                                                                                                                                    gTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGa                                                        <  1:609336/228‑1 (MQ=255)
                                                                                                                                                       taGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATgg                                                     <  1:352483/228‑1 (MQ=255)
                                                                                                                                                            cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGGTTTTTTTTG‑ATTAATTTTAAAATAAAAACAAAGATCCGATAAAAATAAAAAAGCGTTTCAATTGATGTGGTTTTTATTCCTTTTATTTATTTAAATAATGTTTTTCTTTGTTTTCATCAAAACTTAt                                                                                                          >  2:252084/1‑170 (MQ=255)
                                                                                                                                                            cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGGTTTTTTTTG‑ATCAATTTCAAAATAAAAACAAAGATCCGATAAAAAAAAAACAACGTTTCAATTTATGTGGGTTTGATCACTTTTATTGATTAAATGAATGTTTATCTTCGTTTCCTTCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTTAAAAAAAGTCGGCAGGGGATGGtaaa                                                 >  2:531453/1‑226 (MQ=255)
                                                                                                                                                             gCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCaa                                                                                                                                                                                                                     <  1:671954/62‑1 (MQ=255)
                                                                                                                                                             gCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATTATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATg                                                                                                         <  1:252474/170‑1 (MQ=255)
                                                                                                                                                             gCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtg                                                                                                                                                                      <  1:588208/109‑1 (MQ=255)
                                                                                                                                                             gCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATgg                                                     <  1:532311/222‑1 (MQ=255)
                                                                                                                                                               atatCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGGTTTTTTTTG‑GTCAATTTTAAAAAAAAAAAAAAGAACCGATAAAAATAAAACAGCGGTTTaa                                                                                                                                                                            >  2:215480/1‑101 (MQ=255)
                                                                                                                                                                tatCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAAt                                                                                                                                                                           <  1:215810/101‑1 (MQ=255)
                                                                                                                                                                       ttAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGat                                                                                                                                                                                                 <  1:336791/72‑1 (MQ=255)
                                                                                                                                                                        tAAAAAATGAACCATCGCCAACGGCGGT‑GGGTTTTTTTTG‑ATTAAATTCAAAAAAAAAAAAAAGATCCGATAAAAATAAAACAGCGGTTCAATTGATGTGGGTTTGATCACTTTTATTGATTAAATGAAAGTCTATTTTCGTTTTCATCAAAACTGATGCTCCATTGAGGAATTAAGCCACAGTCCTTTAAAAAATTCCGATAGGGTATGGGAATGcctgttactct                                     >  2:246853/1‑217 (MQ=255)
                                                                                                                                                                         aaaaaaTGAACCATCACCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGt                                                          <  1:247232/205‑1 (MQ=255)
                                                                                                                                                                                                            ttGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCg  <  1:363051/227‑1 (MQ=255)
                                                                                                                                                                                                              |                                                                                                                                                                                                                                 
TGGATGTGTGGGAACATGCTTACTACCTGAAATTCCAGAACCGCCGTCCGGACTACATTAAAGAGTTCTGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGT‑GGTTTTTTTGTG‑ATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCG  >  NC_000913/4101309‑4101738

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