Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,184,605	C→T	14.3%	intergenic (‑221/‑191)	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,605	0	C	T	14.3%	119.4 / 11.0	49	intergenic (‑221/‑191)	mltF/purL	membrane‑bound lytic murein transglycosylase F/phosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base C (23/19);  new base T (4/3);  total (27/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 = 8.44e-01

AGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGAAT  >  NC_002947/1184468‑1184753                                                                                                                                          |                                                                                                                                                      aGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAAt                                                                                                                                            <  1:211414/149‑1 (MQ=255)        aTTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTc                                                                                                                                       >  6:389870/1‑147 (MQ=255)           gcagcaCGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGCGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGcc                                                                                                                                  <  7:315532/149‑1 (MQ=255)              tCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGa                                                                                                                               >  4:200502/1‑149 (MQ=255)                aGGGGTGCCAGGGGGGCGGCTTGAGATACTCAGCGCGTGTGAGTTCGAGCGCCGCCATGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCTAGCCGGag                                                                                                                              <  2:91265/148‑2 (MQ=25)                   gCTGCCAGGGGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATACTCAGCGCCTATGAGCTCGCGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCt                                                                                                                          <  5:26161/149‑1 (MQ=255)                       ccAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑gcgg                                                                                                                     <  3:368897/149‑1 (MQ=255)                           gTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCt                                                                                                                 <  5:389870/149‑1 (MQ=255)                                               aGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGa                                                                                             >  8:230450/1‑149 (MQ=255)                                                              gAGCGCCGCCAGGTGTCCGCCGTAAGATATTCAGCGCCTATGAGAACGAGCGCCGAGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGATCGATCTCGACTCCAGCAGGACACCTCTGataa                                                                                 >  3:122033/1‑146 (MQ=38)                                                                   ccgccAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAcccc                                                                         <  7:121947/149‑1 (MQ=255)                                                                   ccgccAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAcccc                                                                         >  5:291217/1‑149 (MQ=255)                                                                      ccAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGaa                                                                      <  2:122709/149‑1 (MQ=255)                                                                       cAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAg                                                                     >  7:18745/1‑149 (MQ=255)                                                                        aGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGt                                                                    <  6:291217/149‑1 (MQ=255)                                                                          gtgtCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTcc                                                                  <  3:200502/149‑1 (MQ=255)                                                                                gCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTcc                                                            >  7:98817/1‑149 (MQ=255)                                                                                  cTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGaa                                                                      >  2:328162/1‑137 (MQ=255)                                                                                  cTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGaa                                                                      <  1:328162/137‑1 (MQ=255)                                                                                  cTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTccgc                                                          <  7:230450/149‑1 (MQ=255)                                                                                     aaGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                                                                                       >  6:9427/1‑54 (MQ=255)                                                                                     aaGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                                                                                       <  5:9427/54‑1 (MQ=255)                                                                                       gATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAATGTGTCGAGCCGGACACCA‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAGCTGACCCCGACGACCGGTTCCGCCGccc                                                     >  4:285143/1‑149 (MQ=38)                                                                                       gATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGccc                                                     >  3:225578/1‑149 (MQ=255)                                                                                       cacatTCTGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGccc                                                     >  7:200171/4‑149 (MQ=25)                                                                                        atatTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaa                                                                                   <  7:290143/118‑1 (MQ=255)                                                                                        atatTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaa                                                                                   >  8:290143/1‑118 (MQ=255)                                                                                           ttCTGCGGCTATGTCATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACGT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAGCCCGAAGTCCGGTTCCGCCGCCCggga                                                 <  4:122033/149‑3 (MQ=38)                                                                                             cAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGTCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAg                                                <  8:45469/148‑1 (MQ=255)                                                                                                cgcCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGt                                            >  8:120417/1‑149 (MQ=255)                                                                                                   cTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGc                                         <  4:225578/149‑1 (MQ=255)                                                                                                   cTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCAGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGc                                         <  1:273280/149‑1 (MQ=255)                                                                                                      tGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCt                                        >  7:47746/1‑147 (MQ=255)                                                                                                      tGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTa                                      <  8:47746/149‑1 (MQ=255)                                                                                                         gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTagag                                   >  7:49348/1‑149 (MQ=255)                                                                                                         gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTagag                                   >  4:23466/1‑149 (MQ=255)                                                                                                           tCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTagag                                   >  5:58612/1‑147 (MQ=255)                                                                                                            cGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg                                >  4:18282/1‑149 (MQ=255)                                                                                                            cGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg                                >  1:70334/1‑149 (MQ=255)                                                                                                             gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaaa                                                                                  >  2:308227/1‑98 (MQ=255)                                                                                                             gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaaa                                                                                  <  1:308227/98‑1 (MQ=255)                                                                                                             gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAgg                               >  2:209784/1‑149 (MQ=255)                                                                                                              aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg                                >  5:376597/1‑147 (MQ=255)                                                                                                              aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCt                                                                                                                 <  5:28932/66‑1 (MQ=255)                                                                                                              aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCt                                                                                                                 >  6:28932/1‑66 (MQ=255)                                                                                                                      ccgcgcggcgCTCGATCTCATAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAg                      <  8:200171/145‑1 (MQ=255)                                                                                                                                   gATCTCATAGGCGCTGAATATCTTAAGGCGGACACCTGGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGACTTAGAGCAGGTGCCGAAAGCCATCGAATTa           <  4:355644/148‑1 (MQ=25)                                                                                                                                   gATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCACGGCAGGCGTGGCTTAGAGCAGGTGCCGAACGCCATCGAATTAgg         >  7:331429/1‑149 (MQ=255)                                                                                                                                   gATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGAACTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGAGCCGAAAGCCATCGAATTAg          >  1:91212/1‑148 (MQ=255)                                                                                                                                     tctcACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGc        >  1:262485/1‑148 (MQ=255)                                                                                                                                     tctcACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGc        >  3:331574/1‑148 (MQ=255)                                                                                                                                        cacaGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAg     >  1:292936/1‑148 (MQ=255)                                                                                                                                         acaGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACct                                                                                        <  4:379907/64‑1 (MQ=255)                                                                                                                                         acaGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTa      <  6:376597/146‑1 (MQ=255)                                                                                                                                         acaGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGaa   >  7:54506/1‑149 (MQ=255)                                                                                                                                          caGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGAAt  <  8:49348/149‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                      AGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGAAT  >  NC_002947/1184468‑1184753

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