Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	3,827,923	G→T	R30S (CGT→AGT)	ptxS ←	2‑ketogluconate utilization repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,827,923	0	G	T	100.0%	56.4 / NA	17	R30S (CGT→AGT)	ptxS	2‑ketogluconate utilization repressor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (4/13);  total (4/13)

GGCCGCGCTTTAGGCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAACCGGTTT  >  NC_002947/3827806‑3828052                                                                                                                      |                                                                                                                                  ggCCGCGCTTTAGGCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACTGGAGACGGTGGCTTTGGa                                                                                                                 >  1:389544/1‑136 (MQ=255)      ggcTTTAGGCCGCGAGCCATCTGTTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACTGGAGACGGTGGCTTTGGAAACGc                                                                                                            <  2:387159/135‑1 (MQ=255)      cgcTTTAGGCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACTGGAGACGGTGGCTTTGGAAACGc                                                                                                            <  2:305547/136‑1 (MQ=255)                cgcgAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACTGGAGACGGTGGCTTTGGAAACGCCAGCAACAcg                                                                                                  <  2:537597/136‑1 (MQ=255)                cgcgAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACTGGAGACGGTGGCTTTGGAAACGCCAGCAACAcg                                                                                                  <  2:452130/136‑1 (MQ=255)                       cATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACTGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACt                                                                                           <  1:35123/136‑1 (MQ=255)                                ggggCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACTGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATg                                                                                  >  2:564404/1‑136 (MQ=255)                                   gCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACTGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTg                                                                               >  2:548013/1‑136 (MQ=255)                                         gCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACTGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGt                                                                         <  2:313626/136‑1 (MQ=255)                                         gCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACTGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGt                                                                         <  2:372079/136‑1 (MQ=255)                                         gCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACTGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGt                                                                         <  2:516469/136‑1 (MQ=255)                                         gCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACTGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGt                                                                         <  1:522002/136‑1 (MQ=255)                                         gCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACTGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGt                                                                         <  1:412595/136‑1 (MQ=255)                                                   tCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACTGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTAt                                                               <  1:548013/136‑1 (MQ=255)                                                           ctctTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACTGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGgc                                                       <  2:389544/136‑1 (MQ=255)                                                                                                         ccACCGATGTAACTGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGaaa         >  2:478373/1‑136 (MQ=255)                                                                                                                tGTAACTGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAACCGGttt  <  1:524951/136‑1 (MQ=255)                                                                                                                      |                                                                                                                                  GGCCGCGCTTTAGGCCGCGAGCCATCTGGTTGGGGCGATAGCCGAGGCGTTCGATGACCTCTTCCAGCCGCTTGGCCGTGGCTTCGGCCAGCAATTGGCGGTCGCCACCGATGTAACGGGAGACGGTGGCTTTGGAAACGCCAGCAACACGCGCCACTTCACTGATGGTGACCCGTTCGCGGGTATGGGCAGGCGCGTCGGTCATGGGTTTTCCTTATGATTGTTTTAATAGTGATGAAACCGGTTT  >  NC_002947/3827806‑3828052

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