Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	104,352	A→G	E124G (GAG→GGG)	ftsA →	ATP‑binding cell division FtsK recruitment protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	104,352	0	A	G	100.0%	234.9 / NA	65	E124G (GAG→GGG)	ftsA	ATP‑binding cell division FtsK recruitment protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/32);  total (33/32)

CGCTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGAGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATAT  >  NC_000913/104226‑104481                                                                                                                               |                                                                                                                                  cgttttCTGGTGAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATATCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGtgct                                                                                                                       >  1:967223/4‑139 (MQ=255) cgcTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGc                                                                                                                                 >  2:156166/1‑129 (MQ=255)       cTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCAtgt                                                                                                                 <  1:121057/139‑1 (MQ=255)                acaTCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCa                                                                                                        >  1:1072913/1‑139 (MQ=255)                  aTCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                         <  2:506488/120‑1 (MQ=255)                  aTCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                         >  1:506488/1‑120 (MQ=255)                        tGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATg                                                                                                <  1:452749/139‑1 (MQ=255)                               aTGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGa                                                                                         <  2:1072913/139‑1 (MQ=255)                                tGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGa                                                                                         <  1:929172/138‑1 (MQ=255)                                tGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGa                                                                                         >  2:929172/1‑138 (MQ=255)                                           tGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAgg                                                                             <  1:369834/139‑1 (MQ=255)                                            ggTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                         <  2:453015/94‑1 (MQ=255)                                            ggTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                         >  1:453015/1‑94 (MQ=255)                                            ggTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAggg                                                                            >  1:1565346/1‑139 (MQ=255)                                                ccTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTAt                                                                                     >  2:918598/1‑126 (MQ=255)                                                ccTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTAt                                                                                     <  1:918598/126‑1 (MQ=255)                                                ccTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGaa                                                                               <  1:325965/132‑1 (MQ=255)                                                ccTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGaa                                                                               >  2:325965/1‑132 (MQ=255)                                                      tCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATc                                                                  <  1:842092/139‑1 (MQ=255)                                                       cTGAAGAAGAAGTGACGCGAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATc                                                                         >  1:926311/1‑131 (MQ=255)                                                       cTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATc                                                                         <  2:926311/131‑1 (MQ=255)                                                          aagaagaagTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGaa                                                                               >  2:823769/1‑122 (MQ=255)                                                          aagaagaagTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGaa                                                                               <  1:823769/122‑1 (MQ=255)                                                               gaagTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                         >  2:144560/1‑75 (MQ=255)                                                               gaagTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                         <  1:144560/75‑1 (MQ=255)                                                               gaagTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGa                                                                                                               >  1:1664453/1‑85 (MQ=255)                                                               gaagTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGa                                                                                                               <  2:1664453/85‑1 (MQ=255)                                                                     aCGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                         <  2:868986/69‑1 (MQ=255)                                                                     aCGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                         >  1:868986/1‑69 (MQ=255)                                                                     aCGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGa                                                                                                               >  2:388098/1‑79 (MQ=255)                                                                     aCGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGa                                                                                                               <  1:388098/79‑1 (MQ=255)                                                                               aTGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAggg                                                                            <  1:94286/104‑1 (MQ=255)                                                                               aTGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAggg                                                                            >  2:94286/1‑104 (MQ=255)                                                                                 gtgGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                         >  1:726983/1‑57 (MQ=255)                                                                                 gtgGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                         <  2:726983/57‑1 (MQ=255)                                                                                 gtgGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGc                                       >  2:704395/1‑139 (MQ=255)                                                                                    gAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCg                                                                                                                             <  1:1091600/50‑1 (MQ=255)                                                                                    gAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATAGGCATCg                                                                                                                             >  2:1091600/1‑50 (MQ=255)                                                                                     aaaaCGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCAtgt                                                                                                                 >  2:723298/1‑61 (MQ=255)                                                                                     aaaaCGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCAtgt                                                                                                                 <  1:723298/61‑1 (MQ=255)                                                                                         cgtcgtCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGa                                                                                >  2:106375/1‑90 (MQ=255)                                                                                         cgtcgtCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGa                                                                                <  1:106375/90‑1 (MQ=255)                                                                                         cgtcgtCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCggg                                    <  1:968848/134‑1 (MQ=255)                                                                                         cgtcgtCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAgg                                    >  2:968848/1‑134 (MQ=255)                                                                                         cgtcgtCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCaaaa                               >  1:240356/1‑139 (MQ=255)                                                                                            cgtcCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                         >  2:374358/1‑46 (MQ=255)                                                                                            cgtcCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                         <  1:374358/46‑1 (MQ=255)                                                                                               ccATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCAc                         <  2:1565346/139‑1 (MQ=255)                                                                                                 aTACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCt                       <  2:967223/139‑1 (MQ=255)                                                                                                   aCCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGa                     >  1:100646/1‑139 (MQ=255)                                                                                                    ccGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAgg                                    >  2:1590753/1‑123 (MQ=255)                                                                                                    ccGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAgg                                    <  1:1590753/123‑1 (MQ=255)                                                                                                       cgAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCAc                         <  2:625252/131‑1 (MQ=255)                                                                                                       cgAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCAc                         >  1:625252/1‑131 (MQ=255)                                                                                                       cgAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCAc                         >  1:1514233/1‑131 (MQ=255)                                                                                                       cgAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCAc                         <  2:1514233/131‑1 (MQ=255)                                                                                                         aaaTCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGa                                                                     >  1:488398/1‑85 (MQ=255)                                                                                                         aaaTCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGa                                                                     <  2:488398/85‑1 (MQ=255)                                                                                                                    tgtgCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGTAGGCAAAAGTGCACCTGATCACATGt             >  2:514771/1‑130 (MQ=255)                                                                                                                    tgtgCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGTAGGCAAAAGTGCACCTGATCACATGt             <  1:514771/130‑1 (MQ=255)                                                                                                                    tgtgCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATcaca                >  2:976549/1‑127 (MQ=255)                                                                                                                    tgtgCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATcaca                <  1:976549/127‑1 (MQ=255)                                                                                                                      tgCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGa                     <  1:341625/120‑1 (MQ=255)                                                                                                                      tgCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGa                     >  2:341625/1‑120 (MQ=255)                                                                                                                      tgCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGatat  >  2:478341/1‑139 (MQ=255)                                                                                                                               |                                                                                                                                  CGCTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGAGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATAT  >  NC_000913/104226‑104481

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