Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	855,864	A→G	8.0%	T239A (ACC→GCC) ‡	murI →	glutamate racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	855,864	0	A	G	8.0%	134.4 / 2.8	50	T239A (ACC→GCC) ‡	murI	glutamate racemase
Reads supporting (aligned to +/- strand):  ref base A (22/24);  new base G (2/2);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

CCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGA  >  NC_002947/855716‑856012                                                                                                                                                     |                                                                                                                                                     cccTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGa                                                                                                                                                      >  8:300255/1‑149 (MQ=255)   cTTTCCTTCGCCCAGTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGAc                                                                                                                                                     >  7:305101/1‑148 (MQ=255)            gCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGa                                                                                                                                                      <  3:55697/138‑1 (MQ=255)            gCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGa                                                                                                                                                      >  4:55697/1‑138 (MQ=255)             cccATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGAcc                                                                                                                                          >  1:38276/1‑149 (MQ=255)              caaTTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGAcc                                                                                                                                          >  3:161499/3‑148 (MQ=255)                    cgggCCGGCATGGTCCCGGATGACGTGGCCCTCATCGACACCGGCGCTGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTTCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGc                                                                                                                                   <  6:351644/147‑1 (MQ=255)                       gCCGGCATGGTCCCGGATGATGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCtt                                                                                                                                <  5:169818/149‑1 (MQ=255)                       gCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGGTCAGCGTTGGTCCAGAAGcg                                                                                                                                  <  5:52977/147‑2 (MQ=255)                         cGGCATGGTCCCGGATGACGAGGCGATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCTACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTaa                                                                                                                              <  2:50866/149‑1 (MQ=255)                                    cGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTa                                                                                                                   >  5:125835/1‑149 (MQ=255)                                            gTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACctgtgct                                                                                                           >  3:237887/1‑149 (MQ=255)                                              ggCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGTCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACctgtgctgt                                                                                                         >  1:244621/1‑149 (MQ=255)                                               gCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGACCTGTTGGCCGAAGGCCCGTCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACctgtgctgtg                                                                                                        <  7:188230/149‑1 (MQ=255)                                                 catcatCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTgg                                                                                                       >  4:129795/1‑148 (MQ=255)                                                  atcatcGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCa                                                                                                     <  6:125835/149‑1 (MQ=255)                                                     atcGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGGTCAGCGCTGGTCCAGAAGcgt                                                                                                                                 >  4:205602/1‑116 (MQ=39)                                                     atcGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGGTCAGCGCTGGTCCAGAAGcgt                                                                                                                                 <  3:205602/118‑3 (MQ=39)                                                     atcGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATaa                                                                                                  <  2:38276/149‑1 (MQ=255)                                                      tcGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCAGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAg                                                                                                 >  2:174799/1‑149 (MQ=255)                                                         acacCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTcc                                                                                              >  8:74710/1‑149 (MQ=255)                                                           acCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTgg                                                                                                       >  5:153331/1‑138 (MQ=255)                                                           acCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTgg                                                                                                       <  6:153331/138‑1 (MQ=255)                                                            ccGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCcgac                                                                                           <  7:300255/149‑1 (MQ=255)                                                             cGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCcgacg                                                                                          <  4:161499/149‑1 (MQ=255)                                                              ggCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCcgacga                                                                                         >  2:18704/1‑149 (MQ=255)                                                                           cgcgTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCttt                                                                            <  6:139909/149‑1 (MQ=255)                                                                                   ttGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGtt                                                                       <  7:332044/146‑1 (MQ=255)                                                                                   ttGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGtt                                                                       >  8:332044/1‑146 (MQ=255)                                                                                   ttGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTtgt                                                                     >  7:102208/1‑148 (MQ=255)                                                                                     gCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGa                                                                   <  3:129795/148‑1 (MQ=255)                                                                                     gCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGa                                                                   >  1:265418/1‑148 (MQ=255)                                                                                       aaCGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGaaaa                                                                <  5:290398/149‑1 (MQ=255)                                                                                         cGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGAGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGaaaaa                                                               >  7:234631/1‑148 (MQ=255)                                                                                           gCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAAc                                                            >  2:48207/1‑149 (MQ=255)                                                                                                 gggTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGaa                                                       <  8:305101/148‑1 (MQ=255)                                                                                                               tGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTtg                                                                      <  4:82729/119‑1 (MQ=255)                                                                                                               tGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTtg                                                                      >  3:82729/1‑119 (MQ=255)                                                                                                                 ttGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACt                                          <  2:244621/145‑1 (MQ=255)                                                                                                                    gCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTc                                   <  1:48207/149‑1 (MQ=255)                                                                                                                           gCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAGAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCt                            <  1:174799/149‑1 (MQ=255)                                                                                                                              cggccggcGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCt                         <  8:326820/149‑1 (MQ=255)                                                                                                                                gccggcGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAATAAGGCTGAACTATCGTTCCACTTCTGCttt                       <  7:74710/149‑1 (MQ=255)                                                                                                                                   ggcGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTa                    >  6:107334/1‑149 (MQ=255)                                                                                                                                      gcTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCg                 <  5:107334/149‑1 (MQ=255)                                                                                                                                        tGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGGACTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCt               >  1:93764/1‑149 (MQ=255)                                                                                                                                           aagCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTTTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTg            <  1:149373/147‑1 (MQ=255)                                                                                                                                             gCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGcc          >  6:336310/1‑149 (MQ=255)                                                                                                                                             gCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGcc          <  6:267878/149‑1 (MQ=255)                                                                                                                                               ttCTGGGTCAGCGCTGGTCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCt         >  1:12442/1‑148 (MQ=255)                                                                                                                                                 cTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTga       <  8:102208/148‑1 (MQ=255)                                                                                                                                                  tGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAACTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTAACGTTCAACTTCTGCTATCTACCGCTTTGCCt         >  4:76333/1‑145 (MQ=255)                                                                                                                                                    gACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTACGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGAt    >  1:99345/1‑148 (MQ=255)                                                                                                                                                     aCCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGGCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGa  <  4:237887/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGA  >  NC_002947/855716‑856012

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