Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	933,387	C→T	A55V (GCC→GTC)	ftsK →	DNA translocase at septal ring sorting daughter chromsomes

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	933,387	0	C	T	91.4%	205.7 / 6.3	70	A55V (GCC→GTC)	ftsK	DNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base T (34/30);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCA  >  NC_000913/933291‑933479                                                                                                 |                                                                                             tGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTg                                                                                           <  2:362260/100‑1 (MQ=255) tGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTg                                                                                           <  2:1223371/100‑1 (MQ=255)   gAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTggn                                                                                         <  2:1879259/100‑2 (MQ=255)    aGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCa                                                                                        >  1:548782/1‑100 (MQ=255)    aGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCa                                                                                        >  1:1819466/1‑100 (MQ=255)       cGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGa                                                                                     >  2:1412083/1‑100 (MQ=255)       cGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGa                                                                                     >  1:365259/1‑100 (MQ=255)        gTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGaa                                                                                    <  2:1351275/100‑1 (MQ=255)         ttGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAAc                                                                                   <  1:992692/100‑1 (MQ=255)          tGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAAcc                                                                                  <  2:1825333/100‑1 (MQ=255)              gATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTAtn                                                                              <  1:954962/100‑2 (MQ=255)              gATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATc                                                                              <  2:1769730/100‑1 (MQ=255)               aTCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATcc                                                                             <  1:994055/100‑1 (MQ=255)                 ccTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCAt                                                                           >  2:224972/1‑100 (MQ=255)                 ccTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCAt                                                                           >  1:472189/1‑100 (MQ=255)                 ccTTATTGTCCTGTTTGCCGTCTGGGTGATGGCTGCCTTACGAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCAt                                                                           <  2:750586/100‑1 (MQ=255)                     aTTGTCCTGTTTGCCGTCTGGTTGAGGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAAtt                                                                       <  1:1616042/100‑1 (MQ=255)                       tGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTa                                                                     <  2:1816464/100‑1 (MQ=255)                        gTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACTCCAGCTGTTCGCAAACGGTCTGGCTTGAACCTTGCGATAAttt                                                                      >  2:771136/1‑98 (MQ=255)                        gTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAg                                                                    >  1:500924/1‑100 (MQ=255)                         tCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAgg                                                                   <  1:131438/100‑1 (MQ=255)                          ccTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGt                                                                  >  1:867102/1‑100 (MQ=255)                                  ccGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGcc                                                          >  2:485665/1‑100 (MQ=255)                                   cGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGccc                                                         >  1:265603/1‑100 (MQ=255)                                       tGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTg                                                     >  2:238532/1‑100 (MQ=255)                                        ggTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGc                                                    >  2:699564/1‑100 (MQ=255)                                        ggTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGc                                                    >  1:1857291/1‑100 (MQ=255)                                        ggTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGc                                                    >  1:1472556/1‑100 (MQ=255)                                        ggTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGc                                                    >  2:584163/1‑100 (MQ=255)                                         gTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCg                                                   >  2:913319/1‑100 (MQ=255)                                         gTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCg                                                   <  2:395088/100‑1 (MQ=255)                                          ttGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGt                                                  <  2:1546171/100‑1 (MQ=255)                                               ggCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTg                                             >  2:638813/1‑100 (MQ=255)                                                 cTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGc                                           <  2:1083310/100‑1 (MQ=255)                                                 cTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGc                                           <  2:1863034/100‑1 (MQ=255)                                                  tctcTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCa                                          <  2:1713367/97‑1 (MQ=255)                                                  tGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGcn                                          <  1:1388054/100‑2 (MQ=255)                                                  tGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCa                                          <  2:1194346/100‑1 (MQ=255)                                                  tGCCTTACTAAGCTTTAACCCGTCGGGCCGCAGCTGGTCGCGAACGGCCTGGCATGAACCTATGCATAATTTAGGTGGGATGCCCGGTGCGTGGGTGGc                                           <  1:956375/99‑1 (MQ=255)                                                   gCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAg                                         <  2:1793462/100‑1 (MQ=255)                                                   gCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAg                                         >  2:1210276/1‑100 (MQ=255)                                                    ccTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGa                                        >  2:4500/1‑100 (MQ=255)                                                     cTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCaaga                                       >  1:1713356/1‑97 (MQ=255)                                                       tACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGAATGCCCGGTGCGTGGTTGGCAGATAc                                     <  2:1186964/100‑1 (MQ=255)                                                          tAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCt                                  >  1:1058456/1‑100 (MQ=255)                                                               tttAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCt                             >  2:1340852/1‑100 (MQ=255)                                                                ttAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCtt                            >  2:518781/1‑100 (MQ=255)                                                                 tAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCttt                           >  2:241359/1‑100 (MQ=255)                                                                  aaCCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTa                          >  1:639882/1‑100 (MQ=255)                                                                   aCCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAt                         >  2:801743/1‑100 (MQ=255)                                                                    cccTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAACTTAGGGGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCttt                           >  1:37306/1‑97 (MQ=255)                                                                     ccTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAttt                       >  2:529502/1‑100 (MQ=255)                                                                      cTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGTTGCGTGGTTGGCAGATACGCTGTTCTTTAtttn                      <  1:1172798/100‑2 (MQ=255)                                                                      cTTCAGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAtttt                      >  2:873998/1‑100 (MQ=255)                                                                       ttCGGACACCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAAGTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAttttt                     >  2:1820454/1‑100 (MQ=255)                                                                         cGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTgg                   <  1:1820443/100‑1 (MQ=255)                                                                          ggACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGc                  <  1:854235/100‑1 (MQ=255)                                                                           gACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCg                 >  1:64611/1‑100 (MQ=255)                                                                            aCCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGt                <  2:1543824/100‑1 (MQ=255)                                                                            aCCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGt                <  2:329976/100‑1 (MQ=255)                                                                                 aGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATg            <  1:1506672/99‑1 (MQ=255)                                                                                  gCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGc          <  1:4500/100‑1 (MQ=255)                                                                                  cgatctCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGc          <  2:1245525/95‑1 (MQ=255)                                                                                       tCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGACGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTaca     <  1:1240725/100‑1 (MQ=255)                                                                                        cGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTacac    >  1:1865589/1‑100 (MQ=255)                                                                                        cGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGGTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGcagatcg    >  1:1245521/1‑94 (MQ=255)                                                                                         gCAAACGGTCTGGCATGAACCTATCCATAATTTATGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACAcc   <  2:257940/100‑1 (MQ=255)                                                                                         gCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACAcc   <  1:119632/100‑1 (MQ=255)                                                                                          cAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGTCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCa  >  2:1014308/1‑100 (MQ=255)                                                                                          cAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGTCTTacac    >  1:776147/1‑98 (MQ=255)                                                                                                 |                                                                                             TGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACACCA  >  NC_000913/933291‑933479

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