Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,808,278	A→G	7.4%	D110G (GAC→GGC)	kdsA‑I →	2‑dehydro‑3‑deoxyphosphooctonate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,808,278	0	A	G	7.4%	147.3 / 2.0	54	D110G (GAC→GGC)	kdsA‑I	2‑dehydro‑3‑deoxyphosphooctonate aldolase
Reads supporting (aligned to +/- strand):  ref base A (29/21);  new base G (2/2);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

AACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGCTCAGTTCCTCGCACCTCAGGAAATGAAACACATCCTGGCCAAGTGCGAAGAGG  >  NC_002947/1808133‑1808421                                                                                                                                                   |                                                                                                                                                 aaCCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGGcat                                                                                                                                               >  6:434429/1‑149 (MQ=255) aaCCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGGcat                                                                                                                                               >  2:63425/1‑149 (MQ=255) aaCCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGAcat                                                                                                                                               >  5:359233/1‑149 (MQ=255)                  gTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCtt                                                                                                                             >  7:282121/1‑149 (MQ=255)                    ccTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCgcgacag                                                                                                                           >  5:8924/1‑143 (MQ=255)                     cTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTc                                                                                                                           <  6:359233/148‑1 (MQ=255)                     cTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTc                                                                                                                           <  5:395782/148‑1 (MQ=255)                      tACCGTGGCCCTGGCATGGAAGAGGGGATGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCt                                                                                                                         >  7:273590/1‑149 (MQ=255)                           tGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGAATCCGACATAATCCAGCT‑GCCGGCTTTCCTGTccc                                                                                                                    >  3:276366/1‑147 (MQ=255)                            ggCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTcgcg                                                                                                                   <  1:210195/149‑1 (MQ=255)                                 tGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGa                                                                                                              >  2:87587/1‑149 (MQ=255)                                  ggCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGAc                                                                                                             >  7:174341/1‑149 (MQ=255)                                   gCATGGAAGAGGGGCTGAAGATCTTCGAAGTGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGa                                                                                                              <  7:289172/147‑1 (MQ=255)                                       ggAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACatcatc                                                                                                                                           <  6:19528/115‑1 (MQ=255)                                       ggAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACatcatc                                                                                                                                           >  5:19528/1‑115 (MQ=255)                                        gAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGAcc                                                                                                        <  2:147187/148‑1 (MQ=255)                                                     aGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGtggccat                                                                                           <  4:240463/148‑1 (MQ=255)                                                       aTCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGtggccatggc                                                                                        <  8:273590/149‑1 (MQ=255)                                                            cGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCaaaa                                                                                   >  1:171515/1‑149 (MQ=255)                                                             gAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCaaaa                                                                                   <  1:87587/148‑1 (MQ=255)                                                              aaGAGATCAAGCGCACCTTCAACGTGCCAGGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAAc                                                                                  >  1:247810/1‑149 (MQ=255)                                                                  gaTCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCg                                                                             >  1:110529/1‑149 (MQ=255)                                                                   aTCAAGCGCACCTTCAACGTGCC‑GGTGATCAGCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCgt                                                                            >  5:83459/1‑149 (MQ=255)                                                                    tCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCgtg                                                                           <  8:282121/149‑1 (MQ=255)                                                                     cAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCgtgg                                                                          >  3:244854/1‑149 (MQ=255)                                                                      aaGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCgtgg                                                                          <  3:321801/148‑1 (MQ=255)                                                                             ccTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGGCATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAAc                                                                  <  1:63425/149‑1 (MQ=255)                                                                             ccTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCCTCGCCAAGGTATGCGGCATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAAc                                                                  <  5:434429/149‑1 (MQ=255)                                                                                        cc‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGc                                                       >  2:108992/1‑149 (MQ=255)                                                                                           ggTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCaag                                                            >  8:349331/1‑142 (MQ=255)                                                                                           ggTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCaag                                                            <  7:349331/142‑1 (MQ=255)                                                                                            gTGATAACCGATTTGCACGAGACTTAAAAGGCCGAGCCCGTCGACAAGTTATGCGACATCATACAGCT‑GCCGGCTTTCCTGTCGAGGCAGACTGACCTGGTTGTCGACATTGCCAAAAACGGCGTGGTCATCAACATCAACAAGGCt                                                      >  2:115416/1‑147 (MQ=255)                                                                                                 cACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCGACATCAAGAAGGCTCAGTTcc                                               >  6:73168/1‑149 (MQ=255)                                                                                                 cACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGCTCAGTTcc                                               >  2:313107/1‑149 (MQ=255)                                                                                                           cACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGCTCAGTTCCTCGCACCt                                       >  7:364643/1‑147 (MQ=255)                                                                                                            aCGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGCCGTGGTGATCAACATCAAGAAGGCTCAGTTCCTCGCACCTCAg                                    >  1:197107/1‑149 (MQ=255)                                                                                                                 ccTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGCTCAGTTCCTCGCACCTTAGGAAAt                               >  5:204789/1‑149 (MQ=255)                                                                                                                  cTTACCAGGCCGAGCCCGTCGCCAAGCTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTGGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGCTCAGTTCCTCGCACCTCAGGAAAt                               <  6:142452/148‑1 (MQ=255)                                                                                                                   ttACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCa                                                              >  7:362161/1‑116 (MQ=255)                                                                                                                   ttACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCa                                                              <  8:362161/116‑1 (MQ=255)                                                                                                                      ccAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGCTCAGTTCCTCGCACCTCAGGAAATGaaa                           <  2:247810/148‑1 (MQ=255)                                                                                                                       cAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGCTCAGTTCCTCGCACCTCAGGAAATGAAaca                         >  1:271295/1‑149 (MQ=255)                                                                                                                                cccGTCGCCAAGGTATGCGACATCATCCAGCTGGCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGCTCAGTTCCTCGCACCTCAGGAAATGAAACACATCCTgg                 >  1:273011/1‑149 (MQ=255)                                                                                                                                cccGTCGCCAAGGTATGCGACATCATCCAGCT‑TCCGGCTTTCCTTTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGCTCAGTTCCTCGCACCTCAGGAAATGAAACACATCCTGGc                <  4:276366/149‑1 (MQ=255)                                                                                                                                      gCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGCTCAGTTCCTCGCACCTCAGGAAATGAAACACATCCTGGCCAAGt           <  2:273011/148‑1 (MQ=255)                                                                                                                                        cAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGCTCAGTTCCTCGCACCTCAGGAAATGAAACACATCCTGACCAAGTGc         >  3:184246/1‑148 (MQ=255)                                                                                                                                            gTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCaa                                                             <  4:171852/92‑1 (MQ=255)                                                                                                                                            gTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCaa                                                             >  3:171852/1‑92 (MQ=255)                                                                                                                                            gTATGCGAAATCATCCAGCT‑GACGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACAGGAGTTGTGATCAACATAAAGAATGCTCAGTTACTCGAACCTCAGGAAATGAACCACATACTGGACAAGTGCg        >  8:349962/1‑145 (MQ=255)                                                                                                                                              aTGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGCTCAGTTCCTCGCACCTCAGGAAATGAAACACATCCTGGCCAAGTGCGaa      >  3:154358/1‑145 (MQ=255)                                                                                                                                              aTGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGCTCAGTTCCTCGCACCTCAGGAAATGAAACACATCCTGGCCAAGTGCGaa      <  4:154358/145‑1 (MQ=255)                                                                                                                                               tGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGCTCAGTTCCTCGCACCTCAGGAAATGAAACACATCCTGGCCAAGTGCGAAGAgg  <  8:364643/148‑1 (MQ=255)                                                                                                                                                 cGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGCTCAGTTCCTCGCACCTc                                      <  8:129852/110‑1 (MQ=255)                                                                                                                                                 cGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGCTCAGTTCCTCGCACCTc                                      >  7:129852/1‑110 (MQ=255)                                                                                                                                                   |                                                                                                                                                 AACCGTTCGTCGGTAACGTCCTACCGTGGCCCTGGCATGGAAGAGGGGCTGAAGATCTTCGAAGAGATCAAGCGCACCTTCAACGTGCC‑GGTGATCACCGATGTGCACGAGCCTTACCAGGCCGAGCCCGTCGCCAAGGTATGCGACATCATCCAGCT‑GCCGGCTTTCCTGTCGCGGCAGACTGACCTGGTGGTGGCCATGGCCAAAACCGGCGTGGTGATCAACATCAAGAAGGCTCAGTTCCTCGCACCTCAGGAAATGAAACACATCCTGGCCAAGTGCGAAGAGG  >  NC_002947/1808133‑1808421

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