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%	204.7 / NA	63	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 (33/30);  total (33/30)

AATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGC  >  NC_000913/3398740‑3399010                                                                                                                                        |                                                                                                                                        aaTACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTAt                                                                                                                                      >  2:872072/1‑139 (MQ=255)             tGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGgc                                                                                                                          >  1:488551/1‑139 (MQ=255)             tGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAAGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGgc                                                                                                                          >  1:1043430/1‑139 (MQ=255)                   acacCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCAcgc                                                                                                                    >  1:549349/1‑139 (MQ=255)                           tttGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCgg                                                                                                            <  2:1043430/139‑1 (MQ=255)                                  ggAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGca                                                                                                     <  1:1224511/139‑1 (MQ=255)                                                 tATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcggcggc                                                                                      <  1:872072/139‑1 (MQ=255)                                                   tCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGATCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTccc                                                                                                                             <  2:1689489/98‑1 (MQ=255)                                                   tCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTccc                                                                                                                             >  1:1689489/1‑98 (MQ=255)                                                    ccGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTa                                                                                   <  2:549349/139‑1 (MQ=255)                                                     cGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAg                                                                                  >  1:1377802/1‑139 (MQ=255)                                                     cGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAg                                                                                  >  2:773666/1‑139 (MQ=255)                                                      ggCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGc                                                                                 <  1:773666/139‑1 (MQ=255)                                                       gCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGcc                                                                                >  1:1517017/1‑139 (MQ=255)                                                         cAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCt                                                                              >  2:907121/1‑139 (MQ=255)                                                            ggCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTcc                                                                                                                              >  1:660706/1‑88 (MQ=255)                                                            ggCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTcc                                                                                                                              <  2:660706/88‑1 (MQ=255)                                                                 gATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGgag                                                                      >  1:152657/1‑139 (MQ=255)                                                                  aTTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGgaga                                                                     <  2:152657/139‑1 (MQ=255)                                                                  aTTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGgaga                                                                     <  2:1377802/139‑1 (MQ=255)                                                                     tGCAGCAACTGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCg                                                                  <  1:1285489/139‑1 (MQ=255)                                                                          caACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTag                                                             >  2:1668632/1‑139 (MQ=255)                                                                                cgcAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGccc                                                                               >  1:1287717/1‑115 (MQ=255)                                                                                cgcAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGccc                                                                               <  2:1287717/115‑1 (MQ=255)                                                                                    aTGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGagag                                                                    <  2:196416/122‑1 (MQ=255)                                                                                    aTGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGagag                                                                    >  1:196416/1‑122 (MQ=255)                                                                                        ggAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCACACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTc                                               >  2:913243/1‑139 (MQ=255)                                                                                             gagagCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGgcgc                                          <  1:907121/139‑1 (MQ=255)                                                                                              agagCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCACACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCgca                                         <  1:913243/139‑1 (MQ=255)                                                                                                 gCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCgcagca                                      >  1:80930/1‑139 (MQ=255)                                                                                                 gCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCgcagca                                      >  1:820879/1‑139 (MQ=255)                                                                                                   cAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGagag                                                                    <  1:748035/107‑1 (MQ=255)                                                                                                   cAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGagag                                                                    >  2:748035/1‑107 (MQ=255)                                                                                                        tACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcggcggc                                                                                      <  1:1518881/84‑1 (MQ=255)                                                                                                        tACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcggcggc                                                                                      >  2:1518881/1‑84 (MQ=255)                                                                                                          cccAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTa                             <  2:80930/139‑1 (MQ=255)                                                                                                           ccAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGgc                                            <  2:1167340/123‑1 (MQ=255)                                                                                                           ccAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCAgc                                            >  1:1167340/1‑123 (MQ=255)                                                                                                                    cTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCgg                                                                                                            >  1:141473/1‑50 (MQ=255)                                                                                                                    cTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCgg                                                                                                            <  2:141473/50‑1 (MQ=255)                                                                                                                    cTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCtgtg                                                 <  1:571467/109‑1 (MQ=255)                                                                                                                    cTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCtgtg                                                 >  2:571467/1‑109 (MQ=255)                                                                                                                       gCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTg                <  2:2841/139‑1 (MQ=255)                                                                                                                          aCGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCg             >  1:516191/1‑139 (MQ=255)                                                                                                                            gATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAgg                                                                                                               >  1:671615/1‑39 (MQ=255)                                                                                                                            gATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAgg                                                                                                               <  2:671615/39‑1 (MQ=38)                                                                                                                               aGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGt        <  2:199533/139‑1 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTagcag                                                          >  1:1690313/1‑87 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTagcag                                                          <  2:1690313/87‑1 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTg                >  1:141976/1‑129 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTg                <  2:141976/129‑1 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCgg         <  1:220759/136‑1 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCgg         >  2:220759/1‑136 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg      >  2:307186/1‑139 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg      >  1:907812/1‑139 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGC‑‑CGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAgg                                    <  2:571800/107‑1 (MQ=255)                                                                                                                                 cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGC‑‑CGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAgg                                    >  1:571800/1‑107 (MQ=255)                                                                                                                                  tCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCt                          <  2:1733769/118‑1 (MQ=255)                                                                                                                                  tCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCt                          >  1:1733769/1‑118 (MQ=255)                                                                                                                                     ccATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGcac                                                                                                    <  1:479067/41‑1 (MQ=255)                                                                                                                                     ccATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCTGCGTAGCCCTTTGTGGAGCGCGATTAGCAGCAGGCTGGGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtgctgc  >  2:1496771/1‑139 (MQ=255)                                                                                                                                     ccATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCtgtg                                                 <  2:668435/92‑1 (MQ=255)                                                                                                                                     ccATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCtgtg                                                 >  1:668435/1‑92 (MQ=255)                                                                                                                                        |                                                                                                                                        AATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGC  >  NC_000913/3398740‑3399010

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: ‑