Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,185,401	G→T	7.5%	R353R (CGC→CGA)	aceK ←	isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,185,401	0	G	T	7.5%	166.5 / 4.1	53	R353R (CGC→CGA)	aceK	isocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base G (26/23);  new base T (2/2);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

GCGCCACCTCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAGGTGG  >  NC_002947/5185259‑5185542                                                                                                                                               |                                                                                                                                              gcgcCACCTCCATCAACTCGGCCAGGCAATCCGGCTCTAAATTGCTGCGCGTAAAGCTAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATAC‑TGTCGATCACCGTGGCGCGATCCa                                                                                                                                        <  7:52449/149‑1 (MQ=255)        cTCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGATCGCGCCACGGtga                                                                                                                                  <  3:233710/149‑3 (MQ=255)         tCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGtttt                                                                                                                                 >  1:350626/1‑149 (MQ=255)            aGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGa                                                                                                                              <  8:344960/149‑1 (MQ=255)             gCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGAt                                                                                                                             <  7:65034/149‑1 (MQ=255)               aaCTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATgg                                                                                                                           >  4:91177/1‑149 (MQ=255)                 cTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGc                                                                                                                          >  1:133729/1‑148 (MQ=255)                 cTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCg                                                                                                                         <  3:91177/149‑1 (MQ=255)                 cTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGATCGCGCCACGGtga                                                                                                                                  <  7:313226/140‑3 (MQ=255)                 cTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGATCGCGCCACGGtga                                                                                                                                  >  8:313226/1‑138 (MQ=255)                      ccAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAg                                                                                                                    <  2:59200/149‑1 (MQ=255)                              tCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCt                                                                                                             <  3:136139/148‑1 (MQ=255)                              tCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCt                                                                                                             >  4:136139/1‑148 (MQ=255)                                cGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTtg                                                                                                           <  7:23187/148‑1 (MQ=255)                                    tCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTgatgat                                                                                                      >  5:322031/1‑149 (MQ=255)                                      gAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCt                                                                                                    >  1:237182/1‑149 (MQ=255)                                          ttGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGAt                                                                                                                             <  2:235508/120‑1 (MQ=255)                                          ttGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGAt                                                                                                                             >  1:235508/1‑120 (MQ=255)                                               gctcgGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTAGACACTCTTCCACAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTAGATGATCTTGAACACCg                                                                                           >  5:259148/4‑149 (MQ=255)                                                       aGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAg                                                                                   <  2:350626/149‑1 (MQ=255)                                                           aaaaTCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCgg                                                                               <  2:237182/149‑1 (MQ=255)                                                                 ggCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCGACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAg                                                                          <  2:133729/148‑1 (MQ=255)                                                                    aaaCTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACCTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTg                                                                       >  4:309051/1‑148 (MQ=255)                                                                      aCTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCg                                                                                >  5:152202/1‑137 (MQ=255)                                                                      aCTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCg                                                                                <  6:152202/137‑1 (MQ=255)                                                                        tCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACAGTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGt                                                                        <  3:144295/143‑1 (MQ=255)                                                                             ggTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCa                                                             <  5:353562/149‑1 (MQ=255)                                                                                 tCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGAc                                                         >  3:208564/1‑149 (MQ=255)                                                                                   ggCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCa                                                       >  6:248104/1‑149 (MQ=255)                                                                                           gCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCg                                                >  6:238721/1‑148 (MQ=255)                                                                                           gCCCCACCAGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCg                                                >  7:173568/1‑148 (MQ=255)                                                                                             cccACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGcccc                                                 >  5:36539/1‑145 (MQ=255)                                                                                             cccACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGcccc                                                 <  6:36539/145‑1 (MQ=255)                                                                                              ccACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAg                                                                                   >  8:313221/1‑110 (MQ=255)                                                                                              ccACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAg                                                                                   <  7:313221/110‑1 (MQ=255)                                                                                              ccACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGccc                                                  <  6:129870/143‑1 (MQ=255)                                                                                              ccACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGccc                                                  >  5:129870/1‑143 (MQ=255)                                                                                                    ggTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGt                                      <  5:248104/149‑1 (MQ=255)                                                                                                        gACACTCTTCACCAGCCGATACATGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCACGCACCCCGGGTGcc                                   >  7:236929/1‑148 (MQ=255)                                                                                                          cacTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCAt                                 <  4:208564/148‑1 (MQ=255)                                                                                                                  aCCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGc                         <  8:173568/148‑1 (MQ=255)                                                                                                                    cAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACCCGCTCATGACCATGCCCCGCACCACTGGTGCCATGACAAAGCGGt                      >  8:328395/1‑149 (MQ=255)                                                                                                                    cAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGt                      >  2:24654/1‑149 (MQ=255)                                                                                                                     aGCCGATACCTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCTGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGtc                     <  3:309051/149‑1 (MQ=255)                                                                                                                      gCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCAt                                 >  6:373914/1‑136 (MQ=255)                                                                                                                      gCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCAt                                 <  5:373914/136‑1 (MQ=255)                                                                                                                           tACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCAct               >  1:231783/1‑149 (MQ=255)                                                                                                                              ttGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGCACCCGCTCATGACCATGCCCCGCAACCCTGGTGCCCTGACAACGCGGGCATCCctgct            >  8:201040/1‑149 (MQ=255)                                                                                                                              ttGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCActgct            >  3:369011/1‑149 (MQ=255)                                                                                                                                gTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGc                                                                <  7:125397/95‑1 (MQ=255)                                                                                                                                gTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGc                                                                >  8:125397/1‑95 (MQ=255)                                                                                                                                      aaCCGTGGATCGCGCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAggt    >  6:171017/2‑149 (MQ=255)                                                                                                                                        ccGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAggtgg  >  6:193477/1‑149 (MQ=255)                                                                                                                                               gcgATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCg                        <  3:206934/120‑1 (MQ=255)                                                                                                                                               gcgATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCg                        >  4:206934/1‑120 (MQ=255)                                                                                                                                               |                                                                                                                                              GCGCCACCTCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAGGTGG  >  NC_002947/5185259‑5185542

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑