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%	82.1 / NA	25	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 (17/8);  total (17/8)

TCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCA‑CAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGA  >  NC_002947/352922‑353200                                                                                                                                                  |                                                                                                                                       tCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACAAATTTCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGAACGTGTCTGCTTTGCAGCCCTTCGCggg                                                                                                                                      <  6:138554/148‑1 (MQ=255)               ccAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTccc                                                                                                                         >  8:56123/1‑147 (MQ=255)                 aGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCca‑ca                                                                                                                     <  3:81785/148‑1 (MQ=255)                           tGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCTCGCGGGATCGTCGTTCACTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCA‑CAGGTTTGCTga                                                                                                           <  2:7705/148‑1 (MQ=255)                           tGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCA‑CAGGTTTGCTga                                                                                                           >  1:121915/1‑148 (MQ=255)                                  cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTCCCA‑CAGGTTTGCTGAGctctct                                                                                                    >  7:45655/1‑148 (MQ=255)                                             aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTCCCA‑CACGTTTGCTGAGATCACAGTTGGAGCggg                                                                                         >  8:23956/1‑148 (MQ=38)                                                ttGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCA‑CAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGt                                                                                      <  8:163182/148‑1 (MQ=255)                                                   gagAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCA‑CCGGGTTGCTGAGGTCTCTGTGGGCTCGGGCggggc                                                                                   >  6:181203/3‑144 (MQ=38)                                                        tACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCA‑CAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCcgc                                                                                >  8:142046/1‑146 (MQ=255)                                                        tACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCACCCACCCGGTTTG‑GGAGATATCAGTGGGGGCGGGCGGGCCcgcca                                                                              >  4:18161/1‑146 (MQ=255)                                                         acacAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCCCGCCCC‑CACGTTTGCTGAGATCGCTGGGGGAGCGGGCGTGCCCGCCaa                                                                             >  6:183479/1‑148 (MQ=255)                                                                 atgtaGGCCTTGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCA‑CAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcgc                                                                      >  6:95621/5‑147 (MQ=255)                                                                           aGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCA‑CAGGGTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGccc                                                            >  4:75038/1‑147 (MQ=255)                                                                                  cATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCA‑CACGGTTGCTGAGGTCTGAGTGGGAGCGGGCGTGCCCCCCCAGGGCCGCGCAGCGGGCCCACTTGCt                                                    >  4:140173/1‑148 (MQ=37)                                                                                   aTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCA‑CAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGc                                                     <  6:56854/146‑1 (MQ=255)                                                                                       ctgcgcgcgTGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCA‑CAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAAc                                               <  8:88144/145‑1 (MQ=255)                                                                                        ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTCCCA‑CAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCaaac                                              >  8:185347/1‑146 (MQ=255)                                                                                           gcgcgGGATCGGCGTGCGCTGATCCATTGACGGGGGCTGCTTTGCAGCCCTGCGCGGGCACGCCCGCTCCCA‑CAGGTTTGCTGAGATCTGTGTGGGTGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTg                                           <  4:122009/148‑1 (MQ=38)                                                                                            cgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCGCTCCCA‑CAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGCAGGGCCGCGAAGCGGCCCCAAATTCTTCAACCGTgg                                          >  2:113872/1‑148 (MQ=38)                                                                                                      gcgtTCGCTGAGCCATTGACCGGGGCGGCTTTGCAGCCCTGCGCGGGCACGCCCGCTCCCA‑CAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGa                                <  3:19305/147‑1 (MQ=255)                                                                                                                  ccATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCA‑CAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTgg                     >  6:32379/1‑147 (MQ=255)                                                                                                                    aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCA‑CAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGCAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGt                    >  8:103414/1‑146 (MQ=255)                                                                                                                         gaagggCTGCAAAGCAGCCCTTCGCGGGCCCGCCCCCAACCA‑CAGGTTAGCTGAGATGTATGTGGGAGCGGGCGTGCCCGCGCAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGGGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTc             >  6:93017/4‑148 (MQ=25)                                                                                                                                   tttGCAGCCCTTCGCGGGCCCGCCCGCTCCCA‑CAGGTTTGCTGAGATCTCTGTGGGCGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGGCCCAATTGCTTCAACCCTGGTTCTTTCCGAGCAGCGCGGGGTAAAGCGCGGTATccc     >  2:26576/1‑146 (MQ=38)                                                                                                                                    ttGCAGCCCTGCGCGGGCACGCCCGCTCCCA‑CTGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCAAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa  <  2:91271/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                       TCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCA‑CAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGA  >  NC_002947/352922‑353200

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