Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,552,889	G→A	D180N (GAC→AAC)	amiA →	N‑acetylmuramoyl‑L‑alanine amidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,552,889	0	G	A	100.0%	142.9 / NA	39	D180N (GAC→AAC)	amiA	N‑acetylmuramoyl‑L‑alanine amidase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/17);  total (22/17)

GCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACGCTCGG  >  NC_000913/2552814‑2552965                                                                            |                                                                             gCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAg                                                                         >  2:450128/1‑81 (MQ=255)  cAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAg                                                                         >  2:450133/1‑80 (MQ=255)   aaGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGTAGTTGCCGGTAAAAAGGCGACTAACGAGGa                                                                       <  2:450085/81‑1 (MQ=255)    aGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGAt                                                                      <  2:450086/81‑1 (MQ=255)          gcAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTa                                                                >  2:450120/1‑81 (MQ=255)           cAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTAt                                                               <  1:450087/81‑1 (MQ=255)             aTGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTg                                                             >  1:450126/1‑81 (MQ=255)               ggCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGc                                                            >  2:450123/1‑80 (MQ=255)               ggCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGc                                                            >  1:450118/1‑80 (MQ=255)                gCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGcaa                                                          >  2:450131/1‑81 (MQ=255)                gCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGcaa                                                          <  2:450088/81‑1 (MQ=255)                            tCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGtt                                               <  2:450089/80‑1 (MQ=255)                               gAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGAt                                           >  2:450143/1‑81 (MQ=255)                                 aCGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCt                                         >  1:450124/1‑81 (MQ=255)                                 aCGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCt                                         >  1:450115/1‑81 (MQ=255)                                     gAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTg                                     <  1:450090/81‑1 (MQ=255)                                       aaaCCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCa                                   <  2:450091/81‑1 (MQ=255)                                        aaCCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCa                                   >  2:450125/1‑80 (MQ=255)                                        aaCCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCa                                   >  2:450121/1‑80 (MQ=255)                                          ccGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAAc                                <  1:450093/81‑1 (MQ=255)                                          ccGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAAc                                <  1:450092/81‑1 (MQ=255)                                               ccGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATa                           <  2:450094/81‑1 (MQ=255)                                                cGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATAc                          >  2:450122/1‑81 (MQ=255)                                                cGATGAAGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATAc                          <  1:450095/81‑1 (MQ=255)                                                     aaGTTGCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCAtt                      >  1:450135/1‑80 (MQ=255)                                                         tGCCGGTAAAAAGGCCACTAACAATGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTaaaaa                 >  1:450163/1‑81 (MQ=255)                                                          gCCGGTAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTaaaaa                 >  1:450127/1‑80 (MQ=255)                                                            cGGTAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAg              >  2:450130/1‑81 (MQ=255)                                                              gTAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTc            <  1:450096/81‑1 (MQ=255)                                                               tAAAAAGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCt           <  1:450097/81‑1 (MQ=255)                                                                 aaaaGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTg          <  2:450098/80‑1 (MQ=255)                                                                 aaaaGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGa         >  2:450138/1‑81 (MQ=255)                                                                  aaaGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGAc        <  1:450099/81‑1 (MQ=255)                                                                   aaGGCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGAc        >  1:450179/1‑80 (MQ=255)                                                                      gCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACgctc    <  2:450100/81‑1 (MQ=255)                                                                      gCGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACgctc    <  1:450101/81‑1 (MQ=255)                                                                       cGACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACgctcg   >  1:450153/1‑81 (MQ=255)                                                                        gACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACgctcgg  >  2:450132/1‑81 (MQ=255)                                                                        gACTAACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACgctcgg  >  1:450134/1‑81 (MQ=255)                                                                            |                                                                             GCAAGTAGCGCAATGGCGAAATACCTGTCTGAACGCGAAAACCGCGCCGATGAAGTTGCCGGTAAAAAGGCGACTGACAAGGATCACCTATTGCAACAAGTGCTGTTTGATCTGGTGCAAACAGATACCATTAAAAATAGTCTGACGCTCGG  >  NC_000913/2552814‑2552965

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