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%	86.3 / NA	28	*478* (TAG→TAA)	dacB	D‑alanyl‑D‑alanine carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/13);  total (15/13)

GAGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCG‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTG  >  NC_000913/3330264‑3330518                                                                                                                                      |                                                                                                                           gagcgGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCaaa                                                                                                                       >  1:735485/1‑139 (MQ=255)     gggCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGa                                                                                                                   >  1:146480/1‑139 (MQ=255)     gggCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGa                                                                                                                   <  1:221814/139‑1 (MQ=255)                       gtgCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCg                                                                                                 >  2:222135/1‑139 (MQ=255)                              atCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTc                                                                                          <  2:269949/139‑1 (MQ=255)                                  ttCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAATCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGat                                                                                      >  1:670513/1‑139 (MQ=255)                                       gCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGGGCGATTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTg                                                                                 >  1:331251/1‑139 (MQ=255)                                                  gAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAATAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTgcgc                                                                      <  2:574644/139‑1 (MQ=255)                                                       tGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGc                                                                                                   >  2:468205/1‑105 (MQ=255)                                                       tGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGc                                                                                                   <  1:468205/105‑1 (MQ=255)                                                       tGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGt                                                                 >  1:53607/1‑139 (MQ=255)                                                        gCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGcac                                                                                                        <  2:93512/99‑1 (MQ=255)                                                        gCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGcac                                                                                                        >  1:93512/1‑99 (MQ=255)                                                                          cGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa                                              <  1:222135/139‑1 (MQ=255)                                                                                 cGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCg                                                                                                 <  2:425037/81‑1 (MQ=255)                                                                                 cGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCg                                                                                                 >  1:425037/1‑81 (MQ=255)                                                                                    tAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTAt                                    <  1:738852/139‑1 (MQ=255)                                                                                       tgcgtttttGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCg                                 >  2:691992/6‑139 (MQ=255)                                                                                         cgtttttGAAAGCCGTTTGTATAAAGATAGTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGAGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGc                               >  1:514929/4‑139 (MQ=255)                                                                                                            tataAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt            <  2:767872/139‑1 (MQ=255)                                                                                                                      tttATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTg  >  1:650578/1‑139 (MQ=255)                                                                                                                      tttATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTg  <  2:146480/139‑1 (MQ=255)                                                                                                                            aGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCg                                 >  1:539202/1‑102 (MQ=255)                                                                                                                            aGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCg                                 <  2:539202/102‑1 (MQ=255)                                                                                                                              aaCAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCg                                 <  2:100596/100‑1 (MQ=255)                                                                                                                              aaCAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCg                                 >  1:100596/1‑100 (MQ=255)                                                                                                                              aaCAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATG‑CGGATGCGGCGTAAACGCCTTATc                                   >  1:766549/1‑97 (MQ=255)                                                                                                                              aaCAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATG‑CGGATGCGGCGTAAACGCCTTATc                                   <  2:766549/97‑1 (MQ=255)                                                                                                                                      |                                                                                                                           GAGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCG‑TTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTG  >  NC_000913/3330264‑3330518

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