Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,269,485	G→A	L24L (CTG→CTA)	aphA →	acid phosphatase/phosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,269,485	0	G	A	100.0%	41.3 / NA	13	L24L (CTG→CTA)	aphA	acid phosphatase/phosphotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/1);  total (12/1)

AATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTAC‑TTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGA  >  NC_000913/4269276‑4269729                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      aaTATAGTCATCATGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGACAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGATGACCTAGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGc                                                                                                                                                                                                              >  2:26675/1‑251 (MQ=255)                                 cGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTAGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTggg                                                                                                                                                                              >  1:75727/1‑251 (MQ=255)                                                                            gATGGTTTTTGCAACCAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTAGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGcc                                                                                                                                   >  1:117215/1‑251 (MQ=255)                                                                                                tcACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTAGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGat                                                                                                                >  1:218144/1‑250 (MQ=255)                                                                                                             tttCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACCCAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTAGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGaggtaacggg                                                                                                  >  1:25354/1‑242 (MQ=255)                                                                                                                     tACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTAGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTAC‑ttttt                                                                                          <  2:152886/250‑1 (MQ=255)                                                                                                                      aCTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTAGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTAC‑TTTTTTc                                                                                        >  1:1553/1‑251 (MQ=255)                                                                                                                           cTATTAAGGGAAAAATATGCGCAAGATCACACCGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTAGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATTTCGATGATACGGTACCTTTTTTcaat                                                                                     >  2:228336/1‑247 (MQ=255)                                                                                                                                               gcAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTAGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGACACGGTAC‑TTTTTTCCAGTCCGGGCTTCTGGCGCGGCaa                                                                >  2:52808/1‑250 (MQ=255)                                                                                                                                                 aaGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTAGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTAC‑TTTTTTCCCGTCCGGGCTTCTGGCGCGGCaaaa                                                              >  2:89352/1‑250 (MQ=255)                                                                                                                                                  aGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTACACAGTTCCGCTGTTGCCCTAGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATTGATGACACGGTAC‑TTTTTTCCAGTCCGGGCTTCTGGCGCGGCaaaaa                                                             >  2:64946/1‑250 (MQ=255)                                                                                                                                                                               ttATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTAGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTAC‑TTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCt                               >  1:202781/1‑251 (MQ=255)                                                                                                                                                                                                             gCCCTAGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGACACGGTAC‑TTTTTTTCAGTCCGGGCTTTTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATTCTGTGTTCTGGGAAAAAATGa  >  1:75542/1‑250 (MQ=255)                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      AATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTAC‑TTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGA  >  NC_000913/4269276‑4269729

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