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

AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAG  >  NC_000913/3781058‑3781242
                                                                                                |                                                                                          
cgatctccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGa                                                                                       <  2:1708417/94‑1 (MQ=255)
 ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaa                                                                                      >  2:891503/1‑100 (MQ=255)
   ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                                    >  2:1472163/1‑100 (MQ=255)
   ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                                    >  2:1028415/1‑100 (MQ=255)
           cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCAc                                                                            >  1:309349/1‑100 (MQ=255)
           cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCAc                                                                            >  2:309349/1‑100 (MQ=255)
               gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCt                                                                        >  1:205159/1‑100 (MQ=255)
               gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCt                                                                        >  2: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                       >  1:206711/1‑100 (MQ=255)
                                                                ccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                       >  2:206711/1‑100 (MQ=255)
                                                                ccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                       >  2:2085417/1‑100 (MQ=255)
                                                                   cGCCCCGCGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa                    <  2:590876/100‑1 (MQ=255)
                                                                   cGCCCCGCGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa                    <  1:590876/100‑1 (MQ=255)
                                                                    gCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt                   >  2:1645597/1‑100 (MQ=255)
                                                                     ccccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                  >  2:2232980/1‑100 (MQ=255)
                                                                     ccccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                  >  2:2009404/1‑100 (MQ=255)
                                                                     ccccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                  >  2:2232914/1‑100 (MQ=255)
                                                                        cgccgggggggggggCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                 <  2:410264/94‑1 (MQ=255)
                                                                        cgccgggggggggggCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                 <  1:410264/94‑1 (MQ=255)
                                                                            gggggggggggGCGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt             <  2:1834530/98‑1 (MQ=255)
                                                                             ggggAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt           >  2:1911163/1‑100 (MQ=255)
                                                                               ggggggggCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc         <  2:2629/95‑1 (MQ=255)
                                                                               ggggggggCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc         <  1:2629/95‑1 (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       >  1:440124/1‑100 (MQ=255)
                                                                                 agagGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa       >  2:440124/1‑100 (MQ=255)
                                                                                    gggCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc    >  1:357517/1‑100 (MQ=255)
                                                                                    gggCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc    >  2:357517/1‑100 (MQ=255)
                                                                                    gggCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc    >  1:709498/1‑100 (MQ=255)
                                                                                    gggCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc    >  2:709498/1‑100 (MQ=255)
                                                                                      gCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg  >  2:1420041/1‑100 (MQ=255)
                                                                                                |                                                                                          
AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAG  >  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: