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%	159.9 / NA	48	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/28);  total (20/28)

GAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGA  >  NC_002947/352926‑353200                                                                                                                                              |                                                                                                                                        gAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑c                                                                                                                                   <  5:180287/147‑1 (MQ=255)     ggCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGc                                                                                                                             <  6:265661/149‑1 (MQ=255)        aTGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCc‑cccccccgcc                                                                                                                          >  6:29780/1‑144 (MQ=255)         cgCCAGCACCGCGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGCTCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTccc                                                                                                                         >  1:43241/2‑149 (MQ=39)            cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCca                                                                                                                        >  5:85242/1‑147 (MQ=255)             aGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑g                                                                                                                    <  4:158116/149‑1 (MQ=255)              gcgcGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑gg                                                                                                                   >  1:242421/1‑149 (MQ=255)               ccaggTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑g                                                                                                                    <  7:126888/144‑1 (MQ=255)                  ggTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTg                                                                                                               <  3:225641/149‑1 (MQ=255)                  ggTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTg                                                                                                               <  3:304726/149‑1 (MQ=255)                     gCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCGCCCCCCCCCCCA‑GGttttgg                                                                                                             >  1:368104/1‑146 (MQ=255)                     gCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTg                                                                                                            >  7:180335/1‑149 (MQ=255)                      cTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTga                                                                                                           <  3:10368/149‑1 (MQ=255)                       tGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCt                                                                                                             >  5:137701/1‑146 (MQ=255)                        ggcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCt                                                                                                             >  2:114811/1‑145 (MQ=255)                         gcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTgaga                                                                                                         >  6:117160/1‑148 (MQ=255)                                  aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTggg                                                                                               <  2:43241/149‑1 (MQ=255)                                        cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCgg                                                                                          <  8:131814/148‑1 (MQ=255)                                         aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCggg                                                                                         >  1:130633/1‑148 (MQ=255)                                         aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCggg                                                                                         <  2:130633/148‑1 (MQ=255)                                         aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                                        <  1:114811/149‑1 (MQ=255)                                                    tACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGccc                                                                                  <  6:296641/144‑1 (MQ=255)                                                    tACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGccc                                                                                  >  5:296641/1‑144 (MQ=255)                                                    tACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGGTTGCTGAGATATCAGTGGGAGCGGGCGTTCCcgcg                                                                               >  3:152270/1‑147 (MQ=38)                                                     acacAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CCCCCCCTACCCCACGGTTTG‑TGAGATAGATGTGGGAGCGGGGGTGGCCGCGaac                                                                            >  1:342308/1‑148 (MQ=255)                                                           tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCG‑TGCCCGcgaag                                                                                                                         >  4:247970/1‑95 (MQ=37)                                                           tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCG‑TGCCCGcgaag                                                                                                                         <  3:247970/99‑5 (MQ=37)                                                             ccGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTgg                                                                                                <  5:123759/121‑1 (MQ=255)                                                             ccGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTgg                                                                                                >  6:123759/1‑121 (MQ=255)                                                               gtaGGCCAAGGCCTAATTGGCTTCGCGGGTTATGGTGTTGCGCTTATCATTTCACCGGGCGTGCTTTGCTGCCGTGCGGGGGCC‑CCCCCGCTCCCACT‑GGTTTGCTTAGTTCTTTTTGGTTGCGGGCGTGCCCGCGAAGGGCCGCGAAg                                                                  <  8:180335/147‑1 (MQ=255)                                                                    ccAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGcc                                                             >  5:271635/1‑149 (MQ=255)                                                                     cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAg                                                                            <  5:162952/133‑1 (MQ=255)                                                                     cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGGTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAg                                                                            >  6:162952/1‑133 (MQ=255)                                                                      aaGGCCTACATGGCTTCGCGCGGGATCGGCGTTCGCTGATCCATTGACCGGGGCTGCTTGGCAGCCCTGCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGcccc                                                           <  3:237270/149‑1 (MQ=255)                                                                       aGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCa                                                          <  6:137701/149‑1 (MQ=255)                                                                         gCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTGCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAAt                                                        <  8:216243/149‑1 (MQ=255)                                                                           cTACATGGCTTCGCGCGGGAGCGTCGTCCGCTGATCCTTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACT‑GGTTTGCTGAGAGCTCTGTGGTAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTg                                                      <  2:242421/149‑1 (MQ=38)                                                                               aTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTc                                                  <  6:85242/149‑1 (MQ=255)                                                                                  ggttCGCGGGTGTTCGTCGTTCCGTCATCATTTCACGGGGGGTGCTTTCCTCCCGTTCGCGGCCC‑CGCCCGGTCCCACC‑GGTTTGCTTAGATGTCTGGGGTCGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCaa                                                <  2:342308/146‑1 (MQ=12)                                                                                                            agccATTGCCCGGGGCGGGTTTCCCGCCCGGCGCGGGCA‑CGCCCGCCCCCCCA‑GGTTGCCGGCCCTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTgg                     <  5:29780/147‑1 (MQ=21)                                                                                                                   gACCGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCt              >  5:381627/1‑149 (MQ=255)                                                                                                                     ccGGGGCTGCTTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg            <  6:271635/149‑1 (MQ=255)                                                                                                                       ggggCTGCTTTGCAGCCGTTCGCGGGCA‑CGCCCGCTCCCACT‑GGTTTGCTTTGTTCTCTGTGGGTGGGGGGGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGt          <  4:152270/149‑1 (MQ=38)                                                                                                                          gctgctTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAAc                                               <  1:200597/109‑1 (MQ=255)                                                                                                                              ctTTGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGCAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcccc    >  8:9360/1‑148 (MQ=255)                                                                                                                                 tGCAGCCCTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa  <  5:117160/147‑1 (MQ=255)                                                                                                                                        cTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGGTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGt          <  6:123124/132‑1 (MQ=255)                                                                                                                                        cTTCGCGGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGt          >  5:123124/1‑132 (MQ=255)                                                                                                                                              |                                                                                                                                        GAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCA‑CGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGA  >  NC_002947/352926‑353200

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