Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,123,177	T→C	T55A (ACG→GCG)	ftsN ←	essential cell division protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,123,177	0	T	C	100.0%	196.6 / NA	52	T55A (ACG→GCG)	ftsN	essential cell division protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/31);  total (21/31)

GCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTGAGGTGCTCTTTTTTCGCCGCGAAGGTGCCGGTTGGCTGCGGCGTACAT  >  NC_000913/4123029‑4123323                                                                                                                                                     |                                                                                                                                                   gCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCaa                                                                                                                                                  <  2:941955/151‑1 (MQ=255) gCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCaa                                                                                                                                                  <  1:957720/151‑1 (MQ=255) gCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCaa                                                                                                                                                  <  1:605849/151‑1 (MQ=255)     aaGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGa                                                                                                                                               >  1:809789/1‑150 (MQ=255)       ggTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGt                                                                                                                                            <  1:421007/151‑1 (MQ=255)          tCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTAc                                                                                                                                          >  1:278219/1‑150 (MQ=255)             gtgtGCGCACGCACTCCCGGCTTGCGGCTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTGCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGac                                                                                                                                      <  2:350446/151‑1 (MQ=255)                   gcacgcacTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGa                                                                                                                                <  1:403679/151‑1 (MQ=255)                        cacTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATaa                                                                                                                             >  1:769660/1‑149 (MQ=255)                         acTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAgg                                                                                                                          >  2:664115/1‑151 (MQ=255)                             ccGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTcac                                                                                                                      <  2:960932/151‑1 (MQ=255)                              cGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGCAGTCCGTTTCCGGTCACTTTCTGGCTTTGCGGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTcac                                                                                                                      >  1:199552/1‑150 (MQ=255)                               ggCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACaa                                                                                                                    <  1:37805/151‑1 (MQ=255)                                 cTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTcaca                                                                                                                     >  2:853798/1‑148 (MQ=255)                                   ggCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAAc                                                                                                                <  1:614169/151‑1 (MQ=255)                                      gACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAAcggc                                                                                                             <  1:612013/151‑1 (MQ=255)                                           tCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGgca                                                                                                         <  2:204145/150‑1 (MQ=255)                                                  ctTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGc                                                                                                 >  1:426609/1‑151 (MQ=255)                                                      aaTGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGAcc                                                                                             >  2:710850/1‑151 (MQ=255)                                                          tAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGGGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGaac                                                                                             >  2:246957/1‑145 (MQ=255)                                                           aGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAg                                                                                         <  2:445611/150‑1 (MQ=255)                                                                cAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGgag                                                                                   <  2:141611/151‑1 (MQ=255)                                                                 aGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCggg                                                                                     >  1:240515/1‑148 (MQ=255)                                                                   cGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGaa                                                                                 >  2:1048344/1‑150 (MQ=255)                                                                     ttcttcTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAAcc                                                                              >  1:160080/1‑151 (MQ=255)                                                                     ttcttcTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAAc                                                                               <  2:676072/150‑1 (MQ=255)                                                                         tctGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCgca                                                                           <  2:308464/150‑1 (MQ=255)                                                                         tctGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCgc                                                                            >  1:820891/1‑149 (MQ=255)                                                                         tctGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCgc                                                                            >  2:354365/1‑149 (MQ=255)                                                                            ggTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCgcaggca                                                                       <  1:818319/151‑1 (MQ=255)                                                                            ggTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCgcagg                                                                         <  2:278195/149‑1 (MQ=255)                                                                                 tggtggTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGAtt                                                                   <  1:594250/150‑1 (MQ=255)                                                                                     ggtAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGtt                                                              <  1:387929/151‑1 (MQ=255)                                                                                            cGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCttttt                                                       >  1:113622/1‑151 (MQ=255)                                                                                             gTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTc                                                      <  1:952199/151‑1 (MQ=255)                                                                                             gTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTc                                                      <  2:530310/151‑1 (MQ=255)                                                                                                 ccGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGt                                                   <  2:552658/150‑1 (MQ=255)                                                                                                 ccGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTg                                                  >  1:796667/1‑151 (MQ=255)                                                                                                       aCTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTGAGGTGc                                            >  1:388833/1‑151 (MQ=255)                                                                                                            cTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTGAGGTGCTCttt                                       >  1:685653/1‑151 (MQ=255)                                                                                                                  ttGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTGAGGTGCTCTTTTTTc                                   <  1:628239/149‑1 (MQ=255)                                                                                                                    gCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTGAGGTGCTCTTTTTTcgccg                               <  1:664149/151‑1 (MQ=255)                                                                                                                     cAGCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTGAGGTGCTCTTTTTTcgccgc                              >  2:642009/1‑151 (MQ=255)                                                                                                                       gCGTCTCGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTGAGGTGCTCTTTTTTCGCCGCGa                            <  2:15491/151‑1 (MQ=255)                                                                                                                             cGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTGAGGTGCTCTTTTTTCGCCGCGAAgg                         <  1:238342/148‑1 (MQ=255)                                                                                                                             cGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTGAGGTGCTCTTTTTTCGCCGCGAAGGTg                       <  2:581723/150‑1 (MQ=255)                                                                                                                             cGGACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTGAGGTGCTCTTTTTTCGCCGCGAAGGTGc                      >  1:180351/1‑151 (MQ=255)                                                                                                                              ggACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTGAGGTGCTCTTTTTTCGCCGCGAAGGTGcc                     <  1:184569/151‑1 (MQ=255)                                                                                                                               gACTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTGAGGTGCTCTTTTTTCGCCGCGAAGGTGCCg                    <  2:681498/151‑1 (MQ=255)                                                                                                                                aCTCTTCTTTCTTGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTGAGGTGCTCTTTTTTCGCCGCGAAGGTGCCgg                   <  2:820845/151‑1 (MQ=255)                                                                                                                                         tcttGTGATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTGAGGTGCTCTTTTTTCGCCGCGAAGGTGCCGGTTGGCTgc           <  1:177803/150‑1 (MQ=255)                                                                                                                                               gATGCGCAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTGAGGTGCTCTTTTTTCGCCGCGAAGGTGCCGGTTGGCTGCGGCGTAc    >  1:656853/1‑151 (MQ=255)                                                                                                                                                  gcgcAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTGAGGTGCTCTTTTTTCGCCGCGAAGGTGCCGGTTGGCTGCGGCGTACAt  <  2:747185/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCAGGCAGATTTCGTTGCTTTTTCCGTGAGGTGCTCTTTTTTCGCCGCGAAGGTGCCGGTTGGCTGCGGCGTACAT  >  NC_000913/4123029‑4123323

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