Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,545,816	C→A	S15* (TCG→TAG)	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,545,816	0	C	A	100.0%	92.5 / NA	28	S15* (TCG→TAG)	murQ	N‑acetylmuramic acid 6‑phosphate etherase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (12/16);  total (12/16)

AAGCCTGTCAATTCGGGATATATGATTCCATACATCTGCAATTATTTGATCGTAAATAGTAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTCGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGATAAAACCGTACCGCTTGCCGTTGAGCGCGTACTGCCGGATATCGCCGCGGCGATCGATGTTATCCACGCCCA  >  NC_000913/2545700‑2545945                                                                                                                     |                                                                                                                                  aaGCCTGTCAATTCGGGATATATGATTCCATACATCTGCAATTATTTGATCGTAAATAGTAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGa                                                                                                            <  2:404847/140‑1 (MQ=255)     cTGTCAATTCGGGATATATGATTCCATACATCTGCAATTATTTGATCGTAAATAGTAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCt                                                                                                        >  1:298557/1‑140 (MQ=255)         cAATTCGGGATATATGATTCCATACATCTGCAATTATTTGATCGTAAATAGTAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGaa                                                                                                    <  2:143464/140‑1 (MQ=255)               gggATATATGATTCCATACATCTGCAATTATTTGATCGTAAATAGTAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGc                                                                                              <  1:52831/140‑1 (MQ=255)                ggATATATGATTCCATACATCTGCAATTATTTGATCGTAAATAGTAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGcc                                                                                             <  1:407609/140‑1 (MQ=255)                         aTTCCATACATCTGCAATTATTTGATCGTAAATAGTAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCa                                                                                    >  2:156861/1‑140 (MQ=255)                         aTTCCATACATCTGCAATTATTTGATCGTAAATAGTAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCa                                                                                    <  2:390401/140‑1 (MQ=255)                             cATACATCTGCAATTATTTGATCGTAAATAGTAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGa                                                                                >  1:52883/1‑140 (MQ=255)                              aTACATCTGCAATTATTTGATCGTAAATAGTCAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTTGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGAt                                                                               >  2:337095/1‑140 (MQ=255)                              aTACATCTGCAATTATTTGATCGTAAATAGTAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGAt                                                                               <  2:266476/140‑1 (MQ=255)                              aTACATCTGCAATTATTTGATCGTAAATAGTAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGAt                                                                               <  1:427548/140‑1 (MQ=255)                              aTACATCTGCAATTATTTGATCGTAAATAGTAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGAt                                                                               <  1:72196/140‑1 (MQ=255)                               tACATCTGCAATTATTTGATCGTAAATAGTAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATg                                                                              >  1:352842/1‑140 (MQ=255)                                    cTGCAATTATTTGATCGTAAATAGTAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGAt                                                                         >  2:245516/1‑140 (MQ=255)                                             tttGATCGTAAATAGTAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGATAAAACCGTa                                                                <  2:52883/140‑1 (MQ=255)                                                            tAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGATAAAACCGTACCGCTTGCCGTTGAg                                                 >  2:139026/1‑140 (MQ=255)                                                            tAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGATAAAACCGTACCGCTTGCCGTTGAg                                                 >  1:143234/1‑140 (MQ=255)                                                                gTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGATAAAACCGTACCGCTTGCCGTTGAGcgcg                                             <  2:42848/140‑1 (MQ=255)                                                                gTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGATAAAACCGTACCGCTTGCCGTTGAGcgcg                                             <  1:170003/140‑1 (MQ=255)                                                                gTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGATAAAACCGTACCGCTTGCCGTTGAGcgcg                                             <  2:143234/140‑1 (MQ=255)                                                                    ccaccGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGATAAAACCGTACCGCTTGCCGTTGAGCGCGTACt                                         >  1:333787/1‑140 (MQ=255)                                                                          aTGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGATAAAACCGTACCGCTTGCCGTTGAGCGCGTACTGCCGGa                                   >  2:297697/1‑140 (MQ=255)                                                                                    aaaaGATGATTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGATAAAACCGTACCGCTTGCCGTTGAGCGCGTACTGCCGGATATCGCcgcg                         <  1:245516/140‑1 (MQ=255)                                                                                         atgatTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGATAAAACCGTACCGCTTGCCGTTGAGCGCGTACTGCCGGATATCGCCGCGGcgat                    <  2:418032/140‑1 (MQ=255)                                                                                         atgatTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGATAAAACCGTACCGCTTGCCGTTGAGCGCGTACTGCCGGATATCGCCGCGGcgat                    <  2:9571/140‑1 (MQ=255)                                                                                            atTACTGAAGGCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGATAAAACCGTACCGCTTGCCGTTGAGCGCGTACTGCCGGATATCGCCGCGGcgatcga                 >  2:54531/1‑140 (MQ=255)                                                                                                     ggCTCGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGATAAAACCGTACCGCTTGCCGTTGAGCGCGTACTGCCGGATATCGCCGCGGCGATCGATGTTATCCa        >  1:262640/1‑140 (MQ=255)                                                                                                         cGAACACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGATAAAACCGTACCGCTTGCCGTTGAGCGCGTACTGCCGGATATCGCCGCGGCGATCGATGTTATCCACGcc    <  1:312144/140‑1 (MQ=255)                                                                                                           aaCACCGCCTAGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGATAAAACCGTACCGCTTGCCGTTGAGCGCGTACTGCCGGATATCGCCGCGGCGATCGATGTTATCCACGCCCa  >  1:265700/1‑140 (MQ=255)                                                                                                                     |                                                                                                                                  AAGCCTGTCAATTCGGGATATATGATTCCATACATCTGCAATTATTTGATCGTAAATAGTAAGGTCACCACCGATGCAATTTGAAAAGATGATTACTGAAGGCTCGAACACCGCCTCGGCTGAAATTGACCGCGTATCGACGCTGGAAATGTGCCGGATTATCAACGATGAAGATAAAACCGTACCGCTTGCCGTTGAGCGCGTACTGCCGGATATCGCCGCGGCGATCGATGTTATCCACGCCCA  >  NC_000913/2545700‑2545945

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