Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	93,668	A→G	E168G (GAG→GGG)	murE →	UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate:meso‑ diaminopimelate ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	93,668	0	A	G	100.0%	131.1 / NA	36	E168G (GAG→GGG)	murE	UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate:meso‑ diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/20);  total (16/20)

ACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCA  >  NC_000913/93526‑93809                                                                                                                                               |                                                                                                                                              aCCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCgg                                                                                                                                       >  2:494035/1‑151 (MQ=255)  ccACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCggg                                                                                                                                      <  2:719784/151‑1 (MQ=255)     cGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCt                                                                                                                                   <  2:725298/151‑1 (MQ=255)           cccAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGa                                                                                                                             <  1:386424/151‑1 (MQ=255)            ccAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCtggggg                                                                                                                              >  2:683284/1‑149 (MQ=255)                  gTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAgg                                                                                                                        <  2:178691/149‑1 (MQ=255)                       cgcAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGGGGATCAGGGcgcg                                                                                                                   >  1:151961/1‑149 (MQ=255)                            tGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGGGGATCAGGGGGCGACGtttt                                                                                                            >  1:447195/1‑151 (MQ=255)                                       gCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATgg                                                                                                   <  1:152053/149‑1 (MQ=255)                                       gCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATgg                                                                                                   >  2:712277/1‑149 (MQ=255)                                       gCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGaa                                                                                                 <  1:563053/151‑1 (MQ=255)                                           ggCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGtt                                                                                              <  1:736588/150‑1 (MQ=255)                                                   cAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCa                                                                                      <  2:597414/150‑1 (MQ=255)                                                       gcgGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACggg                                                                                  >  2:507925/1‑150 (MQ=255)                                                        cgGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCt                                                                                <  1:637152/151‑1 (MQ=255)                                                         ggTAATGGGCACCGTTGGTAGCGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCt                                                                                <  2:679091/150‑1 (MQ=255)                                                          gTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGGGGGTCAGGGCGCGACGGTTTGCGCCATGGAAGTTTCCTCCCCCGGGCTgg                                                                              >  1:599690/1‑151 (MQ=255)                                                            aaTGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTa                                                                            <  2:677092/151‑1 (MQ=255)                                                                              cGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGGACAGCACCGTgggg                                                               >  1:551578/1‑146 (MQ=255)                                                                              cGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGGGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACCGCACCGTGTggcgg                                                            >  2:313100/1‑149 (MQ=255)                                                                                    gctgGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGGGGATCAGGGCGCGGAGTTTTGCGCAATTGAAGTTTCCTCCCCCGGGCTGGTACCGCACCGTGTGGCGGCCTTgga                                                     >  1:724455/1‑148 (MQ=255)                                                                                           aaaGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGGGGATCAGGGCGCGACGTTTTGCGCAATGGAAGGTTCCTCCCACGGGCTGGGACAGCACCGTGTGGCGGCATTGAAATTTgcgg                                             >  1:336798/1‑151 (MQ=255)                                                                                                    ccGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCt                                    <  1:386895/151‑1 (MQ=255)                                                                                                    ccGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGGACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTc                                     >  2:321486/1‑150 (MQ=255)                                                                                                        cAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCttt                                  >  2:508402/1‑149 (MQ=255)                                                                                                           aaaaTACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCaa                             <  1:546111/151‑1 (MQ=255)                                                                                                           aaaaTACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCaa                             <  1:101960/151‑1 (MQ=255)                                                                                                           aaaaTACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCa                              >  2:13253/1‑150 (MQ=255)                                                                                                               tACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACtt                          >  2:417546/1‑150 (MQ=255)                                                                                                                 cAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAg                       <  1:521232/151‑1 (MQ=255)                                                                                                                     cGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCcgc                   <  2:151994/151‑1 (MQ=255)                                                                                                                          cGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATc               <  2:737359/150‑1 (MQ=255)                                                                                                                          cGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATc               >  2:752980/1‑150 (MQ=255)                                                                                                                               gTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCGGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTg         <  1:447777/151‑1 (MQ=255)                                                                                                                                       tCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCa  <  1:582168/150‑1 (MQ=255)                                                                                                                                       tCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCa  <  1:412672/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                              ACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCA  >  NC_000913/93526‑93809

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