Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	828,160	T→A	7.4%	A208A (GCT→GCA)	ppkB →	polyphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	828,160	0	T	A	7.4%	151.6 / 3.4	54	A208A (GCT→GCA)	ppkB	polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (27/23);  new base A (2/2);  total (29/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 = 6.99e-01

GATGAGCAGACCCGGCGCCTGCAGGACCGCATCAACGACGGCCGCAAGCTGTGGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATAACCCGCGACCAGCAGGTGAAGCTGCCCAAACGCGGCAAGATCGGCAAGT  >  NC_002947/828017‑828305                                                                                                                                                  |                                                                                                                                                  gATGAGCAGACCCGGCGCCTGCAGGACCGCATCAACGACGGCCGCAAGCTGTGGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGc                                                                                                                                              <  7:213406/149‑1 (MQ=255)       cAGACCCGGCGCCTGCAGGACCGCATCAACGACGGCCGCAAGCTGTGGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTg                                                                                                                                        >  4:123565/1‑149 (MQ=255)              gtcgcCTGCAGGACCGCATCAACGTCGGCCGCAAGCTGTGTAAGTTGTCGCCGATTG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATg                                                                                                                                 <  4:38007/147‑1 (MQ=255)              ggCGCCTGCAGGACCGCATCAACGACGGCCGCAAGCTGTGGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATg                                                                                                                                 >  8:291971/1‑149 (MQ=255)                 gcCTGCAGGACCGCATCAACGACGGCCGGAAGCTGTGGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCAGCTGGGATGAATACACCAAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATggg                                                                                                                               >  1:172056/1‑147 (MQ=255)                    tGCAGGACCGCATCAACGACGGCCGCAAGCTGTGGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAAt                                                                                                                           >  6:280568/1‑149 (MQ=255)                      cAGGACCGCATCAACGACGGCCGCAAGCTGTGGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAAtt                                                                                                                          <  3:123565/148‑1 (MQ=255)                         gACCGCATCAACGACGGCCGCAAGCTGTGGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCaa                                                                                                                      >  6:79361/1‑149 (MQ=255)                            cGCATCAACGACGGCCGCAAGCTGTGGAAGTTGTCGCCGATGGCTTATATAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAAc                                                                                                                     <  1:102254/149‑1 (MQ=255)                              cATCAACGACGGCCGCAAGCTGTGGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGAc                                                                                                                  >  4:68411/1‑148 (MQ=255)                                    cgacgGCCGCAAGCTGTGGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTc                                                                                                           <  3:68411/149‑1 (MQ=255)                                     gacgGCCGCAAGCTGTGGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTcg                                                                                                          >  1:154250/1‑149 (MQ=255)                                     gacgGCCGCAAGCTGTGGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATAAACCCAAGCGCGTGACGAAATGTTCCATGCGTCGGATACCGCTTGGGCACCGTGTTTAATGGCCAATTCCAACGACAAACGTcg                                                                                                          >  6:84398/1‑149 (MQ=255)                                        ggCCGCAAGCTGTGGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGGTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCTAGTCCAACGACAACCGTCgcgc                                                                                                       >  2:264574/1‑149 (MQ=255)                                          ccGCAAGCTGTGGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGccc                                                                                                     >  4:76661/1‑149 (MQ=255)                                                 cTGTGGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGaa                                                                                              >  6:75162/1‑149 (MQ=255)                                                 cTGTGGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGaa                                                                                              <  5:299729/149‑1 (MQ=255)                                                  tgtgGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGaa                                                                                              >  4:130218/1‑148 (MQ=255)                                                    tgGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAata                                                                                            <  2:300142/148‑1 (MQ=255)                                                         gTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATAtcatca                                                                                      <  7:291971/149‑1 (MQ=255)                                                          ttGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAg                                                                                     <  5:79361/149‑1 (MQ=255)                                                             tCGCCGATGG‑ACCT‑CAAGTAGTATAGCCGCTGGGATGAATACACCCAAGCGAGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGTCCAATTCCAACGACAAACGTCGCGCCCGGCTGAAGATCATCAGccc                                                                                  >  4:62150/1‑148 (MQ=255)                                                              cGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTg                                                                                                                                                 <  7:149049/85‑1 (MQ=255)                                                              cGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTg                                                                                                                                                 >  8:149049/1‑85 (MQ=255)                                                                          cT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCg                                                                      >  5:284686/1‑149 (MQ=255)                                                                               aGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTAtt                                                                  <  3:130218/149‑1 (MQ=255)                                                                                 gcGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTcc                                                                <  7:220829/148‑1 (MQ=255)                                                                                  cGTATAGCCGCTGGTATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTcc                                                                <  4:4074/148‑1 (MQ=255)                                                                                     ataGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTa                                                             >  1:115802/1‑148 (MQ=255)                                                                                      taGCCGCTGGTATGAATACACCCAAGCGCTTGACGAAATGTTCCACGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACCTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTAc                                                            <  5:84398/148‑1 (MQ=255)                                                                                      taGCCGCTGGGATGAATACACCCATGCGCGTGACGAAATGTTGCAGGCGTCGGATACGGCTTGGGCACCGTGGTTCATGGCAATTTCAAACTAAAAACGTGGAGCCAGGCTGATTTTCTTCTGCCACCTGATATGCCGTATTccagac                                                            >  8:286256/1‑144 (MQ=255)                                                                                      taGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTAc                                                            <  8:43710/148‑1 (MQ=255)                                                                                         ccGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAg                                                         >  7:40990/1‑148 (MQ=255)                                                                                          cGCTGGGATTAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGa                                                       <  6:284686/149‑1 (MQ=255)                                                                                              gggATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCATGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACAt                                                    >  3:173132/1‑148 (MQ=255)                                                                                              gggATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCATGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATa                                                   >  3:176916/1‑149 (MQ=255)                                                                                               ggATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATa                                                   >  5:286791/1‑148 (MQ=255)                                                                                                    aaTACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATAACCcgc                                             >  1:46413/1‑149 (MQ=255)                                                                                                           ccAAGCGAGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGTCACCGTGGTTCATGGACAATTCCAACGACAAACGTCGCGCCCGGCTTAATATCATCAGCCACCTGCTCATCCGTATTCCCTACAAGGACATAACCCGCGACcagc                                      <  1:239633/149‑1 (MQ=255)                                                                                                             aaGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATAACCCGCGACcagcag                                    >  4:194347/1‑149 (MQ=255)                                                                                                               gcgcGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATAACCCGCGACCAGCAGGt                                  <  5:280568/149‑1 (MQ=255)                                                                                                                        gAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATAACCCGCGACCAGCAGGTGAAGCTGcc                         >  4:105820/1‑149 (MQ=255)                                                                                                                        gAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATAACCCGCGACCAGCAGGTGAAGCTGcc                         >  3:270573/1‑149 (MQ=255)                                                                                                                         aaaTGTTCCAGGCGTCGGATACCGCATGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATAACCCGCGACCAGCAGGTGAAGCTGccc                        <  4:173132/149‑1 (MQ=255)                                                                                                                          aaTGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATAACCCGCGACCAGCAGGTGAAGCTGccc                        >  7:84651/1‑148 (MQ=255)                                                                                                                          aaTGTTCCAGGCGTCGGATACCGCATGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATAACCCGCGACCAGCAGGTGAAGCTGccc                        <  4:176916/148‑1 (MQ=255)                                                                                                                                 cAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATAACCCGCGACCAGCAGGTGAAGCTGCCCAAACGCgg                >  5:147156/1‑149 (MQ=255)                                                                                                                                    gCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATAACCCGCGACCAGCAGGTGAAGCTGCCCAAACGCGGCaa             <  5:75162/149‑1 (MQ=255)                                                                                                                                          gATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATAACCCGCGACCAGCAGGTGAAGCTGCCCAAACGCGGCAAGATCgg       >  6:97933/1‑149 (MQ=255)                                                                                                                                          gATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATAACCCGCGACCAGCAGGTGAAGCTGCCCAAACGCGGCAAGATCgg       >  7:85771/1‑149 (MQ=255)                                                                                                                                          gATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATAACCCGCGACCAGCAGGTGAAGCTGCCCAAACGCGGCAAGATCgg       <  2:154250/149‑1 (MQ=255)                                                                                                                                              ccGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATAACCCGCGACCAGCAGGTGAAGCTGCCCAAACGCGGCAAGATCGGCaa    <  2:46413/148‑1 (MQ=255)                                                                                                                                              cagCTTTGGCAACGTGGGTCATGGCCAAATCCAACGACAAAAGTCGCGCCCGGATTAAAATCATCACCCACCTCCTCAGCCGTATTCCCTACAAGGACATAACCCGCGACCAGCAGGTGAAGCTGCCCAAACGCGGCAAGATCGGCaa    <  7:286256/146‑1 (MQ=255)                                                                                                                                               cGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATAACCCGCGACCAGCAGGTGAAGCTGCCCAAACGCGGCAAGATCGGCAAGt  <  6:286791/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  GATGAGCAGACCCGGCGCCTGCAGGACCGCATCAACGACGGCCGCAAGCTGTGGAAGTTGTCGCCGATGG‑ACCT‑CAAGTCGTATAGCCGCTGGGATGAATACACCCAAGCGCGTGACGAAATGTTCCAGGCGTCGGATACCGCTTGGGCACCGTGGTTCATGGCCAATTCCAACGACAAACGTCGCGCCCGGCTGAATATCATCAGCCACCTGCTCAGCCGTATTCCCTACAAGGACATAACCCGCGACCAGCAGGTGAAGCTGCCCAAACGCGGCAAGATCGGCAAGT  >  NC_002947/828017‑828305

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