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	97.5%	137.0 / ‑2.5	40	coding (139/1221 nt)	nagC	N‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (18/21);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

ATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGC‑GGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCATACTTT  >  NC_000913/701307‑701599                                                                                                                                                       |                                                                                                                                                  aTGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGC‑g                                                                                                                                                  >  2:553795/1‑151 (MQ=255)     cgacCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGC‑gggg                                                                                                                                               <  2:365078/150‑1 (MQ=255)             ccGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTg                                                                                                                                       >  2:735132/1‑151 (MQ=255)             ccGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTg                                                                                                                                       <  2:633756/151‑1 (MQ=255)             ccGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTg                                                                                                                                       <  2:454020/151‑1 (MQ=255)             ccGTACGCCGATTGCGTGGAAATTGCGGGTTGCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTg                                                                                                                                       <  1:241245/151‑1 (MQ=255)                   gCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATGAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCttt                                                                                                                                  >  2:497341/1‑150 (MQ=255)                         tGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCgg                                                                                                                            <  1:231739/150‑1 (MQ=255)                          gCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGc                                                                                                                           <  1:553743/150‑1 (MQ=255)                           cGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCaa                                                                                                                         <  2:741975/151‑1 (MQ=255)                                     gCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCcg                                                                                                               <  1:1430/151‑1 (MQ=255)                                            tCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑Att                                                                                                       >  1:120275/1‑151 (MQ=255)                                             cGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTg                                                                                                      <  2:341533/151‑1 (MQ=255)                                                 gACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGccc                                                                                                  <  2:86766/151‑1 (MQ=255)                                                    gATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGt                                                                                                >  1:394564/1‑150 (MQ=255)                                                          gaTAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTgg                                                                                           >  1:160218/1‑149 (MQ=255)                                                           aTAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCa                                                                                        <  2:530916/151‑1 (MQ=255)                                                                gcggcgGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCa                                                                                    <  1:696404/150‑1 (MQ=255)                                                                       cccccGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATaa                                                                            <  2:644428/151‑1 (MQ=255)                                                                       cccccGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATaa                                                                            <  1:300940/151‑1 (MQ=255)                                                                         cccGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAAc                                                                          <  2:144073/151‑1 (MQ=255)                                                                           cGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAAcc                                                                         >  2:170194/1‑150 (MQ=255)                                                                              tGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAAccgcc                                                                      >  1:727517/1‑150 (MQ=255)                                                                                    ccTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAGGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGt                                                              <  2:31815/151‑1 (MQ=255)                                                                                       gctgATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGGAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGGCAATCAGGCGATAAAACGCCGCG‑CTGTTa                                                            >  1:294282/1‑150 (MQ=255)                                                                                          gATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGc                                                         <  1:32536/150‑1 (MQ=255)                                                                                           aTCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGCTg                                                       >  1:677353/1‑151 (MQ=255)                                                                                           aTCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGCTg                                                       <  2:562242/151‑1 (MQ=255)                                                                                                 ttcttTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGCTGtttt                                                   >  1:14578/1‑149 (MQ=255)                                                                                                        aTCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGCTGTTTTACGAGATCa                                          <  2:142089/151‑1 (MQ=255)                                                                                                         tCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGCTGTTTTACGAGATc                                           >  2:490393/1‑149 (MQ=255)                                                                                                         tCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGCTGTTTTACGAGATCa                                          >  1:187936/1‑150 (MQ=255)                                                                                                           aGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGC‑GGGGGGCAACTGGGTTTTGTCCGCCATTTGAAACCGCGAGAATGGGCCCTAACGGGCCAACAAGGGGATAAA‑CGCCGCGCCTGTTAA‑CTGGTTTTAGAGAttca                                         >  2:702031/1‑145 (MQ=255)                                                                                                            gCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGCTGTTTTACGAGATCAACAt                                      >  2:314841/1‑151 (MQ=255)                                                                                                             cccGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGCTGTTTTACGAGATCAACAtt                                     <  2:387819/151‑1 (MQ=255)                                                                                                               cgcgTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGCTGTTTTACGAGATCAACATTAc                                   <  2:487883/151‑1 (MQ=255)                                                                                                                  gTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGCTGTTTTACGAGATCAACATTACCt                                 >  2:376008/1‑150 (MQ=255)                                                                                                                       aTAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTg                            >  2:318414/1‑150 (MQ=255)                                                                                                                         aaGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGc                         >  2:190069/1‑151 (MQ=255)                                                                                                                         aaGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGc                         <  1:744475/151‑1 (MQ=255)                                                                                                                                       ttttGGTTACGCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGgtg           >  2:50206/1‑151 (MQ=255)                                                                                                                                                 gCTGGCGGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGCTGGTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCATACttt  >  1:309065/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                  ATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGC‑GGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAG‑ATTGGCCCGTACTGGTCAATCAGGCGATAAACCGCCGCG‑CTGTTAAGCTGTTTTACGAGATCAACATTACCTATCTGAGCTTGTCCGCCTGGTGTCATACTTT  >  NC_000913/701307‑701599

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