Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	93,440	A→G	E92G (GAG→GGG)	murE →	UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	93,440	0	A	G	89.5%	52.7 / ‑3.7	19	E92G (GAG→GGG)	murE	UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (12/5);  total (12/7)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGGGCGTCGATATATCCCGCAGGCGATAGCGCAAGGTGTGGCTGCCATTATTGCAGAGGCGAAAGATGAGGCGACCGATGGTGAAATCCGTGAAATGCACGGCGTACCGGTCATCTATCTCAGCCAGCTCAACGAGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAAT  >  NC_000913/93304‑93575                                                                                                                                         |                                                                                                                                        gACGGGCGTCGATATATCCCGCAGGCGATAGCGCAAGGTGTGGCTGCCATTATTGCAGAGGCGAAAGATGAGGCGACCGATGGTGAAATCCGTGAAATGCACGGCGTACCGGTCATCTATCTCAGCCAGCTCAACGGGCg                                                                                                                                      >  2:483238/1‑140 (MQ=255) gACGGGCGTCGATATATCCCGCAGGCGATAGCGCAAGGTGTGGCTGCCATTATTGCAGAGGCGAAAGATGAGGCGACCGATGGTGAAATCCGTGAAATGCACGGCGTACCGGTCATCTATCTCAGCCAGCTCAACGGGCg                                                                                                                                      <  1:483238/140‑1 (MQ=255)      gCGTCGATATATCCCGCAGGCGATAGCGCAAGGTGTGGCTGCCATTATTGCAGAGGCGAAAGATGAGGCGACCGATGGTGAAATCCGTGAAATGCACGGCGTACCGGTCATCTATCTCAGCCAGCTCAACGGGCGTTTAt                                                                                                                                 >  2:457100/1‑140 (MQ=255)                    cGCAGGCGATAGCGCAAGGTGTGGCTGCCATTATTGCAGAGGCGAAAGATGAGGCGACCGATGGTGAAATCCGTGAAATGCACGGCGTACCGGTCATCTATCTCAGCCAGCTCAACGGGCGTTTATCTGCACTGGCGGGc                                                                                                                   >  2:644634/1‑140 (MQ=255)                                   aaGGTGTGGCTGCCATTATTGCAGAGGCGAAAGATGAGGCGACCGATGGTGAAATCCGTGAAATGCACGGCGTACCGGTCATCTATCTCAGCCAGCTCAACGGGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGaa                                                                                                    <  2:374471/140‑1 (MQ=255)                                                 attattGCAGAGGCGAAAGATGAGGCGACCGATGGTGAAATCCGTGAAATGCACGGCGTACCGGTCATCTATCTCAGCCAGCTCAACGGGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAAttt                                                                                      >  2:51328/1‑140 (MQ=255)                                                             gCGAAAGATGAGGCGACCGATGGTGAAATCCGTGAAATGCACGGCGTACCGGTCATCTATCTCAGCCAGCTCAACGGGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTggg                                                                          >  2:99/1‑140 (MQ=255)                                                                       aGGCGACCGATGGTGAAATCCGTGAAATGCACGGCGTACCGGTCATCTATCTCAGCCAGCTCAACGGGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGc                                                                >  2:444687/1‑140 (MQ=255)                                                                             ccGATGGTGAAATCCGTGAAATGCACGGCGTACCGGTCATCTATCTCAGCCAGCTCAACGGGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAAc                                                          >  2:120353/1‑140 (MQ=255)                                                                                        aTCCGTGAAATGCACGGCGTACCGGTCATCTATCTCAGCCAGCTCAACGGGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCAc                                               <  1:99/140‑1 (MQ=255)                                                                                                         cGTACCGGTCATCTATCTCAGCCAGCTCAACGGGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTgg                              >  2:142538/1‑140 (MQ=255)                                                                                                          gTACCGGTCATCTATCTCAGCCAGCTCAACGGGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGgc                             <  1:457100/140‑1 (MQ=255)                                                                                                            aCCGGTCATCTATCTAAGCCAGTTCCACGAG‑‑‑‑TATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTg                       <  1:169248/140‑1 (MQ=255)                                                                                                                gTCATCTATCTCAGCCAGCTCAACGGGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTg                       <  1:444687/140‑1 (MQ=255)                                                                                                                   atctatctCAGCCAGCTCAACGGGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGAAAAACCACGACTGCCCCTATGTTGGCTCAGTGGAg                    >  2:310860/1‑140 (MQ=255)                                                                                                                        tctcAGCCAGCTCAACGGGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGGAACGGGCACCAACGGCAAAAACACGACTACCCCGCTGTTGGCGCAGTGGAGCCAAc               >  2:558873/1‑140 (MQ=255)                                                                                                                                agcTCAACGGGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCCTGGGCGTAACGGGCACCCACCGCAAAACCACGACCACCCCGCTGTTTGCGCAGGGGAGCCCACTCCTTggg       >  2:579928/1‑139 (MQ=255)                                                                                                                                  cTCAACGGGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGa     >  2:634727/1‑140 (MQ=255)                                                                                                                                     aaCGAGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAAt  <  1:329530/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        GACGGGCGTCGATATATCCCGCAGGCGATAGCGCAAGGTGTGGCTGCCATTATTGCAGAGGCGAAAGATGAGGCGACCGATGGTGAAATCCGTGAAATGCACGGCGTACCGGTCATCTATCTCAGCCAGCTCAACGAGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAAT  >  NC_000913/93304‑93575

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑