Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,326,450	0	T	G	21.1%	109.5 / 19.3	52	G114G (GGT→GGG)	atoB	acetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (2/39);  new base G (11/0);  total (13/39)
Fisher's exact test for biased strand distribution p-value = 1.29e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTAT  >  NC_000913/2326312‑2326576                                                                                                                                           |                                                                                                                                ggCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTa                                                                                                                                <  2:416700/140‑1 (MQ=255)   cACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑Tg                                                                                                                             <  2:1176636/140‑1 (MQ=255)           aaagGCGGGCTGGCAGAAACGGTGTGCGGATTCTCGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAatat                                                                                                                     <  1:222922/140‑1 (MQ=255)                  ggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGtttt                                                                                                              >  1:355556/1‑139 (MQ=255)                    cTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGc                                                                                                            <  1:1060563/140‑1 (MQ=255)                     tGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGcc                                                                                                           <  1:137490/140‑1 (MQ=255)                     tGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGcc                                                                                                           <  2:694649/140‑1 (MQ=255)                      ggCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGccc                                                                                                          <  1:959925/140‑1 (MQ=255)                      ggCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGccc                                                                                                          <  1:746767/140‑1 (MQ=255)                      ggCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGccc                                                                                                          <  2:254798/140‑1 (MQ=255)                       gCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGcccc                                                                                                         <  2:1120936/140‑1 (MQ=255)                                 tgtgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGG‑TGGAAAAAATGAGTTTAGCCCCCCACTTACTc                                                                                               >  2:1058579/1‑140 (MQ=255)                                   tgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGa                                                                                             <  1:186968/140‑1 (MQ=255)                                   tgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGa                                                                                             <  1:589441/140‑1 (MQ=255)                                     cGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATg                                                                                           <  1:630791/140‑1 (MQ=255)                                           cACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGA‑TGGGAAAAATGGGTTTAGCCCCCCACTTACTCGATGGaaaaa                                                                                     >  2:1346996/1‑139 (MQ=255)                                            aCGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGc                                                                                    <  1:1058579/140‑1 (MQ=255)                                             cGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCa                                                                                   <  2:1331303/140‑1 (MQ=255)                                             cGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCa                                                                                   <  1:1101794/140‑1 (MQ=255)                                                 cAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGCGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGct                                                                               <  2:489474/140‑1 (MQ=255)                                                  aaTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGctc                                                                              <  1:94564/140‑1 (MQ=255)                                                  aaTAAAGTATGTGGCTCGGGTCTTAAAAGTGTGGCGCTTGTCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGctc                                                                              >  1:647018/1‑140 (MQ=255)                                                       aGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGtt                                                                         <  1:506992/140‑1 (MQ=255)                                                         tATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTAt                                                                       <  1:9820/140‑1 (MQ=255)                                                          aTGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACg                                                                                 >  1:884786/1‑129 (MQ=255)                                                          aTGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACg                                                                                 <  2:884786/129‑1 (MQ=255)                                                              ggTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGG‑TGGAAAATATGAGGTTTGCCCCCCACTTACTCGATGCAAAAAGACGCTCTGGTTATTGTCt                                                                  >  2:950008/1‑140 (MQ=255)                                                              ggTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCt                                                                  >  1:1352692/1‑140 (MQ=255)                                                                ttCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGA‑TGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTg                                                                >  1:452120/1‑140 (MQ=255)                                                                   gggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGA‑TGGGAAAAATGAGGTTAGCCCCCCCCTTACTCCAGGCAAAAACACACCCTGGTTATCGTCtttggg                                                             >  2:1386373/1‑136 (MQ=255)                                                                    ggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGgaga                                                            <  2:1352692/140‑1 (MQ=255)                                                                     gTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGTATTGTGGCGGGGGGGATTGGAAAAAAGGGGTTTTGGCCCCT‑CTTCCTCGCTACCAAAACACCCGTTCGTTATAGGCTTGgaggc                                                           >  2:1400688/1‑138 (MQ=255)                                                                         tAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGAc                                                       <  2:266071/140‑1 (MQ=255)                                                                          aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGA‑TGGAAAAAATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGGGACGGACa                                                      >  1:391115/1‑140 (MQ=255)                                                                           aaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAg                                                     <  2:608007/140‑1 (MQ=255)                                                                                gtgGCGCTTGTCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTa                                                <  2:647018/140‑1 (MQ=255)                                                                                gtgGCACTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTa                                                <  1:1400688/140‑1 (MQ=255)                                                                                 tgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTAt                                               >  1:844306/1‑140 (MQ=255)                                                                                 tgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGA‑TGGAAAAAATGAGTTTAGCCCCCTCCTTACTCCAAGGAAAAACGCGCCCTGGTTATCGTCTTGGAGACGGACAGGTTTAt                                               >  2:1141852/1‑140 (MQ=255)                                                                                 tgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGA‑TGGAAAAAATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATTGTCTTGGGGACGGACAGGTTTAt                                               >  2:745896/1‑140 (MQ=255)                                                                                    cgcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGA‑TGGGAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGGTTATGAc                                            >  2:1133310/1‑140 (MQ=255)                                                                                               cAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTgcg                                 <  1:512864/140‑1 (MQ=255)                                                                                                 ggCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGcgcg                               <  2:844306/140‑1 (MQ=255)                                                                                                       tCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGA‑TGGGAAAAATGAGTTTAGCCCCCCACTTACTCCATGCAAAAACACGCTCTGGTTATCGTCTTGGGGACGGACAGGGTTTTGAACTAATCCTTCGCGATGGcc                         >  2:969722/1‑140 (MQ=255)                                                                                                               gTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTgc                 <  1:745896/140‑1 (MQ=255)                                                                                                               gTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTgc                 <  1:1255524/140‑1 (MQ=255)                                                                                                                tCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTgcg                <  2:355556/140‑1 (MQ=255)                                                                                                                tCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTgcg                <  1:1141852/140‑1 (MQ=255)                                                                                                                tCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTgca                <  1:407920/140‑2 (MQ=255)                                                                                                                    gcgcAGAGCATTGTGGCGGGGGGGA‑TGGAAAAAATGAGGTTAGCCCCCCCCTTACTCGATGCAAAAACACGCCCTGGTTATCGTCTTTGAGACGGACAGGGTTATTACGGAAAACCGCGCGGTGGCCCGGTGTGcaccac            >  2:407920/1‑140 (MQ=255)                                                                                                                     cgcAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCAcc           <  1:969722/140‑1 (MQ=255)                                                                                                                     cgcAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCAcc           <  1:368034/140‑1 (MQ=255)                                                                                                                       cAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCa         <  2:1161924/140‑1 (MQ=255)                                                                                                                        agagCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCAt        <  2:686572/140‑1 (MQ=255)                                                                                                                           gCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGt     <  2:352853/140‑1 (MQ=255)                                                                                                                            cATTGTGGCGGGGGGGA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGtt    >  1:1315362/1‑140 (MQ=255)                                                                                                                              ttGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTAt  <  1:950008/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                GGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTA‑TGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTAT  >  NC_000913/2326312‑2326576

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑