Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	666,076	T→C	I80V (ATC→GTC)	mrdB ←	cell wall shape‑determining protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	666,076	0	T	C	100.0%	217.8 / NA	55	I80V (ATC→GTC)	mrdB	cell wall shape‑determining protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/29);  total (26/29)

GATAAAGCGCGCAACCATCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCAGGGCGCT  >  NC_000913/665930‑666223                                                                                                                                                   |                                                                                                                                                    gATAAAGCGCGCAACCATCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGAcata                                                                                                                                                  <  1:1029809/150‑1 (MQ=255)   tAAAGCGCGCAACCATCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATag                                                                                                                                                 >  1:277770/1‑149 (MQ=255)    aaaGCGCGCAACCATCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATagag                                                                                                                                               <  2:712717/150‑1 (MQ=255)     aaGCGCGCAACCATCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATag                                                                                                                                                 >  1:451544/1‑147 (MQ=255)             aaCCATCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGcc                                                                                                                                     >  2:54704/1‑151 (MQ=255)             aaCCATCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGc                                                                                                                                      <  1:243346/150‑1 (MQ=255)               ccATCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGccc                                                                                                                                    >  2:868423/1‑150 (MQ=255)                 aTCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGc                                                                                                                                 >  1:612485/1‑151 (MQ=255)                   cAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCt                                                                                                                               >  2:702341/1‑151 (MQ=255)                    aGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCtt                                                                                                                              <  1:309131/151‑1 (MQ=255)                    aGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCtt                                                                                                                              <  2:483696/151‑1 (MQ=255)                      tGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTc                                                                                                                             >  1:299749/1‑150 (MQ=255)                       ggTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCa                                                                                                                            <  1:318575/150‑1 (MQ=255)                           cGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAAc                                                                                                                       <  2:152249/151‑1 (MQ=255)                            ggCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAAc                                                                                                                       >  1:1199891/1‑150 (MQ=255)                             gCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAAcgc                                                                                                                     <  2:405857/151‑1 (MQ=255)                                aTTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAAcgcg                                                                                                                    <  2:260580/149‑1 (MQ=255)                                 ttttGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTgg                                                                                                                 >  2:587241/1‑151 (MQ=255)                                    tGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTggag                                                                                                               <  2:1110200/150‑1 (MQ=255)                                    tGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTggag                                                                                                               <  1:287219/150‑1 (MQ=255)                                            tCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTgcg                                                                                                      >  2:878255/1‑151 (MQ=255)                                                      gAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATcacc                                                                                              <  2:1156876/149‑1 (MQ=255)                                                           cGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCatga                                                                                       <  2:505128/151‑1 (MQ=255)                                                             aaCAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCatgatg                                                                                     <  1:735857/151‑1 (MQ=255)                                                                aaTACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGAcc                                                                                  <  2:469168/151‑1 (MQ=255)                                                                   aCCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGAcc                                                                                  >  2:860693/1‑148 (MQ=255)                                                                    ccGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGCGGAATTTGCGCCATCACCACCAGGATGACCAgcc                                                                              >  2:232995/1‑149 (MQ=255)                                                                    ccGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGAc                                                                              <  2:522272/151‑1 (MQ=255)                                                                        ggTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCAt                                                                          <  2:743343/151‑1 (MQ=255)                                                                         gTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATc                                                                         <  2:683659/151‑1 (MQ=255)                                                                         gTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATc                                                                         <  1:924121/151‑1 (MQ=255)                                                                            cagccagcGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATcgcg                                                                      <  2:648839/151‑1 (MQ=255)                                                                                      tGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCg                                                             <  1:78096/150‑1 (MQ=255)                                                                                       gAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGAt                                                           >  2:671366/1‑151 (MQ=255)                                                                                           aCCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTa                                                       <  2:380864/151‑1 (MQ=255)                                                                                                      tGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCcatca                                             <  2:415392/150‑1 (MQ=255)                                                                                                       ggCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATg                                           <  2:1089281/151‑1 (MQ=255)                                                                                                       ggCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATg                                           <  2:952398/151‑1 (MQ=255)                                                                                                              aaaGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAatat                                    >  2:723930/1‑151 (MQ=255)                                                                                                              aaaGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAatat                                    >  1:1083749/1‑151 (MQ=255)                                                                                                              aaaGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAatat                                    <  1:587612/151‑1 (MQ=255)                                                                                                              aaaGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAatat                                    >  2:1132319/1‑151 (MQ=255)                                                                                                               aaGCATCTACCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATc                                   >  1:705869/1‑151 (MQ=255)                                                                                                                       aCCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGAcc                             <  1:311755/149‑1 (MQ=255)                                                                                                                       aCCGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGc                           >  1:78183/1‑151 (MQ=255)                                                                                                                         cGCCACCAGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTg                         <  2:984323/151‑1 (MQ=255)                                                                                                                            caccaGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGgcg                      >  2:574094/1‑151 (MQ=255)                                                                                                                             accaGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGgcgc                     >  2:534861/1‑151 (MQ=255)                                                                                                                               caGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTc                   >  2:1181046/1‑151 (MQ=255)                                                                                                                               caGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTc                   >  2:806327/1‑151 (MQ=255)                                                                                                                                aGCAAAATAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTcc                  >  2:438716/1‑151 (MQ=255)                                                                                                                                   aaaaTAATACAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTcc                  >  2:789620/1‑148 (MQ=255)                                                                                                                                        aataCAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAAcc          <  2:57408/151‑1 (MQ=255)                                                                                                                                         ataCAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCa         <  2:618488/151‑1 (MQ=255)                                                                                                                                          taCAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAAcc          >  2:28000/1‑149 (MQ=255)                                                                                                                                           aCAGATGACATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCAgg       >  2:1027912/1‑151 (MQ=255)                                                                                                                                              gatgacATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCAGGgcg    <  2:443042/151‑1 (MQ=255)                                                                                                                                                tgacATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCAGGGCGCt  >  2:68039/1‑151 (MQ=255)                                                                                                                                                tgacATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCAGGGCGCt  >  1:850149/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GATAAAGCGCGCAACCATCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAACAATACCGAGGTCCAGCCAGCGTTGAGCACCTTTAGAGATGGCACCGAAAGCATCTACCGCCACCAGCAAAATAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCAGGGCGCT  >  NC_000913/665930‑666223

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