Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,329,743	G→T	D261Y (GAT→TAT)	dacB →	peptidoglycan DD‑endopeptidase DacB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,329,743	0	G	T	100.0%	88.0 / NA	31	D261Y (GAT→TAT)	dacB	peptidoglycan DD‑endopeptidase DacB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (13/18);  total (13/18)

TGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGG‑TTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGC  >  NC_000913/3329502‑3329983                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  tGCCCCAAAGCCTGGGGATATGGCTTTTATTCGCGTTGCACCTTATTACCCCGTTACGTTGTTCAGCCAGGTACCCCCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGACCTGGAGGTTGGGCAAGGCGCCTTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTa                                                                                                                                                                                                                                           <  2:304009/250‑1 (MQ=255)     ccAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGCTGGAGCCAGCTATGCCGGTGCACTTCTGAAATATGAGTTAAAAc                                                                                                                                                                                                                                       >  2:98035/1‑250 (MQ=255)        aaGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGc                                                                                                                                                                                                                                   <  1:285969/251‑1 (MQ=255)                   tatGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCt                                                                                                                                                                                                                        >  1:77441/1‑251 (MQ=255)                     tGGCTTTGATTCGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACTCTCCCCCGTGGTGCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGATCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCGGAAATATGAGTTAAAACAGGCGGGTATCACCTg                                                                                                                                                                                                                       <  2:133919/250‑1 (MQ=255)                            tataCGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCggagatccctcgacccgggc                                                                                                                                                                                                               >  2:274248/1‑233 (MQ=255)                                    tGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACAatgctg                                                                                                                                                                                                       >  1:269612/1‑251 (MQ=255)                                      gCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTgcg                                                                                                                                                                                                     >  1:311375/1‑251 (MQ=255)                                                                            cACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGG‑TTAACGAACCTGGAACGGTAGTTGcc                                                                                                                                                               >  1:19373/1‑250 (MQ=255)                                                                                 ccccccGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGG‑TTAACGAACCTGGAACGGTAGTTGCCAGTaaa                                                                                                                                                         <  1:245372/250‑1 (MQ=255)                                                                                  cccccGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCAACAGGGGCGGAACACTGCTGCGCCAGACTCAGGGTAAACGAAAC‑GGAACGGTATTTGcccgcaaa                                                                                                                                                         >  1:76165/1‑244 (MQ=255)                                                                                             cTGCCGAAGCGCAATACTGCGAACTGGATGTAGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATTAACTGGCGCGGAACACTGCTTCGCCAGACTCAGG‑TTAACGAACCTGGGACGGTAGTTGctagcaaccatctggccccc                                                                                                                                             >  1:271582/1‑232 (MQ=255)                                                                                              tGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGGCAGACTCAGG‑TTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGc                                                                                                                                            >  1:946/1‑251 (MQ=255)                                                                                               gcagAAGCGCAATACTGCGAACTGGATGGGGTGCCAGGCGTCCTGAACCGCTTTACTCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTGGGCTTTTGCCGTGCAGGATGGAGCCGGCTATGCCGTTGCAATTCTGGGATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGG‑TTAACGGACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCg                                                                                                                                             <  1:311490/246‑1 (MQ=255)                                                                                               gCCGAAGCTCAATACTGCGAACTGTAGGTGGTGCCAGGCGACCTGAACCGCTTTACGCGGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTGTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGTCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACGCAGG‑TGAACGAACCGGGAACGGGAGTTGCCAGTAAACAGGCGGcccc                                                                                                                                              <  2:311375/248‑1 (MQ=255)                                                                                                   aaGCGCAAGACGGTGAACTGGATGTGGTGCCAGGCGACGGGAGCCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGGGCCAGCGATGCCGGTGCAATTCTGAAATATGAGTTAAAACATGCGGGGAGCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGG‑TTAATGGACCTGGAACGGGAGTTGCCAGTAAACAGTCGGCCCCGCTGCAc                                                                                                                                       <  1:158705/251‑1 (MQ=255)                                                                                                        aaaTACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGG‑TTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCt                                                                                                                                  <  1:237861/250‑1 (MQ=255)                                                                                                                     tGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTaa                                                                                                                                                                                                                                          <  2:176777/135‑1 (MQ=255)                                                                                                                     tGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTaa                                                                                                                                                                                                                                          >  1:176777/1‑135 (MQ=255)                                                                                                                        agggggtgCCAGGCGACCGCAACCGCGTTACGCTGAGGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGGGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATAGGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGG‑TTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCt                                                                                                                   <  2:177588/246‑1 (MQ=255)                                                                                                                            ggtgCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGG‑TTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGaaaa                                                                                                              <  1:98035/251‑1 (MQ=255)                                                                                                                               gCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGG‑TTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCTCTGCACGATCTGCTTAAGATTATGCTGAAAAAg                                                                                                            <  1:245354/250‑1 (MQ=255)                                                                                                                                                    tgcGCTGACGGTATGCCTGCAACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACTCTGCTGCGCCAGACTCAGG‑TTAACGAACCTGGAACGGTAGTTGCCAGTAAACATGCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACACCATGATCGCCGAt                                                                                       <  2:251467/248‑1 (MQ=255)                                                                                                                                                               gATGCCTGCCACAACGTTCTGAGTCTCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGGCAGCTATGCCGGTGCAATGCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTTCTGCGCCAGACTCAGG‑TTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGTCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGTTACGGTTTTCCGc                                                                           <  2:58793/251‑1 (MQ=255)                                                                                                                                                                  gcctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGG‑TTAACGATCCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGGCGGACAACATGATCGCCGATACGGTTTTCCGCATg                                                                        <  1:260757/251‑1 (MQ=255)                                                                                                                                                                                                 ggCTTTTGCCGTGCAGGATGGAGCCAGCTATGCAGGTGCAATTCTGAAATCTGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCCCCAGACTCAGG‑GTAACGAGACTGGAACGGTAATTGACAGCAAACCGTCGGCCCCGCTGCACAATCTGCTTAAGATTATACTGAAAAAGTCGGACAACATGATAGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATTGGCCTggg                                          >  1:204577/1‑249 (MQ=255)                                                                                                                                                                                                    ttttGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAAACTCAGG‑TTCAAGAACCTGGAACGGTAGGTGCCAGGAAAAAGTCGGGCCCGCTGCCCGATCCGCTTAAGATTTTGCTGAAAAAATCGGACAACATGATCGCCGATAAGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAAcct                                      >  1:221525/1‑249 (MQ=255)                                                                                                                                                                                                               aGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAACCAGGCGGGGATAACCTGGAGCGGAACAATGATGCGCCAGAATCAGG‑GTAACGAACCTGGAACGGTGGATGCCAGTCAACCATCCGCCCCGCTGCACGATCTGGTTAAGATTATGCTGAAAAAGACGGACAACATGAACGCCCATACGGGTTTCCGCATGATGGGCCATACGGCCTTCAATGTGCCTGGGACATgggcggacggg                           >  2:223801/1‑242 (MQ=255)                                                                                                                                                                                                                     gAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGTTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGG‑TTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCTGAc                     <  2:76165/251‑1 (MQ=255)                                                                                                                                                                                                                                    gTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGG‑TTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATc      <  1:134073/251‑1 (MQ=255)                                                                                                                                                                                                                                        aaTTATGAAATAGGAGTTAAAACAGGGGGATAGCACCTGGAGCGCAACATGGCTGCGCTGGACTCAGG‑TTAACGAACCTGGAACGGTAGTTGCCAGTAAACATGCGGCCCCTCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACTTGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTgc  <  2:132757/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  TGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGG‑TTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGC  >  NC_000913/3329502‑3329983

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