Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,326,4500TG45.8% 62.4 / 13.2 48G114G (GGT→GGG) atoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (4/22);  new base G (21/1);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.09e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAA  >  NC_000913/2326303‑2326598
                                                                                                                                                   |                                                                                                                                                    
ggCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATg                                                                                                                                                   <  2:405833/151‑1 (MQ=255)
 gcgcGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATgg                                                                                                                                                  <  1:1051991/151‑1 (MQ=255)
 gcgcGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATgg                                                                                                                                                  <  2:918582/151‑1 (MQ=255)
        aGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGgaa                                                                                                                                               >  2:1203901/1‑147 (MQ=255)
          gCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGTCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATg                                                                                                                                          >  1:824963/1‑150 (MQ=255)
           cACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAg                                                                                                                                        <  1:269026/151‑1 (MQ=255)
               gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGtt                                                                                                                                      <  1:227995/149‑1 (MQ=255)
                ttAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGgtt                                                                                                                                     >  2:617808/1‑149 (MQ=255)
                 tAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAg                                                                                                                                   <  2:119686/150‑1 (MQ=255)
                  aaaaaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGcc                                                                                                                                 >  1:963741/1‑151 (MQ=255)
                    aaaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGaccc                                                                                                                               >  1:941788/1‑151 (MQ=255)
                       gCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCt                                                                                                                             >  1:27122/1‑150 (MQ=255)
                              ggCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtacttact                                                                                                                      <  1:302374/150‑1 (MQ=255)
                               gCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGTAAAAATGGGGTTTTCCCCCTCCTTtctc                                                                                                                     >  2:325074/1‑150 (MQ=255)
                                  gAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGAAGGGAAATATGAGGTTAGCCCCCCACTTACTCGAt                                                                                                                  >  2:1098513/1‑150 (MQ=255)
                                          gtgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGGTTTGCCCCCCACTTTCTCGATGGaaaa                                                                                                            >  2:24004/1‑148 (MQ=255)
                                            gCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCAc                                                                                                        <  1:226657/150‑1 (MQ=255)
                                                      ggTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTTGTTAt                                                                                             >  1:610605/1‑151 (MQ=255)
                                                         cAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGt                                                                                          <  1:505889/151‑1 (MQ=255)
                                                             aaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAAAACGCTCTTGTTATCGTCtt                                                                                       >  1:1185179/1‑150 (MQ=255)
                                                              aaGTATGTGGTTCGGGTCTTAAAAGTGAGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGCTGGGAAATATGAGGTTAGCCCCCTATTTATTCGCTTCAAAACCCCGCTCTCGTTATCGTCtt                                                                                       >  1:715943/1‑149 (MQ=255)
                                                                    gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTTACCCCCCACTTACTCGGTGCAAAAAAACGCTCTGGTTATCGTCCTTGAGACgg                                                                               >  2:535517/1‑151 (MQ=255)
                                                                        ttCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCTAAATACTCGCTGCAAAAACACGGTCCAGCTATCg                                                                                           >  1:856070/1‑135 (MQ=255)
                                                                           gggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGtt                                                                        <  1:650063/151‑1 (MQ=255)
                                                                            ggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGttt                                                                       <  2:696138/151‑1 (MQ=255)
                                                                              tCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGGCCCCTACTTACTCGATGCAAAAGCACGCCCTGGGTATCGTCTTTGAGACGGACAGGgtttt                                                                     >  2:558116/1‑149 (MQ=255)
                                                                                 tAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGAAGGGAAAAATGAGGTTAGCCCCCCACTTACTCGATGGAAAAGCACGCCCTGGGTATTGTCTTGGAGACGGACAGGttt                                                                       >  1:1128305/1‑146 (MQ=255)
                                                                                  aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGAc                                                                  <  1:314065/150‑1 (MQ=255)
                                                                                    aaGTGTGGCGCTTGCCGCGCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTa                                                               <  2:120280/151‑1 (MQ=255)
                                                                                     aGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGTTGGAAAATATTAGGTTTACCCCCCACTTACTCGATGAAAAAGCACCCCCTGGTTCACGTCCTGGGGACCGACAGGTTTATGACGt                                                                >  1:1046934/1‑149 (MQ=255)
                                                                                            cgcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTgcg                                                       <  2:979293/151‑1 (MQ=255)
                                                                                            cgcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTgcg                                                       <  2:1252556/151‑1 (MQ=255)
                                                                                               ttGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGGAAATATGAGTTTAGCCCCCCACTTTCTCGATGCAAAAGCACGCCCTGGGTATTGTCTTGGGGACGGGCCGGGTTATGAAGTAAACCcgcg                                                       >  2:1014554/1‑148 (MQ=255)
                                                                                                tGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTGGGGCAAATATGAGTTTAACCCCCTACTTTTTCCGTGCAAAAACAAGCTCTGGTTATCGTATTGGAGGCAGACAGGCTTATGACGTAATACTACGGGAt                                                   >  1:955909/1‑151 (MQ=255)
                                                                                                tGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGcgcg                                                     <  1:340945/149‑1 (MQ=255)
                                                                                                       cAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGa                                            <  1:1264730/151‑1 (MQ=255)
                                                                                                       cAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGAGGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGGTTATGACGTAATCCTGCGCGATGGCCt                                              >  1:1266056/1‑149 (MQ=255)
                                                                                                        aGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAAAGAGTTTAGTCCCCTACTTAGTTGATGGCAAAAAACGCTCTTGCTATCGGCTTGGAGGCGGAGAGGATTTTGATCTAATTCTGTGCGGTGGCCt                                              >  1:656071/1‑148 (MQ=255)
                                                                                                          gCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGCCCCCTACTTACTCGCTGCAAAAACACGCTCTGGTTATCGTCTTGGGGACGGACAGGGTTATGACGTAATCCTGCGCGATGGCCTg                                             >  1:927039/1‑147 (MQ=255)
                                                                                                                     aggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGTGTTTAGCCCCCCACCTACTCGCTGGAAAAAAACGCTCTGGGTATCGTCTTGGAGAAGGACAGGGTTATGACGTAATCCCGCGCGATGGGCCGATGGGCGCCAccc                                >  2:1195753/1‑149 (MQ=255)
                                                                                                                     aggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCAt                              <  1:690868/151‑1 (MQ=255)
                                                                                                                            gcgcAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATc                       <  1:928276/151‑1 (MQ=255)
                                                                                                                              gcAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCCACTTTCTCGATGCAAAAGCACGCTCTGGGTATCGTCTTGGAGACGGACAGGGTTATGAAGGAATCCTGCGCGATGGGCTGATGTGCGCCACCCATgggtttta                      >  2:251255/1‑144 (MQ=255)
                                                                                                                               cAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCata                    <  2:358268/151‑1 (MQ=255)
                                                                                                                                  agCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCTCCCATGGTTATCATATg                  <  1:1048789/150‑1 (MQ=255)
                                                                                                                                   gCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATggg                <  2:637521/151‑1 (MQ=255)
                                                                                                                                             ggggggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccgcc      >  1:269824/1‑151 (MQ=255)
                                                                                                                                                 ggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCTATTCAAAAGCACGCTCTGGTTATCTTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGaaa  <  2:361378/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAA  >  NC_000913/2326303‑2326598

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