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%	165.0 / NA	52	P182L (CCG→CTG)	anmK	anhydro‑N‑acetylmuramic acid kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/27);  total (25/27)

TTTTGCAGCAGTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCGGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCCCAAGG  >  NC_000913/1718926‑1719187                                                                                                                                     |                                                                                                                                  ttttGCAGCAGTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGcc                                                                                                                             <  2:1222379/139‑1 (MQ=255) ttttGCAGCAGTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGcc                                                                                                                             <  1:330305/139‑1 (MQ=255)                             ccgcccgTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                                 <  2:1046829/139‑1 (MQ=255)                             ccgcccgTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                                 <  2:1020147/139‑1 (MQ=255)                               gcccgTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATgc                                                                                               >  1:672541/1‑139 (MQ=255)                                 ccgTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATgccg                                                                                             >  2:424905/1‑139 (MQ=255)                                    tGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCa                                                                                          >  2:511864/1‑139 (MQ=255)                                                     tATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATtcg                                                                         <  2:672541/139‑1 (MQ=255)                                                            aGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGcc                                                                                                                             >  2:79244/1‑80 (MQ=255)                                                            aGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGcc                                                                                                                             <  1:79244/80‑1 (MQ=255)                                                            aGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                   <  2:971195/90‑1 (MQ=255)                                                            aGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCa                                                                                          <  2:898427/115‑1 (MQ=255)                                                            aGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCa                                                                                          >  1:898427/1‑115 (MQ=255)                                                            aGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCg                                                                  >  1:208303/1‑139 (MQ=255)                                                            aGGATTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                   >  1:971195/1‑90 (MQ=255)                                                                  aCCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTggg                                                            >  1:1385495/1‑139 (MQ=255)                                                                   ccGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGt                                                           <  1:1031992/139‑1 (MQ=255)                                                                        cTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                                 >  1:250343/1‑96 (MQ=255)                                                                        cTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                                 <  2:250343/96‑1 (MQ=255)                                                                         tGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGcc                                                                                                                             >  2:90886/1‑67 (MQ=255)                                                                         tGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGcc                                                                                                                             <  1:90886/67‑1 (MQ=255)                                                                               cAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAgcg                                               <  1:511864/139‑1 (MQ=255)                                                                                       ggCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATgg                                       >  2:520769/1‑139 (MQ=255)                                                                                        gCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGa                                      <  2:208303/139‑1 (MQ=255)                                                                                            ccatcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGc                                  >  2:1198045/1‑139 (MQ=255)                                                                                            ccatcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGc                                  <  1:1198045/139‑1 (MQ=255)                                                                                                       ttACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCgg                                                                 <  2:109650/97‑1 (MQ=255)                                                                                                       ttACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCgg                                                                 >  1:109650/1‑97 (MQ=255)                                                                                                       ttACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCg                       >  2:55930/1‑139 (MQ=255)                                                                                                        tACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCgg                                                                 <  1:690687/96‑1 (MQ=255)                                                                                                        tACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCgg                                                                 >  2:690687/1‑96 (MQ=255)                                                                                                              gACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAatgtt                                                                                     <  2:934681/70‑4 (MQ=255)                                                                                                              gACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAatgtt                                                                                     >  1:934681/1‑67 (MQ=255)                                                                                                                ccGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCtt              >  2:1133024/1‑139 (MQ=255)                                                                                                                 cGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                                 <  1:1256790/55‑1 (MQ=255)                                                                                                                 cGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                                 >  2:1256790/1‑55 (MQ=255)                                                                                                                 cGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCt                                                                    >  1:969013/1‑84 (MQ=255)                                                                                                                 cGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCt                                                                    <  2:969013/84‑1 (MQ=255)                                                                                                                  ggTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGt            <  1:1133024/139‑1 (MQ=255)                                                                                                                   gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCcag                                                    >  2:667310/1‑98 (MQ=255)                                                                                                                   gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCcag                                                    <  1:667310/98‑1 (MQ=255)                                                                                                                     aTCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGtcct         <  2:1385495/139‑1 (MQ=255)                                                                                                                       cGTAGCCCCCACCCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGTAGGTACCAGCGGTGCGCCTTGtcctcc       <  1:926369/139‑1 (MQ=255)                                                                                                                        gTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTccc      <  1:55930/139‑1 (MQ=255)                                                                                                                         tAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCa                                 <  2:310024/111‑1 (MQ=255)                                                                                                                         tAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCa                                 >  1:310024/1‑111 (MQ=255)                                                                                                                         tAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCCCa     >  1:1257964/1‑139 (MQ=255)                                                                                                                            cccccAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCCCAAgg  >  1:412859/1‑139 (MQ=255)                                                                                                                              cccAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGcc                <  2:343040/123‑1 (MQ=255)                                                                                                                              cccAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGcc                >  1:343040/1‑123 (MQ=255)                                                                                                                                 aaCCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                                 >  1:1078561/1‑39 (MQ=255)                                                                                                                                 aaCCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                                 <  2:1078561/39‑1 (MQ=255)                                                                                                                                     |                                                                                                                                  TTTTGCAGCAGTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCGGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCCCAAGG  >  NC_000913/1718926‑1719187

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