Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,326,450	0	T	G	53.3%	33.5 / 46.7	31	G114G (GGT→GGG)	atoB	acetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/13);  new base G (16/0);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 1.17e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.37e-05

GGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCG  >  NC_000913/2326303‑2326595                                                                                                                                                    |                                                                                                                                                  ggCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATg                                                                                                                                                <  2:340662/151‑1 (MQ=255)      gTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGaaaa                                                                                                                                           <  2:543326/151‑1 (MQ=255)         aGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAata                                                                                                                                         >  2:617345/1‑150 (MQ=255)             aCTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCTGGGGGTATGGAAAATATGAGt                                                                                                                                    <  1:47908/151‑1 (MQ=255)               tGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGttt                                                                                                                                  >  2:392321/1‑151 (MQ=255)                gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTa                                                                                                                                 <  2:401405/151‑1 (MQ=255)                   aaaaaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGc                                                                                                                               <  2:565981/150‑1 (MQ=255)                         cGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtac                                                                                                                        <  2:131085/151‑1 (MQ=255)                               ggCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtacttact                                                                                                                   <  1:561120/150‑1 (MQ=255)                                     aaCGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGa                                                                                                                >  2:339540/1‑147 (MQ=255)                                                 aTTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGctc                                                                                                 >  2:113916/1‑150 (MQ=255)                                                    cACGGTCAATAGAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGGCAGGCGCAGAGCATTGTGGCGGCGGGCCTGAAAAAAAAGAGGATCGCCCCCTACTCTCACGATGCaaaaagaagctttgg                                                                                              >  1:554867/1‑139 (MQ=255)                                                    cACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTgg                                                                                              >  2:599906/1‑150 (MQ=255)                                                     aCGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGtt                                                                                            <  1:27641/151‑1 (MQ=255)                                                            aTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCCACCTTCTCGATGCCAAAGCAAGCTCTGGTTATCGTCt                                                                                     >  1:579794/1‑151 (MQ=255)                                                              aaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCtt                                                                                    >  2:80196/1‑150 (MQ=255)                                                                   aTGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATTGAAAATATGAGCTTACCCCCCTAATTCCTCGATGCAAAAGCACGCTCTGGTTATCGTCTTTGAGAc                                                                              >  2:287615/1‑151 (MQ=255)                                                                         ttCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATTAGTTTAGCCCCCCACTTTTTCGATGCAAAAGGCCGCCCTGGTTATCGTCTTGGAGACg                                                                             >  1:638908/1‑146 (MQ=255)                                                                              gTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCGGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTa                                                                   <  2:218467/151‑1 (MQ=255)                                                                                  tAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGa                                                                >  1:396109/1‑150 (MQ=255)                                                                                     aaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATTGGTTTTGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATTGTCTTGGAGACGGACAGGGTTTTGAc                                                               >  2:573982/1‑148 (MQ=255)                                                                                             cgcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATAAGAGGTTAGACCCCCACTTACTCGATGCAAAAGGACGCTCTGGTTATTGTCTTGGAGAAGGACAGGGTTTTGACGTAATcc                                                        >  2:412820/1‑147 (MQ=255)                                                                                               cTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTAGCCCCCCACTTACTCGATGCAAAAAGACGCTCTGGGTATCGGCTTGGAGACGGACAAGGTTATGACGTAATcc                                                        >  1:409247/1‑145 (MQ=255)                                                                                                     gccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAGGTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGGGACGGACAGGGTTATGACGTAATCCTGCGCGATGGcc                                            >  2:276846/1‑151 (MQ=255)                                                                                                         aGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGAt                                        <  2:662688/151‑1 (MQ=255)                                                                                                                     caggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCa                            <  2:654849/151‑1 (MQ=255)                                                                                                                       ggTCAGGCGCAGAGCATTGTGGCCGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCTCCCAt                           <  1:381516/150‑1 (MQ=255)                                                                                                                              cgcAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCa                   <  2:389707/151‑1 (MQ=255)                                                                                                                                cAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGCTGCAAAAACACGCTCTGGTTATCGTCTTTGGGACGGACAGGGTTATTACGTAATCCTGCGCGATGGGCTGATGTGCGGCCCCCCTGGTTTTCat                  >  1:373019/1‑150 (MQ=255)                                                                                                                                      aTTGTGGCGGGGGGGGTGGAAAATATGAGGGTAGCCCCCTCCTTAATCGATGC‑AAAGGACGCTCTGGTTATCGTCTTGGGGACGGACAGGGTTATGACGTAATCCTGCGCGATGGGCTGGTGGGGGCCCCCCCTGGGTTTTATTTggggg           >  1:75094/1‑149 (MQ=255)                                                                                                                                               gggggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccgccg  >  2:531153/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCG  >  NC_000913/2326303‑2326595

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