Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	353,067	A→G	intergenic (+34/+84)	hisF → / ← cbcV	imidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	353,067	0	A	G	100.0%	40.7 / NA	14	intergenic (+34/+84)	hisF/cbcV	imidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (8/6);  total (8/6)

GAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCG  >  NC_002947/352926‑353190                                                                                                                                              |                                                                                                                             gAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACg                                                                                                                        >  5:59675/1‑148 (MQ=255)      gCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCcccccccc                                                                                                                   >  1:163359/1‑146 (MQ=255)                       tGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCAAGCCCGCTCCCACA‑GGTTTGc                                                                                                    >  3:126355/1‑145 (MQ=255)                         gcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCCCCCCCCCACGGTTTgtggga                                                                                                >  1:76729/1‑143 (MQ=255)                                          cTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                              <  1:53571/148‑1 (MQ=255)                                                            catgtaGGCCTTGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGCTCTCTGTGGGCGCGGGCGTGCCCGCGAAGGGCcgc                                                            >  3:75681/6‑148 (MQ=38)                                                              cGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCAACCACA‑CGTTTGCTGACATCACTGTGGGAACGGGCGTTCCCGCCAAGGGCCCCGa                                                          >  2:125570/1‑148 (MQ=255)                                                                     cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGcc                                                   >  3:70612/1‑148 (MQ=255)                                                                        ggCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCGCTCCCACA‑GGTTTGCTGAGCTCTATGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCa                                                >  8:5332/1‑148 (MQ=255)                                                                                tGGCTTCGCGCGGGAGCGTCGTGCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTGCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTc                                        <  1:2384/148‑1 (MQ=255)                                                                                                  gcgtGCGCTTATCCATTGCCCGGGGCTGCTTTGCAGCCGTGCGCGGGCACGCCCGCGCCCACC‑GGTTTGCGGAGATCTCTGTGGGAGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGa                      <  2:32057/147‑1 (MQ=37)                                                                                                                   gagcGGGGCTGCTTGGCAGCCGTGGGCGGGCACGCCCGAGCCCACA‑GTTTTGCTGATATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAg      <  1:125570/144‑1 (MQ=37)                                                                                                                    aCCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGc     <  6:59675/147‑1 (MQ=255)                                                                                                                      cGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg  <  4:38865/148‑1 (MQ=255)                                                                                                                                              |                                                                                                                             GAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCG  >  NC_002947/352926‑353190

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: ‑