Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,184,574	A→G	10.9%	intergenic (‑184/‑228)	mltF ← / → purL	membrane‑bound lytic murein transglycosylase F/phosphoribosylformylglycinamidine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,184,574	0	A	G	10.9%	118.6 / 5.3	46	intergenic (‑184/‑228)	mltF/purL	membrane‑bound lytic murein transglycosylase F/phosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base G (3/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-01

ATCTCGACAACGGTTAGCTTGTCTGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCT  >  NC_002947/1184426‑1184721                                                                                                                                                     |                                                                                                                                                    atCTCGACAACGGTTAGCTTGTCTGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGAGGTGGGGGGGCCGGCTAGAGATGGTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTGAAGATGTTCAGCGCCTAt                                                                                                                                                    >  1:133215/1‑150 (MQ=37)     cGACAACGGTTAGCTTGTCTGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATgaga                                                                                                                                                >  2:441089/1‑150 (MQ=255)                        tGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCccgcgcgg                                                                                                                             <  2:1829643/150‑1 (MQ=255)                             gCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                        <  1:441089/150‑1 (MQ=255)                             gCGGGTTGTGTGGTTGGGGGAGGGGGGAGTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                        <  2:743946/150‑1 (MQ=255)                                ggTTGTGTGGTTGGGGGAGTGTGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                     <  2:814725/150‑1 (MQ=255)                                ggTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                                     <  2:993700/150‑1 (MQ=255)                                       tgGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCG‑CAGGTGTCCGGCTCGACACATTCTGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAgg                                                                                                             >  1:301678/1‑150 (MQ=14)                                            gggggAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCa                                                                                                         >  2:1696032/1‑150 (MQ=255)                                              gggAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGACCTCACAGGCGCAGa                                                                                                       >  2:1063479/1‑150 (MQ=255)                                               ggAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGGTCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                        <  2:922805/132‑1 (MQ=37)                                               ggAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGGTCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                        >  1:922805/1‑132 (MQ=37)                                                  ggggggATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAAtgt                                                                                                   <  2:1107911/150‑1 (MQ=255)                                                   gggggATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAAtgtg                                                                                                  <  1:1696032/150‑1 (MQ=255)                                                     gggATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGNGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTc                                                                                                >  1:1633531/1‑150 (MQ=255)                                                           gcagcaGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCg                                                                                          <  1:746921/150‑1 (MQ=255)                                                             aTCAGGTCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCCCCTTAAGCTATTCAGCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGa                                                                                        <  2:1503947/150‑1 (MQ=37)                                                             aTCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                      <  2:1391152/120‑1 (MQ=25)                                                             aTCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                      >  1:1391152/1‑120 (MQ=25)                                                                 gggCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCgcgggc                                                                                                                               <  2:972433/107‑1 (MQ=255)                                                                 gggCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCgcgggc                                                                                                                               >  1:972433/1‑107 (MQ=255)                                                                   gCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTg                                                                                  >  1:1028076/1‑150 (MQ=255)                                                                          ggTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCt                                                                           <  2:1770646/150‑1 (MQ=255)                                                                           gTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCACCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTc                                                                          <  2:533101/150‑1 (MQ=255)                                                                                         atatTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGGTCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAg                                                            >  2:1282268/1‑150 (MQ=255)                                                                                          tatTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGc                                                           <  2:1633531/150‑1 (MQ=255)                                                                                           atTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGcc                                                          <  2:1028076/150‑1 (MQ=255)                                                                                                       tgAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGgc                                                                                                            >  1:373715/1‑88 (MQ=255)                                                                                                       tgAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGgc                                                                                                            <  2:373715/88‑1 (MQ=255)                                                                                                        gagaTCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGGTCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaa                                             <  1:1282268/150‑1 (MQ=255)                                                                                                                gcgcCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAcc                                     >  1:1817061/1‑150 (MQ=255)                                                                                                                  gccgccAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAcccc                                   >  1:1473095/1‑150 (MQ=255)                                                                                                                    cgccAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGa                                 >  2:1830737/1‑150 (MQ=255)                                                                                                                     gccAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGGTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGaa                                >  1:1608244/1‑150 (MQ=255)                                                                                                                     gccAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGaa                                >  1:640973/1‑150 (MQ=255)                                                                                                                      ccAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAg                               >  1:1142884/1‑150 (MQ=255)                                                                                                                               cGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGaca                                                                                      <  2:1202413/86‑1 (MQ=255)                                                                                                                               cGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGaca                                                                                      >  1:1202413/1‑86 (MQ=255)                                                                                                                                gCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGgcgc                                                                            >  2:551770/1‑95 (MQ=255)                                                                                                                                gCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGgcgc                                                                            <  1:551770/95‑1 (MQ=255)                                                                                                                                        atatTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATgtgt                                                                                                 <  2:733594/66‑1 (MQ=255)                                                                                                                                        atatTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATgtgt                                                                                                 >  1:733594/1‑66 (MQ=255)                                                                                                                                                 gcCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTgg    >  2:314888/1‑150 (MQ=255)                                                                                                                                                 gcCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTgg    >  1:1204677/1‑150 (MQ=255)                                                                                                                                                   cTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCt  <  2:1204677/150‑1 (MQ=255)                                                                                                                                                   cTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCt  >  2:1070500/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                    ATCTCGACAACGGTTAGCTTGTCTGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCT  >  NC_002947/1184426‑1184721

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