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%	122.2 / NA	37	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 (20/17);  total (20/17)

GGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGC‑GGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCC  >  NC_002947/352925‑353197                                                                                                                                               |                                                                                                                                     ggAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACg                                                                                                                                >  8:1562/1‑149 (MQ=255)      ggCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGc                                                                                                                           >  8:169638/1‑149 (MQ=255)       gCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCt                                                                                                                          <  5:128390/149‑1 (MQ=255)       gCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCt                                                                                                                          <  7:269037/149‑1 (MQ=255)             cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTcccc                                                                                                                      >  1:115219/1‑146 (MQ=255)                 gcgGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGtt                                                                                                               <  1:87738/149‑1 (MQ=255)                   gatGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCCCTCCCACACGGttt                                                                                                              >  1:138020/3‑149 (MQ=39)                     tGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCt                                                                                                           <  7:1562/149‑1 (MQ=255)                            ggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGAtc                                                                                                     >  7:13255/1‑148 (MQ=255)                               cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCtg                                                                                                 >  1:25266/1‑149 (MQ=255)                                aGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGAtc                                                                                                     <  3:107074/144‑1 (MQ=255)                                aGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGAtc                                                                                                     >  4:107074/1‑144 (MQ=255)                                 gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCtgtg                                                                                               <  1:96317/149‑1 (MQ=255)                                  cATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTgg                                                                                              >  4:79349/1‑149 (MQ=255)                                  cATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTgg                                                                                              <  7:169638/149‑1 (MQ=255)                                      ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACC‑GGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATGTGTGTGGGAGCg                                                                                         >  3:9079/1‑149 (MQ=255)                                        ccACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGa                                                                                            <  5:89099/145‑1 (MQ=255)                                          aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                                      <  8:13255/149‑1 (MQ=255)                                                        acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTgag                                                                                                        <  3:209410/117‑1 (MQ=255)                                                        acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTgag                                                                                                        >  4:209410/1‑117 (MQ=255)                                                          aGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCATTCGCGGGCACGCCCGCTCCCCCA‑GGGTTGCTGAGATGTGTGTGGGAGCGGGCGTGCCCGCGAAGGGcc                                                                      >  7:290811/1‑149 (MQ=39)                                                              ccGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTACCCCA‑GGGTTGCTGGGGTGTCTGGGGGAGCGGGGGGGGCCGCGCACGGCcgcgc                                                                  >  1:131648/1‑148 (MQ=37)                                                                 gAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCCCCCCCGCCCCCA‑GGGTTTGGGAGGTAAGTGTGGGAGCCGGGGGGCCCCCGCAGGGCcccgcacc                                                               >  4:23620/1‑145 (MQ=255)                                                                    gCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGC‑GGc                                                           >  7:150078/1‑149 (MQ=255)                                                                      cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGC‑GGccc                                                         >  8:178406/1‑149 (MQ=255)                                                                            cTAAATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGC‑GGCCCCAATTg                                                   <  2:91639/149‑1 (MQ=255)                                                                               cATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGC‑GGCCCCAATTGCtt                                                <  2:18947/149‑1 (MQ=255)                                                                                aTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGC‑GGCCCCAATTGCTTc                                               <  2:25266/149‑1 (MQ=255)                                                                                     ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCGCCCCCA‑GGGTTGGTGAGCTCTGTGTGGGGGCGGGCGGGGCCGCGCAGGGCCCCGCAGCGGGCCCCAATTG‑TTCAAccc                                          >  1:279063/1‑148 (MQ=25)                                                                                        gcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGC‑GGCCCCAATTGCTTCAACCGTg                                        <  3:79349/148‑1 (MQ=255)                                                                                           cgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGC‑GGCCCCAATTGCTTCAACCGTGGTTc                                    >  3:171982/1‑149 (MQ=255)                                                                                              gATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCAACCCCA‑GGTTTGCTGAGAACTCTGTGGGAGCGGGCCTGCCCCCGAAGGGCCCCCAATC‑GGCCCCAATTGCGTAAAACCTGGTTCTTg                                 >  4:108526/1‑149 (MQ=255)                                                                                                             aTCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGC‑GGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                   <  7:178406/148‑1 (MQ=255)                                                                                                                 aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCAACCCCA‑GGTTTGGGGAGGTAGCTGTGGGGGCGGGCGGTGCCGCCCAGGGCCCCGGAGC‑GGCCCCAAATGGGTGAACCGTGGGTAGTGCCGGGGCGCGCGTGGTa                >  2:85155/1‑147 (MQ=255)                                                                                                                     aCCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGCAGC‑GGCCCCCATTGCTTCAACCGTGGTTCTTGCCGAGCGGCGCGTGGTAAAGCTc          >  7:230728/1‑149 (MQ=38)                                                                                                                      ccGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGC‑GGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg         <  5:122249/149‑1 (MQ=255)                                                                                                                             tgctTTGCAGCCCTGCGCGGCCACGCCCGCTCCCACA‑GTTTTGCTGAGAGCTCTGTGGGAGGGGGGGTGCCCGCGAAGGGCCGCGAAGC‑GGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATccc  <  3:23620/149‑1 (MQ=38)                                                                                                                                               |                                                                                                                                     GGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGC‑GGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCC  >  NC_002947/352925‑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: ‑