Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,184,588	G→C	44.0%	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	44.0%	11.2 / 53.2	41	intergenic (‑204/‑208)	mltF/purL	membrane‑bound lytic murein transglycosylase F/phosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base G (17/6);  new base C (10/8);  total (27/14)
Fisher's exact test for biased strand distribution p-value = 3.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

GGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCG  >  NC_002947/1184444‑1184729                                                                                                                                                 |                                                                                                                                                tgTAGCGGGTTGGGTGGTTGGGGGAGGGGTGATTGCATCAGGGCTGCGAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCccgcgcg                                                                                                                                              <  8:353477/147‑1 (MQ=37)  gtcgCGGGTTGTGTGGTGGGGGGAGGGGTGATTGCATCAGGGCTGCCAGGTGGGCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCgcgggc                                                                                                                                               <  4:73028/143‑1 (MQ=255)     gCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGg‑cgc                                                                                                                                         <  1:312904/148‑1 (MQ=255)       gggTTGGGTGGGTGGGGGAGGGGGGATGGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGG‑CGCTCg                                                                                                                                      <  3:318174/149‑1 (MQ=37)                      gggAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGACATATAGCGCCGCGCGGGCGG‑CGCTCg                                                                                                                                      >  3:397725/1‑134 (MQ=255)                          ggggggATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGG‑CGCt                                                                                                                                        >  8:8426/1‑128 (MQ=255)                          ggggggATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGG‑CGCt                                                                                                                                        <  7:8426/128‑1 (MQ=255)                            ggggACTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCTGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCCGCTCGATCTCATAGGCGCTGAAtat                                                                                                                  <  1:5799/149‑3 (MQ=21)                             gggATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGggcgg‑c                                                                                                                                           >  6:37669/1‑122 (MQ=255)                             gggATTGCATCAGGGCGGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGggcgg‑c                                                                                                                                           <  5:37669/122‑1 (MQ=255)                                                   gTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGG‑CGCTCGAt                                                                                                                                    >  8:47625/1‑107 (MQ=255)                                                   gTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGG‑CGCTCGAt                                                                                                                                    <  7:47625/107‑1 (MQ=255)                                                           cTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAg                                                                                                                             >  2:131993/1‑106 (MQ=255)                                                           cTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAg                                                                                                                             <  1:131993/106‑1 (MQ=255)                                                                                         cgccgcCAGGTGTCCTCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCTCTCGATCTCGTCTCAAGCCGGACACCTCGGAAAACTGccc                                                   >  6:115187/1‑149 (MQ=255)                                                                                         cgccgcCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAcc                                                   >  6:120943/1‑149 (MQ=255)                                                                                                  gtgtCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTcc                                          >  1:24815/1‑149 (MQ=255)                                                                                                         ccTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTccg                                   <  2:24815/149‑1 (MQ=255)                                                                                                                     cAGCGCCTATGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAgg                       >  2:256655/1‑149 (MQ=255)                                                                                                                      aGCGCCTATGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCg                                                >  2:289824/1‑123 (MQ=255)                                                                                                                      aGCGCCTATGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGCCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCg                                                <  1:289824/123‑1 (MQ=255)                                                                                                                        cgcCTGTGAGATCGAGCGCCGCCCGCGCGG‑CGCTCGATCTCATAGGCGCTGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGt                    >  6:374454/1‑149 (MQ=25)                                                                                                                          ccTATGAGATCGAGCGCCGTGCGGGCGG‑CGCTCGGTCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTgg                  >  1:55964/1‑149 (MQ=255)                                                                                                                           cTATGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGc                 >  4:133182/1‑149 (MQ=255)                                                                                                                            tATGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTAGATCTCGACTCAACCCGGACACCTCAGAAAACTGACGCCGACGTCCCGATCCGCCGCCCGGCAGGCGTGGcg                >  4:291788/1‑148 (MQ=38)                                                                                                                             aTGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCt                >  6:168783/1‑148 (MQ=255)                                                                                                                              tGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGt                                            >  4:148174/1‑119 (MQ=255)                                                                                                                              tGAGATCGAGCGCCGCCCGCGCGG‑CGCTCGATCTCATAGGCGCTGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCtt               <  5:374454/148‑1 (MQ=255)                                                                                                                                agaTCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTa              >  6:93084/1‑147 (MQ=255)                                                                                                                                 gaTCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAcc                                                   <  5:79020/109‑1 (MQ=255)                                                                                                                                 gaTCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAcc                                                   >  6:79020/1‑109 (MQ=255)                                                                                                                                 gaTCGAGCGCCGCCCGCGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGc                 >  7:16858/1‑143 (MQ=255)                                                                                                                                 gaTCGAGCGCCGCCCGCGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGc                 <  8:16858/143‑1 (MQ=255)                                                                                                                                  aTCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTagag           >  1:60432/1‑148 (MQ=255)                                                                                                                                   tCGAGCGCCGCCCGCGCGG‑CGCTCGATCTCATAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTg                   <  7:60474/139‑1 (MQ=255)                                                                                                                                   tCGAGCGCCGCCCGCGCGG‑CGCTCGATCTCATAGACGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTg                   >  8:60474/1‑139 (MQ=255)                                                                                                                                    cGAGCGCCGCCCGCGCGG‑CGCTCGATCTCATAGGCGCTGAATATCTTAAGGCGGACACCTGGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTgg                  <  3:339863/140‑1 (MQ=17)                                                                                                                                    cGAGCGCCGCCCGCGCGG‑CGCTCGATCTCATAGGCGCTGAATATCTTAAGGCGGACACCTGGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTgg                  >  4:339863/1‑140 (MQ=17)                                                                                                                                     gAGCGCCGCCCGCGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAgg       >  2:103632/1‑149 (MQ=255)                                                                                                                                      aGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGt      >  7:86730/1‑149 (MQ=255)                                                                                                                                      aGCGCCGCCCGCGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGCCACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGt      >  6:59460/1‑149 (MQ=255)                                                                                                                                       gcgcCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTgg                  >  8:13138/1‑136 (MQ=255)                                                                                                                                       gcgcCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTgg                  <  7:13138/136‑1 (MQ=255)                                                                                                                                       gcgcCGCCCGCGCGT‑CGCTCGATCTCACAGTCGCAGAATGAGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGCCACCGCTGAACACAGACCCCGAAGACCGGTTCCGCCGCCCGGCAGGCGTTGCTGAGAGAAGGTg     >  8:332121/1‑149 (MQ=255)                                                                                                                                           cgcccgCGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCg  >  4:147364/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                GGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCG  >  NC_002947/1184444‑1184729

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