Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,279,468	0	C	G	45.5%	70.5 / 11.2	44	G330G (GGC→GGG)	ghxP	guanine/hypoxanthine permease, high affinity; guanine/hypoxanthine:H+ symporter
Reads supporting (aligned to +/- strand):  ref base C (3/21);  new base G (20/0);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.02e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATG‑CAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAACGTGGCGAAAATC  >  NC_000913/4279327‑4279618                                                                                                                                                |                                                                                                                                                          gAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGG‑GGG‑GGGGGGGaaaaa                                                                                                                                                   >  2:499445/1‑149 (MQ=255) gAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGC‑GGG‑GGGGGGGaaaa                                                                                                                                                    >  1:197399/1‑149 (MQ=255)       gctgGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGttt                                                                                                                                             <  1:379706/150‑1 (MQ=255)       gctgGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGtt                                                                                                                                              <  2:100196/149‑1 (MQ=255)         tgGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGGGGGGAAACAGGgttt                                                                                                                                            >  1:566582/1‑148 (MQ=255)          ggATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGAGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGGCGGA‑GGG‑GGGGGGGGGAACCggg                                                                                                                                               >  1:296996/1‑144 (MQ=255)                aGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATc                                                                                                                                   <  2:221728/151‑1 (MQ=255)                    gggCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGC‑GGCGGGGGGGGAAACCCGGTTGGAGGGccttacc                                                                                                                                 >  1:366551/1‑144 (MQ=255)                      gCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGtt                                                                                                                           <  2:766821/151‑1 (MQ=255)                          atcatcGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGGAAAAACGGTTTGTCGGCCATC‑A‑CCgcttggg                                                                                                                       >  2:802741/1‑145 (MQ=255)                           tcatcGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGGGGGAAAAACGGGTTGACGGGTTTC‑A‑CCGGTGGCg                                                                                                                       >  2:350666/1‑150 (MQ=255)                               cGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTg                                                                                                                  <  2:609402/151‑1 (MQ=255)                                      gggAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTtcctcc                                                                                                           <  2:783788/151‑1 (MQ=255)                                      gggAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTtcctcc                                                                                                           <  1:55023/151‑1 (MQ=255)                                        gAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTg                                                                                                         <  2:970704/151‑1 (MQ=255)                                            gCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTc                                                                                                     <  1:568213/151‑1 (MQ=255)                                             cactgaccactgacTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGCGGGGGGGGAAAACCGGTTTG‑CGGGTATT‑A‑CCCTTGGGGGGCTGTTTCTCCTGGtttt                                                                                                    >  2:699430/1‑149 (MQ=255)                                               ctgaccactgacTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGCGGGGGGGGAAAACCGGTTT‑ACGGCTCTC‑A‑CCCTTGGCGGGCTGGTTCTCCTGGttt                                                                                                     >  2:374864/1‑147 (MQ=255)                                                  accactgacTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGAAGGC‑GGG‑GGGGGGGGAAAACGGGTTGAAGGGGCTC‑A‑CCCTTTGGGTGGTGTTCCTCCTGATTTTGtt                                                                                                 >  2:424335/1‑149 (MQ=255)                                                    cactgacTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCctc                                                                                             <  2:621870/151‑1 (MQ=255)                                                         acTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGC‑GGG‑GGGGGGGAAAAACGGGTTGAGGGGTATC‑A‑CCCGTGGCGGGCTGTTTCTCCTGATTTTTTTTCTCTCTcc                                                                                        >  2:139527/1‑151 (MQ=255)                                                          cTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAAGGC‑GGG‑GGGGGGGAAAAACGGGTTTGCGGGTATC‑A‑CCCTTTGGGGGCTGTTTCTCCTGGTTTTGTTTccctct                                                                                          >  2:114361/1‑148 (MQ=255)                                                              aTGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTTAGGC‑GGC‑CGGGGGTTAAACCCGGTTTACCGGTAAT‑A‑CCCCTTGCGGGCTTTTTCTTCCGATTCTGTTTctcttttccct                                                                                     >  2:660153/1‑142 (MQ=255)                                                              aTGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGCGGGGGGGAAAAACCGGTTT‑ACGGCCATC‑A‑CCCGTGGCGGGCTGTTTCTCCTGATTTTTTTTCTCTCTcc                                                                                        >  2:491353/1‑146 (MQ=255)                                                              aTGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACCGG‑GGG‑GGGGGGGAAAAACGGGTTGTCGGCTATC‑T‑ACCCTGGGGCGTGGTTCCTCCCGATTCTTTTTct                                                                                              >  1:793194/1‑140 (MQ=255)                                                                gAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGctctc                                                                                  <  2:869277/150‑1 (MQ=255)                                                                    agcGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGGAGGC‑GGCGGGGGGGAAAAACCGGTTTTGCGGC‑ATC‑C‑CCCGTGGGGGGCTGTTCCTCCCCATTTTCTTTCTCTCCCCCCTCTCttttc                                                                             >  1:960858/1‑148 (MQ=255)                                                                           tctcTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGC‑GGCGGGGGGGAAAAACCGGTTTG‑CGGCTATC‑A‑CCCTTTGGGTGCTGTTTCTCCCGATTTTTTTTCTCTCTCCCCCCTCtt                                                                                >  2:72416/1‑141 (MQ=255)                                                                             tctGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGGCGGG‑GGC‑GGGGGGGAAAAAAGGGTTTAAGGGGATT‑A‑CCCGTGGGGGGGTGGTTCTCCCGGTTTCGTTTTTCTCTCCGCTCTCTTTcc                                                                             >  1:796151/1‑142 (MQ=255)                                                                                        gggtgCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATG‑CAAc                                                         <  2:210348/150‑1 (MQ=255)                                                                                          gtgCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATG‑CAACGGc                                                      <  1:704666/151‑1 (MQ=255)                                                                                                  ccGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATG‑CAACGGCTCCGgcgc                                              <  2:341247/151‑1 (MQ=255)                                                                                                  ccGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGC‑GGC‑GGGGGGGGAAACCGGGTTTGCGGCCTTC‑ACCCGTTGG‑GGGCCGTTTCTCCCCGTTTTTTTTCTCTCTCCCCCCTCTTTCCTCCTTCCGGGGGGTG‑CAACGGcccc                                                   >  1:376458/1‑146 (MQ=255)                                                                                                     gCAGCGGTATACATCGAGTCTGCGGCGGGGACGGC‑GGC‑GGGGGGGAAAAACGGGTTTGCGGGTATCAA‑CCGTTGG‑GTGCTGGTTCTCCTTATTTTTTTTCCCTCTCCGCCCTCTTTCCTCCCTTCGGGGGATG‑CAACCGCCCCGgcg                                               >  2:1015287/1‑147 (MQ=255)                                                                                                      cAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATG‑CAACGGCTCCGGCGCTGAt                                          <  2:701821/151‑1 (MQ=255)                                                                                                        gCGGTATACATCGAGTCTGCGGCGGGGACGGG‑GGG‑GGGCGGTAAAAACGGGTTTACGGCTATC‑A‑CCGTTGGCGGGCTGTTCCTCCCGATTCTGTTTCTCTCTCCCCCCTCCTTCCTCGTTCCCGGGGATG‑CAACGGCTCCCGCGCTGgtgt                                        >  2:651777/1‑151 (MQ=255)                                                                                                           gTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCCGCTTTGCGGCCTTC‑A‑CCCCTGGCGGGCTGTTTCCCCCGCTTCTTTTTCTCTCCCCCCTCCCTTCCCCCCTTCCCGGGGGTG‑CCACCGCCCCCGCGCTGGTGTAc                                      >  1:307227/1‑150 (MQ=255)                                                                                                             ataCATCGAGTCTGCGGCGGGGACGGCGGGG‑GGGGGGAAAAACCGGTTTG‑CGGCCGTC‑ACCCGTTGG‑GGGCTGGTTCTTCTGGTTTTGTTTCCCTCTCCCCCCTCCTTCCTCTTTCCCGGGGGTGCCAACGGccccccccttgtttaccct                                    >  1:653760/1‑136 (MQ=255)                                                                                                                   cGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATG‑CAACGGCTCCGGCGCTGATGTACGTTGGCct                              <  2:15274/150‑1 (MQ=255)                                                                                                                          gcggcgGGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATG‑CAACGGCTCCGGCGCTGATGTACGTTGGCCtgctgatg                       <  1:504613/150‑1 (MQ=255)                                                                                                                            ggcggGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATG‑CAACGGCTCCGGCGCTGATGTACGTTGGCCtgctgatgctg                    <  2:331984/151‑1 (MQ=255)                                                                                                                                 gTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATG‑CAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCaa               <  1:147088/151‑1 (MQ=255)                                                                                                                                    cggc‑ggc‑ggGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATG‑CAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAACGt            <  1:809858/151‑1 (MQ=255)                                                                                                                                     ggc‑ggc‑ggGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATG‑CAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAACGTg           <  1:582213/151‑1 (MQ=255)                                                                                                                                           c‑gggcggTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATG‑CAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAACGTGGCGa       >  1:437764/1‑150 (MQ=255)                                                                                                                                             gggcggTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATG‑CAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAACGTGGGGaaa     <  1:390422/151‑1 (MQ=255)                                                                                                                                               gcggTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATG‑CAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAACGTGGCGAAAAt   >  1:681778/1‑151 (MQ=255)                                                                                                                                                cggTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATG‑CAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAACGTGGCGAAAATc  >  2:718839/1‑151 (MQ=255)                                                                                                                                                cggTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATG‑CAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAACGTGGCGAAAATc  <  2:201224/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                          GAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGC‑GGC‑GGGCGGTAAAACCGGTTTGACGGCTATC‑A‑CCGTTGGCGTGCTGTTCCTCCTGATTCTGTTCCTCTCTCCGCTCTCTTACCTCGTTCCGGGGTATG‑CAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAACGTGGCGAAAATC  >  NC_000913/4279327‑4279618

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