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%	80.2 / NA	26	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 (11/15);  total (11/15)

TCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCC  >  NC_002947/352922‑353197                                                                                                                                                  |                                                                                                                                    tCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCggg                                                                                                                                   <  3:68392/148‑1 (MQ=255)                 tggaagatGCTGGCCGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCcaca                                                                                                                   >  2:39901/8‑148 (MQ=255)                 aGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCcaca                                                                                                                   >  3:158777/1‑148 (MQ=255)                         gCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCt                                                                                                          <  2:18978/148‑1 (MQ=255)                         gCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCt                                                                                                          <  2:2473/148‑1 (MQ=255)                                 ccAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGAtctc                                                                                                  >  7:144005/1‑148 (MQ=255)                                  cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATctct                                                                                                 <  4:176037/148‑1 (MQ=255)                                      aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTgg                                                                                             <  4:124734/148‑1 (MQ=255)                                              cTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                                     <  2:127498/148‑1 (MQ=255)                                                  ggCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGt                                                                                   <  8:144005/146‑1 (MQ=255)                                                          cacaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGACCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCA‑GGGTTGGTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGCAg                                                                         >  1:134711/1‑148 (MQ=38)                                                             aGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCAGCCCCCTACCACACGGTTTG‑GGTGATCTAGGTGGGAGCGGGCGTGCCCGCGGAgggg                                                                      >  1:165975/1‑147 (MQ=255)                                                                  cGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCCCTCCCCCC‑CGGTTGGGGGGCTATCTGTGGGGGCGGGGGGGGCCGCGCCCGGCcgcgc                                                                 >  1:172902/1‑147 (MQ=255)                                                                         cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTCCCACA‑CGTTTGCTGAGAAATCTGTGGGAGCCGGCGTGCCCGCGAAGGGCCGCGCAGCGGcc                                                          >  6:134468/1‑148 (MQ=38)                                                                                   aTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCtt                                                <  8:92487/148‑1 (MQ=255)                                                                                        ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGc                                                                                                           >  5:47617/1‑84 (MQ=255)                                                                                        ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAAcc                                           >  2:16909/1‑148 (MQ=255)                                                                                         tCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCt                                                                                                          <  6:47617/84‑1 (MQ=255)                                                                                                         ttCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagca                          <  2:163722/148‑1 (MQ=255)                                                                                                                 tCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                   <  5:39754/147‑1 (MQ=255)                                                                                                                 tCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                   <  1:16909/147‑1 (MQ=255)                                                                                                                     ttGACCGGGGCTGCTTTGCTGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaaa              <  5:92640/148‑1 (MQ=255)                                                                                                                          cGGGGCTGCTTTGCAGCCCTGCGCGGGCACGCCCGCCCCCACC‑GGTTTGCTGAGAGCGCGGGGGGAGCGGGCGGGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg         <  1:40763/148‑1 (MQ=38)                                                                                                                          cGGGGCGGCTTTGCAGCCCGTCGCGGGCACGCCCGCTCCCACC‑TGTTTGCTGAGATCTCTGTGGGAGCGGGCGGGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg         <  5:175505/148‑1 (MQ=38)                                                                                                                             ggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTa      <  2:134711/148‑1 (MQ=255)                                                                                                                                   tttGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGCTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCCACCGTGGGTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATccc  >  6:81895/1‑146 (MQ=255)                                                                                                                                   tttGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATccc  >  2:32562/1‑146 (MQ=255)                                                                                                                                                  |                                                                                                                                    TCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCC  >  NC_002947/352922‑353197

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