Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,326,450	0	T	G	38.5%	96.6 / 10.6	52	G114G (GGT→GGG)	atoB	acetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/31);  new base G (19/1);  total (20/32)
Fisher's exact test for biased strand distribution p-value = 5.09e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAA  >  NC_000913/2326301‑2326598                                                                                                                                                      |                                                                                                                                                       ccGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGt                                                                                                                                                        <  1:1071752/150‑1 (MQ=255)    gcgcGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATgg                                                                                                                                                    <  1:404901/151‑1 (MQ=255)        gTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGaaaa                                                                                                                                                <  1:862718/151‑1 (MQ=255)          cAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAat                                                                                                                                               <  2:3887/150‑1 (MQ=255)          cAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGgaaa                                                                                                                                                >  1:965312/1‑149 (MQ=255)           aGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGgaaa                                                                                                                                                >  2:150603/1‑148 (MQ=255)                 tGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGA‑Ggtt                                                                                                                                      >  2:444185/1‑151 (MQ=255)                  gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑Gttt                                                                                                                                      <  1:159480/150‑1 (MQ=255)                  gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑Gtt                                                                                                                                       <  2:857290/149‑1 (MQ=255)                   ttAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGA‑Gttt                                                                                                                                      >  1:188929/1‑149 (MQ=255)                                  gCAGAAACGGTGTGTGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCg                                                                                                                     <  2:101704/151‑1 (MQ=255)                                   cAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAATATGA‑GGTTTGCCCCCcactt                                                                                                                          >  2:540556/1‑145 (MQ=255)                                    aGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGAt                                                                                                                   <  1:66285/151‑1 (MQ=255)                                     gAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATg                                                                                                                  >  1:628241/1‑151 (MQ=255)                                          ggTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGA‑GTTTAGCCCCCCACTTACTCGATGCaaa                                                                                                              >  1:104590/1‑150 (MQ=255)                                                cGGATTGACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCACTTGCCGCCCAGGCCAATCAGGCAGGTCAGGCGCAGAGAATTGTGGCGTGGGGTATAGAAAATATTA‑GTTTACCCCCCTACATACTCGATGCaaaatcccg                                                                                                        <  2:498229/150‑6 (MQ=255)                                                cGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGc                                                                                                       <  1:1199387/151‑1 (MQ=255)                                                 ggATTCACGGTCAATGAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCATGTCAGGCGCAGAGCATTGTGGCGGAGGGTATCGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGct                                                                                                      <  2:312759/151‑1 (MQ=255)                                                     tCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTg                                                                                                   <  1:100540/150‑1 (MQ=255)                                                       aCGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGtt                                                                                                <  1:689201/151‑1 (MQ=255)                                                          gTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATc                                                                                             <  1:988354/151‑1 (MQ=255)                                                             aaTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑c                                                                                         <  1:679615/151‑1 (MQ=255)                                                                 aaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑Ctt                                                                                       <  2:141196/149‑1 (MQ=255)                                                                   gTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGgag                                                                                   <  2:701858/151‑1 (MQ=255)                                                                     aTGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGgaga                                                                                  <  1:296143/150‑1 (MQ=255)                                                                       gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGGAGACgg                                                                               <  1:1176296/151‑1 (MQ=255)                                                                       gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGGAGACg                                                                                <  1:421908/150‑1 (MQ=255)                                                                             cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGGAGACGGACAGGt                                                                         <  1:1203349/151‑1 (MQ=255)                                                                                  cTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGA‑GGTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGGAGACGGACAGGgtt                                                                       >  1:841070/1‑148 (MQ=255)                                                                                  cTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGATCAGGCGCAGAACATTGTGGGGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGGAGATGGACAGGTTTa                                                                      <  2:1089568/149‑1 (MQ=255)                                                                                    tAAAAGTGTGGCGCTTGCCGGCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTGGGAGACGGACAGGTTTATGa                                                                   <  2:1274889/150‑1 (MQ=255)                                                                                    tAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGG‑GTTTAGCCCCCTACTTCCTCGATGCAAAAGCACACCCTGGTTATCGT‑CTTGGAGAAGGACAGGGTTAt                                                                     >  1:455956/1‑148 (MQ=255)                                                                                     aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGTGTATGGAAAATATGA‑GTTTAGCCGCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGGAGTCGAACAGGTTTATgcc                                                                  <  2:434495/150‑3 (MQ=255)                                                                                        aGTGAGGCGCTTGCCGCCCAGGCCATTCAGGCAGGCCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGGTTTAG‑CCCCTACTTAATCGATGCAAAAAGACGCTCTGGTTATCGT‑GTTGGAGAAGGACAGGGTTATTACGTa                                                               >  2:733602/1‑150 (MQ=255)                                                                                         gtgtGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGGAGACGGACAGGTTTATGACGTAAt                                                             <  1:407518/151‑1 (MQ=255)                                                                                              gctcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATAG‑GTTTAGCCCCCCACTTACTCGATGCCAAAGGACGCTCTGGTTCTTAT‑CTTGGGGACGGAAAGGGTTTTGACGTa                                                               >  2:242173/4‑144 (MQ=255)                                                                                                   tGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGGAGACGGACAGGTTTATGACGTAATCCTGcgcg                                                     <  1:1015279/149‑1 (MQ=255)                                                                                                   tGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGGAGACGGACAGGTTTATGACGTAATCCTGcgcg                                                     <  2:210654/149‑1 (MQ=255)                                                                                                   tGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGA‑GTTTAGCCCCCTACTTTCTCGGTGCAAAAACACGCTCTGGTTTTCGT‑CTTGGAGACGGACAGGGTTACCACGTaaaactgggc                                                      >  1:402477/1‑140 (MQ=255)                                                                                                       gccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGA‑GGTTTGCCCCCTACTTTCTCGATGCAAAAGCACGCTCTGGTTATTGT‑CTTTGAGACGGACAGGTTTATGACGTAATCCTGCGCGATgg                                                 >  1:850236/1‑149 (MQ=255)                                                                                                         ccAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTg                                             <  1:885283/151‑1 (MQ=255)                                                                                                             gCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGAtgt                                         <  2:1256959/151‑1 (MQ=255)                                                                                                              ccATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTTGGGACGGACAGGTTTATGACGGAATCCTTCGCGCTGGCCTGAtgtg                                        >  2:1270677/1‑151 (MQ=255)                                                                                                               cATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGA‑GGTTAGCCCCCCACTTACTCGATGCAAAAACACGCTCTGGGTATCGT‑CTTGGAGACGGACAGGGTTATGACGTAATCCTGCGCGATGGCCTGAtgtg                                        >  1:703519/1‑150 (MQ=255)                                                                                                               cATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGA‑GGTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGT‑CTTTGAGACGGACAGGGTTATGACGTAATCCTGCGCGATGGCCTGAtgtg                                        >  1:10016/1‑150 (MQ=255)                                                                                                                         ggTCAGGCGCAGAGCATTGTGGCGGTGGGGGCGGAAAAAATGG‑ACTTAGCCCCCCACTTTATCCATGCAAAAGCACCCGCCGGGTTTCGA‑TCTGGCAACGGGCAAGGTTTTTAAGTCAACCTGCGCGATAGCCTGAtggtcggcaacac                               >  2:751769/1‑138 (MQ=255)                                                                                                                            cAGGCGCAGTGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTAGTTATGCGATGCAAAAGCACGCTCTGGTTATCGT‑GTTGTAGACGGACAGGTTTATGAGGTAAACCTGCGCGATGCCCTGATGTgcgcgtcgcctgg                            <  2:874802/149‑10 (MQ=255)                                                                                                                              ggCGCAGAGCATTGGGGCGAGGGGGATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCGCGCTCTGGTTATCGG‑CTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGGCGGATGTGCGCCAGCCATGGTTa                         <  2:188856/150‑1 (MQ=255)                                                                                                                                   agagCATTGTGGCGGGGGGGATGGAAAATATTA‑GTTTAGGCCCCTACTTAATCCATGCAAAAGCACGCTCTGGTTATTGTGCTTGG‑GAAGGGCACGGTTATAACCGTATTCTGCGCGATGGGCTGATGTGCGCCCCCCCAGGTTTTc                       >  2:622834/1‑147 (MQ=255)                                                                                                                                       cATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATTGT‑CTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATgggg               <  1:676564/151‑1 (MQ=255)                                                                                                                                       cATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATgggg               <  1:1103833/151‑1 (MQ=255)                                                                                                                                        aTTGTGGCGGGGGGGATGGAAAAAATGA‑GTTTAGGCCCCTACTTACTCGCTGCAAAAAGAAGCTCTGGTTATTGG‑CTTGGAGAAGGGCAGGGTTTTTAAGTAATCCCGCGCGATGGCCTGATGTGCGCCACCCATGGGTATCATATggggg              >  2:1089746/1‑150 (MQ=255)                                                                                                                                                  ggggTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccgccg     >  2:916482/1‑150 (MQ=255)                                                                                                                                                     ttATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGaaa  >  2:738875/2‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                       CCGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGA‑GTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGT‑CTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAA  >  NC_000913/2326301‑2326598

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