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%	122.1 / NA	39	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 (22/17);  total (22/17)

ACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTG  >  NC_000913/3398746‑3399006                                                                                                                                  |                                                                                                                                    acaaGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGGGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATCATTGCCCt                                                                                                                            >  2:1109448/1‑139 (MQ=255)  caaGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTAt                                                                                                                                  >  1:690305/1‑132 (MQ=255)  caaGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTAt                                                                                                                                  <  2:690305/132‑1 (MQ=255)      aTGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCgg                                                                                                                       <  1:79369/139‑1 (MQ=255)         aGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGgcgc                                                                                                                    >  2:168664/1‑139 (MQ=255)            aaCACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTAGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCAcg                                                                                                                 >  1:237846/1‑139 (MQ=255)               acCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACgcgc                                                                                                              >  1:433422/1‑139 (MQ=255)                cccAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCa                                                                                                             >  2:537223/1‑139 (MQ=255)                  cAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAgg                                                                                                           <  1:537223/139‑1 (MQ=255)                  cAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTAGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAgg                                                                                                           <  2:237846/139‑1 (MQ=255)                   aGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGc                                                                                                          >  2:619889/1‑139 (MQ=255)                   aGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGc                                                                                                          >  2:1073296/1‑139 (MQ=255)                       tGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGtt                                                                                                      >  1:1014530/1‑139 (MQ=255)                             gAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGcac                                                                                                <  2:718549/139‑1 (MQ=255)                                        ggTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcggc                                                                                     >  2:1187771/1‑139 (MQ=255)                                              ccGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTa                                                                               <  1:168664/139‑1 (MQ=255)                                                 gCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGcc                                                                                                                              <  1:787009/89‑1 (MQ=255)                                                 gCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGcc                                                                                                                              >  2:787009/1‑89 (MQ=255)                                                  ccAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCa                                                                                                                   <  1:739290/99‑1 (MQ=255)                                                  ccAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCa                                                                                                                   >  2:739290/1‑99 (MQ=255)                                                          tGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGga                                                                   <  2:1014530/139‑1 (MQ=255)                                                           gATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGgag                                                                  >  2:685044/1‑139 (MQ=255)                                                                gcagcaACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGa                                                             >  1:478299/1‑139 (MQ=255)                                                                                           gCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTtgt                                                                      <  1:722936/103‑1 (MQ=255)                                                                                           gCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTtgt                                                                      >  2:722936/1‑103 (MQ=255)                                                                                           gCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCgcagca                                  <  2:478299/139‑1 (MQ=255)                                                                                                          gTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGca                                                                                                 <  2:787650/61‑1 (MQ=255)                                                                                                          gTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGca                                                                                                 >  1:787650/1‑61 (MQ=255)                                                                                                                  cTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTagcag                                                      >  2:499730/1‑96 (MQ=255)                                                                                                                  cTACGATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTagcag                                                      <  1:499730/96‑1 (MQ=255)                                                                                                                      gATAGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGcc       >  2:1291427/1‑139 (MQ=255)                                                                                                                         aGCTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGt    <  1:1187771/139‑1 (MQ=255)                                                                                                                           cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTa                                                                               >  2:604315/1‑62 (MQ=255)                                                                                                                           cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTa                                                                               <  1:604315/62‑1 (MQ=255)                                                                                                                           cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAgg                                                  <  2:552109/91‑1 (MQ=255)                                                                                                                           cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAgg                                                  >  1:552109/1‑91 (MQ=255)                                                                                                                           cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCt             <  2:3361/128‑1 (MQ=255)                                                                                                                           cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCt             >  1:3361/1‑128 (MQ=255)                                                                                                                           cTCGCCATTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtg  >  1:73908/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                    ACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTG  >  NC_000913/3398746‑3399006

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