Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1520AG35.3% 23.5 / 18.7 35intergenic (+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 (19/3);  new base G (0/12);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 8.30e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAG  >  NC_000913/3781053‑3781242
                                                                                                    |                                                                                           
gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGa                                                                                             >  2:628097/1‑100 (MQ=255)
gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGa                                                                                             >  1:772573/1‑100 (MQ=255)
 ggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGAg                                                                                            >  1:1054429/1‑100 (MQ=255)
   agagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGAggg                                                                                          >  2:1276903/1‑100 (MQ=255)
         gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGAGGGGA‑AAAg                                                                                   >  1:324739/1‑100 (MQ=255)
           ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGAGGGGA‑AAAg                                                                                   >  1:56353/1‑98 (MQ=255)
                  gcgcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGAGGGGA‑AAAgcgcgc                                                                              >  2:1170300/1‑91 (MQ=255)
                           gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGAGGGGA‑AAAGGCCGCACTGCTCACCAAt                                                                 >  1:36169/1‑100 (MQ=255)
                             aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCTCTCGCCCCTCC‑GGGGAGAGGGTCGGGGGGAGGGGA‑AAAGGCCGCACTGCTCACCAAttt                                                               >  2:876875/1‑100 (MQ=255)
                              cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGAGGGGA‑AAAGGCCGCACTGCTCACCAAtttt                                                              >  1:991129/1‑100 (MQ=255)
                                 ccTCCCCCTGAAACGACCCCGCCCTTTTTTTTCCTCCGGCGCCCCCCG‑GGGGGGGGGGCGGGGGTGAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTg                                                           <  2:140304/100‑1 (MQ=255)
                                        ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGAGGGGA‑AAAGGCCGCACTGCTAACCAATTTTCTGCTATCCt                                                    >  1:217398/1‑100 (MQ=255)
                                         tGAAACGACACCGCACTCTTTTTTTCCCCCTCGCCCCTCC‑GGGGGGGGGGCCGGGGTGGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                    <  1:131666/99‑1 (MQ=39)
                                                 cacCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCa                                           >  1:356135/1‑100 (MQ=255)
                                                   ccGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                         >  1:1071380/1‑100 (MQ=255)
                                                   ccGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                         >  2:327081/1‑100 (MQ=255)
                                                      ccccttttttttCCCCCCCCCCCCCCCGGGGGAGGGGGCGGGGGGGGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa                                       <  1:661941/95‑1 (MQ=255)
                                                          cttctgtccccccTCGCCCGCCG‑GGGGAGAGGGCCGGGGGGGGGGAA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                  <  2:1058830/90‑1 (MQ=255)
                                                           tttttttCCCCCCCGCCCCGCG‑GGGGGGA‑GGCCGGCGGGGGGGGATAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                                   <  2:878621/98‑1 (MQ=255)
                                                           ttttctccccccTCGCCCCGCG‑GGGGGGGGGCCCGGGGTGGGGGAA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgg                                     <  1:319480/87‑1 (MQ=255)
                                                             tttttCTTTCTCTCCCCTCC‑GGGGAGAGGGCCGGGGTGAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                               >  2:535040/1‑100 (MQ=255)
                                                             ttttcCCCCCCCCCCCCCCG‑GGGGGGGGGGCCGGGGGGGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                               <  1:329310/100‑1 (MQ=255)
                                                              ttttCTCCCTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAga                              >  2:590272/1‑99 (MQ=255)
                                                              ttttCCCCCCCTCCCTCCG‑GGGGGGAGGGGGGGGGGGCGGGCA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                 <  1:590276/97‑1 (MQ=255)
                                                                  ctccctccCCCCCCG‑GGGGGGGGGGCCGGGGTGAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATg                           <  2:772568/99‑1 (MQ=255)
                                                                   gcccTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                         >  1:720620/2‑100 (MQ=255)
                                                                    cccTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                        >  2:316459/1‑100 (MQ=255)
                                                                    cccTCCCCCCCCG‑GGGGGGAGGGCGGGGGTGGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                        <  2:94804/100‑1 (MQ=255)
                                                                     cccccccccTCG‑GGGGAGGGGCCGGGGGTGGGGGAA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                       <  1:215031/95‑1 (MQ=255)
                                                                          ccccTCG‑GGGGGGGGGGCGGGGGGGGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa                    <  1:1076698/98‑1 (MQ=255)
                                                                             ccccgggggAGGGGGCCGGGGGGAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                <  1:1170311/95‑1 (MQ=255)
                                                                              tCC‑GGGGAGAGGGCCGGGGTGAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc              >  1:75142/1‑100 (MQ=255)
                                                                                   gggAGAGGGCCGGGGTGAGGGAA‑AAAGGCCGCACTGCGCACCAATTGTCTGCTATCCTGCCACCGCACTAAGGGGGCAGTATGCTAAACATCGTAGTTTa          >  2:1147101/1‑100 (MQ=255)
                                                                                     gaggGGGCGGGGGGGGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg        <  2:991120/100‑1 (MQ=255)
                                                                                         gggCGGGGGGGGGGGAA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc    <  1:862393/100‑1 (MQ=255)
                                                                                         ggccgggggTGGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg        <  2:937911/91‑1 (MQ=255)
                                                                                           gCCGGGGTGAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg  >  1:532/1‑100 (MQ=255)
                                                                                                    |                                                                                           
GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCCGGGGTGAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAG  >  NC_000913/3781053‑3781242

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