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%	57.0 / NA	18	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 (7/11);  total (7/11)

GCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTA  >  NC_002947/352935‑353193                                                                                                                                     |                                                                                                                               gCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGCTCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTccc                                                                                                                 >  5:169891/1‑148 (MQ=255)         ggTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGttt                                                                                                         <  1:66319/148‑1 (MQ=255)              tGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGc                                                                                                       >  6:68807/1‑145 (MQ=255)                         aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                         <  3:8519/148‑1 (MQ=255)                               cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGa                                                                                       <  3:144647/144‑1 (MQ=255)                                 cTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                                 <  5:68807/148‑1 (MQ=255)                                        gAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTCCCCCAGGGTTGGTGAGCTCGCGGTGGGCGCGGGCGGGCCcg                                                                          >  7:68534/1‑148 (MQ=37)                                             cacaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAg                                                                     >  6:116811/1‑148 (MQ=255)                                              acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGGTGAGCTCTCTGTGGGAGCGGGCGTGCCCGCGAAgg                                                                    >  8:140426/1‑148 (MQ=255)                                                      ggcggCCAAGGCCTACTGGGCTGCGGGGGGGATCGGCGTCCGCTGACCCATTGCCGGGGGCGGCTTTGCCGCCGTTCGGGGGCCCGCCGGCTCCCACTGGTTTGCTGCTCGCTCTGGGGGAGGGGGCGTGCCCGCGAAGGGCCGCGaa                                                            <  6:145504/145‑1 (MQ=255)                                                      ggAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTTCGGAGATATCGGTGGGAGCGGGCGGGCCCGCCAACGGCcgcgca                                                            >  4:21397/1‑146 (MQ=37)                                                          gCCAAGTCCTACATGGCTGCGCGCGGTAGCGTCGTGCGCTGATCAATGGCCCGGGGCTGCTTTGCTGCCCTTCGCGGTCATGCCCGCTCCCACAGGTTTGCTTAGATCTCTGGGGGAGCGGGCGGGCCCGCGAAGGGCCGCGAAGCgg                                                        <  6:117468/148‑1 (MQ=255)                                                                      aTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCGCCCACAGGATTGCAGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCCCCAAGCGGCCCCAAATGca                                             >  6:72284/1‑146 (MQ=37)                                                                               cgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTgg                                   <  6:71379/148‑1 (MQ=255)                                                                                      gcggggtTCGATGATCGATTGCCGGGGGCGGGTTGGCAGCCCTGCGCGGGCACGCCCGCGCCCCCTGTTTTGCTGTGTGTGCTGTGGGTGCGGGGGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGc                            <  5:116811/143‑1 (MQ=19)                                                                                                     ccATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGt             <  2:58957/148‑1 (MQ=255)                                                                                                       atggACCGGGGCTGCTTTGCAGCCCGTCGCGGGCACGCCCGCTGCCACAGGTTTGGTGTGATCTCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTa            <  2:44373/144‑1 (MQ=38)                                                                                                                ggCTGCTTTGCAGCCCTGCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTa  <  8:68534/148‑1 (MQ=255)                                                                                                                                     |                                                                                                                               GCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTA  >  NC_002947/352935‑353193

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