Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1520AG27.9% 99.2 / 26.6 61intergenic (+135/‑63)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (37/7);  new base G (0/17);  total (37/24)
Fisher's exact test for biased strand distribution p-value = 6.45e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAG  >  NC_000913/3781058‑3781242
                                                                                                |                                                                                          
cgatctccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGa                                                                                       <  2:1708417/94‑1 (MQ=255)
 ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGaa                                                                                      >  2:891503/1‑100 (MQ=255)
   ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGaaaa                                                                                    >  2:1472163/1‑100 (MQ=255)
   ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGaaaa                                                                                    >  2:1028415/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                                                                        >  1:205159/1‑100 (MQ=255)
               gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCt                                                                        >  2: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                       >  1:206711/1‑100 (MQ=255)
                                                                ccTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                       >  2:206711/1‑100 (MQ=255)
                                                                ccTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                       >  2:2085417/1‑100 (MQ=255)
                                                                   cGCCCCGCG‑GGGGGGGGGGCG‑GGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa                    <  2:590876/100‑1 (MQ=255)
                                                                   cGCCCCGCG‑GGGGGGGGGGCG‑GGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa                    <  1:590876/100‑1 (MQ=255)
                                                                    gCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt                   >  2:1645597/1‑100 (MQ=255)
                                                                     ccccTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                  >  2:2232980/1‑100 (MQ=255)
                                                                     ccccTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                  >  2:2009404/1‑100 (MQ=255)
                                                                     ccccTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                  >  2:2232914/1‑100 (MQ=255)
                                                                        cgccgggggggggggCG‑GGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                 <  2:410264/94‑1 (MQ=255)
                                                                        cgccgggggggggggCG‑GGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                 <  1:410264/94‑1 (MQ=255)
                                                                            gggggggggggGC‑GGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt             <  2:1834530/98‑1 (MQ=255)
                                                                             ggggAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt           >  2:1911163/1‑100 (MQ=255)
                                                                               ggggggggCG‑GGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc         <  2:2629/95‑1 (MQ=255)
                                                                               ggggggggCG‑GGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc         <  1:2629/95‑1 (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       >  1:440124/1‑100 (MQ=255)
                                                                                 agagGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa       >  2:440124/1‑100 (MQ=255)
                                                                                    gggCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc    >  1:357517/1‑100 (MQ=255)
                                                                                    gggCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc    >  2:357517/1‑100 (MQ=255)
                                                                                    gggCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc    >  1:709498/1‑100 (MQ=255)
                                                                                    gggCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc    >  2:709498/1‑100 (MQ=255)
                                                                                      gCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg  >  2:1420041/1‑100 (MQ=255)
                                                                                                |                                                                                          
AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAG  >  NC_000913/3781058‑3781242

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: ‑