Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	855,865	C→T	8.0%	T239I (ACC→ATC) ‡	murI →	glutamate racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	855,865	0	C	T	8.0%	140.1 / 2.2	50	T239I (ACC→ATC) ‡	murI	glutamate racemase
Reads supporting (aligned to +/- strand):  ref base C (21/25);  new base T (2/2);  total (23/27)
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.54e-01

CTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGA  >  NC_002947/855718‑856012                                                                                                                                                    |                                                                                                                                                    cTTTCCTTCGCCCAGTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGAc                                                                                                                                                     >  7:305101/1‑148 (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                                                                                                                                 <  3:205602/118‑3 (MQ=39)                                                   atcGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGGTCAGCGCTGGTCCAGAAGcgt                                                                                                                                 >  4:205602/1‑116 (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                                                                       >  8:332044/1‑146 (MQ=255)                                                                                 ttGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGtt                                                                       <  7:332044/146‑1 (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)                                                                                                                                                    |                                                                                                                                                    CTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTGCCAATGATCTGTTGGCCGAAGGCCCGGCCGGCGCTGCACGCTTCTGGACCAGCGCTGACCCAGAAGCCTTAAGAAAAATCCTACCTGTGCTGTGGCATAAGTCCGACGATGTGCAAAGCTTTGCGTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGA  >  NC_002947/855718‑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: ‑