Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,602,371	C→T	G126E (GGG→GAG)	ftsE ←	cell division ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,602,371	0	C	T	100.0%	198.5 / NA	50	G126E (GGG→GAG)	ftsE	cell division ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/22);  total (28/22)

GCCTTCCGACAGCGCGTCGTCCAGGTTACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTGGCGGCGC  >  NC_000913/3602223‑3602517                                                                                                                                                     |                                                                                                                                                   gCCTTCCGACAGCGCGTCGTCCAGGTTACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCg                                                                                                                                                  <  1:115882/151‑1 (MQ=255)   cTTCCGACAGCGCGTCGTCCAGGTTACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGAc                                                                                                                                                >  2:389287/1‑151 (MQ=255)   cTTCCGACAGCGCGTCGTCCAGGTTACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGAc                                                                                                                                                >  1:420808/1‑151 (MQ=255)     tCCGACAGCGCGTCGTCCAGGTTACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACt                                                                                                                                               <  1:299598/150‑1 (MQ=255)      ccGACAGCGCGTCGTCCAGGTTACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACttt                                                                                                                                             <  1:599459/151‑1 (MQ=255)           aGCGCGTCGTCCAGGTTACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATcc                                                                                                                                         >  2:169152/1‑150 (MQ=255)           aGCGCGTCGTCCAGGTTACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCa                                                                                                                                        <  1:79345/151‑1 (MQ=255)                              cAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGcg                                                                                                                     <  1:366047/151‑1 (MQ=255)                              cAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGcg                                                                                                                     <  1:4011/151‑1 (MQ=255)                               aGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGcg                                                                                                                     >  2:370490/1‑150 (MQ=255)                                      tCGTCCGCCAGCAGTACCGCGGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAata                                                                                                              <  1:150601/150‑1 (MQ=255)                                            gCCAGCAGTACCGCGGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCAc                                                                                                       <  2:29716/151‑1 (MQ=255)                                             ccAGCAGTACCGCGGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCAcc                                                                                                      <  2:287551/151‑1 (MQ=255)                                                  agTACCGCGGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTgg                                                                                                 >  2:377280/1‑151 (MQ=255)                                                      ccGCGGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCAcc                                                                                             >  1:509829/1‑151 (MQ=255)                                                         cgGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGc                                                                                          >  1:426663/1‑151 (MQ=255)                                                          gggCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCg                                                                                         >  2:595580/1‑151 (MQ=255)                                                               tGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATaa                                                                                     <  1:608897/150‑1 (MQ=255)                                                                  tCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCa                                                                                  <  1:108917/150‑1 (MQ=255)                                                                     ccaccGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCgg                                                                              >  1:543446/1‑151 (MQ=255)                                                                                  aaTGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTa                                                                 <  2:213614/151‑1 (MQ=255)                                                                                     gCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCg                                                              >  2:62495/1‑151 (MQ=255)                                                                                        aaCACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGt                                                             >  2:15422/1‑149 (MQ=255)                                                                                             gCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGt                                                      <  2:317692/151‑1 (MQ=255)                                                                                             gCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGt                                                      >  2:283616/1‑151 (MQ=255)                                                                                                 ttgttCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACgg                                                  >  1:415658/1‑151 (MQ=255)                                                                                                      cACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCa                                              <  2:168680/150‑1 (MQ=255)                                                                                                        ccGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATc                                            >  1:406837/1‑150 (MQ=255)                                                                                                        ccGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCa                                           <  1:389590/151‑1 (MQ=255)                                                                                                            ccGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTAc                                                    <  1:591928/138‑1 (MQ=255)                                                                                                            ccGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTAc                                                    >  2:591913/1‑138 (MQ=255)                                                                                                            ccGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGt                                         >  2:173694/1‑149 (MQ=255)                                                                                                             cGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGa                                      <  2:124163/151‑1 (MQ=255)                                                                                                              gAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGt                                         >  1:215646/1‑147 (MQ=255)                                                                                                               aaaGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATg                                    >  1:354234/1‑151 (MQ=255)                                                                                                                      aaTAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATcc                             >  2:396738/1‑151 (MQ=255)                                                                                                                        tAGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTg                           >  2:451133/1‑151 (MQ=255)                                                                                                                         aGGGAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTgg                          <  1:173699/151‑1 (MQ=255)                                                                                                                            gAAGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGaa                        >  2:222694/1‑150 (MQ=255)                                                                                                                             aaGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGaaa                       >  2:70360/1‑150 (MQ=255)                                                                                                                              aGTTCTTCGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAAt                     >  2:632475/1‑151 (MQ=255)                                                                                                                                     cGCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGcc               >  2:135130/1‑150 (MQ=255)                                                                                                                                      gCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATg                 <  2:487484/147‑1 (MQ=255)                                                                                                                                      gCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATg                 >  1:487494/1‑147 (MQ=255)                                                                                                                                      gCTTTGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCa              >  1:105019/1‑150 (MQ=255)                                                                                                                                         ttGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATc           <  1:83924/150‑1 (MQ=255)                                                                                                                                         ttGTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCt          >  2:53383/1‑151 (MQ=255)                                                                                                                                           gTCCAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTg         <  2:420801/150‑1 (MQ=255)                                                                                                                                              cAGTAGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTggcgg     >  1:32296/1‑151 (MQ=255)                                                                                                                                                  aGCTCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTGGCGgcgc  <  1:236764/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GCCTTCCGACAGCGCGTCGTCCAGGTTACCAGTCGGTTCGTCCGCCAGCAGTACCGCGGGCTTGTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTGGCGGCGC  >  NC_000913/3602223‑3602517

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