Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,330,396	G→A	*478* (TAG→TAA)	dacB →	D‑alanyl‑D‑alanine carboxypeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,330,396	0	G	A	100.0%	168.4 / NA	52	*478* (TAG→TAA)	dacB	D‑alanyl‑D‑alanine carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/26);  total (26/26)

AGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCG‑‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTA  >  NC_000913/3330265‑3330531                                                                                                                                      |                                                                                                                                        agcgGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCaaaa                                                                                                                                   >  2:1434995/1‑139 (MQ=255)    gggCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGa                                                                                                                                <  1:1434995/139‑1 (MQ=255)                         cAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCgggg                                                                                                           <  1:720842/139‑1 (MQ=255)                                           gCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAAccc                                                                                                                           <  2:656326/105‑1 (MQ=255)                                           gCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAAccc                                                                                                                           >  1:656326/1‑105 (MQ=255)                                                  aaCCTGCGGGTCAAGTTAATTGCCGTATTCCTTTAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCt                                                                                  <  1:990059/139‑1 (MQ=255)                                                           aTCAGCGTAATCGCCGTATTCCGTTAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGca                                                                                                                      <  2:1269678/94‑1 (MQ=255)                                                           aTCAGCGTAATCGCCGTATTCCGTTAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGca                                                                                                                      >  1:1269678/1‑94 (MQ=255)                                                                    aTCGCCGTATTCCGTTAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATg                                                                                                                  >  2:20965/1‑89 (MQ=255)                                                                    aTCGCCGTATTCCGTTAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATg                                                                                                                  <  1:20965/89‑1 (MQ=255)                                                                     tCGCCGTATTCCGTTAGTGCGTTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                <  2:31379/139‑1 (MQ=39)                                                                          gTATTCCGTTAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                >  2:508859/1‑133 (MQ=255)                                                                          gTATTCCGTTAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                <  1:508859/133‑1 (MQ=255)                                                                          gTATTCCGTTAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaa                                                          >  1:472039/1‑139 (MQ=255)                                                                              tCCGTTAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGc                                                      <  2:1085735/139‑1 (MQ=255)                                                                                  ttAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACGTGGCCGGGGCTTCAGat                                                                                                   >  2:1016653/1‑90 (MQ=255)                                                                                  ttAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGcc                                                                                                               <  1:432262/78‑1 (MQ=255)                                                                                  ttAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGcc                                                                                                               >  2:432262/1‑78 (MQ=255)                                                                                  ttAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGat                                                                                                   <  1:1016653/90‑1 (MQ=255)                                                                                  ttAGTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTa                                                  >  2:455079/1‑139 (MQ=255)                                                                                     gTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGcc                                                                                                               >  2:1273803/1‑75 (MQ=255)                                                                                     gTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGcc                                                                                                               <  1:1273803/75‑1 (MQ=255)                                                                                     gTGCGT‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATcc                                               >  1:1599197/1‑139 (MQ=255)                                                                                       tgcgtttttGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATgc                                                                 >  2:1147951/6‑120 (MQ=255)                                                                                       tgcgtttttGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATgc                                                                 <  1:1147951/115‑1 (MQ=255)                                                                                           tttttGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATg                                                                  >  2:474767/2‑115 (MQ=255)                                                                                           tttttGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATg                                                                  <  1:474767/114‑1 (MQ=255)                                                                                             tttGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAAccc                                                                                                                           <  1:407427/56‑1 (MQ=255)                                                                                             tttGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAAccc                                                                                                                           >  2:407427/1‑56 (MQ=255)                                                                                                      cGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAg                                                                                                     >  1:522753/1‑69 (MQ=255)                                                                                                      cGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAg                                                                                                     <  2:522753/69‑1 (MQ=255)                                                                                                       gTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                              <  1:455079/139‑1 (MQ=255)                                                                                                                 aGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCg                                                                                                              <  1:1113429/49‑1 (MQ=255)                                                                                                                 aGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                              >  2:409951/1‑129 (MQ=255)                                                                                                                 aGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                              <  1:409951/129‑1 (MQ=255)                                                                                                                 aGATATTTATCAGAACAATAAATCAAAAAGAAACCCCGGCACATGGCCg                                                                                                              >  2:1113429/1‑49 (MQ=255)                                                                                                                      tttATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGat                                                                                                   <  1:1101051/55‑1 (MQ=255)                                                                                                                      tttATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGat                                                                                                   >  2:1101051/1‑55 (MQ=255)                                                                                                                      tttATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCg                                              <  1:908622/108‑1 (MQ=255)                                                                                                                      tttATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCg                                              >  2:908622/1‑108 (MQ=255)                                                                                                                      tttATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                         <  2:1178743/129‑1 (MQ=255)                                                                                                                      tttATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                         >  1:1178743/1‑129 (MQ=255)                                                                                                                      tttATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTg               >  2:18543/1‑139 (MQ=255)                                                                                                                      tttATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTg               >  2:89664/1‑139 (MQ=255)                                                                                                                           cAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATc                                                <  2:1667116/101‑1 (MQ=255)                                                                                                                           cAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATc                                                >  1:1667116/1‑101 (MQ=255)                                                                                                                            aGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtc         <  1:18543/139‑1 (MQ=255)                                                                                                                             gAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGt                                                                                       >  1:701099/1‑60 (MQ=255)                                                                                                                             gAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGt                                                                                       <  2:701099/60‑1 (MQ=255)                                                                                                                               aCAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCt                                                                                  <  2:1413545/63‑1 (MQ=255)                                                                                                                               aCAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCt                                                                                  >  1:1413545/1‑63 (MQ=255)                                                                                                                                   ttAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTa  >  1:1103672/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                        AGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCG‑‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTA  >  NC_000913/3330265‑3330531

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