Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,326,450	0	T	G	58.3%	35.2 / 67.1	48	G114G (GGT→GGG)	atoB	acetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/19);  new base G (28/0);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 1.73e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-10

GTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAA  >  NC_000913/2326308‑2326599                                                                                                                                               |                                                                                                                                                       gTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGaaaa                                                                                                                                                <  1:237749/151‑1 (MQ=255) gTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGgaa                                                                                                                                                 >  1:190973/1‑150 (MQ=255)    aGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGgaaa                                                                                                                                                >  1:1019173/1‑148 (MQ=255)      gCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAT‑ATGa                                                                                                                                          <  1:134238/151‑1 (MQ=255)           gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAT‑ATGAGtt                                                                                                                                       <  2:1111913/149‑1 (MQ=255)             tAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAT‑ATGAGTTTAg                                                                                                                                    <  1:415541/150‑1 (MQ=255)                 aaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAA‑ATGAGGTTATgcccc                                                                                                                               >  2:1118458/1‑151 (MQ=255)                      ggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAT‑AAGAGGTTAGcccccccctt                                                                                                                          >  2:1056240/1‑146 (MQ=255)                         tGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATAATGGGTTTGAcccccacct                                                                                                                           >  2:389707/1‑144 (MQ=255)                           gCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCg                                                                                                                     <  1:906158/151‑1 (MQ=255)                             aGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAA‑ATGAGGTTAGCCCCCCACTTTCTCg                                                                                                                     >  2:789328/1‑149 (MQ=255)                                 aCGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGAGGGAAAA‑ATGAGATTAGTCCCCcaccttct                                                                                                                       >  2:1044385/1‑138 (MQ=255)                                  cGGTGAGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAGAAT‑ATGAGGTTAGCCCCCTCCTTACTTGactca                                                                                                                >  2:721133/1‑145 (MQ=255)                                      gtgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCaaaa                                                                                                             >  1:219820/1‑148 (MQ=255)                                               cACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAT‑ATGAGTGTAGCCCCCTAATTACTCTATGCaaaaagaagcgctgg                                                                                                  >  2:124156/1‑139 (MQ=255)                                                     cAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAT‑ATGAGGTTTGCCCCCTACTTACTCCATGCAAAAACACGCTCTGGtt                                                                                                >  2:532831/1‑146 (MQ=255)                                                      aaTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAT‑ATGAGGTTTGGCCCCCACTTACTCGATGCAAAAACACGCCCTGGTTATCGt                                                                                           >  1:372683/1‑150 (MQ=255)                                                         aaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCtt                                                                                        <  1:111454/150‑1 (MQ=255)                                                             tATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAT‑ATTAGTTTTGCCCCCTACTTTCTCGATGCAAAAGAACCCTCTGGGTATCGTCtt                                                                                        >  1:981211/1‑146 (MQ=255)                                                                gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACgg                                                                                <  2:353900/151‑1 (MQ=255)                                                                  ggTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAT‑ATGAGGTTAGCCCCCCACTTACTCGATGCAAAAACACGCCCTGGGTATCGTCTTTGAGACgg                                                                                >  1:24346/1‑149 (MQ=255)                                                                   gTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGGTATCGTCTTGGAGACgg                                                                                >  1:275548/1‑148 (MQ=255)                                                                      cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAT‑ATGAGTTTAGGCCCCCACTTACTCGATGCAAAAACACGCTCTGGGTATCGTCTTGGAGACGGACAgg                                                                           >  1:394882/1‑150 (MQ=255)                                                                      cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAA‑AAGAGGTTTGCCCCCCACTTACTCGATGCAAAAACACGCTCTGGGTATCGTCTTTGGGACGGACa                                                                             >  2:889473/1‑148 (MQ=255)                                                                          tCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAT‑ATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTTGGGACGGACAGGGTTAt                                                                      >  1:621873/1‑151 (MQ=255)                                                                              aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAT‑ATGGGTTTATGCCCCCACTTTCTCGATGCAAAAACACGCCCTGGGTATCGTCTTTGGGGAGGGCAGGGTTAt                                                                      >  2:188112/1‑147 (MQ=255)                                                                              aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAT‑ATGAGGTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTTGAGACGGACAGGgtt                                                                        >  1:617187/1‑145 (MQ=255)                                                                                aaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAT‑ATGAGGTTAGCCCCCCCCTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGGGACGGACAGGGTTATGACGt                                                                 >  1:122816/1‑150 (MQ=255)                                                                                       gcgcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAT‑AAGAGGTTACCCCCCTCATTACTCCATGCAAAAGCAAGCTCTGGGTATCGTCTTGGGGACGGACAGGGTTATGAAGTAATCCTgc                                                         >  2:535433/1‑151 (MQ=255)                                                                                          cTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAT‑ATGAGGTTAGCCCCCCACTTACTTGATGCAAAAACACGCTCTGGTTATCGTCTTGGGGAAGGGCAGGGTTATGACGTAATACTgcg                                                        >  1:783747/1‑149 (MQ=255)                                                                                            tGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGcgcg                                                      <  1:336053/149‑1 (MQ=255)                                                                                            tGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAT‑ATGAGGTTTGCCCCCCACTTACTCGATGCAAAAACACGCTCTGGTTATAGTTTTGGAGACGGACAGGGTTTTTAAGTAATCCcgcg                                                        >  2:1144443/1‑147 (MQ=255)                                                                                              ccgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATgg                                                  <  1:1145291/151‑1 (MQ=255)                                                                                               cgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAT‑ATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGCCCTGGTTATCGTCTTGGGGACGGACAGGGTTATGACGTAATCCTGCGCGATgg                                                  >  1:342799/1‑150 (MQ=255)                                                                                                 cccAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAA‑ATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACGCTCTGGGTATCGTCTTTGAGACGGACAGGGTTATGACGTAATCCTGCGCGATGGcc                                                >  1:409286/1‑150 (MQ=255)                                                                                                  ccAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAG‑ATGAGGTTATGACCCCACTTACTCGGTGGCAAACCACGCCCCCGTTATCGTCTTGGGGGCGGAAAGGGTTTTTACGGAATCTTACGCGGTGGCCTg                                              >  1:300209/1‑151 (MQ=255)                                                                                                        cATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAAT‑ATGAGTTTTGACCCCTCCTTACTCGATGCAAAAGCCCCCTCTGGGTATCGTCTTTGACACGGACAGGTTTATGACGTAAACCTACGCGATGGCCTGAtgtg                                         >  1:1050321/1‑150 (MQ=255)                                                                                                            caggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCa                                    <  2:143908/151‑1 (MQ=255)                                                                                                              ggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAA‑ATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGCTCTGGTTATTGTCTTGGGGACGGACAGGGTTATTACGTAATCCTGCGCGGTGGCCTGATGTGCGCCAcc                                  >  2:108315/1‑151 (MQ=255)                                                                                                                caggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCa                                <  1:1150557/151‑1 (MQ=255)                                                                                                                  ggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCTCCCATg                              <  1:373282/151‑1 (MQ=255)                                                                                                                      aGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTa                          <  2:607242/151‑1 (MQ=255)                                                                                                                                   tgtgGCGGGGGGTATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGAt              <  2:1119505/150‑1 (MQ=255)                                                                                                                                   tgtgGCGGGGGGTATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGAt              <  1:797471/150‑1 (MQ=255)                                                                                                                                      ggCGGGGGGTATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAcc          <  1:128161/151‑1 (MQ=255)                                                                                                                                      ggCGGGGGGTATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAc           <  2:844799/150‑1 (MQ=255)                                                                                                                                         ggggggTATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccgcc       <  2:142678/151‑1 (MQ=255)                                                                                                                                          gggggTATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccgccg      <  2:190964/151‑1 (MQ=255)                                                                                                                                               tATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGaaaa  <  2:783707/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                       GTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAT‑ATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAA  >  NC_000913/2326308‑2326599

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