Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,662,224	G→A	8.0%	A110T (GCT→ACT)	mmgF →	2‑methylisocitrate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,662,224	0	G	A	8.0%	143.7 / 2.6	50	A110T (GCT→ACT)	mmgF	2‑methylisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base G (21/25);  new base A (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 = 4.10e-01

GGCATCACTGGCCTGGACGACGTGCTCACTGACGTGCGTCGCATCACCGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGCAGCGATGACAGCTTCGTGATCATGGC  >  NC_002947/2662077‑2662369                                                                                                                                                    |                                                                                                                                                  ggCATCACTGGCCTGGACGACGTGCTCACTGACGTGCGTCGCATCACCGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGgcgc                                                                                                                                                  >  4:92899/1‑149 (MQ=255)           gCCTGGACGTCGAGCTCACTGACGTGCGTCGCATCACCGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGccat                                                                                                                                         <  8:330891/148‑1 (MQ=255)           gCCTGGACGACGTGCTCACTGACGTGTGTCGCATCACCGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGccat                                                                                                                                         <  6:275783/148‑1 (MQ=255)           gCCTGGACGACGTGCTCACTGACGTGCGTCGCATCACCGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGccatc                                                                                                                                        <  8:276784/149‑1 (MQ=255)            ccTGGACGACGTGCTCACTGACGTGCGTCGCATCACCGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGccatcc                                                                                                                                       <  3:92899/149‑1 (MQ=255)                    acgTGCTCACTGACGTGCGTCGCATCACCGACGTGGGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAg                                                                                                                               <  4:280943/149‑1 (MQ=255)                        gCTCACTGACGTGCGTCGCATCACCGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGAc                                                                                                                            >  3:151507/1‑148 (MQ=255)                          tCACTGACGTGCGTCGCATCACCGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAg                                                                                                                         <  8:240692/149‑1 (MQ=255)                          tCACTGACGTGCGTCGCATCACCGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTAGGCGCTGCCGCCATCCATATCGAGGACCAg                                                                                                                         >  6:355250/1‑149 (MQ=255)                            aCTGACGTGCGTTGCATCACCGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGt                                                                                                                       >  2:127614/1‑149 (MQ=255)                            aCTGACGTGCGTCGCATCACCGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAg                                                                                                                         <  4:234851/147‑1 (MQ=255)                              tGACGTGCGTCGCATCACCGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCACTGCCGCCATCCATATCGAGGACCAGGTTg                                                                                                                     >  6:305720/1‑149 (MQ=255)                                aCGTGCGTCGCATCACCGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGgc                                                                                                                   <  5:241006/149‑1 (MQ=255)                                     cGTCGCATCACCGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGcc                                                                                                                >  5:58185/1‑147 (MQ=255)                                       tCGCATCACCGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGc                                                                                                            >  2:162946/1‑149 (MQ=255)                                           aTCACCGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTg                                                                                                        <  6:58185/149‑1 (MQ=255)                                               ccGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGc                                                                                                    >  3:103252/1‑149 (MQ=255)                                                cGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGc                                                                                                    >  4:76962/1‑148 (MQ=255)                                                  aCGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCAc                                                                                                 <  6:19142/149‑1 (MQ=255)                                                                  gctgctGGTGGATGTCGACACCGGCTTCTGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGaaa                                                                                 <  5:355250/149‑1 (MQ=255)                                                                      ctgGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGt                                                                             <  3:76962/149‑1 (MQ=255)                                                                       tggtggATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGt                                                                             >  7:304265/1‑148 (MQ=255)                                                                         gtggATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTc                                                                            <  4:151507/147‑1 (MQ=255)                                                                              tGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTcagccag                                                                      >  4:278850/1‑148 (MQ=255)                                                                               gTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTcagccag                                                                      >  6:110451/1‑147 (MQ=255)                                                                                tCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCcagcag                                                                   <  5:294299/149‑1 (MQ=255)                                                                                tCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTAGCCCGTACCGTGCGCTCCATGACCAAGGTCGGCGCAGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCcagcag                                                                   <  5:289268/149‑1 (MQ=255)                                                                                 cGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCcagcag                                                                   >  1:313326/1‑148 (MQ=255)                                                                                  gACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAgg                                                                  >  2:196913/1‑148 (MQ=255)                                                                                  gACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCACTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTcagcc                                                                        >  2:93947/1‑142 (MQ=255)                                                                                  gACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCACTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTcagcc                                                                        <  1:93947/142‑1 (MQ=255)                                                                                       cGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATgg                                                            >  6:373989/1‑149 (MQ=255)                                                                                          cttcggcgcttcggcTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCg                                                         <  1:127614/149‑1 (MQ=255)                                                                                               gcgcttcggcTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCg                                                     >  8:9925/1‑148 (MQ=255)                                                                                                  cttcggcTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATc                                                 <  1:162946/149‑1 (MQ=255)                                                                                                    tcggcTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATAGTCAGCCAGCAGGAAATGGTCGACCGTATCa                                                >  7:271840/1‑148 (MQ=255)                                                                                                     cggcTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCaa                                               <  3:96/148‑1 (MQ=255)                                                                                                      ggcTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCaaa                                              <  2:313326/148‑1 (MQ=255)                                                                                                         tttCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCATCAGGAAATGGTCGCCCGTATCAAAGCTg                                          >  5:289620/1‑149 (MQ=255)                                                                                                         tttCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTg                                          >  3:255568/1‑149 (MQ=255)                                                                                                            cAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGcc                                        <  4:103252/148‑1 (MQ=255)                                                                                                                  tGCCCGTACCGTGCGCTAGATGATCAAGTTCGGCACTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGAc                                  <  5:305720/148‑1 (MQ=255)                                                                                                                   gCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTgg                                                                                                                    <  3:60357/65‑1 (MQ=255)                                                                                                                   gCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTgg                                                                                                                    >  4:60357/1‑65 (MQ=255)                                                                                                                             tGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGCAGCGAt                      <  8:304265/149‑1 (MQ=255)                                                                                                                                gcTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGCAGCGATGAc                   >  2:306660/1‑149 (MQ=255)                                                                                                                                         gCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGCAGCGATGACAGCTTCGt           >  2:171805/1‑148 (MQ=255)                                                                                                                                          cAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGCAGCGATGACAGCTTCGTGa         <  5:110451/149‑1 (MQ=255)                                                                                                                                                 ggCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGCAGCGATGACAGCTTCGTGATCATgg   <  1:171805/148‑1 (MQ=255)                                                                                                                                                 ggCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGCAGCGATGACAGCTTCGTGATCATGgc  >  2:87200/1‑149 (MQ=255)                                                                                                                                                 ggCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGCAGCGATGACAGCTTCGTGATCATGgc  <  2:152161/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GGCATCACTGGCCTGGACGACGTGCTCACTGACGTGCGTCGCATCACCGACGTGTGCGACCTGCCGCTGCTGGTGGATGTCGACACCGGCTTCGGCGCTTCGGCTTTCAACGTTGCCCGTACCGTGCGCTCGATGAGCAAGTTCGGCGCTGCCGCCATCCATATCGAGGACCAGGTTGGCGCCAAGCGTTGCGGCCACCGTCCCAACAAGGAAATTGTCAGCCAGCAGGAAATGGTCGACCGTATCAAAGCTGCCGTCGACGCCCGCAGCGATGACAGCTTCGTGATCATGGC  >  NC_002947/2662077‑2662369

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