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%	94.1 / NA	30	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 (12/18);  total (12/18)

ATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCG  >  NC_002947/352933‑353190                                                                                                                                       |                                                                                                                            aTGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTATGGCGACTACACAGTGCCCGAGGCCAACGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGCGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTc                                                                                                                >  8:82720/1‑148 (MQ=255)   gCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCGccc                                                                                                              >  1:163361/1‑148 (MQ=255)          cgGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGtt                                                                                                       <  6:13539/148‑1 (MQ=255)            gTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTg                                                                                                     <  5:83612/148‑1 (MQ=255)                tGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                 <  3:110322/148‑1 (MQ=255)                  gcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGGCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTgaga                                                                                               <  1:72288/148‑1 (MQ=255)                   cggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTgag                                                                                                <  1:8075/146‑1 (MQ=255)                    ggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCCCACCCCCAGGGTTGcggcggtc                                                                                             >  8:147649/1‑141 (MQ=255)                       cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTCCCCCACGGTTGCTGAGATatct                                                                                          >  3:136927/1‑148 (MQ=39)                       cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATctct                                                                                          <  1:191422/148‑1 (MQ=255)                         gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgt                                                                                        <  6:72212/148‑1 (MQ=255)                              ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAg                                                                                   >  1:164183/1‑148 (MQ=255)                               tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTACCACAGGGTTGCTGAGCTCTCTGTGGGAGc                                                                                  >  3:149388/1‑148 (MQ=39)                                  aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCggg                                                                               <  6:82421/148‑1 (MQ=255)                                   cTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                              <  2:164183/148‑1 (MQ=255)                                              acacAGTGCCGGAGGCCAAGGCCTACATGGCTGCGCGCGGGAGCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCGTCGCGGGCACGCCCGCTCCCACAGGGTTGCTGAGATCTCTGGGGGAGCGGGCGTGCCCGCGaa                                                                   <  5:79219/148‑1 (MQ=38)                                                acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCACCCCCACGTTTGCTGAGATCTCGGTGGGGGCGGGCGTGCCCCCCAAgg                                                                 >  2:157110/1‑148 (MQ=255)                                                             ccAAGGCCTACATGGCTTGGCGCGGTATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTGCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGc                                                    <  8:119877/148‑1 (MQ=255)                                                              cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGGGCCCGCGCACGGCCGCGCAGCGGcc                                                   >  3:91518/1‑148 (MQ=38)                                                              cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTACCCCACGGTTGCTGAGATCTCTGTGGGAGCCGGCGGGCCCGCGAACGGCCCCGAAGCgggc                                                   >  5:138905/1‑146 (MQ=255)                                                                   ccTACATGGCTGGGCGCGGGAGGGTCGTGCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACTGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAAt                                              <  4:102614/148‑1 (MQ=37)                                                                        aTGGCTTCGCGCGGGAGCGGCGTGCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTGCGCGGGCACGCCCGCGCCCACCGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCtt                                         <  5:77588/148‑1 (MQ=38)                                                                           gCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCCTGTTTGCTGAGCTGTCTGGGGGAGCGGGCGTGCCCGCGCAGGGCCGCGCAGCGGCCCCACTTGCGTcga                                      >  7:111185/1‑146 (MQ=37)                                                                                 cgcgGTAGCGGCGTCCGCGTACCCATTGCCCGGGGCGGCTTTCCCCCCCTGCGCGGGCACGCCCGCTCCCACAGGTTTGCTGATATCTGTGTGGGTGGGGGGGGGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTgg                                <  2:138217/148‑1 (MQ=255)                                                                                           tcgtTCGCTGATCCTTTGACCGGGGGGGGTTTGCCGCCGTGCGCGGGAACGCCGGCTCCCACTGGTTGGCTGAGATCTCTGTGGGTGGGGGGGGGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGa                      <  2:156622/148‑1 (MQ=255)                                                                                           tcgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGCAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGa                      >  8:159348/1‑148 (MQ=255)                                                                                                        cATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCACCCACAGGTTTGCTGAGCTATCTGTGGGAGCGGGCGTGCCCGCGCAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGt          >  8:34928/1‑147 (MQ=38)                                                                                                          ttGACGGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaaa       <  4:136927/148‑1 (MQ=255)                                                                                                               cGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg  <  2:21957/148‑1 (MQ=255)                                                                                                               cGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg  <  8:106217/148‑1 (MQ=255)                                                                                                                                       |                                                                                                                            ATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCG  >  NC_002947/352933‑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: ‑