Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,184,605	C→T	17.9%	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	17.9%	87.5 / 9.7	39	intergenic (‑221/‑191)	mltF/purL	membrane‑bound lytic murein transglycosylase F/phosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base C (23/9);  new base T (3/4);  total (26/13)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGAA  >  NC_002947/1184471‑1184752                                                                                                                                        |                                                                                                                                                     ggggACTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCTGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCCGCTCGATCTCATAGGCGCTGAAtat                                                                                                                                         <  1:5799/149‑3 (MQ=21)                                cTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAg                                                                                                                                                    <  1:131993/106‑1 (MQ=255)                                cTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAg                                                                                                                                                    >  2:131993/1‑106 (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                                        <  8:16858/143‑1 (MQ=255)                                                                                                      gaTCGAGCGCCGCCCGCGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGc                                        >  7:16858/1‑143 (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                                         >  4:339863/1‑140 (MQ=17)                                                                                                         cGAGCGCCGCCCGCGCGG‑CGCTCGATCTCATAGGCGCTGAATATCTTAAGGCGGACACCTGGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTgg                                         <  3:339863/140‑1 (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)                                                                                                                         gg‑cgCTCGATCTCACTGGCGCAGACTGTGTCGAGCCGGAGACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCTGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCg               <  7:332121/149‑1 (MQ=255)                                                                                                                               tCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGc                           >  8:77542/1‑132 (MQ=255)                                                                                                                               tCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGc                           <  7:77542/132‑1 (MQ=255)                                                                                                                                    ctcACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg                               <  3:173823/123‑1 (MQ=255)                                                                                                                                    ctcACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg                               >  4:173823/1‑123 (MQ=255)                                                                                                                                       acaGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg                               >  2:91975/1‑120 (MQ=255)                                                                                                                                       acaGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg                               <  1:91975/120‑1 (MQ=255)                                                                                                                                       acaGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTa     <  5:168783/146‑1 (MQ=255)                                                                                                                                       acaGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGaa  >  2:154347/1‑149 (MQ=255)                                                                                                                                       acaGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGa   <  1:256655/148‑1 (MQ=255)                                                                                                                                        |                                                                                                                                                     GGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGG‑CGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCT‑GCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGAA  >  NC_002947/1184471‑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: ‑