Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,678,723	A→C	G48G (GGT→GGG)	glyA ←	serine hydroxymethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,678,723	0	A	C	100.0%	49.0 / NA	16	G48G (GGT→GGG)	glyA	serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (7/9);  total (7/9)

AGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTT  >  NZ_CP009273/2678637‑2678799                                                                                       |                                                                             aGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGAATAACCTTCAGCATATTTGTTGGTCAGCTGAGAcccc                                                                           >  2:471867/1‑90 (MQ=255)             atatCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCAtt                                                               >  2:214396/1‑90 (MQ=255)                   aCATACCCCCAACCGCCGTAGTAGCGTTTGCCCGAATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTAcgcgc                                                         <  1:471867/90‑1 (MQ=255)                     aTACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCgg                                                       >  1:189256/1‑90 (MQ=255)                       aCTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGc                                                     >  2:389676/1‑90 (MQ=255)                                            gTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGa                                >  1:189815/1‑90 (MQ=255)                                                     gATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGa                       <  1:101445/90‑1 (MQ=255)                                                      aTAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGAt                      <  1:438686/90‑1 (MQ=255)                                                       tAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGAtg                     <  2:487896/90‑1 (MQ=255)                                                              cAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCgg                                                       <  2:5440/49‑1 (MQ=255)                                                              cAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCgg                                                       >  1:5440/1‑49 (MQ=255)                                                                     tttGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACg       <  2:189815/90‑1 (MQ=255)                                                                     tttGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACg       <  1:46774/90‑1 (MQ=255)                                                                      ttgttgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGt      <  1:476445/90‑1 (MQ=255)                                                                          tgGTCAGCTGAGACCCCTGCGCCTGCATTAcgcgc                                                         >  1:48612/1‑35 (MQ=38)                                                                          tgGTCAGCTGAGACCCCTGCGCCTGCATTAcgcgc                                                         <  2:48612/35‑1 (MQ=38)                                                                          tgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACtt  >  2:364983/1‑90 (MQ=255)                                                                                       |                                                                             AGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTT  >  NZ_CP009273/2678637‑2678799

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