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

GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAG  >  NC_000913/3781053‑3781242
                                                                                                    |                                                                                           
gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                                                             >  2:628097/1‑100 (MQ=255)
gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                                                             >  1:772573/1‑100 (MQ=255)
 ggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                                                                                            >  1:1054429/1‑100 (MQ=255)
   agagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                                                                          >  2:1276903/1‑100 (MQ=255)
         gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAg                                                                                   >  1:324739/1‑100 (MQ=255)
           ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAg                                                                                   >  1:56353/1‑98 (MQ=255)
                  gcgcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAgcgcgc                                                                              >  2:1170300/1‑91 (MQ=255)
                           gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAt                                                                 >  1:36169/1‑100 (MQ=255)
                             aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCTCTCGCCCCTCCGGGGAGAGGGTCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAAttt                                                               >  2:876875/1‑100 (MQ=255)
                              cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                                                              >  1:991129/1‑100 (MQ=255)
                                 ccTCCCCCTGAAACGACCCCGCCCTTTTTTTTCCTCCGGCGCCCCCCGGGGGGGGGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                           <  2:140304/100‑1 (MQ=255)
                                        ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTAACCAATTTTCTGCTATCCt                                                    >  1:217398/1‑100 (MQ=255)
                                         tGAAACGACACCGCACTCTTTTTTTCCCCCTCGCCCCTCCGGGGGGGGGGCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                    <  1:131666/99‑1 (MQ=39)
                                                 cacCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCa                                           >  1:356135/1‑100 (MQ=255)
                                                   ccGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                         >  1:1071380/1‑100 (MQ=255)
                                                   ccGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                         >  2:327081/1‑100 (MQ=255)
                                                      ccccttttttttCCCCCCCCCCCCCCCGGGGGAGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa                                       <  1:661941/95‑1 (MQ=255)
                                                          cttctgtccccccTCGCCCGCCGGGGGAGAGGGCCGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                  <  2:1058830/90‑1 (MQ=255)
                                                           tttttttCCCCCCCGCCCCGCGGGGGGGAGGCCGGCGGGGGGGGATAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                                   <  2:878621/98‑1 (MQ=255)
                                                           ttttctccccccTCGCCCCGCGGGGGGGGGGCCCGGGGTGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgg                                     <  1:319480/87‑1 (MQ=255)
                                                             tttttCTTTCTCTCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                               >  2:535040/1‑100 (MQ=255)
                                                             ttttcCCCCCCCCCCCCCCGGGGGGGGGGGCCGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                               <  1:329310/100‑1 (MQ=255)
                                                              ttttCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAga                              >  2:590272/1‑99 (MQ=255)
                                                              ttttCCCCCCCTCCCTCCGGGGGGGAGGGGGGGGGGGCGGGCAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                 <  1:590276/97‑1 (MQ=255)
                                                                  ctccctccCCCCCCGGGGGGGGGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATg                           <  2:772568/99‑1 (MQ=255)
                                                                   gcccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                         >  1:720620/2‑100 (MQ=255)
                                                                    cccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                        >  2:316459/1‑100 (MQ=255)
                                                                    cccTCCCCCCCCGGGGGGGAGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                        <  2:94804/100‑1 (MQ=255)
                                                                     cccccccccTCGGGGGAGGGGCCGGGGGTGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                       <  1:215031/95‑1 (MQ=255)
                                                                          ccccTCGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa                    <  1:1076698/98‑1 (MQ=255)
                                                                             ccccgggggAGGGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                <  1:1170311/95‑1 (MQ=255)
                                                                              tCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc              >  1:75142/1‑100 (MQ=255)
                                                                                   gggAGAGGGCCGGGGTGAGGGAAAAAGGCCGCACTGCGCACCAATTGTCTGCTATCCTGCCACCGCACTAAGGGGGCAGTATGCTAAACATCGTAGTTTa          >  2:1147101/1‑100 (MQ=255)
                                                                                     gaggGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg        <  2:991120/100‑1 (MQ=255)
                                                                                         gggCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc    <  1:862393/100‑1 (MQ=255)
                                                                                         ggccgggggTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg        <  2:937911/91‑1 (MQ=255)
                                                                                           gCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg  >  1:532/1‑100 (MQ=255)
                                                                                                    |                                                                                           
GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAG  >  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: