Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	453,291	G→A	L85L (CTC→CTT)	ampG ←	muropeptide transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	453,291	0	G	A	100.0%	130.4 / NA	38	L85L (CTC→CTT)	ampG	muropeptide transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/18);  total (20/18)

AAGACTATATCCTGGGAAGCAGAGCAAAAAGCGATCACCACTGCCAGCGCCGCCATCCAGCGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGA  >  NC_000913/453163‑453424                                                                                                                                 |                                                                                                                                      aaGACTATATCCTGGGAAGCAGAGCAAAAAGCGATCACCACTGCCAGCGCCGCCATCCAGCGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGcgccgccg                                                                                                                   >  3:114976/1‑149 (MQ=255)                      gagCAAAAAGCGATCACCACTGCCAGCGCCGCCATCCAGCGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAag                                                                                                                                     >  3:84645/1‑110 (MQ=255)                      gagCAAAAAGCGATCACCACTGCCAGCGCCGCCATCCAGCGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAag                                                                                                                                     <  4:84645/110‑1 (MQ=255)                      gagCAAAAAGCGATCACCACTGCCAGCGCCGCCATCCAGCGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCgtg                                                                                                >  2:33934/1‑147 (MQ=255)                       agCAAAAAGCGATCACCACTGCCAGCGCCGCCATCCAGCGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAg                                                                                             >  8:80066/1‑149 (MQ=255)                            aaaGCGATCACCACTGCCAGCGCCGCCATCCAGCGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTc                                                                                        >  4:3792/1‑149 (MQ=255)                             aaGCGATCACCACTGCCAGCGCCGCCATCCAGCGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTc                                                                                        >  5:86777/1‑148 (MQ=255)                               gCGATCACCACTGCCAGCGCCGCCATCCAGCGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCAt                                                                                     <  1:33934/149‑1 (MQ=255)                                        acTGCCAGCGCCGCCATCCAGCGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGa                                                                            >  7:112656/1‑149 (MQ=255)                                                    gccATCCAGCGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAg                                                                                                     >  1:101897/1‑112 (MQ=255)                                                    gccATCCAGCGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAg                                                                                                     <  2:101897/112‑1 (MQ=255)                                                     ccATCCAGCGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAAttt                                                                >  1:78667/1‑148 (MQ=255)                                                     ccATCCAGCGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTa                                                               <  4:114976/149‑1 (MQ=255)                                                        tCCAGCGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTaaaa                                                            <  2:78667/149‑1 (MQ=255)                                                           aGCGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACg                                                         >  5:49917/1‑149 (MQ=255)                                                                  ttGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCt                                                  <  6:42454/149‑1 (MQ=255)                                                                    gggTGCATGGTTCGAGAAAACCAATAGCCCCAATGCCGCCTAATAACAGGATTTGCGTGGCATCCTGCCACCCGCGCCGCCGCCCAAAAAATGGAGGGGTGAAGCGGTCAAGCAGCGGTTCCCAGAGGATTTTAAAAACGTAAGCCTgg                                                >  2:81704/1‑149 (MQ=255)                                                                              ttCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCa                                         >  1:32492/1‑146 (MQ=255)                                                                                  agaAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCag                                                                        >  4:94821/1‑111 (MQ=255)                                                                                  agaAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCag                                                                        <  3:94821/111‑1 (MQ=255)                                                                                   gaAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTaaaa                                                            <  1:112091/122‑1 (MQ=255)                                                                                   gaAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTaaaa                                                            >  2:112091/1‑122 (MQ=255)                                                                                     aaaCCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGagaaga                                >  1:78116/1‑148 (MQ=255)                                                                                       aCCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAAc                             <  2:78116/149‑1 (MQ=255)                                                                                        cccATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTg                                                 >  3:90916/1‑128 (MQ=255)                                                                                        cccATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTg                                                 <  4:90916/128‑1 (MQ=255)                                                                                        cccATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAAcc                            <  3:3792/149‑1 (MQ=255)                                                                                                 gcAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGtttt                   <  7:80066/149‑1 (MQ=255)                                                                                                 gcAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGACGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGtttt                   <  6:86777/149‑1 (MQ=255)                                                                                                  cAATGGCGACTTATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTg                  <  8:112656/149‑1 (MQ=255)                                                                                                  cAATGGCGACTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTg                  <  6:49917/149‑1 (MQ=255)                                                                                                           cTAATAACAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGATGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGata         <  1:81704/149‑1 (MQ=255)                                                                                                                aaCAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATtctc    >  8:85020/1‑149 (MQ=255)                                                                                                                  cAGGATTTGCGTGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGa  >  1:14822/1‑149 (MQ=255)                                                                                                                             tGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGtttt                   >  6:31863/1‑121 (MQ=255)                                                                                                                             tGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGtttt                   <  5:31863/121‑1 (MQ=255)                                                                                                                             tGGCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGAt              >  3:65073/1‑126 (MQ=255)                                                                                                                              ggCAAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAAttt                                                                <  4:110929/75‑1 (MQ=255)                                                                                                                                 |                                                                                                                                      AAGACTATATCCTGGGAAGCAGAGCAAAAAGCGATCACCACTGCCAGCGCCGCCATCCAGCGGAGTTGGGTGCCTGGTTCGAGAAAACCCATCGCCGCAATGGCGACTAATAACAGGATTTGCGTGGCGAGCAGCCAACCGCGCCGCCGCCCAAAAAATGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGA  >  NC_000913/453163‑453424

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