Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1400AG31.2% 67.7 / 24.7 49intergenic (+123/‑75)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (27/6);  new base G (0/15);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 4.96e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGA  >  NC_000913/3781043‑3781237
                                                                                                  |                                                                                                  
ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAgg                                                                                                  <  2:731113/100‑1 (MQ=255)
ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAgg                                                                                                  <  1:731113/100‑1 (MQ=255)
               cgatctccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGa                                                                                  <  2:1708417/94‑1 (MQ=255)
                ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGaa                                                                                 >  2:891503/1‑100 (MQ=255)
                  ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGaaaa                                                                               >  2:1028415/1‑100 (MQ=255)
                  ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGaaaa                                                                               >  2:1472163/1‑100 (MQ=255)
                          cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCAc                                                                       >  1:309349/1‑100 (MQ=255)
                          cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCAc                                                                       >  2:309349/1‑100 (MQ=255)
                              gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCt                                                                   >  2:205159/1‑100 (MQ=255)
                              gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCt                                                                   >  1:205159/1‑100 (MQ=255)
                               cTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTc                                                                  >  2:428726/1‑100 (MQ=255)
                               cTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTc                                                                  >  1:428726/1‑100 (MQ=255)
                                    tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCaa                                                             >  2:780748/1‑100 (MQ=255)
                                    tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCaa                                                             >  1:780748/1‑100 (MQ=255)
                                    tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTATTTCTCCCTCGCCCCTCC‑GGGGAGCGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCa                                                              >  2:1810096/1‑99 (MQ=255)
                                    tgcccccccTCACCCGGAAAGGCCACCCCCCTCTTTTTTCCCCCCTCCCCCCGCG‑GGGGAGGGGGCC‑GGGGTGGGGGGAAAAGGCCGCACTGCTCAcc                                                               <  2:167412/92‑1 (MQ=255)
                                    tgcccccccTCACCCGGAAAGGCCACCCCCCTCTTTTTTCCCCCCTCCCCCCGCG‑GGGGAGGGGGCC‑GGGGTGGGGGGAAAAGGCCGCACTGCTCAcc                                                               <  1:167412/92‑1 (MQ=255)
                                       cctccctccCTCTGAAACGCCCCCGCACTTTTTTTTTCTCCCCCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAttt                                                          <  2:2296283/99‑1 (MQ=38)
                                        cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                                                         >  2:1576898/1‑100 (MQ=255)
                                          cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                       >  2:1820111/1‑100 (MQ=255)
                                           cccccctcgaaaaCGACCCGCCCCTCTTTTTTTCCCCCGC‑CCCCCCC‑GGGGGGGGGGCCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                      <  2:263559/90‑1 (MQ=255)
                                           cccccctcgaaaaCGACCCGCCCCTCTTTTTTTCCCCCGC‑CCCCCCC‑GGGGGGGGGGCCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                      <  1:263559/90‑1 (MQ=255)
                                           ccTCACTCGGAAACGACACCGCCCTTTTTTTTCCTCCCCCCCCCCTCC‑GGGGAGGGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                      <  2:1766470/100‑1 (MQ=255)
                                            cccACTCTAAAACGCCCCCGC‑CTTTTTTTTTCCCCCCCCCCCCCCCGGGGGGGGGGGCG‑GGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                     <  2:1265468/98‑1 (MQ=255)
                                                  tgaaaaCGCCCCCGCCCTCTTTTTTTCCCCCCCCCCCCCCC‑GGGGAGGGGCCC‑GGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                <  2:1708221/96‑1 (MQ=255)
                                                   tGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                              >  2:1910780/1‑100 (MQ=255)
                                                   tGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                              >  2:922238/1‑100 (MQ=255)
                                                   gaaaaCGACCCCCCCCCTTTTTTTCCCCCCCCGCCCCGCG‑GGGGAGGGGGCC‑GGGGTGGGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                <  2:2099558/96‑1 (MQ=255)
                                                     aaaCGACCCCGCACTCTTTTTTCCCCCCTCCCCCCTCG‑GGGGAGGGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                            <  2:1460971/100‑1 (MQ=38)
                                                          acacCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGc                                       >  2:2190603/1‑100 (MQ=255)
                                                                  ttttttttttCCCCCTCCCCCCCCG‑GGGGAGGGGGCG‑GGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                               <  2:491764/97‑1 (MQ=255)
                                                                  ttttttttttCCCCCTCCCCCCCCG‑GGGGAGGGGGCG‑GGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                               <  1:491764/97‑1 (MQ=255)
                                                                  ctctTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                               >  2:1890201/1‑100 (MQ=255)
                                                                   cctTTTTTTTCCCCCCCCCCCCCC‑GGGGGGGGGGCCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgg                                <  2:2195760/98‑1 (MQ=255)
                                                                           tctcCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCCCACTAAGGGGGCAGTATg                      >  1:249971/1‑100 (MQ=255)
                                                                           tctcCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCCCACTAAGGGGGCAGTATg                      >  2:249971/1‑100 (MQ=255)
                                                                              cccgccccccGCG‑GGGGGGGGGGCG‑GGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                       <  2:2093499/90‑1 (MQ=255)
                                                                               ccTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                  >  2:206711/1‑100 (MQ=255)
                                                                               ccTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                  >  2:2085417/1‑100 (MQ=255)
                                                                               ccTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                  >  1:206711/1‑100 (MQ=255)
                                                                                  cGCCCCGCG‑GGGGGGGGGGCG‑GGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa               <  1:590876/100‑1 (MQ=255)
                                                                                  cGCCCCGCG‑GGGGGGGGGGCG‑GGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa               <  2:590876/100‑1 (MQ=255)
                                                                                   gCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt              >  2:1645597/1‑100 (MQ=255)
                                                                                    ccccTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc             >  2:2009404/1‑100 (MQ=255)
                                                                                    ccccTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc             >  2:2232914/1‑100 (MQ=255)
                                                                                    ccccTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc             >  2:2232980/1‑100 (MQ=255)
                                                                                       cgccgggggggggggCG‑GGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg            <  1:410264/94‑1 (MQ=255)
                                                                                       cgccgggggggggggCG‑GGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg            <  2:410264/94‑1 (MQ=255)
                                                                                           gggggggggggGC‑GGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt        <  2:1834530/98‑1 (MQ=255)
                                                                                            ggggAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt      >  2:1911163/1‑100 (MQ=255)
                                                                                              ggAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc    <  1:734727/100‑1 (MQ=255)
                                                                                              ggAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc    <  2:734727/100‑1 (MQ=255)
                                                                                                agagGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGTTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa  >  2:313686/1‑100 (MQ=255)
                                                                                                agagGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGTTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa  >  1:313686/1‑100 (MQ=255)
                                                                                                agagGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa  >  2:440124/1‑100 (MQ=255)
                                                                                                agagGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa  >  1:440124/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
CTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGA  >  NC_000913/3781043‑3781237

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