Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,398,875	C→T	M1M (GTG→ATG) † *368* (TAG→TAA)	mreD ← mreC ←	cell wall structural complex MreBCD transmembrane component MreD cell wall structural complex MreBCD transmembrane component MreC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,398,875	0	C	T	100.0%	196.4 / NA	61	M1M (GTG→ATG) † *368* (TAG→TAA)	mreD mreC	cell wall structural complex MreBCD transmembrane component MreD cell wall structural complex MreBCD transmembrane component MreC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/27);  total (34/27)

CAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGC  >  NC_000913/3398747‑3399010                                                                                                                                 |                                                                                                                                        caaGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTc                                                                                                                               >  2:125059/1‑139 (MQ=255)   aGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTg                                                                                                                                    <  1:1147275/132‑1 (MQ=255)   aGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTg                                                                                                                                    >  2:1147275/1‑132 (MQ=255)    gATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCg                                                                                                                            >  2:362869/1‑139 (MQ=255)    gATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCg                                                                                                                            >  2:1391230/1‑139 (MQ=255)    gATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCg                                                                                                                            >  2:1236310/1‑139 (MQ=255)         gTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCa                                                                                                                       >  1:1163765/1‑139 (MQ=255)                                           aTCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCgg                                                                                                            <  1:1119010/116‑1 (MQ=255)                                           aTCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCgg                                                                                                            >  2:1119010/1‑116 (MQ=255)                                             ccGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTa                                                                                   <  2:1163765/139‑1 (MQ=255)                                             ccGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGACGTa                                                                                   <  1:362869/139‑1 (MQ=255)                                                gCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCgg                                                                                                            >  2:561541/1‑111 (MQ=255)                                                gCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCgg                                                                                                            <  1:561541/111‑1 (MQ=255)                                                  cAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCt                                                                              <  2:1162869/139‑1 (MQ=255)                                                     ggCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGc                                                                                 <  2:287335/133‑1 (MQ=255)                                                     ggCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGc                                                                                 >  1:287335/1‑133 (MQ=255)                                                           aTTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAgg                                                                                                               <  2:485325/97‑1 (MQ=255)                                                           aTTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAgg                                                                                                               >  1:485325/1‑97 (MQ=255)                                                               gcagcaACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAgg                                                                                                               <  1:732692/93‑1 (MQ=255)                                                               gcagcaACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAgg                                                                                                               >  2:732692/1‑93 (MQ=255)                                                               gcagcaACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGa                                                                 >  2:657322/1‑139 (MQ=255)                                                               gcagcaACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCACTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGa                                                                 >  1:1504080/1‑139 (MQ=255)                                                                        gcgcAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcgg                                                                                          >  2:1406037/1‑105 (MQ=255)                                                                        gcgcAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcgg                                                                                          <  1:1406037/105‑1 (MQ=255)                                                                        gcgcAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAgcagca                                                        >  2:148100/1‑139 (MQ=255)                                                                        gcgcAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAgcagca                                                        <  1:125059/139‑1 (MQ=255)                                                                             aTGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCtg                                                   >  2:1214885/1‑139 (MQ=255)                                                                              tGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCtgt                                                  >  1:1325122/1‑139 (MQ=255)                                                                              tGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCtgt                                                  >  1:740165/1‑139 (MQ=255)                                                                                aGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGgcg                                                                                                                         >  1:1361273/1‑66 (MQ=255)                                                                                aGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGgcg                                                                                                                         <  2:1361273/66‑1 (MQ=255)                                                                                        gagCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCa                                                                                                                 >  2:1365813/1‑66 (MQ=255)                                                                                        gagCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCa                                                                                                                 <  1:1365813/66‑1 (MQ=255)                                                                                          gCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGAt                                                                <  1:792718/113‑1 (MQ=255)                                                                                          gCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGAt                                                                >  2:792718/1‑113 (MQ=255)                                                                                           ccAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGcagcag                                     >  2:524480/1‑139 (MQ=255)                                                                                             aGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGgc                                   <  1:1391230/139‑1 (MQ=255)                                                                                              gATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGgcg                                  <  1:657322/139‑1 (MQ=255)                                                                                               aTTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGgcgc                                 <  1:1236310/139‑1 (MQ=255)                                                                                                ttACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGTAGGCGCa                                <  1:1214885/139‑1 (MQ=255)                                                                                                  accatgtgacccTAGCTTCGATAGCTCGCCATTATTGCCCACCCGGCGCAGGGGCAGGAGGTTGAGCACCACTTTGCAGTGACGTGGCCCTTTGTGGAGAGCGATTAGCAGCAGGATGTGTCGGCGCAGCAGGCACAGt                              <  2:848102/130‑1 (MQ=255)                                                                                                   ccccGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTa                             <  2:304400/139‑1 (MQ=255)                                                                                                        cGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGc                        >  2:971154/1‑139 (MQ=255)                                                                                                         gTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTa                             <  1:28082/133‑1 (MQ=255)                                                                                                         gTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTa                             >  2:28082/1‑133 (MQ=255)                                                                                                             cTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCATTTTgcggcggc                                                                                      <  1:567196/72‑1 (MQ=255)                                                                                                             cTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcggcggc                                                                                      >  2:567196/1‑72 (MQ=255)                                                                                                                gCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTg                <  2:1504080/139‑1 (MQ=255)                                                                                                                     gATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAgg                                    >  2:996982/1‑114 (MQ=255)                                                                                                                     gATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAgg                                    <  1:996982/114‑1 (MQ=255)                                                                                                                        ggCTCGCCATGATTGCCCTCCCGGCGAACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGt        <  2:1140072/138‑1 (MQ=255)                                                                                                                          cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAg                                                                                                       >  2:1526788/1‑42 (MQ=255)                                                                                                                          cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAg                                                                                                       <  1:1526788/42‑1 (MQ=255)                                                                                                                          cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCa                                >  1:1166751/1‑113 (MQ=255)                                                                                                                          cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCa                                <  2:1166751/113‑1 (MQ=255)                                                                                                                          cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTa                             >  2:1431348/1‑116 (MQ=255)                                                                                                                          cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTa                             <  1:1431348/116‑1 (MQ=255)                                                                                                                          cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg      >  2:1242353/1‑139 (MQ=255)                                                                                                                          cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg      >  1:460366/1‑139 (MQ=255)                                                                                                                          cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCAGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg      >  2:967380/1‑139 (MQ=255)                                                                                                                           tCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtgc     >  1:216278/1‑139 (MQ=255)                                                                                                                              ccATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtgctgc  >  2:1531877/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                        CAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGC  >  NC_000913/3398747‑3399010

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