Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	701,378	+T	coding (217/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,378	1	.	T	83.8%	84.3 / 12.7	38	coding (217/1221 nt)	nagC	N‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base . (4/2);  new base T (14/17);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCGGGTAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑C‑CCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGG  >  NC_000913/701243‑701501                                                                                                                                              |                                                                                                                            gCAGCGGGTAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑GTG‑C‑C‑CCCCGt                                                                                                                       >  2:665661/1‑143 (MQ=255)         tAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCct                                                                                                               >  1:281746/1‑143 (MQ=255)              ttgttctGCCACCCCTTTGCTGCGCTGATCACACAGCGTGATGGTGGCGTCATGACGCCCGCTCCGTACGCCGCTTGCGTCGTCA‑TTGCAGGTTTCGGTGACGATGGATATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGa                                                                                                          <  2:714238/140‑1 (MQ=255)                  tctGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCaa                                                                                                      <  1:372895/143‑1 (MQ=255)                      ccAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAActtc                                                                                                  >  1:258898/1‑143 (MQ=255)                            ttttGCTGCTGAGATCAAACAGAGCGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGAt                                                                                            <  1:3929/142‑1 (MQ=255)                            cTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGAt                                                                                            <  2:304622/143‑1 (MQ=255)                             tttGCTGCTGAGATCAAACAGAGTGAAGCTAGCGTCATGTCGACCAAGCCGTACGCCGATTGCGAGGAAAGTTGC‑GGATCGGAGGACGAGGGAGATAGCGCGAC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGAGCAACGTCTGTGATc                                                                                           >  2:398454/1‑143 (MQ=255)                                  tgctgAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGccc                                                                                      >  1:213733/1‑143 (MQ=255)                                        gaTCAAACTGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTc                                                                                <  1:708555/143‑1 (MQ=255)                                         aTCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGACA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCTCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTCCg                                                                               >  1:322690/1‑143 (MQ=255)                                         aTCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCg                                                                               >  1:107721/1‑143 (MQ=255)                                          tCAAACAGAGTGATGGTGGCGTCATGACGACCAAGACGTAAGCCGATTGGGTGGAAA‑TTGCGGGTTTCGGTGAAGATGGAGATAGCGCGGC‑G‑G‑CAT‑CCCGGTGGAGGCCTGCTGATCAACTTCTTTGCTCAGCGCGCGCTCGa                                                                              >  2:557497/1‑143 (MQ=255)                                              aCAGAGGGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTGTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAg                                                                          <  1:35673/143‑1 (MQ=255)                                                   gtggtggtgGCGTCGTGAGGGCCAAGCTGTGGGGCGATTGGGTGGAGA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGAc                                                                     <  2:458038/143‑1 (MQ=255)                                                    tggtggtggGGTCAGGCCGGCGCAGGCGTACGCGGAGTGCGTGGAAA‑GTGCGGGTTTCGGTGACGAGGGAGATAGCGGGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGatg                                                                    <  2:650250/140‑3 (MQ=255)                                                    tGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑GCC‑C‑CCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATACGCTGAcg                                                                    >  1:583283/1‑143 (MQ=255)                                                      aTGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGGTCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGAcgcg                                                                  >  2:407501/1‑143 (MQ=255)                                                               tCATGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTgg                                                         >  2:94261/1‑143 (MQ=255)                                                                cATGACGGCCAAGCTGTACGCCGATGGTGTTGAGA‑GTGCGGGTTGCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGCTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGt                                                        <  2:41970/143‑1 (MQ=255)                                                                 aTGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGAC‑G‑G‑C‑CTCCCGGTGGAGGCCGGCTGATCAACTACTTTGATAAGCCCCGGTACGATAAGACTACGCGAAAATCTGGtg                                                       >  2:611474/1‑142 (MQ=255)                                                                  tGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGG‑G‑C‑T‑CCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTGCGGTAGGCTGACGCGTCATTTTGGTTa                                                      >  2:39472/1‑143 (MQ=255)                                                                  tGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTa                                                      >  1:22494/1‑143 (MQ=255)                                                                   gacgacCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTAc                                                     <  1:477632/143‑1 (MQ=255)                                                                        cccaGCCGGCCGCCGACTGCGTGGCAA‑GTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTgg                                                <  2:419853/140‑1 (MQ=255)                                                                          aaGCCGTACGCCGATTGCGGGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑C‑CCCGGTGGAGGCCTGCTGATCAACTTCTTTGAGCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCgg                                             <  1:214642/143‑1 (MQ=255)                                                                           aGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCgg                                             <  1:667587/143‑1 (MQ=255)                                                                                 aCGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCaa                                       >  1:331442/1‑143 (MQ=255)                                                                                  cGCCGATGGCGTGGAAG‑TTGCGGGTTTCGGGGACGGTGGAGATAGCGCGGC‑G‑G‑CGT‑CCCGGTGGAGGCCCGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAg                                      <  1:399891/143‑1 (MQ=255)                                                                                      gATTGCGTGGGAA‑TTGCGGGTTGCGGTGGCGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATGAGCGCGCGTTGGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTgg                                  <  2:538698/143‑1 (MQ=255)                                                                                            gTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTCCGCTGGCGGGGGCAAGATGGCTTTGc                            >  2:110514/1‑143 (MQ=255)                                                                                               gaca‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCg                         <  1:236452/140‑1 (MQ=255)                                                                                                 acgtggCGGGTTCCGGTGACGTTGGAGATAGCGCGGC‑C‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCGTTGATCAGCCCGCGTTCGACAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCACCCTGGCTGGGCTCgg                        <  2:525112/138‑1 (MQ=255)                                                                                                      gCGGGTTTCGGGGACGATGGAGATAGCGCGGC‑G‑G‑C‑CGCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATc                   <  2:8962/143‑1 (MQ=255)                                                                                                                 tGACGGTGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCGGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCga        <  2:258898/143‑1 (MQ=255)                                                                                                                     gATGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTGCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTGCGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGGtt    >  2:633926/1‑143 (MQ=255)                                                                                                                      aTGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTg   <  2:107721/143‑1 (MQ=255)                                                                                                                       tGGAGATAGCGCGGC‑G‑G‑C‑CTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATggg  >  2:495192/1‑143 (MQ=255)                                                                                                                                              |                                                                                                                            GCAGCGGGTAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAA‑TTGCGGGTTTCGGTGACGATGGAGATAGCGCGGC‑G‑G‑C‑C‑CCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGG  >  NC_000913/701243‑701501

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