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

CCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCGGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCC  >  NC_000913/1718957‑1719182                                                                                                      |                                                                                                                             cccgTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAATCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATgcc                                                                                         >  1:1028549/1‑139 (MQ=255)   cgTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATgccgc                                                                                       <  1:813041/139‑1 (MQ=255)     tGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGa                                                                                                        >  1:1251232/1‑120 (MQ=255)     tGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGa                                                                                                        <  2:1251232/120‑1 (MQ=255)                      tATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGaa                                                                                                                  >  1:1448829/1‑93 (MQ=255)                      tATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGaa                                                                                                                  <  2:1448829/93‑1 (MQ=255)                                         cTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGcc                                                 >  2:698164/1‑139 (MQ=255)                                          tGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCAc                                                                          <  2:502962/113‑1 (MQ=255)                                          tGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCAc                                                                          >  1:502962/1‑113 (MQ=255)                                          tGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGcc                                                 <  1:1032245/138‑1 (MQ=255)                                          tGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGcc                                                 >  2:1032245/1‑138 (MQ=255)                                            aCGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTg                                                                                                         >  1:1230424/1‑80 (MQ=255)                                            aCGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTg                                                                                                         <  2:1230424/80‑1 (MQ=255)                                                      cAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCatgat                                    >  2:658499/1‑139 (MQ=255)                                                            tccatcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACg                              >  2:746724/1‑139 (MQ=255)                                                              catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCaa                                                                                    <  1:1559525/83‑1 (MQ=255)                                                              catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCaa                                                                                    >  2:1559525/1‑83 (MQ=255)                                                              catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATtcgtcg                                                                 >  1:316457/1‑102 (MQ=255)                                                              catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATtcgtcg                                                                 <  2:316457/102‑1 (MQ=255)                                                              catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCa                            >  2:547060/1‑139 (MQ=255)                                                              catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCa                            >  2:111094/1‑139 (MQ=255)                                                                  aGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTa                        >  2:399767/1‑139 (MQ=255)                                                                   gCATGTTATCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTAc                       <  2:1028549/139‑1 (MQ=255)                                                                    cATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAg                                                                                                           <  1:177895/54‑1 (MQ=255)                                                                    cATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAg                                                                                                           >  2:177895/1‑54 (MQ=255)                                                                    cATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTAcc                      <  1:111094/139‑1 (MQ=255)                                                                         tACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCgg                 >  2:1517643/1‑139 (MQ=255)                                                                          aCCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGt                >  1:1330884/1‑139 (MQ=255)                                                                            cAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTgc              <  2:1330884/139‑1 (MQ=255)                                                                                aCCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                            >  1:1654593/1‑57 (MQ=255)                                                                                aCCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                            <  2:1654593/57‑1 (MQ=255)                                                                                  cGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGa                                                                                                               <  2:218532/36‑1 (MQ=255)                                                                                  cGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGa                                                                                                               >  1:218532/1‑36 (MQ=255)                                                                                     tATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTggg                                                       >  2:141813/1‑89 (MQ=255)                                                                                     tATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTggg                                                       <  1:141813/89‑1 (MQ=255)                                                                                        cGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGtcctcc  <  1:658499/139‑1 (MQ=255)                                                                                             cccccAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                            <  2:200436/44‑1 (MQ=255)                                                                                             cccccAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                            >  1:200436/1‑44 (MQ=255)                                                                                                  aaCCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTg                                                         <  2:533633/74‑1 (MQ=255)                                                                                                  aaCCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTg                                                         >  1:533633/1‑74 (MQ=255)                                                                                                      |                                                                                                                             CCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCGGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCC  >  NC_000913/1718957‑1719182

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