Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,184,588	G→C	13.5%	intergenic (‑204/‑208)	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,588	0	G	C	13.5%	139.6 / 13.7	52	intergenic (‑204/‑208)	mltF/purL	membrane‑bound lytic murein transglycosylase F/phosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base G (23/22);  new base C (5/2);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.60e-01

TCTGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGG  >  NC_002947/1184441‑1184724                                                                                                                                                    |                                                                                                                                         tCTGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCgcgg                                                                                                                                          <  7:127940/148‑1 (MQ=255)    ggTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCgcgggcgg                                                                                                                                      >  4:368897/1‑149 (MQ=255)     gTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCccgcgcgg                                                                                                                                      <  8:50011/148‑1 (MQ=255)     gTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCGGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGggcggc                                                                                                                                     >  1:122709/1‑149 (MQ=255)         cGGGTTGTCTGGTGGGGGGAGTGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                 <  4:6947/149‑1 (MQ=255)             tggtgggGATGGGGGCGGTGGGAGTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                             <  4:68572/147‑1 (MQ=37)                tgtgGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtctc                                                                                                                           >  6:26161/1‑148 (MQ=255)                tgtgGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                          >  2:273280/1‑149 (MQ=255)                  tgGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCGGCGTTAAGATATACAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCGCGCtcggccc                                                                                                                        >  8:315532/1‑144 (MQ=38)                            aGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAAt                                                                                                              <  1:211414/149‑1 (MQ=255)                                   aTTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTc                                                                                                         >  6:389870/1‑147 (MQ=255)                                      gcagcaCGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGCGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGcc                                                                                                    <  7:315532/149‑1 (MQ=255)                                        aTCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtctc                                                                                                                           <  2:77626/124‑1 (MQ=255)                                        aTCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtctc                                                                                                                           >  1:77626/1‑124 (MQ=255)                                         tCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGa                                                                                                 >  4:200502/1‑149 (MQ=255)                                           aGGGGTGCCAGGGGGGCGGCTTGAGATACTCAGCGCGTGTGAGTTCGAGCGCCGCCATGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCTAGCCGGag                                                                                                <  2:91265/148‑2 (MQ=25)                                              gCTGCCAGGGGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATACTCAGCGCCTATGAGCTCGCGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCt                                                                                            <  5:26161/149‑1 (MQ=255)                                                  ccAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTgcgg                                                                                        <  3:368897/149‑1 (MQ=255)                                                      gTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCt                                                                                    <  5:389870/149‑1 (MQ=255)                                                                          aGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGa                                                                >  8:230450/1‑149 (MQ=255)                                                                                         gAGCGCCGCCAGGTGTCCGCCGTAAGATATTCAGCGCCTATGAGAACGAGCGCCGAGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGATCGATCTCGACTCCAGCAGGACACCTCTGataa                                                    >  3:122033/1‑146 (MQ=38)                                                                                              ccgccAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAcccc                                            >  5:291217/1‑149 (MQ=255)                                                                                              ccgccAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAcccc                                            <  7:121947/149‑1 (MQ=255)                                                                                                 ccAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGaa                                         <  2:122709/149‑1 (MQ=255)                                                                                                  cAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAg                                        >  7:18745/1‑149 (MQ=255)                                                                                                   aGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGt                                       <  6:291217/149‑1 (MQ=255)                                                                                                     gtgtCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTcc                                     <  3:200502/149‑1 (MQ=255)                                                                                                           gCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTcc                               >  7:98817/1‑149 (MQ=255)                                                                                                             cTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGaa                                         >  2:328162/1‑137 (MQ=255)                                                                                                             cTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGaa                                         <  1:328162/137‑1 (MQ=255)                                                                                                             cTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTccgc                             <  7:230450/149‑1 (MQ=255)                                                                                                                aaGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                                                         >  6:9427/1‑54 (MQ=255)                                                                                                                aaGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                                                         <  5:9427/54‑1 (MQ=255)                                                                                                                  gATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAATGTGTCGAGCCGGACACCAGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAGCTGACCCCGACGACCGGTTCCGCCGccc                        >  4:285143/1‑149 (MQ=38)                                                                                                                  gATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGccc                        >  3:225578/1‑149 (MQ=255)                                                                                                                  cacatTCTGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGccc                        >  7:200171/4‑149 (MQ=25)                                                                                                                   atatTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaa                                                      >  8:290143/1‑118 (MQ=255)                                                                                                                   atatTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaa                                                      <  7:290143/118‑1 (MQ=255)                                                                                                                      ttCTGCGGCTATGTCATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACGTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAGCCCGAAGTCCGGTTCCGCCGCCCggga                    <  4:122033/149‑3 (MQ=38)                                                                                                                        cAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGTCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAg                   <  8:45469/148‑1 (MQ=255)                                                                                                                           cgcCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGt               >  8:120417/1‑149 (MQ=255)                                                                                                                              cTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGc            <  4:225578/149‑1 (MQ=255)                                                                                                                              cTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCAGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGc            <  1:273280/149‑1 (MQ=255)                                                                                                                                 tGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCt           >  7:47746/1‑147 (MQ=255)                                                                                                                                 tGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTa         <  8:47746/149‑1 (MQ=255)                                                                                                                                    gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTagag      >  4:23466/1‑149 (MQ=255)                                                                                                                                    gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTagag      >  7:49348/1‑149 (MQ=255)                                                                                                                                      tCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTagag      >  5:58612/1‑147 (MQ=255)                                                                                                                                       cGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg   >  4:18282/1‑149 (MQ=255)                                                                                                                                       cGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg   >  1:70334/1‑149 (MQ=255)                                                                                                                                        gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaaa                                                     <  1:308227/98‑1 (MQ=255)                                                                                                                                        gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaaa                                                     >  2:308227/1‑98 (MQ=255)                                                                                                                                        gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAgg  >  2:209784/1‑149 (MQ=255)                                                                                                                                         aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg   >  5:376597/1‑147 (MQ=255)                                                                                                                                         aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCt                                                                                    <  5:28932/66‑1 (MQ=255)                                                                                                                                         aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCt                                                                                    >  6:28932/1‑66 (MQ=255)                                                                                                                                                    |                                                                                                                                         TCTGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGG  >  NC_002947/1184441‑1184724

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