Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	701,461	(G)5→6	coding (134/1221 nt)	nagC ←	N‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	701,456	1	.	G	95.7%	158.3 / ‑0.5	47	coding (139/1221 nt)	nagC	N‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (24/21);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGC‑GG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCATACTTTCT  >  NC_000913/701316‑701601                                                                                                                                               |                                                                                                                                                   aaGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGC‑gg‑ggg                                                                                                                                              >  2:685623/1‑146 (MQ=255)              aTTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCCAGCTGGCTTTGCTCg                                                                                                                             >  2:774674/1‑151 (MQ=255)                     ggAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCt                                                                                                                      <  2:836544/151‑1 (MQ=255)                     ggAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCt                                                                                                                      >  2:295981/1‑151 (MQ=255)                      gAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTg                                                                                                                     <  2:792998/151‑1 (MQ=255)                           tGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATcc                                                                                                                >  2:245040/1‑151 (MQ=255)                            gCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATcc                                                                                                                <  1:68276/150‑1 (MQ=255)                                 tttCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCcgc                                                                                                              >  1:645816/1‑147 (MQ=255)                                    cGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGAtt                                                                                                        >  1:166529/1‑150 (MQ=255)                                    cGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGAtt                                                                                                        >  2:368792/1‑150 (MQ=255)                                           gATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGt                                                                                                 >  1:157934/1‑150 (MQ=255)                                            aTGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTAc                                                                                               <  1:579030/151‑1 (MQ=255)                                            aTGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTAc                                                                                               <  1:61627/151‑1 (MQ=255)                                                 gaTAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTc                                                                                          <  1:469934/151‑1 (MQ=255)                                                 gaTAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCACGAGATTGGCCCGTACTgg                                                                                            >  2:530547/1‑149 (MQ=255)                                                      cgcgGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATc                                                                                      >  2:590907/1‑150 (MQ=255)                                                            ggcCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGAt                                                                               >  1:836186/1‑151 (MQ=255)                                                                                    cAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCCCGCTGTTAAGCTGtt                                                       >  1:101724/1‑151 (MQ=255)                                                                                     aaCTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCgaga                                                                                                          >  2:526921/1‑99 (MQ=255)                                                                                     aaCTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCgaga                                                                                                          <  1:526883/99‑1 (MQ=255)                                                                                     aaCTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGttt                                                      >  2:684543/1‑151 (MQ=255)                                                                                          cttTGATCAGCGCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTAc                                                   <  1:754323/149‑1 (MQ=255)                                                                                          cttTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGGCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACg                                                  >  2:715318/1‑150 (MQ=255)                                                                                             tGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACgaga                                               <  1:279666/150‑1 (MQ=255)                                                                                                tCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCa                                            <  1:367259/150‑1 (MQ=255)                                                                                                   gCCCGCTTTCGTTAGGCTGACGCGTAATTTTGGTTACGCTGGGG‑GG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAACCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATc                                             <  1:453591/146‑1 (MQ=255)                                                                                                   gCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACAt                                        >  1:367095/1‑151 (MQ=255)                                                                                                    cccGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACAtt                                       <  2:695394/151‑1 (MQ=255)                                                                                                      cgcgTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTAc                                     >  2:439454/1‑151 (MQ=255)                                                                                                        cgTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCt                                   >  2:304919/1‑151 (MQ=255)                                                                                                         gTTCGATAAGCTGACGGGTAATTTTGGTTACCCT‑GGC‑GGTGGGCAAGCTGGCTTTGTTGGGCAATTTGAATCCGCCAAATTGGCCCGTACTGGTCATTCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTTCCt                                   <  2:614523/150‑1 (MQ=255)                                                                                                         gTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGGTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCt                                   >  1:421933/1‑150 (MQ=255)                                                                                                          ttCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTAt                                 >  2:822691/1‑151 (MQ=255)                                                                                                           tCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTAt                                 <  2:716782/150‑1 (MQ=255)                                                                                                              aTAAGCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGa                             <  2:607037/151‑1 (MQ=255)                                                                                                                  gCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCtt                         <  2:661332/151‑1 (MQ=255)                                                                                                                  gCTGACGCGTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCtt                         >  2:213033/1‑151 (MQ=255)                                                                                                                       cgcgTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTcc                     <  2:376966/150‑1 (MQ=255)                                                                                                                       cgcgTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCg                    <  1:400908/151‑1 (MQ=255)                                                                                                                         cgTAATTTTGGTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGGTAAGCTGTTTTACGAGATCAACCTTACCTATCTGAGCTTGTCCGcc                  >  1:322223/1‑151 (MQ=255)                                                                                                                                  ggTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCa          <  1:598743/150‑1 (MQ=255)                                                                                                                                  ggTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCa          <  1:295956/150‑1 (MQ=255)                                                                                                                                  ggTTACGCT‑GGCGGG‑GGGAAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCa          <  2:838331/150‑1 (MQ=255)                                                                                                                                   gTTACGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCAt         >  1:272172/1‑150 (MQ=255)                                                                                                                                       cGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCATACttt    <  1:310023/151‑1 (MQ=255)                                                                                                                                       cGCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCATACttt    <  2:504685/151‑1 (MQ=255)                                                                                                                                        gCT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGGACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCATACttt    >  2:316737/1‑150 (MQ=255)                                                                                                                                         cT‑GGCGGG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGAGTCATACTTtct  >  2:645723/1‑151 (MQ=255)                                                                                                                                               |                                                                                                                                                   AAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCT‑GGC‑GG‑GGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCGCTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCATACTTTCT  >  NC_000913/701316‑701601

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