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%	108.0 / NA	32	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 (18/14);  total (18/14)

TGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCC  >  NC_002947/352924‑353207                                                                                                                                                |                                                                                                                                             tGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCa                                                                                                                                          >  4:166443/1‑148 (MQ=255)      gggCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGgcccgccc                                                                                                                                     >  1:115969/1‑148 (MQ=255)               aGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCCCCCCCCCCccca                                                                                                                            >  6:115974/1‑148 (MQ=255)               aGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCcaca                                                                                                                            >  1:50397/1‑148 (MQ=255)                      tGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGc                                                                                                                     <  1:96494/148‑1 (MQ=255)                       gCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                                    <  2:50397/148‑1 (MQ=255)                         tGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                                  <  1:68651/148‑1 (MQ=255)                                     tCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCCGGTTTGCTGAGATCGCtgtg                                                                                                        >  7:16793/1‑146 (MQ=255)                                     tCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgtg                                                                                                        >  2:28523/1‑146 (MQ=255)                                        tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCAGTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTCCCCCACGTTTGCTGAGATCACTCTGGGacc                                                                                                   >  5:154939/1‑146 (MQ=38)                                                        cacaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAg                                                                                   >  4:169906/1‑148 (MQ=255)                                                         acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGGGTGGTGAGAGGTCTGTGGGCGCGGGCGGGCCcgcgccgg                                                                                  >  5:161465/1‑144 (MQ=37)                                                           aGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTCCCACAGGGTTGCTGAGATCTCTTTGGGAGCGGGCGTGCCCGCGCAGGGc                                                                                >  6:101772/1‑148 (MQ=39)                                                             tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcg                                                                              >  5:92369/1‑148 (MQ=255)                                                             tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCACCCCCAGGTTTGCTGAGATATCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcg                                                                              >  4:149331/1‑148 (MQ=255)                                                                      ccAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGc                                                                     <  4:51761/148‑1 (MQ=255)                                                                       cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTgaga                                                                                                                >  1:147664/1‑104 (MQ=255)                                                                        aaGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAt                                                                                                               <  2:147664/104‑1 (MQ=255)                                                                           gCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTGGCAGCCCTGCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCaa                                                                <  2:131892/148‑1 (MQ=255)                                                                                  tgggtgcgcgcgGGATGGTCGTTCGCTGATCCTTTGACCGGGGGTGCTTTGCAGCCCGTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTGTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTc                                                         <  4:80946/142‑1 (MQ=37)                                                                                    gCTTCGCGCGGGATCGTCTTGCGCTCAACCATTCACGGGGTCGGCTTTGCAGCCCTGCGCGGGCACGCCCGCTCCCACCGGTTTGCTCCCATCTCTGGGGGTGCGGGGGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTc                                                         <  6:8407/146‑1 (MQ=25)                                                                                              ggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCt                                             <  6:65652/148‑1 (MQ=255)                                                                                               gATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCtt                                            <  6:92369/148‑1 (MQ=255)                                                                                                       ttCGCTGATCCATTGACCGGGTCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCAACCACAGGTTTGCGGAGATGTCTGTGGGAGCGGGCGGGCCCGCGAAGGGCCGCGCAGCGGGCCCAAATGCTTGAACCCGGGTTCTTTCCGagca                                    >  6:177022/1‑148 (MQ=255)                                                                                                          gCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAgcg                                 >  4:52924/1‑148 (MQ=255)                                                                                                                  aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCGCTCCCACAGGTTTTCTTAGATCTCTGTGGGAGCCGGCGTGCCCGCCAACGGCCGCGAACCGGCCCCAATTGCTTAAACCCTGGTTCTTGCCGAGCACCGCCTGGTa                          >  6:147414/1‑147 (MQ=37)                                                                                                                       ccGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTc                    <  4:44811/148‑1 (MQ=255)                                                                                                                        cGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                   <  3:166443/148‑1 (MQ=255)                                                                                                                        cGGGGCTGCTTGGCACCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGACATCGCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                   <  5:177022/148‑1 (MQ=38)                                                                                                                                 tttGCAGCCCTTCGCGGGCACGCCCGCTACCACAGGTTTGCTGAGATATCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAAATGCATCAACCGTGGTTATTTCCGACCAGCGCGTGGTAAACCTAGGTATacc            >  8:29363/1‑146 (MQ=37)                                                                                                                                     cAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAAACGTGGTTCTTGCCGCGCAGCGCGTGGTAAAGCTCGGTATCCCCGCGTa      >  8:105379/1‑148 (MQ=255)                                                                                                                                         ccTGCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATccc  <  2:115969/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                             TGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCC  >  NC_002947/352924‑353207

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