Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	418,582	T→G	9.3%	L121R (CTG→CGG) ‡	PP_0345 →	heptose kinase WapP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	418,582	0	T	G	9.3%	105.3 / 4.5	43	L121R (CTG→CGG) ‡	PP_0345	heptose kinase WapP
Reads supporting (aligned to +/- strand):  ref base T (17/22);  new base G (2/2);  total (19/24)
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.79e-01

GGCCGAATGGCAAAACCTCAAGCAGTTCGCCAAATGGGGTATCCCCACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGAT  >  NC_002947/418433‑418721                                                                                                                                                      |                                                                                                                                            ggCCGAATGGCAAAACCTCAAGCAGTTCGCCAAATGGGGTATCCCCACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCt                                                                                                                                             >  2:530703/1‑150 (MQ=255)   ccGAATGGCAAAACCTCAAGCAGTTCGCCAAATGGGGTATCCCCACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTgg                                                                                                                                           <  1:499501/150‑1 (MQ=255)    cGAATGGCAAAACCTCAAGCAGTTCGCCAAATGGGGTATCCCCACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGc                                                                                                                                          >  1:694621/1‑150 (MQ=255)           cAAAACCTCAAGCAGTTCGCCAAATGGGGTATCCCCACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCAGGGCCGacag                                                                                                                                    <  1:628239/149‑4 (MQ=255)           cAAAACCTCAAGCAGTTCGCCAAATGGGGTATCCCCACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCAGGGCCGacag                                                                                                                                    >  2:628239/1‑146 (MQ=255)            aaaaCCTCAAGCAGTTCGCCAAATGGGGTATCCCCACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAAcaca                                                                                                                                  >  2:170843/1‑150 (MQ=255)            aaaaCCTCAAGCAGTTCGCCAAATGGGGTATCCCCACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCAGGGCCGaca                                                                                                                                     >  2:562714/1‑145 (MQ=255)            aaaaCCTCAAGCAGTTCGCCAAATGGGGTATCCCCACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCAGGGCCGaca                                                                                                                                     <  1:562714/147‑3 (MQ=255)                 cTCAAGCAGTTCGCGAAATGGGGTATCCCCACTGCCGAAGTCGTGGCCTGGGGCCTGGACCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGAt                                                                                                                             <  2:83223/150‑1 (MQ=255)                 cTCAAGCAGTTCGCCAAATGGGGTATCCCCACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGAt                                                                                                                             <  2:582954/150‑1 (MQ=255)                             gCCAAATGGGGTATCCCCACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGcc                                                                                                                 >  1:522467/1‑150 (MQ=255)                               cAAATGGGGTATCCCAACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGTACCGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACCGCCGATCTCCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGa                                                                                                               <  1:450223/150‑1 (MQ=255)                                aaaTGGGGTATCCCCACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGAc                                                                                                              <  1:583255/150‑1 (MQ=255)                                   tGGGGTATCCCCACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACcgc                                                                                                           <  2:463496/150‑1 (MQ=255)                                        tATCCCCACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTg                                                                                                      >  1:223670/1‑150 (MQ=255)                                          tCCCCACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTggg                                                                                                    >  1:6408/1‑150 (MQ=255)                                             ccACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCACCGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTg                                                                                                                                            <  1:699619/107‑1 (MQ=255)                                             ccACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTg                                                                                                                                            >  2:699619/1‑107 (MQ=255)                                               ccTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGAc                                                                                               <  1:27928/149‑1 (MQ=255)                                                      aaGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCa                                                                                        <  1:170843/150‑1 (MQ=255)                                                                   gggcctggAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGgc                                                                           >  2:694802/1‑150 (MQ=255)                                                                             cgcAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATAcc                                                                 <  1:738517/150‑1 (MQ=255)                                                                                                   gTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATg                                                                                           <  1:154390/102‑1 (MQ=255)                                                                                                   gTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATg                                                                                           >  2:154390/1‑102 (MQ=255)                                                                                                        gcgATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGcc                                      >  2:692959/1‑150 (MQ=255)                                                                                                        gcgATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATCCGCGCTTACCACCCTTCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGcc                                      <  1:530703/150‑1 (MQ=255)                                                                                                           atgatCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCAt                                   <  2:522467/150‑1 (MQ=255)                                                                                                               tCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCACATAACg                               >  1:257495/1‑150 (MQ=255)                                                                                                                 aCCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGGTCCTGGGGCGACGGTGTCAGTCGTTAGCGGGCGTGGCTTGCCCGGGTGATGCACGGACACCGCCTGGTCCATAAAGAc                             >  2:180535/1‑150 (MQ=255)                                                                                                                     gTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGa                         >  1:230105/1‑150 (MQ=255)                                                                                                                           tGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGgc                   <  2:223670/150‑1 (MQ=255)                                                                                                                               gCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCaa               >  2:113692/1‑150 (MQ=255)                                                                                                                                cacaCCGAAGACCTGTCGGCCCTGGCCCAACACAACGCTGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCCTGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAAc              <  1:180535/150‑1 (MQ=255)                                                                                                                                 acacCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAAcc             >  1:63189/1‑150 (MQ=255)                                                                                                                                    ccGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACctgc          <  2:694621/150‑1 (MQ=255)                                                                                                                                     cGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACctgct         <  2:230105/150‑1 (MQ=255)                                                                                                                                       aaGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTgg       >  1:525496/1‑150 (MQ=255)                                                                                                                                           ccTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGa   <  1:251554/150‑1 (MQ=255)                                                                                                                                           ccTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGa   <  1:827466/150‑1 (MQ=255)                                                                                                                                            cTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAAc              >  2:685650/1‑138 (MQ=255)                                                                                                                                            cTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAAc              <  1:685650/138‑1 (MQ=255)                                                                                                                                            cTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGAt  <  1:694802/150‑1 (MQ=255)                                                                                                                                                   cccTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATAcc                                                                 >  1:220557/1‑80 (MQ=255)                                                                                                                                                   cccTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATAcc                                                                 <  2:220557/80‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                            GGCCGAATGGCAAAACCTCAAGCAGTTCGCCAAATGGGGTATCCCCACTGCCGAAGTCGTGGCCTGGGGCCTGGAGCGCAACGGCCTGGCCTTTGGCCGTGGCGCGATGATCACCCGTGAACTGCCGCACACCGAAGACCTGTCGGCCCTGGCCGAACACAACGATGCGCGCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGAT  >  NC_002947/418433‑418721

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