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

CTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGA  >  NC_000913/3781043‑3781237
                                                                                                  |                                                                                                  
ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAgg                                                                                                  <  2:731113/100‑1 (MQ=255)
ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAgg                                                                                                  <  1:731113/100‑1 (MQ=255)
               cgatctccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGa                                                                                  <  2:1708417/94‑1 (MQ=255)
                ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaa                                                                                 >  2:891503/1‑100 (MQ=255)
                  ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                               >  2:1028415/1‑100 (MQ=255)
                  ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                               >  2:1472163/1‑100 (MQ=255)
                          cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCAc                                                                       >  1:309349/1‑100 (MQ=255)
                          cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCAc                                                                       >  2:309349/1‑100 (MQ=255)
                              gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCt                                                                   >  2:205159/1‑100 (MQ=255)
                              gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCt                                                                   >  1:205159/1‑100 (MQ=255)
                               cTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTc                                                                  >  2:428726/1‑100 (MQ=255)
                               cTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTc                                                                  >  1:428726/1‑100 (MQ=255)
                                    tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCaa                                                             >  2:780748/1‑100 (MQ=255)
                                    tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCaa                                                             >  1:780748/1‑100 (MQ=255)
                                    tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTATTTCTCCCTCGCCCCTCCGGGGAGCGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCa                                                              >  2:1810096/1‑99 (MQ=255)
                                    tgcccccccTCACCCGGAAAGGCCACCCCCCTCTTTTTTCCCCCCTCCCCCCGCGGGGGAGGGGGCCGGGGTGGGGGGAAAAGGCCGCACTGCTCAcc                                                               <  2:167412/92‑1 (MQ=255)
                                    tgcccccccTCACCCGGAAAGGCCACCCCCCTCTTTTTTCCCCCCTCCCCCCGCGGGGGAGGGGGCCGGGGTGGGGGGAAAAGGCCGCACTGCTCAcc                                                               <  1:167412/92‑1 (MQ=255)
                                       cctccctccCTCTGAAACGCCCCCGCACTTTTTTTTTCTCCCCCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAttt                                                          <  2:2296283/99‑1 (MQ=38)
                                        cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                                                         >  2:1576898/1‑100 (MQ=255)
                                          cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                       >  2:1820111/1‑100 (MQ=255)
                                           cccccctcgaaaaCGACCCGCCCCTCTTTTTTTCCCCCGCCCCCCCCGGGGGGGGGGCCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                      <  2:263559/90‑1 (MQ=255)
                                           cccccctcgaaaaCGACCCGCCCCTCTTTTTTTCCCCCGCCCCCCCCGGGGGGGGGGCCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                      <  1:263559/90‑1 (MQ=255)
                                           ccTCACTCGGAAACGACACCGCCCTTTTTTTTCCTCCCCCCCCCCTCCGGGGAGGGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                      <  2:1766470/100‑1 (MQ=255)
                                            cccACTCTAAAACGCCCCCGCCTTTTTTTTTCCCCCCCCCCCCCCCGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                     <  2:1265468/98‑1 (MQ=255)
                                                  tgaaaaCGCCCCCGCCCTCTTTTTTTCCCCCCCCCCCCCCCGGGGAGGGGCCCGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                <  2:1708221/96‑1 (MQ=255)
                                                   tGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                              >  2:1910780/1‑100 (MQ=255)
                                                   tGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                              >  2:922238/1‑100 (MQ=255)
                                                   gaaaaCGACCCCCCCCCTTTTTTTCCCCCCCCGCCCCGCGGGGGAGGGGGCCGGGGTGGGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                <  2:2099558/96‑1 (MQ=255)
                                                     aaaCGACCCCGCACTCTTTTTTCCCCCCTCCCCCCTCGGGGGAGGGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                            <  2:1460971/100‑1 (MQ=38)
                                                          acacCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGc                                       >  2:2190603/1‑100 (MQ=255)
                                                                  ttttttttttCCCCCTCCCCCCCCGGGGGAGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                               <  2:491764/97‑1 (MQ=255)
                                                                  ttttttttttCCCCCTCCCCCCCCGGGGGAGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                               <  1:491764/97‑1 (MQ=255)
                                                                  ctctTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                               >  2:1890201/1‑100 (MQ=255)
                                                                   cctTTTTTTTCCCCCCCCCCCCCCGGGGGGGGGGCCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgg                                <  2:2195760/98‑1 (MQ=255)
                                                                           tctcCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCCCACTAAGGGGGCAGTATg                      >  1:249971/1‑100 (MQ=255)
                                                                           tctcCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCCCACTAAGGGGGCAGTATg                      >  2:249971/1‑100 (MQ=255)
                                                                              cccgccccccGCGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                       <  2:2093499/90‑1 (MQ=255)
                                                                               ccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                  >  2:206711/1‑100 (MQ=255)
                                                                               ccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                  >  2:2085417/1‑100 (MQ=255)
                                                                               ccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                  >  1:206711/1‑100 (MQ=255)
                                                                                  cGCCCCGCGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa               <  1:590876/100‑1 (MQ=255)
                                                                                  cGCCCCGCGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa               <  2:590876/100‑1 (MQ=255)
                                                                                   gCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt              >  2:1645597/1‑100 (MQ=255)
                                                                                    ccccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc             >  2:2009404/1‑100 (MQ=255)
                                                                                    ccccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc             >  2:2232914/1‑100 (MQ=255)
                                                                                    ccccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc             >  2:2232980/1‑100 (MQ=255)
                                                                                       cgccgggggggggggCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg            <  1:410264/94‑1 (MQ=255)
                                                                                       cgccgggggggggggCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg            <  2:410264/94‑1 (MQ=255)
                                                                                           gggggggggggGCGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt        <  2:1834530/98‑1 (MQ=255)
                                                                                            ggggAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt      >  2:1911163/1‑100 (MQ=255)
                                                                                              ggAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc    <  1:734727/100‑1 (MQ=255)
                                                                                              ggAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc    <  2:734727/100‑1 (MQ=255)
                                                                                                agagGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGTTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa  >  2:313686/1‑100 (MQ=255)
                                                                                                agagGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGTTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa  >  1:313686/1‑100 (MQ=255)
                                                                                                agagGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa  >  2:440124/1‑100 (MQ=255)
                                                                                                agagGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa  >  1:440124/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
CTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGA  >  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: