Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,719,058	G→A	P182L (CCG→CTG)	anmK ←	anhydro‑N‑acetylmuramic acid kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,719,058	0	G	A	100.0%	210.0 / NA	65	P182L (CCG→CTG)	anmK	anhydro‑N‑acetylmuramic acid kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/31);  total (34/31)

GTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCGGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCCCA  >  NC_000913/1718936‑1719184                                                                                                                           |                                                                                                                               gTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                <  1:193198/139‑1 (MQ=255)                tACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                >  2:1510666/1‑124 (MQ=255)                tACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                <  1:1510666/124‑1 (MQ=255)                       ccgTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAAGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATgccg                                                                                          >  2:1080294/1‑139 (MQ=255)                          tGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATg                                                                                                                  <  1:1117650/112‑1 (MQ=255)                          tGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATg                                                                                                                  >  2:1117650/1‑112 (MQ=255)                          tGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCa                                                                                       >  2:509843/1‑139 (MQ=255)                          tGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCa                                                                                       >  2:138975/1‑139 (MQ=255)                                   ggCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGc                                                                              >  1:1507928/1‑139 (MQ=255)                                      aTCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCa                                                                                       <  2:908000/127‑1 (MQ=255)                                      aTCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCa                                                                                       >  1:908000/1‑127 (MQ=255)                                         tttATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                <  1:900469/99‑1 (MQ=255)                                         tttATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                >  2:900469/1‑99 (MQ=255)                                              cGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATtcgtcg                                                                   <  1:1080294/139‑1 (MQ=255)                                               gTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATgc                                                                                            <  1:868150/113‑1 (MQ=255)                                               gTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATgc                                                                                            >  2:868150/1‑113 (MQ=255)                                                  aGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATgccgc                                                                                         >  2:351682/1‑113 (MQ=255)                                                  aGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATgccgc                                                                                         <  1:351682/113‑1 (MQ=255)                                                  aGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCg                                                               >  1:1390737/1‑139 (MQ=255)                                                   ggTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCgg                                                              >  2:432234/1‑139 (MQ=255)                                                       tACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTgg                                                          >  1:965094/1‑139 (MQ=255)                                                        aCCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGcc                                                                                                                          <  2:1308678/74‑1 (MQ=255)                                                        aCCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGcc                                                                                                                          >  1:1308678/1‑74 (MQ=255)                                                          cGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCa                                                                                       <  1:359662/107‑1 (MQ=255)                                                          cGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCa                                                                                       >  2:359662/1‑107 (MQ=255)                                                               tGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                <  2:1134430/77‑1 (MQ=255)                                                               tGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                >  1:1134430/1‑77 (MQ=255)                                                               tGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCca                                                  >  1:779886/1‑139 (MQ=255)                                                                      aGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAgcgc                                           >  2:54614/1‑139 (MQ=255)                                                                       gATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCa                                          <  1:509843/139‑1 (MQ=255)                                                                             ggCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATgg                                    >  2:1455407/1‑139 (MQ=255)                                                                                 tccatcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACg                                <  1:1455407/139‑1 (MQ=255)                                                                                  ccatcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAAt                                                                                                                   <  2:302434/55‑1 (MQ=255)                                                                                  ccatcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAAt                                                                                                                   >  1:302434/1‑55 (MQ=255)                                                                                  ccatcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCa                                          >  1:239455/1‑128 (MQ=255)                                                                                  ccatcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCa                                          <  2:239455/128‑1 (MQ=255)                                                                                  ccatcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGc                               <  1:432234/139‑1 (MQ=255)                                                                                   catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                >  2:114496/1‑57 (MQ=255)                                                                                   catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                <  1:114496/57‑1 (MQ=255)                                                                                   catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATgc                                                                                            <  2:1199276/77‑1 (MQ=255)                                                                                   catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATgc                                                                                            >  1:1199276/1‑77 (MQ=255)                                                                                   catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAgcg                                            >  1:1477413/1‑125 (MQ=255)                                                                                   catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAgcg                                            <  2:1477413/125‑1 (MQ=255)                                                                                   catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGa                                   <  2:968220/134‑1 (MQ=255)                                                                                   catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGa                                   >  1:968220/1‑134 (MQ=255)                                                                                   catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCa                              <  2:779886/139‑1 (MQ=255)                                                                                             ttACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCg                    <  2:965094/139‑1 (MQ=255)                                                                                                       cGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                              <  1:748465/55‑1 (MQ=255)                                                                                                       cGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                              >  2:748465/1‑55 (MQ=255)                                                                                                        ggTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGt         <  2:1390737/139‑1 (MQ=255)                                                                                                        ggTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGt         <  2:839114/139‑1 (MQ=255)                                                                                                         gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTg                                                                                                   >  2:267880/1‑48 (MQ=255)                                                                                                         gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTg                                                                                                   <  1:267880/48‑1 (MQ=255)                                                                                                         gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGtc        >  1:666512/1‑139 (MQ=255)                                                                                                         gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGtc        >  2:575205/1‑139 (MQ=255)                                                                                                         gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGA‑TGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCgg                                                              <  1:1427727/84‑1 (MQ=255)                                                                                                         gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGA‑TGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCgg                                                              >  2:1427727/1‑84 (MQ=255)                                                                                                          tATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                              <  2:133857/52‑1 (MQ=255)                                                                                                          tATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                              >  1:133857/1‑52 (MQ=255)                                                                                                             cGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGtcctcc    <  1:54614/139‑1 (MQ=255)                                                                                                               tAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTAcc                        >  2:1231078/1‑117 (MQ=255)                                                                                                               tAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTAcc                        <  1:1231078/117‑1 (MQ=255)                                                                                                               tAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCCCa  >  1:822086/1‑139 (MQ=255)                                                                                                                        aCCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGa                                                      <  1:31492/78‑1 (MQ=255)                                                                                                                        aCCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGa                                                      >  2:31492/1‑78 (MQ=255)                                                                                                                           |                                                                                                                               GTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCGGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCCCA  >  NC_000913/1718936‑1719184

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