Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	97,297	G→A	A71T (GCC→ACC)	murD →	UDP‑N‑acetylmuramoyl‑L‑alanine:D‑glutamate ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	97,297	0	G	A	100.0%	159.3 / NA	49	A71T (GCC→ACC)	murD	UDP‑N‑acetylmuramoyl‑L‑alanine:D‑glutamate ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/25);  total (24/25)

CCTCGCTCGCGGTGTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCGCCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCTGTCGCGAAGCACAAGCACCGATTGTGGCGATTACCGGTTCTA  >  NC_000913/97161‑97426                                                                                                                                         |                                                                                                                                  ccTCGCTCGAGGTGTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTCGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCAcc                                                                                                                                 <  1:867105/139‑1 (MQ=255)      ctcgcGGTGTGACGCCGCGCGTTATGGATAGGCGTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTcc                                                                                                                            <  1:53318/139‑1 (MQ=255)      ctcgcGGTGTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTcc                                                                                                                            <  1:533252/139‑1 (MQ=255)       tcgcGGTGTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTccc                                                                                                                           <  2:409635/139‑1 (MQ=255)          cgGTGTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGt                                                                                                                        >  2:350757/1‑139 (MQ=255)           ggTGTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTa                                                                                                                       <  2:391642/139‑1 (MQ=255)           ggTGTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTa                                                                                                                       <  1:350757/139‑1 (MQ=255)                     cgcgcgTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGgc                                                                                                             >  2:1394031/1‑139 (MQ=255)                       cgcgTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGt                                                                                                                              <  1:1083190/120‑1 (MQ=255)                       cgcgTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGt                                                                                                                              >  2:1083190/1‑120 (MQ=255)                       cgcgTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGt                                                                                                                              <  1:581488/120‑1 (MQ=255)                       cgcgTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGt                                                                                                                              >  2:581488/1‑120 (MQ=255)                       cgcgTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGgcgc                                                                                                           >  1:731688/1‑139 (MQ=255)                         cgTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGcat                                                                                                         <  2:420314/139‑1 (MQ=255)                          gTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGcatc                                                                                                        <  2:97902/139‑1 (MQ=255)                           ttATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGcatcc                                                                                                       <  1:1394031/139‑1 (MQ=255)                                      cgTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGCAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGcat                                                                                                         >  1:383582/1‑126 (MQ=255)                                      cgTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGCAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGcat                                                                                                         <  2:383582/126‑1 (MQ=255)                                        tATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCt                                                                                          >  2:1295598/1‑139 (MQ=255)                                               ccGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTg                                                                                         >  2:658080/1‑133 (MQ=255)                                               ccGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTg                                                                                         <  1:658080/133‑1 (MQ=255)                                                   ctggcctggATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGc                                                                               <  1:340726/139‑1 (MQ=255)                                                   ctggcctggATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGc                                                                               >  2:136236/1‑139 (MQ=255)                                                                      cGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGc                                                            >  2:173409/1‑139 (MQ=255)                                                                        aaGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGa                                                          >  2:715048/1‑139 (MQ=255)                                                                         aGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGca                                                                                                          >  1:1142115/1‑90 (MQ=255)                                                                         aGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGca                                                                                                          <  2:1142115/90‑1 (MQ=255)                                                                                          cGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCTGTcg                                        <  2:1150007/139‑1 (MQ=255)                                                                                           gggCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAgcg                                                                                            >  2:278996/1‑86 (MQ=255)                                                                                           gggCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAgcg                                                                                            <  1:278996/86‑1 (MQ=255)                                                                                                gTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCTGTCGCGaagc                                  >  2:267630/1‑139 (MQ=255)                                                                                                 tCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCTGTCGCGaagca                                 <  1:715048/139‑1 (MQ=255)                                                                                                  cTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTg                                                                                   >  1:1160591/1‑88 (MQ=255)                                                                                                  cTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTg                                                                                   <  2:1160591/88‑1 (MQ=255)                                                                                                     aaTGATGAATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCTGTCGCGaagcacaag                             >  2:751367/1‑139 (MQ=255)                                                                                                           gaATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCTGTCGCGaagca                                 >  1:561653/1‑129 (MQ=255)                                                                                                           gaATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCTGTCGCGaagca                                 <  2:561653/129‑1 (MQ=255)                                                                                                           gaATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCTGTCGCGa                                     <  1:656018/125‑1 (MQ=255)                                                                                                           gaATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCTGTCGCGa                                     >  2:656018/1‑125 (MQ=255)                                                                                                           gaATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCTGTCGCGa                                     <  1:1361595/125‑1 (MQ=255)                                                                                                           gaATGGCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCTGTCGCGa                                     >  2:1361595/1‑125 (MQ=255)                                                                                                                gCTGATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCTGTCGCGAAGCACAAGCACCGATtgtg                  >  1:579875/1‑139 (MQ=255)                                                                                                                   gATGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCTGTCGCGAAGCACAAGCACCGATTGTGGCg               >  2:1089410/1‑139 (MQ=255)                                                                                                                    aTGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCTGTCGCGAAGCACAAGCACCGATtgt                   >  1:294936/1‑134 (MQ=255)                                                                                                                    aTGGCGGCAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCTGTCGCGAAGCACAAGCACCGATtgt                   <  2:294936/134‑1 (MQ=255)                                                                                                                          gcAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTg                                                                                         >  1:519621/1‑58 (MQ=255)                                                                                                                          gcAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTg                                                                                         <  2:519621/58‑1 (MQ=255)                                                                                                                           cAGATCTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCACTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCGGTCGCGAAGCACAAGCACCGATTGTGGCGATTACCgg       >  2:918039/1‑139 (MQ=255)                                                                                                                                cTGATTGTCACCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCTGTCGCGAAGCACAAGCACCGATTGTGGCGATTACCGGTTCTa  <  2:579875/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  CCTCGCTCGCGGTGTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGAAGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCGCCAGTCCCGGTATTGCACTGGCGCATCCATCCTTAAGCGCTGCCGCTGATGCCGGAATCGAAATCGTTGGCGATATCGAGCTGTTCTGTCGCGAAGCACAAGCACCGATTGTGGCGATTACCGGTTCTA  >  NC_000913/97161‑97426

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