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%	46.0 / NA	15	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/7);  total (8/7)

TCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTA  >  NC_002947/352922‑353193                                                                                                                                                  |                                                                                                                               tCTGGAACGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTGCACAGGGCCGGAGGCCAAGGCCGACATGGCTGCGCGCGGGATCGTCGTTCGCTGATCCATTGCCCGGGGCTGCTTTGCAGCCCTTCGCggg                                                                                                                              <  7:111939/148‑1 (MQ=255)       aGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCGTCGCGGGCACGcc                                                                                                                        <  6:41319/148‑1 (MQ=255)          gCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGc                                                                                                                     >  7:107098/1‑148 (MQ=255)            ccgCCAGCACCGCGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTc                                                                                                                   >  1:47406/3‑148 (MQ=39)            aTGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTc                                                                                                                   >  7:134073/1‑148 (MQ=255)                                gcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCGCACCCACAGGTTTGCTGAGAtct                                                                                               >  6:42167/1‑148 (MQ=255)                                                          cacaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCCCCCCCTCCCCCCGGGTTGGGGGGGTCTCTGTGGGTGCGGGCGTGCCcgcgccg                                                                     >  1:40289/1‑145 (MQ=255)                                                                gCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcgc                                                               <  8:134073/148‑1 (MQ=255)                                                                gCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGAGCGTCGTGCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCGGCGCGGGCACGCCCGCTCCCACCGGTTTGCTGAGAGCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcgc                                                               <  5:58705/148‑1 (MQ=38)                                                                                   cgggCTGCGCGCGGGAGCGTCGTTCGCTGAGCCATTGACCGGGGCGGCTTTGCAGCCCTGCGCGGGCACGCCCGCTCCCACAGGTTGGCTGAGATCGCTGTGGGAGCGGGCGGGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCtt                                            <  7:124498/146‑1 (MQ=37)                                                                                                   tcgtcgTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGCGGAGATCGCGGTGGGAGCGGGCGTGGCCGCGAAGGGCCGCGCAGCGGCCCCAATTGCTTCAACCGTGGGTCGTGc                            >  4:62673/1‑148 (MQ=37)                                                                                                            gCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagcagc                    <  6:73824/147‑1 (MQ=255)                                                                                                                     ccgcgaagggCTGCAAAGCAGCCCTTCGCGGGCCCGCCCGCTCCCCCAGGGTTGGTGAGATCTCGGTGGGCGCGGGCGTGGCCGCGAAGGGCCGCGCAGCGGCCCCAATTGCTTCCCCCGTGGTTCTTGCCGAGCAGCGCGGGGGaaa          >  8:64398/8‑148 (MQ=255)                                                                                                                        aCCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCACCCACAAGTTTGCTGAGAAAACTGTGGGAGCGGGCGTGCCCCCGCACGGCCCCCAACCGGCCCCAAATTCTTAAAACCTGGTTATTGACCAGCACCCCCTGGTaaaaca       >  2:119881/1‑145 (MQ=255)                                                                                                                          cGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCGCCCACAGGTTTGCTGAGATCTCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg     <  3:51602/148‑1 (MQ=255)                                                                                                                             ggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGACATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTa  <  4:95079/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                               TCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTA  >  NC_002947/352922‑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: ‑