Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,011,323	+T	intergenic (‑258/+192)	PP_0871 ← / ← SBRG_Ecoli_galM	glycine betaine/carnitine/choline ABC transporter permease/SBRG_Ecoli_galM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,011,323	1	.	T	100.0%	34.8 / NA	12	intergenic (‑258/+192)	PP_0871/SBRG_Ecoli_galM	glycine betaine/carnitine/choline ABC transporter permease/SBRG_Ecoli_galM
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (3/9);  total (3/9)

AAGCGGGTTGCAGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATG‑CTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGT  >  NC_002947/1011204‑1011427                                                                                                                         |                                                                                                         aaGCGGGTTGCAGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCtttt                                                                                           <  1:316044/136‑1 (MQ=255)          tcAGGTGTGTGTAGGTCTCCGTGGGCTTAACGGGACCAATCTTAGGCTGCATTTGTCAATTGTTAATCGCTGGGTTGCGGTTTTTTGGGCGAGAAAACCAGCGACTTGATGTCTCAAATCATCTTTTAATGataat                                                                                  <  1:15272/135‑1 (MQ=255)                                                     tggggtcATTTGGCCAGTTTTAATCTTTGGGTTGCGGTTTTTTGGGCGACAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTAc                                       <  1:408421/130‑1 (MQ=255)                                                                         tAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATa                   >  1:42610/1‑136 (MQ=255)                                                                            ttgttGGGTTGTGGGTTTTTTTGCTATAATAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATATTTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTg                <  2:233693/132‑1 (MQ=255)                                                                            tCGCTGGGTTGCGGGTTTTTGGGCGATAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCATCAGCAATGCTTACGCATAACCATAGCGAAAATAGTg                <  2:357796/136‑1 (MQ=255)                                                                            tCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTg                <  1:187505/136‑1 (MQ=255)                                                                            tCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTg                <  2:394730/136‑1 (MQ=255)                                                                            tCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTg                <  2:42610/136‑1 (MQ=255)                                                                                     tGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAgtgt       >  1:297645/1‑136 (MQ=255)                                                                                          gTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGgt  >  1:95837/1‑136 (MQ=255)                                                                                          gTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGgt  <  2:221650/136‑1 (MQ=255)                                                                                                                         |                                                                                                         AAGCGGGTTGCAGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATG‑CTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGT  >  NC_002947/1011204‑1011427

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