Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,184,612	A→T	16.4%	intergenic (‑228/‑184)	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,612	0	A	T	16.4%	127.5 / 15.4	54	intergenic (‑228/‑184)	mltF/purL	membrane‑bound lytic murein transglycosylase F/phosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base A (20/25);  new base T (5/4);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

GGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGAATGCACGG  >  NC_002947/1184469‑1184759                                                                                                                                                |                                                                                                                                                     ggggggATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGACATTCTGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAAta                                                                                                                                                 <  8:88275/149‑2 (MQ=21)  gggggATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAAtgt                                                                                                                                                <  6:176786/149‑1 (MQ=255)           aaTCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGATATCGAGCGCCGCCAGGTGTACGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCAAAGGCGCAGAATGTGTCGAGcag                                                                                                                                       >  5:151964/2‑147 (MQ=255)                     gCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑gcg                                                                                                                            <  4:55835/149‑1 (MQ=255)                      ccAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑gcgg                                                                                                                           >  7:188095/1‑149 (MQ=255)                      ccAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑gcgg                                                                                                                           >  4:202323/1‑149 (MQ=255)                          gcggCCGGCTTGAGATATTCCGCGCCGGTGAGCTCGAGCGCCGCCAGGTGGCCGCCTTAAGATATGCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCt                                                                                                                       <  5:366711/147‑1 (MQ=38)                            ggCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCg                                                                                                                     <  3:74320/149‑1 (MQ=255)                             gCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGACGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGa                                                                                                                    >  7:294052/1‑149 (MQ=255)                                  cTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGacac                                                                                                                                  <  6:170931/131‑1 (MQ=255)                                  cTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGacac                                                                                                                                  >  5:170931/1‑131 (MQ=255)                                     gagaTATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGtct                                                                                                            >  3:353192/1‑149 (MQ=255)                                     gagaTATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCAATCTCGt                                                                                                              <  6:340372/147‑1 (MQ=255)                                             cAGCGCCTGTTAGATCGAGCGCCGCCAGGTGTCCTCCTTAAGATTTTCAGCGCCTTTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCg                                                                                                     <  4:62605/148‑1 (MQ=255)                                             cAGCGACTGTGAGTGCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATGGCGCGCGGCGCGGGCGTCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCg                                                                                                     <  5:100926/148‑1 (MQ=38)                                              aGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGa                                                                                                   >  5:277138/1‑149 (MQ=255)                                                 gcCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGacac                                                                                                <  4:136538/149‑1 (MQ=255)                                                   cTGTGAGATCGAGCGACTCCAGGTGTCCGCATTAAGATATTCAGCTCCTATGAGATCTAGCGCCGCGCGGGCTGCGCTCGATCTCACAGGCGCATAATGTGTCGAGCCTGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACct                                                                                              <  6:151964/149‑1 (MQ=38)                                                             gAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCATGCCGGACACCTCTGAAAACTg                                                                                    >  7:122805/1‑149 (MQ=255)                                                                     ccAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCg                                                                              >  8:32779/1‑147 (MQ=255)                                                                            tCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGt                                                                     <  6:277138/149‑1 (MQ=255)                                                                                  ttAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCCCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTccgcc                                                               >  3:355176/1‑149 (MQ=255)                                                                                      gATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGccc                                                           <  6:127751/149‑1 (MQ=255)                                                                                           tCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAg                                                      >  4:273664/1‑149 (MQ=255)                                                                                             aGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGc                                                    >  4:270990/1‑149 (MQ=255)                                                                                             aGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGc                                                    >  8:89346/1‑149 (MQ=255)                                                                                               cgcCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGt                                                  <  7:32779/149‑1 (MQ=255)                                                                                                 ccTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTgg                                                <  3:273664/149‑1 (MQ=255)                                                                                                  cTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGc                                               <  3:202323/149‑1 (MQ=255)                                                                                                    gtGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCgg                                                         <  4:125146/136‑1 (MQ=25)                                                                                                    gtGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCgg                                                         >  3:125146/2‑137 (MQ=25)                                                                                                    gtGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATCTTAAGGCGGACACCTGGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCt                                              >  2:229915/2‑149 (MQ=17)                                                                                                      gagaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTa                                            <  8:122805/148‑1 (MQ=255)                                                                                                       agaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCCCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTaga                                          <  4:355176/149‑1 (MQ=255)                                                                                                       agaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATCTTAAGGCGGACACCTGGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTa                                            <  1:229915/148‑1 (MQ=18)                                                                                                        gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTccg                                                                 >  1:225553/1‑125 (MQ=255)                                                                                                        gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTccg                                                                 <  2:225553/125‑1 (MQ=255)                                                                                                          tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCa                                       <  7:89346/149‑1 (MQ=255)                                                                                                          tCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGt                                                  <  7:296532/138‑1 (MQ=37)                                                                                                          tCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGt                                                  >  8:296532/1‑138 (MQ=37)                                                                                                           cGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg                                      >  4:273466/1‑149 (MQ=255)                                                                                                            gAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATCTTAAGGCGGACACCTGGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg                                      >  3:133054/1‑149 (MQ=18)                                                                                                             aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGt                                    <  7:192760/149‑1 (MQ=255)                                                                                                               cgccgcGCTGTCGGAGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGc                                  <  8:323308/149‑1 (MQ=255)                                                                                                                   gcgcgGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCAGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGa                               >  1:275566/1‑148 (MQ=255)                                                                                                                       gggcggCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGcc                          <  5:151131/149‑1 (MQ=255)                                                                                                                             cgcTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGaa                    <  5:47659/149‑1 (MQ=255)                                                                                                                                  gATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAgg               >  4:46172/1‑149 (MQ=255)                                                                                                                                       cacaGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAg           >  1:115460/1‑148 (MQ=255)                                                                                                                                        acaGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAg           <  4:353192/147‑1 (MQ=255)                                                                                                                                        acaGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGaa         >  2:228358/1‑149 (MQ=255)                                                                                                                                        acaGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGaa         >  2:24466/1‑149 (MQ=255)                                                                                                                                            gcgcAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGAATGCa     >  5:6553/1‑149 (MQ=255)                                                                                                                                               cAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGAATGCACgg  <  8:188095/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     GGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGAATGCACGG  >  NC_002947/1184469‑1184759

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