Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,379,933	T→C	D22G (GAC→GGC)	zapE ←	divisome ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,379,933	0	T	C	100.0%	176.7 / NA	50	D22G (GAC→GGC)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/24);  total (26/24)

TGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGTCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTACACTACGCGATAT  >  NC_000913/3379788‑3380079                                                                                                                                                  |                                                                                                                                                   tGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGt                                                                                                                                               >  1:491446/1‑151 (MQ=255)  ggCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGt                                                                                                                                               <  2:723374/150‑1 (MQ=255)  ggCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGt                                                                                                                                               <  2:589869/150‑1 (MQ=255)     gTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTTGCCGGGTTGAt                                                                                                                                           <  2:418886/151‑1 (MQ=255)       aTGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGTGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATgg                                                                                                                                         >  1:500932/1‑151 (MQ=255)       aTGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATgg                                                                                                                                         <  2:255885/151‑1 (MQ=255)         gCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCt                                                                                                                                       <  2:568040/151‑1 (MQ=255)                 tCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCAt                                                                                                                               >  2:457508/1‑151 (MQ=255)                   ttCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCAtt                                                                                                                              >  1:215536/1‑150 (MQ=255)                    tCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTa                                                                                                                             >  2:157858/1‑150 (MQ=255)                    tCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTa                                                                                                                             >  2:5731/1‑150 (MQ=255)                     cgcgTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAg                                                                                                                           <  1:157825/151‑1 (MQ=255)                        gTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAgcg                                                                                                                         >  2:579150/1‑150 (MQ=255)                            aCCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCtt                                                                                                                     <  1:506562/150‑1 (MQ=255)                            aCCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTc                                                                                                                    >  2:376897/1‑151 (MQ=255)                                   aGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTAtt                                                                                                             <  1:85764/151‑1 (MQ=255)                                   aGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTAtt                                                                                                             <  1:211271/151‑1 (MQ=255)                                                   ccATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAAc                                                                                             <  1:702816/151‑1 (MQ=255)                                                         gCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTg                                                                                       <  1:462317/151‑1 (MQ=255)                                                             aCTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGCCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATg                                                                                   >  1:268531/1‑151 (MQ=255)                                                                 gTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCt                                                                               <  2:554494/151‑1 (MQ=255)                                                                          gCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTAt                                                                      <  2:690316/151‑1 (MQ=255)                                                                             ggtggCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTAtt                                                                     <  1:85381/149‑1 (MQ=255)                                                                             ggtggCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTcc                                                                   <  2:774639/151‑1 (MQ=255)                                                                             ggtggCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTc                                                                    <  1:286682/150‑1 (MQ=255)                                                                              gtggCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGTCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCt                                                                  <  2:487132/151‑1 (MQ=255)                                                                              gtggCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCt                                                                  <  1:63729/151‑1 (MQ=255)                                                                                ggCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTg                                                                >  2:253643/1‑151 (MQ=255)                                                                                  cGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGa                                                               >  2:196635/1‑150 (MQ=255)                                                                                        aTTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCg                                                        <  2:220200/151‑1 (MQ=255)                                                                                           gatgaGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTg                                                     >  1:514062/1‑151 (MQ=255)                                                                                              gaGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTgcgc                                                  >  2:438128/1‑151 (MQ=255)                                                                                                gTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAAGTATTGCGATGTTGGGGTAACGCATTGCATGATCTTATTGTTATTCCTTGAATAATCGGt                                                      <  1:563224/145‑1 (MQ=255)                                                                                                             aatTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTg                                   >  2:635868/1‑151 (MQ=255)                                                                                                                  ccAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGa                               >  1:93348/1‑150 (MQ=255)                                                                                                                    aGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGaaa                             >  1:135542/1‑150 (MQ=255)                                                                                                                    aGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGaaa                             >  1:358475/1‑150 (MQ=255)                                                                                                                     ggcggcTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGaaaaa                           <  1:769033/151‑1 (MQ=255)                                                                                                                     ggcggcTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGaaaa                            >  2:151354/1‑150 (MQ=255)                                                                                                                      gcggcTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGaaaaa                           >  2:444437/1‑150 (MQ=255)                                                                                                                      gcggcTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGaaaaa                           >  2:222136/1‑150 (MQ=255)                                                                                                                               ggCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTc                 >  1:672914/1‑151 (MQ=255)                                                                                                                               ggCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTc                 >  2:253357/1‑151 (MQ=255)                                                                                                                                 ccTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCt                >  1:293489/1‑150 (MQ=255)                                                                                                                                 ccTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTa               <  2:363064/151‑1 (MQ=255)                                                                                                                                  ctctTTTTGAACGTCGCCGGGTTGATGGCTGCATTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTac              >  1:756315/1‑151 (MQ=255)                                                                                                                                    ctTTTTGATCTTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTac              <  2:477037/149‑1 (MQ=255)                                                                                                                                      ttttGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTacac            <  2:479825/149‑1 (MQ=255)                                                                                                                                      ttttGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTACACt           >  2:584988/1‑150 (MQ=255)                                                                                                                                           aaCGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTACACTAcgcg      >  1:546018/1‑150 (MQ=255)                                                                                                                                             cgtcgcCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTACACTACGCGata   <  2:748373/151‑1 (MQ=255)                                                                                                                                              gtcgccGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTACACTACGCGatat  <  2:619870/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   TGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGTCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTACACTACGCGATAT  >  NC_000913/3379788‑3380079

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