Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,153,5730TG41.0% 63.0 / 18.7 39Y92D (TAC→GAC) pabC4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base T (2/21);  new base G (16/0);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 4.06e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAA‑CGAGGGGATTAC‑GTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACA  >  NC_000913/1153428‑1153721
                                                                                                                                                 |                                                                                                                                                       
ccAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCGTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACa                                                                                                                                                      <  1:532627/149‑1 (MQ=255)
   gCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGcaca                                                                                                                                                 <  2:864569/151‑1 (MQ=255)
        tACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACAtt                                                                                                                                               >  2:331518/1‑148 (MQ=255)
         aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGa                                                                                                                                             <  2:941621/149‑1 (MQ=255)
          cAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACa                                                                                                                                          <  2:261520/151‑1 (MQ=255)
              aTGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCgg                                                                                                                                      <  2:279975/151‑1 (MQ=255)
                gCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGGGcacaattaaaa                                                                                                                                          >  1:627230/1‑138 (MQ=255)
                gCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGAc                                                                                                                                    <  1:474539/151‑1 (MQ=255)
                   tGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAACACAATTAAAAGCgg                                                                                                                                      >  2:713688/1‑146 (MQ=255)
                   tGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGAcc                                                                                                                                   <  2:81440/149‑1 (MQ=255)
                   tGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGAcc                                                                                                                                   <  2:519729/149‑1 (MQ=255)
                        gCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAAcg                                                                                                                            <  1:1032807/151‑1 (MQ=255)
                           gTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGC‑gg                                                                                                                        <  1:531854/151‑1 (MQ=255)
                                  atTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGC‑GGATtctct                                                                                                                 <  1:195197/151‑1 (MQ=255)
                                                            tGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCACAATGTAACGCGGGCCGGCCACGC‑GGGTTCGTTCCCGTTCCGCGTATCCTGcccctt                                                                                         >  1:515886/1‑144 (MQ=255)
                                                               aCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAAAACACAATGAACACCGGACCCGCCACGC‑GGGTTTTCTCTCCTACGGCTTAt                                                                                                   >  2:139644/1‑136 (MQ=255)
                                                                  agagATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGt                                                                                 <  2:549171/151‑1 (MQ=255)
                                                                    agaTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACATTTCAAAGCGGACCGGCAACGC‑GGGGTCTCTCCCTTTCCGGTTATTCTTCCCATTATGAcccgtt                                                                               >  2:704234/1‑147 (MQ=255)
                                                                      aTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGc                                                                             <  2:646296/151‑1 (MQ=255)
                                                                            aCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCACAATTAACAGCGGACCGGCAACGC‑GGGTTCTCTCCCTTTCGGCTTTTTCTGCACAATACGACCGTTTGCGTaa                                                                         >  2:304980/1‑149 (MQ=255)
                                                                              gCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAGCACAATGAACAACGGACCGGGCACCC‑CGGGTCTCTCCCTTTCGGCGTTTTCTCCACCATACCACCGTTTTCGtga                                                                         >  1:964565/1‑145 (MQ=255)
                                                                              gCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAA‑CGAg                                                                    <  2:758107/151‑1 (MQ=255)
                                                                                   cagcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAA‑CGAGGGGAt                                                               <  1:635133/151‑1 (MQ=255)
                                                                                      cagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAGGACATTTAAAAAAGGGCCGGCGACGC‑GGGGTCTCTCCCTTTCGGCCTTTTCTTCACAATTCGACCGTTTTCGTGA‑CGAggggggt                                                              >  1:65338/1‑146 (MQ=255)
                                                                                       agAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGCAGGGGGGTCAGCACCATTAACAACGGACCCGCAACGC‑GGGGTCTCCCCGTTTCGGCGTTTACTGCCCATTACCATCACTTGCGGAA‑CGAgggg                                                                 >  1:691531/1‑145 (MQ=255)
                                                                                        gAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGTGTAA‑CGAGGGGATTAC‑Gt                                                         <  2:95837/151‑1 (MQ=255)
                                                                                              cAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAACACCATTAACACCGGGCCGGGAAAGCGGGATTTCTTCCTTTCCGGCTT‑TCCTCCCCCTTTTTCCCGTTTGCGTAACCAAGGGGTTTCC‑GTTg                                                       >  2:803439/1‑148 (MQ=255)
                                                                                                 aaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAACACAATTAAAAGCGGGCCGGCAACGC‑CGGGTTTCTCCCTTTCGGGTTATTCTGCCCAATAAGACCGGTTTTGTAA‑CGAGGGGGTTAC‑Gtt                                                        >  2:601940/1‑143 (MQ=255)
                                                                                                          cTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGAACAATGTAAAGCGGACCGGCAAAGC‑GGGGTGTCTCCCTTACGGCTTATCCCTCACATTAAGACCCCTTTCGGAA‑AGAGGGGGGTAC‑GTTGGGGCCAAGcc                                             >  1:237063/1‑145 (MQ=255)
                                                                                                           tGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGGAGAACATTTAACAGCGGACCGGCAAAGC‑GGGGTTTCTCCCTTTCGGCTTTTTCTGCACATTACGACCGGTTTCGTTA‑AGAGGGGGTTAC‑GTTGGGGCTAAGCCCCGGGCg                                      >  2:742261/1‑151 (MQ=255)
                                                                                                             aaaGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGAGGAAAGCCGCATGAAAAACCGGCCGGGCACGC‑GGGGTCTCTCCCTTTCGGCTTTTTCCCCACAATACCACCGTTTTGGTAA‑CGAGGGGGATACGGTTGGGCCTAGCCCCGGTggg                                      >  2:821585/1‑148 (MQ=255)
                                                                                                                gTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAA‑CGAGGGGATTAC‑GTTGGCGCTAAGCCCGGTGCGGCTgg                                 <  1:820378/151‑1 (MQ=255)
                                                                                                                         aGTCGCGGTAGTGGCGGGCGGGGGGACAGCAACTTGAAAAACGGGCCGGCCACGC‑GGGGTCTCTCCCTTTCGGCTTTTTCTTCACCTTACGACCGCTTGCGGGA‑AGAGGGGGGTAC‑GGTGGTGCTAAGCCCGGTGGGGCCgggggg                             >  2:43325/1‑144 (MQ=255)
                                                                                                                            cgcgGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAA‑CGAGGGGATTAC‑GTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCt                      <  2:498620/150‑1 (MQ=255)
                                                                                                                                gTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAA‑CGAGGGGATTAC‑GTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCAtat                 <  2:29134/151‑3 (MQ=255)
                                                                                                                                   gTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAA‑CGAGGGGATTAC‑GTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGc              <  1:416172/151‑1 (MQ=255)
                                                                                                                                            gAGGGTACAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAA‑CGAGGGGATTAC‑GTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTa      >  2:209283/1‑150 (MQ=255)
                                                                                                                                              gggTACAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAA‑CGAGGGGATTAC‑GTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAAc   >  2:997254/1‑151 (MQ=255)
                                                                                                                                               ggTACAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAA‑CGAGGGGATTAC‑GTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACa  <  2:702414/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                       
CCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAA‑CGAGGGGATTAC‑GTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACA  >  NC_000913/1153428‑1153721

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