Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,184,605	C→T	18.6%	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	18.6%	100.2 / 15.3	49	intergenic (‑221/‑191)	mltF/purL	membrane‑bound lytic murein transglycosylase F/phosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base C (16/23);  new base T (5/4);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

GTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGAA  >  NC_002947/1184460‑1184752                                                                                                                                                  |                                                                                                                                                     gTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAgg                                                                                                                                                   >  8:192760/1‑149 (MQ=255)          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                                                                                                                           >  5:170931/1‑131 (MQ=255)                                           cTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGacac                                                                                                                           <  6:170931/131‑1 (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                                             >  8:89346/1‑149 (MQ=255)                                                                                                      aGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGc                                             >  4:270990/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                                           >  8:296532/1‑138 (MQ=37)                                                                                                                   tCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGt                                           <  7:296532/138‑1 (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)                                                                                                                                                  |                                                                                                                                                     GTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGAA  >  NC_002947/1184460‑1184752

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: ‑