Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,546,031	A→G	T87A (ACC→GCC)	murQ →	N‑acetylmuramic acid 6‑phosphate etherase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,546,031	0	A	G	100.0%	206.7 / NA	55	T87A (ACC→GCC)	murQ	N‑acetylmuramic acid 6‑phosphate etherase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/32);  total (23/32)

GTACTGCCGGATATCGCCGCGGCGATCGATGTTATCCACGCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCACCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCGGGAAGGCGGTGTTAATGATCTGAAAAATATTAATTTAACGGCAC  >  NC_000913/2545902‑2546166                                                                                                                                  |                                                                                                                                        gTACTGCCGGATATCGCCGCGGCGATCGATGTTATCCACGCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGt                                                                                                                               <  1:177846/140‑1 (MQ=255) gTACTGCCGGATATCGCCGCGGCGATCGATGTTATCCACGCCCAGGAAAGCGGCGGCGGGCGTCTGATTTCCCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAGTGTCCGCCCGCCTACGGCGt                                                                                                                               <  2:900737/140‑1 (MQ=255)  tACTGCCGGATATCGCCGCGGCGATCGATGTTATCCACGCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGGCCACGGCGTg                                                                                                                              >  2:512802/1‑140 (MQ=255)    cTGCCGGATATCGCCGCGGCGATCGATGTTATCCACGCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGaa                                                                                                                            <  1:433289/140‑1 (MQ=255)       ccGGATATCGCCGCGGCGATCGATGTTATCCGCGCCCAGGTAGGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACACCGGGTCGTCGGGGGATTCGGGAGGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAAcc                                                                                                                         <  2:803009/140‑1 (MQ=255)         ggATATCGCCGCGGCGATCGATGTTATCCACGCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCgg                                                                                                                       <  2:349009/140‑1 (MQ=255)            tatCGCCGCGGCGATCATTGTTATCCACGCCCGGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTc                                                                                                                    <  2:217771/140‑1 (MQ=255)              tCGCCGCGGCGTTTGATGTTATCCCCGCCCTGGATAGCGGTGGTGGGCATCTGATTTACCTCGGTGCGGGAGCATCCGGTCGTCTGGGGATTCTGGATGCCAGCGATTGTCCGCCCGCCTACGGCGTGAAACCGGGTCtg                                                                                                                  <  2:965899/140‑1 (MQ=255)              tCGCCGCGGCGATCGATGTTATCCACGCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCtg                                                                                                                  >  1:1253639/1‑140 (MQ=255)              tCGCCGCGGCGATCGATGTTATCCACGCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCtg                                                                                                                  >  1:678372/1‑140 (MQ=255)               cgccgcGGCGATCGATGTTATCCACGCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCtgg                                                                                                                 >  2:1174123/1‑140 (MQ=255)                 ccgcGGCGATCGATGTTATCCACGCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCtggtg                                                                                                               >  1:1440508/1‑140 (MQ=255)                 ccgcGGCGATCGATGTTATCCACGCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCtggtg                                                                                                               >  2:1418188/1‑140 (MQ=255)                   gcggcgTTCGATGTTATCCAGGCCCAGGTAGGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTGGTTTGGGGATTCTGGGTGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTggtggt                                                                                                             <  2:101130/140‑1 (MQ=255)                         atcgatGTTATCCACGCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGttt                                                                                                       >  1:541877/1‑140 (MQ=255)                              tGTTATCCACGCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTg                                                                                                  <  1:34960/140‑1 (MQ=255)                                  aTCCACGCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTgg                                                                                              >  1:914674/1‑140 (MQ=255)                                   tCCACGCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTggc                                                                                             <  1:134720/140‑1 (MQ=255)                                   tCCACGCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTggc                                                                                             >  2:1014003/1‑140 (MQ=255)                                   cccACGCCCGGGTAGGCGGCGGCGGGCGTCTGATTACCCTCGGTGCGGGACCATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTggc                                                                                             <  2:1204617/139‑1 (MQ=255)                                    ccccGCCGAGGTAAGCGGGGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCGGGTCGTCGGGGGATTCTGGATGCCAGCGAATGCCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTggcg                                                                                            <  2:695901/137‑1 (MQ=255)                                      ccGCCCGGGTAGGCGGCGGGGGGCGTCTGATTTCCCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTggcggc                                                                                          <  2:1108606/139‑1 (MQ=255)                                        gCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGgcggcgg                                                                                        >  1:423175/1‑139 (MQ=255)                                        gCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGa                                                                                        >  2:272957/1‑140 (MQ=255)                                          ccAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGGTTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGGGAAAGCGGGTCGGGTGGGTGGGTTGATTGACAGCGGCTAat                                                                                      >  2:1152675/1‑140 (MQ=255)                                          ccAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAat                                                                                      <  1:369302/140‑1 (MQ=255)                                             ggTAGGCGGCGGCGGGCGTCTGATTNNCCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGCCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATg                                                                                   <  2:74200/140‑1 (MQ=255)                                             ggTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATg                                                                                   <  1:63966/140‑1 (MQ=255)                                             ggTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCGGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATg                                                                                   <  2:619942/140‑1 (MQ=255)                                              gTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGc                                                                                  >  2:1254594/1‑140 (MQ=255)                                                  gcggcggcgGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCAtt                                                                              >  2:474019/1‑140 (MQ=255)                                                  gcggcggcgGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCAtt                                                                              >  2:350998/1‑140 (MQ=255)                                                      cggcggGCGTCTGATTTCCCTCGGTGCGGGAACATCAGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGc                                                                          <  1:1290433/140‑1 (MQ=255)                                                        gcgggcgTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCAc                                                                        >  1:891135/1‑140 (MQ=255)                                                        gcgggcgTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCAc                                                                        >  1:849232/1‑140 (MQ=255)                                                        gcgggcgTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCAc                                                                        >  2:902350/1‑140 (MQ=255)                                                         cgggcgTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCAcg                                                                       >  2:974297/1‑140 (MQ=255)                                                                      tACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGgcgc                                                          >  2:959133/1‑140 (MQ=255)                                                                       aCCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGcgcg                                                         <  1:1014003/140‑1 (MQ=255)                                                                        ccTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCgg                                                        >  1:995811/1‑140 (MQ=255)                                                                             gTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGAt                                                   >  2:738804/1‑140 (MQ=255)                                                                             gTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGCGGAAGGCGCGGAAGAt                                                   >  1:457294/1‑140 (MQ=255)                                                                               gCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAg                                                 <  1:1012799/140‑1 (MQ=255)                                                                                    aaCATCCGGCCGTCGGGGGATTCGGGATGCCCGGGATGGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCggg                                            <  2:245430/140‑1 (MQ=255)                                                                                      cATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCGGGaa                                          <  1:1022705/140‑1 (MQ=255)                                                                                        tCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCGGGAAgg                                        <  2:1307173/140‑1 (MQ=255)                                                                                        tCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCGGGAAgg                                        <  2:541877/140‑1 (MQ=255)                                                                                        tCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCGGGAAgg                                        <  2:914674/140‑1 (MQ=255)                                                                                        tCCGGTCGGCGGGGGATTCTGGCTGCCCGCCAATGCCCCCCCGCCTACGGCGTGAAACCGGGTCTGGTGCTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCGGGAAgg                                        <  1:512802/140‑1 (MQ=255)                                                                                        cccGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCGGGAAgg                                        <  1:207258/139‑1 (MQ=255)                                                                                         ccGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCGGGAAGGc                                       <  1:198733/140‑1 (MQ=255)                                                                                           ggTCGTGTGGGGATTCTGGATGCCAGCGAATGCCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCGGGAAGGCgg                                     <  2:1183389/140‑1 (MQ=255)                                                                                                tcTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCGGGAAGGCGGTGTTa                                <  1:1254594/140‑1 (MQ=255)                                                                                                tcTGGGGATTCTGGATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCGGGAAGGCGGTGTTa                                >  1:105071/1‑140 (MQ=255)                                                                                                           gggttGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCGGGAAGGCGGTGTTAATGATCTGaaa                     <  2:685444/136‑1 (MQ=255)                                                                                                           gggATGCCAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCGGGAAGGCGGTGTTAATGATCTGaaa                     <  1:155147/139‑1 (MQ=255)                                                                                                                 ccAGCGAATGTCCGCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCGGGAAGGCGGTGTTAATGATCTGAAAAATAtt               >  2:978026/1‑140 (MQ=255)                                                                                                                            ccgcccgcCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCGGGAAGGCGGTGTTAATGATCTGAAAAATATTAATTTAACGGc    <  2:849232/140‑1 (MQ=255)                                                                                                                              gCCCGCCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCGGGAAGGCGGTGTTAATGATCTGAAAAATATTAATTTAACGGcac  <  2:105071/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                        GTACTGCCGGATATCGCCGCGGCGATCGATGTTATCCACGCCCAGGTAAGCGGCGGCGGGCGTCTGATTTACCTCGGTGCGGGAACATCCGGTCGTCTGGGGATTCTGGATGCCAGCGAATGTCCGCCCACCTACGGCGTGAAACCGGGTCTGGTGGTTGGTTTGATTGCTGGCGGCGAATATGCCATTCAGCACGCGGTGGAAGGCGCGGAAGATAGCCGGGAAGGCGGTGTTAATGATCTGAAAAATATTAATTTAACGGCAC  >  NC_000913/2545902‑2546166

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