Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,326,450	0	T	G	42.6%	97.5 / 10.0	54	G114G (GGT→GGG)	atoB	acetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (4/27);  new base G (23/0);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 3.23e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAA  >  NC_000913/2326303‑2326597                                                                                                                                                    |                                                                                                                                                    ggCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATg                                                                                                                                                  <  2:791381/151‑1 (MQ=255)    gcgTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGaa                                                                                                                                               <  1:2609/151‑1 (MQ=255)      gTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGaaaa                                                                                                                                             <  2:683922/151‑1 (MQ=255)      gTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGaaa                                                                                                                                              <  1:137215/150‑1 (MQ=255)                gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTa                                                                                                                                   <  2:158211/151‑1 (MQ=255)                  tAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAg                                                                                                                                  <  1:270225/150‑1 (MQ=255)                        gCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGccccc                                                                                                                             >  1:720091/1‑149 (MQ=255)                           ggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtactt                                                                                                                        <  1:351951/151‑1 (MQ=255)                             cTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCcactt                                                                                                                        >  2:324541/1‑149 (MQ=255)                                  aGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGGCCCcttctttc                                                                                                                      >  1:256302/1‑144 (MQ=255)                                    aaaCGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTCCTTACTCGATg                                                                                                                >  2:944009/1‑150 (MQ=255)                                      aCGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGCGGGGATGGGAAAAATGTGGTTAGCCCCCTACTTACTCTCTGCa                                                                                                              >  2:1037851/1‑150 (MQ=255)                                        ggTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTCCTCGAt                                                                                                                 <  2:460209/145‑1 (MQ=255)                                        ggTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTATCCCCCTACTTACTCGATGCaaa                                                                                                            >  2:258457/1‑150 (MQ=255)                                         gtgtGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTTGCCCCCTACTTACTCGATGCaaa                                                                                                            >  1:257367/1‑149 (MQ=255)                                         gtgtGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGACCCCTCCTTACTCGATGCaaaa                                                                                                           >  1:43363/1‑150 (MQ=255)                                          tgtgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGGTTAGGCCCCTCCCTCCTCGATGCa                                                                                                              >  2:1119720/1‑146 (MQ=255)                                             gCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGTCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACg                                                                                                      <  1:1156704/151‑1 (MQ=255)                                                       ggTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGGTAt                                                                                            >  1:619989/1‑151 (MQ=255)                                                        gTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTTCTCGATGCAAAAGCACGCTCTGGTTAt                                                                                            >  2:844338/1‑150 (MQ=255)                                                           aaTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTc                                                                                        <  1:825741/151‑1 (MQ=255)                                                              aaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCtt                                                                                      >  1:159498/1‑150 (MQ=255)                                                                aGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGga                                                                                   <  2:161903/151‑1 (MQ=255)                                                                aGTATGTGGGTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGGGGCAAACATCAATTTGTACCCCTACTTAGTGGCTGGAAAAACAAGTTCTTGtttttgtttt                                                                                      >  2:544930/1‑140 (MQ=255)                                                                    tgtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGGTATCGTCtttg                                                                                    >  1:813936/1‑144 (MQ=255)                                                                       ggTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTGGGGGAAAAATGGGGTTAGCCCCCCACTTTCTTGTTGCAAAAAGACACTCTGGGTTTCCTCtttggg                                                                                  >  1:607582/1‑141 (MQ=255)                                                                        gTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTTGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGGTATCGTCTTGGAGACGGAc                                                                            >  2:1152337/1‑150 (MQ=255)                                                                          tCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGAGACGGACAgg                                                                         >  1:547003/1‑151 (MQ=255)                                                                            gggTCTTAAAAGTGTGGCGCTTGCCGGCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGtt                                                                       <  1:780751/151‑1 (MQ=255)                                                                               tCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGACAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTTGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGAGACGGACAGGgtt                                                                      >  2:327533/1‑149 (MQ=255)                                                                                cTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGCTCTGGGTATCGTCTTTGAGACGGACAGGgtt                                                                      >  1:400553/1‑148 (MQ=255)                                                                                  tAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTAt                                                                    <  2:67422/148‑1 (MQ=255)                                                                                   aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGAc                                                                 <  1:1004334/150‑1 (MQ=255)                                                                                     aaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTa                                                              <  2:859680/151‑1 (MQ=255)                                                                                        tgtgGAGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTGTGGGAAATATGAGGGTAGCCCCCTAATTAGTTGCGGCAAAAACAAGCTCTCGGTATCGTCTTTGGGATGAACAGGGGTATGACGTAAt                                                            >  1:1035584/1‑150 (MQ=255)                                                                                             cgcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTgc                                                       <  1:905071/150‑1 (MQ=255)                                                                                              gcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTgcgc                                                     <  1:266501/151‑1 (MQ=255)                                                                                                  gccgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATg                                                 <  2:275781/151‑1 (MQ=255)                                                                                                    cgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGc                                               <  2:11728/151‑1 (MQ=255)                                                                                                        cAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGa                                           <  1:1136688/151‑1 (MQ=255)                                                                                                                tCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTTCTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGGGACGGACAGGGTTATGACGTAATTCTGCGCGATTGCCTGATGGGCGcc                                   >  1:851738/1‑151 (MQ=255)                                                                                                                 caggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGGTTAGCCCCCCACTTACTCGATGCAAAAGAACGCCCTGGTTATCGTCTTGGAGACGGACAGGGTTATGACGGAATCCTGCGCGCTGGCCTGATGGGGGcc                                   >  2:193656/1‑150 (MQ=255)                                                                                                                    gcaggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCAccc                               <  1:257543/151‑1 (MQ=255)                                                                                                                    gcaggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTTTCGTCTTGGGGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCAccc                               >  2:938305/1‑151 (MQ=255)                                                                                                                        gTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATgg                           <  2:41847/151‑1 (MQ=255)                                                                                                                           aGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTa                        <  2:633091/151‑1 (MQ=255)                                                                                                                              cgcAGAGCATTGTGGCGGGGGGGATTGAAAATATGAGGTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGAAGGACAGGGTTATGACGTAATCCTGCGCGCTGGCCTGATGTGCGCCACCCATgggt                         >  1:977000/1‑145 (MQ=255)                                                                                                                              cgcAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGGTATCGTCTTGGGGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCCGGTGTGCGCCACCCAAGGGTAt                       >  2:361160/1‑149 (MQ=255)                                                                                                                                cAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCata                   <  1:656514/151‑1 (MQ=255)                                                                                                                                         gtgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGAtt           <  2:115866/150‑1 (MQ=255)                                                                                                                                         gtgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTa          <  1:247126/151‑1 (MQ=255)                                                                                                                                         gtgGCGGGGGGGGTGGGAAATATGAGTTTAGCCCCCTACTTACTCGGTGCAAAAACACGCTCTGGGTATTGTCTTTGGGACGGACAGGGTTATGACGGAATCCTGCGCGGTGGCCTGATGTGCGCCACCCCAGGTTATCATATGGGGGtt           >  2:238241/1‑150 (MQ=255)                                                                                                                                           ggCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAcc        >  2:267496/1‑151 (MQ=255)                                                                                                                                                 gggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGaa  <  1:735908/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAA  >  NC_000913/2326303‑2326597

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