Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC61.8% 23.6 / 43.7 34T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (13/0);  new base C (0/21);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 1.08e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.06e-07

TTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCG‑TTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376358‑1376655
                                                                                                                                                       |                                                                                                                                                      
ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCG‑TTAAA‑CCGCTATCAGGATCATATGATCAACAccc                                                                                                                                                     >  2:348919/1‑151 (MQ=255)
                             ggACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCG‑TTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                         >  2:709272/1‑150 (MQ=255)
                             ggACGGCGACCGAAAAAGGGTGAGTGATCCCCGGCGCGTACTGAAACGGGTCGCCGCGCGTACCGGACCC‑GT‑GGTTTATCTCGACCCG‑TAAAC‑CCCCTATCAGGATCATATGACCAACCCCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                          <  2:224798/149‑1 (MQ=255)
                                   cgcccaaaaaaTGGTGAGTGTTCCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGACCC‑GG‑GGTGTATCTCGACCCT‑TTAAC‑CCGCTATCAGGATCATATGATCACCCCCCTCCCCGATCCCCGCCGTTACATCGTCGAAAACGAt                                                                                                                  <  1:624854/144‑1 (MQ=255)
                                          aataTGGGGAGTGATCCCCTGCGCGTCATTGAACAGGTCGCCGCACGTACGGAACCC‑GT‑GGTGTCTCTCTAACCG‑TTAAC‑CCCTGACCGGGCACATATGAACACCCCCCTCGCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                            <  2:357152/147‑1 (MQ=255)
                                              tGGTGAGTGATCCCCTGCGCGTACTGAAACAGGCCCCCGCGCGTCCCGGACCC‑GT‑GGTGTATCTCGACCCT‑TAAAC‑CCGCTATCAGGATCATATAACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                        <  1:78842/150‑1 (MQ=255)
                                               ggTGAGTGTTTCCCGGCGGGTCCGGGAAAAGGTCGCGGCGCGTCCCGGACCG‑GG‑GGTGTATCTC‑AACCGTTTAAC‑CCGCTATAAGGATCAAATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                         <  1:459706/148‑1 (MQ=255)
                                               ggTGAGTGTTCCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGACCC‑GT‑GGTTTATCTCGAACGG‑TTAAC‑CCGCTATCAGAATCATATGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTa                                                                                                      <  1:544204/151‑1 (MQ=255)
                                                          cccTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCG‑TTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                           >  1:16848/1‑151 (MQ=255)
                                                             tGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCG‑TTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                         >  1:293807/1‑150 (MQ=255)
                                                                    aCTGGAACAGGTCGCCGCGCGTACCGGACCC‑GT‑GGTGTATCTCGACCCG‑TTAAC‑CCGCTATCAGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                   <  1:124525/149‑1 (MQ=255)
                                                                        gAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCG‑TTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                              >  1:629251/1‑150 (MQ=255)
                                                                         acaagttCGCCGCGCGTACCGGGCCG‑GT‑GTTTTATCTGGAACCG‑TAAAC‑CCGCTAACAGGATCATATAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                             <  2:701050/144‑1 (MQ=255)
                                                                         aaagggTCGCCGCGCGTCCCGGAACC‑GT‑GTTGTA‑CCCGGACCC‑TTAAACCCGCTATGAGGATCATATGATCACCCCCCTCGCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                            <  1:36554/147‑1 (MQ=255)
                                                                         aaCAGGTCGCCGCGCGTACCGGAACT‑GG‑GGTTTATCTCGAACGG‑TAAAC‑CCGCTATCAGTATCATATGACCACCCCCCCCCCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                              <  2:374052/149‑1 (MQ=255)
                                                                             ggTCGCCGCGCGTCCCGGAACG‑GT‑GTTGTATCTCGACCCT‑TTAAC‑CCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                          <  1:557145/149‑1 (MQ=255)
                                                                                    gcgcgTACCGGAACC‑GT‑GGTGTATCTCGAACCG‑TTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTATATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                   >  1:319051/1‑149 (MQ=255)
                                                                                         tACCGGAACC‑GT‑GGTGTATCTCGAACCG‑TTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                >  1:589186/1‑147 (MQ=255)
                                                                                           ccGGAACCGGT‑GGGTTATC‑CCAACCC‑TTAAACCCGTTAACAGAATCAAAT‑AACCACCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                           <  2:377021/150‑1 (MQ=255)
                                                                                           ccGGAACC‑GT‑GGTGTATCTCGAACCG‑TTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                           >  1:388900/1‑150 (MQ=255)
                                                                                            cGGAACC‑GTGGGTTTTTCTAGACCCG‑TTAAA‑CCGCTACCAGAAGAATATAATCACCCCCCCCGCGTATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                            <  2:57095/149‑1 (MQ=255)
                                                                                                       ggtgTATCTCGAACCG‑TTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                                  >  1:375204/1‑149 (MQ=255)
                                                                                                             cactccacccttaaaa‑CCGCTATGAGCACCATGGTAACAACCCCCTCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                              <  1:192887/135‑1 (MQ=255)
                                                                                                              ctcGACCCG‑TAAAC‑CCGCTATCAGGATCATATGATCAACCCCCCCGCCAATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                               <  2:320944/145‑1 (MQ=255)
                                                                                                                      cttaaaa‑CCCCTATCAGAATCATATAATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                                   <  1:504390/146‑1 (MQ=255)
                                                                                                                        ttAAC‑CCGCTATAAGAATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAc                                 <  2:216548/150‑1 (MQ=255)
                                                                                                                        ttAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                                >  1:62241/1‑151 (MQ=255)
                                                                                                                          aaacccGCTATCAGGATCATATAACCAACCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                               <  1:451119/147‑1 (MQ=255)
                                                                                                                           aa‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                             >  2:321411/1‑151 (MQ=255)
                                                                                                                            acccGCTATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                             <  1:484031/149‑1 (MQ=255)
                                                                                                                              ccGCTAACAGGATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                                <  2:32567/146‑1 (MQ=255)
                                                                                                                                   aTCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGc                       >  2:713319/1‑150 (MQ=255)
                                                                                                                                            aaatGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt               <  2:188835/147‑1 (MQ=255)
                                                                                                                                                  tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat       >  2:347999/1‑151 (MQ=255)
                                                                                                                                                       acCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCgg  >  1:520706/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
TTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GT‑GGTGTATCTCGAACCG‑TTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376358‑1376655

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