Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,326,4500TG52.5% 44.5 / 57.3 40G114G (GGT→GGG) atoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (4/15);  new base G (21/0);  total (25/15)
Fisher's exact test for biased strand distribution p-value = 9.64e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.90e-06

GCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAA  >  NC_000913/2326306‑2326598
                                                                                                                                                |                                                                                                                                                     
gcgTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGaa                                                                                                                                                 <  1:303858/151‑1 (MQ=255)
 cgTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATgg                                                                                                                                                   >  1:77284/1‑148 (MQ=255)
         aCTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGt                                                                                                                                        <  2:659749/151‑1 (MQ=255)
          cTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATTAg                                                                                                                                         >  2:381382/1‑149 (MQ=255)
            gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCAGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGtt                                                                                                                                       <  2:446129/149‑1 (MQ=255)
                     cGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCta                                                                                                                             <  1:160196/150‑1 (MQ=255)
                      gggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGTAAATATGAGGTTAGACCCCtcc                                                                                                                            >  1:186527/1‑148 (MQ=255)
                        gCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAACATGTGTTTATCCCCCca‑cct                                                                                                                         >  2:275153/1‑148 (MQ=255)
                           ggCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTA‑CTTACTc                                                                                                                     >  1:351298/1‑151 (MQ=255)
                                                cACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTGTAGGCCCCTC‑CTTAATCGATGCAAAAGCACTCTCCgg                                                                                                 >  2:202760/1‑150 (MQ=255)
                                                   ggTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAGCACGCTCTGGTTAt                                                                                             <  2:201975/151‑1 (MQ=255)
                                                          aaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTg                                                                                      <  2:473699/151‑1 (MQ=255)
                                                              tATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTCTGGAAAAAATGAGTTTTGCCCCCTA‑CTTACTCGATGGAAAAGAACACCCTGGGTATCGTCTTggg                                                                                    >  2:594279/1‑148 (MQ=255)
                                                              tATGTGGTTCGGGTCTTAAAAGCGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTC‑CTTACTCGATGCAAAAACACGCTCTGGTTATCGTCtt                                                                                       >  2:152405/1‑146 (MQ=255)
                                                                tgtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGAc                                                                                 <  1:206285/150‑1 (MQ=255)
                                                                    gTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTTGCCCCCCA‑CTTTCTCGATGCAAAAACACGCTCTGGGTATCGTCTTGgaggag                                                                                >  2:687486/1‑144 (MQ=255)
                                                                     ttCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAGCACGCCCTGGTTATTGTCTTGGAGACggcc                                                                             >  2:575422/1‑147 (MQ=255)
                                                                       cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGt                                                                         <  2:489031/151‑1 (MQ=255)
                                                                       cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAACACGCTCTGGGTATCGTCTTGGAGACGGACAgg                                                                          >  1:202307/1‑150 (MQ=255)
                                                                        gggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGTGGTTTCCCCCCTACCTTTCT‑GCGGCCCAAGCACGCTCTGGTTATCTTCTTTGGGACgg                                                                               >  2:547684/1‑144 (MQ=255)
                                                                         ggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGGGACGGACCgggt                                                                        >  1:560534/1‑148 (MQ=255)
                                                                               aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACg                                                                 <  2:478231/151‑1 (MQ=255)
                                                                               aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGCTGGAAAATATGAGGTTTGCCCCCCA‑CTTACTCGATGCCAAAACACCCTCTGGTTATCGTCTTGGGGACGGACAGGGTTATGACg                                                                 >  2:625771/1‑151 (MQ=255)
                                                                                aaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTTGAGGCGGACAGGTTTTAGACg                                                                 >  2:683937/1‑150 (MQ=255)
                                                                                            ttGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCCA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGGTTATGACGTAATCCTGcgcg                                                     >  1:459727/1‑150 (MQ=255)
                                                                                             tGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGGAGAAGCACGCTCTGGTTATCGTCTTTGAGAAGGACCgg                                                                          >  2:280175/1‑128 (MQ=255)
                                                                                              gccgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATg                                                  <  1:218041/151‑1 (MQ=255)
                                                                                              gccgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATg                                                  <  2:541406/151‑1 (MQ=255)
                                                                                                 gccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATTAGGTTCACCCCCCC‑CTTTAACGATGTAAAAAACCCCTCCGGTTAATGTATTGGAGACGGACAGGTTTATGACGTAAACCTGCGCGATg                                                  >  2:693895/1‑148 (MQ=255)
                                                                                                  cccAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCCA‑CTTACTCGGTGGCAAAGCACGCCCTGGGTATCGTCTTGGAGACGGACAGGTTTATGACGGAATCCTGCGCGGTGGcc                                               >  2:103173/1‑150 (MQ=255)
                                                                                                    cAGGCCATTCAGGCAGGTCAGGCGCAGAGCCTTGTGGCGGGGGGGCTGGAAAATATTAGCTTTGGCCCCCC‑CTTTATCGATGCAAAAGGACGCTCCGGTTTTCGTCTTGGGGACGGGCAGGGTTATGACGTACTTCATCGCGATGGACt                                              >  2:149377/1‑149 (MQ=255)
                                                                                                      ggCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGTGTTTTGCCCCCTA‑CTTACTCGATGCAAAAGCCCGCCCTGGGTATCGTCTTGGAGACCGACAGGGTTTTGAAGGAAACCCGCGCGATGGGCTGAt                                           >  2:68228/1‑150 (MQ=255)
                                                                                                        ccATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATTGAAAATATGAGGTTAGCCCCCCA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTTGAGACGGACAGGGTTATGACGTAATCCCGCGCGATGGCCTGAtg                                          >  2:77722/1‑149 (MQ=255)
                                                                                                                 caggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGGTTAGCCCCCTA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGGGACGGACAGGTTTTTGACGTAATCCTGCGCGATGGCCTGATGTGCGCCAcccc                               >  2:296294/1‑150 (MQ=255)
                                                                                                                   ggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATcc                                                           <  1:247581/121‑1 (MQ=255)
                                                                                                                   ggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCCA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATcc                                                           >  2:247536/1‑121 (MQ=255)
                                                                                                                       aGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTa                         <  1:16644/151‑1 (MQ=255)
                                                                                                                              gagCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATg                  <  1:80396/151‑1 (MQ=255)
                                                                                                                                     gtgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGAtt            <  2:612975/150‑1 (MQ=255)
                                                                                                                                               gTATGGAAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGaaa  >  2:67782/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTA‑CTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAA  >  NC_000913/2326306‑2326598

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