Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,951,043	T→G	14.7%	K657T (AAG→ACG)	ppk ←	polyphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,951,043	0	T	G	14.7%	85.3 / 11.8	41	K657T (AAG→ACG)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (15/20);  new base G (3/3);  total (18/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 = 3.64e-01

CTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGAATGCCCG  >  NC_002947/5950895‑5951190                                                                                                                                                     |                                                                                                                                                    cTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCtt                                                                                                                                                    <  1:461331/150‑1 (MQ=255)                      gcgcACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtccatcca                                                                                                                               <  1:515532/150‑1 (MQ=255)                          aCGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGTTCCATCCAGTCg                                                                                                                           >  1:783715/1‑150 (MQ=255)                          aCGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCg                                                                                                                           >  1:502183/1‑150 (MQ=255)                                  cccGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGa                                                                                                                   <  2:470416/150‑1 (MQ=255)                                   ccGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGAc                                                                                                                  <  1:778511/150‑1 (MQ=255)                                    cGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACa                                                                                                                 >  2:274291/1‑150 (MQ=255)                                          gCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACACTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGa                                                                                                           >  2:625072/1‑150 (MQ=255)                                            tGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATc                                                                                                         >  2:693854/1‑150 (MQ=255)                                              cAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGCCACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTg                                                                                                       <  1:701367/150‑1 (MQ=255)                                                 ggTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGctc                                                                                                    <  1:140561/150‑1 (MQ=255)                                                     cAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACACTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTcg                                                                                                <  1:625072/150‑1 (MQ=255)                                                             ggTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCAtt                                                                                        <  2:193156/150‑1 (MQ=255)                                                                   gTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCATGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGaa                                                                                  <  2:455517/150‑1 (MQ=255)                                                                    tCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGa                                                                                      <  2:254274/145‑1 (MQ=255)                                                                    tCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGa                                                                                      >  1:254274/1‑145 (MQ=255)                                                                            gTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAacac                                                                         <  2:757907/150‑1 (MQ=255)                                                                              aGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCt                                                                                                        <  1:462731/117‑1 (MQ=255)                                                                              aGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCt                                                                                                        >  2:462731/1‑117 (MQ=255)                                                                                ccTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAg                                                                                                                <  2:472554/107‑1 (MQ=255)                                                                                ccTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAg                                                                                                                >  1:472554/1‑107 (MQ=255)                                                                                 cTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCcgcg                                                                    >  1:842721/1‑150 (MQ=255)                                                                                  ttCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACAAGCGTGTCGAgac                                                                                                                                            <  1:853124/77‑3 (MQ=255)                                                                                  ttCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACAAGCGTGTCGAgac                                                                                                                                            >  2:853124/1‑75 (MQ=255)                                                                                           cttcttCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGTTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAg                                                          <  2:783715/150‑1 (MQ=255)                                                                                           cttcttCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACAAGCGTGTCGAgac                                                                                                                                            <  1:35094/68‑3 (MQ=38)                                                                                           cttcttCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACAAGCGTGTCGAgac                                                                                                                                            >  2:35094/1‑66 (MQ=38)                                                                                                    ccGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACATGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGcc                                                 >  1:150454/1‑150 (MQ=255)                                                                                                    ccGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATc                                                                                                         <  2:345844/94‑1 (MQ=255)                                                                                                    ccGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATc                                                                                                         >  1:345844/1‑94 (MQ=255)                                                                                                    ccGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGcc                                                 >  2:616800/1‑150 (MQ=255)                                                                                                                          tCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCTATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGAtg                           <  2:268419/150‑1 (MQ=255)                                                                                                                             atcgtgAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGAtgttg                        <  2:732406/145‑1 (MQ=255)                                                                                                                               cGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTtgtg                      <  1:693854/150‑1 (MQ=255)                                                                                                                               cGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCACTCGGCGCTGGACAGGGAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTtgtg                      >  1:550855/1‑150 (MQ=255)                                                                                                                                gggAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGc                     <  2:842721/150‑1 (MQ=255)                                                                                                                                      cAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGAtgttgt                       <  2:166973/142‑1 (MQ=255)                                                                                                                                      cAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGAtgttgt                       >  1:166973/1‑142 (MQ=255)                                                                                                                                        aacctcGACAAGCGTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGc                                                     >  1:826083/4‑110 (MQ=255)                                                                                                                                        aacctcGACAAGCGTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGc                                                     <  2:826083/107‑1 (MQ=255)                                                                                                                                          tctcGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCg           <  2:583644/150‑1 (MQ=255)                                                                                                                                                   gCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGCACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCGGGAATGCCCg  >  2:152316/1‑150 (MQ=255)                                                                                                                                                     ttGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGc                                                               >  2:270944/1‑87 (MQ=255)                                                                                                                                                     ttGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGc                                                               <  1:270944/87‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGAATGCCCG  >  NC_002947/5950895‑5951190

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