Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	855,874	A→G	9.6%	D242G (GAC→GGC) ‡	murI →	glutamate racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	855,874	0	A	G	9.6%	135.2 / 5.6	52	D242G (GAC→GGC) ‡	murI	glutamate racemase
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base G (2/3);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

CCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGCG  >  NC_002947/855728‑856016                                                                                                                                                   |                                                                                                                                               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                                                                                                           <  6:153331/138‑1 (MQ=255)                                               acCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTgg                                                                                                           >  5:153331/1‑138 (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)                                                                                                                                     ctgggtcAGCGCTGGTCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTga           <  2:12442/142‑1 (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)                                                                                                                                              cgcTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGCg  >  8:227296/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                               CCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGCG  >  NC_002947/855728‑856016

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