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%	110.8 / NA	34	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/16);  total (18/16)

TGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCG‑AGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCA  >  NC_002947/352924‑353209                                                                                                                                                 |                                                                                                                                                tcgAAGGGCATGCGAGCGCGGTGCGGGGGGCGATCATCTTCCACTTTGGCGCTACTCACAGTGCAGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCa                                                                                                                                             <  1:307419/147‑1 (MQ=255)                    ggTGCTGGCGGCCAGCATCTTCCACTTTGGCG‑AGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTg                                                                                                                        <  1:315056/149‑1 (MQ=255)                                cAGCATCTTCCACTTTGGCG‑AGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCtg                                                                                                            >  7:237202/1‑149 (MQ=255)                                 aGCATCTTCCACTTTGGCG‑AGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGcc                                                                                                                                         >  7:200632/1‑120 (MQ=255)                                 aGCATCTTCCACTTTGGCG‑AGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGcc                                                                                                                                         <  8:200632/120‑1 (MQ=255)                                         ccACTTTGGCG‑AGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCCCTCCCACA‑CGTTTGCTGAGATATCTGTGGGAGCgg                                                                                                   >  5:227515/1‑149 (MQ=39)                                                     aGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCCCTCCCCCACGGTTTG‑GGAGATCTCTGTGGGAGCGGGCGTGCCcgcg                                                                                        >  6:14590/1‑149 (MQ=255)                                                              tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCtg                                                                                                            <  3:158877/120‑1 (MQ=255)                                                              tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCtg                                                                                                            >  4:158877/1‑120 (MQ=255)                                                                 cGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGaa                                                                            <  6:345202/149‑1 (MQ=255)                                                                    cggCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTGCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTGCGCGGGCACGCCCGCTCCCACC‑GGTTTGCTGACATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCg                                                                         <  5:331023/148‑1 (MQ=39)                                                                        cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTcc                                                                                                                                   >  2:63668/1‑88 (MQ=255)                                                                             ccTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCca                                                                                                                                 >  6:379516/1‑85 (MQ=255)                                                                             ccTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCca                                                                                                                                 <  5:379516/85‑1 (MQ=255)                                                                             ccTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTcc                                                                                                                                   <  3:281604/83‑1 (MQ=25)                                                                             ccTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCCGCTTCGCGGCCCTTCGCGGGCACGCCCGCTcc                                                                                                                                   >  4:281604/1‑83 (MQ=25)                                                                                aCATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCt                                                             >  8:342387/1‑149 (MQ=255)                                                                                    ggCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTgg                                                                                                         <  4:316101/101‑1 (MQ=255)                                                                                    ggCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTgg                                                                                                         >  3:316101/1‑101 (MQ=255)                                                                                            gcgGGATCGTCGTTGGCTGATCCATTGACCGTTGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTc                                                                                                                                    <  8:231004/67‑1 (MQ=21)                                                                                            gcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTc                                                                                                                                    >  7:231004/1‑67 (MQ=255)                                                                                               ggATCGGCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCtt                                              <  5:14590/149‑1 (MQ=255)                                                                                                gATCGTCGTTCGCTGATCCTTTCACCGGGGCTGCTTTGCAGCCGTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGGGGGAGGGGGGGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTg                                             <  6:227515/149‑1 (MQ=38)                                                                                                       gtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagc                                       <  5:346068/148‑1 (MQ=255)                                                                                                       gtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGag                                        >  3:116136/1‑147 (MQ=255)                                                                                                                 ccATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCGCCCCCA‑GGTTTGGTGGGAACACGGTGGGAGCGGGCGGGCCCGCGACGGGCCGCGCAGCGGGCCCAAATGCTTAAACCGGGGGTCGTGCCGAGCCGCGCggggg                             >  4:60863/1‑147 (MQ=19)                                                                                                                   aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTa                            >  1:47561/1‑147 (MQ=255)                                                                                                                               tgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATccc              >  6:105801/1‑149 (MQ=255)                                                                                                                                gctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCCGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcccc             >  6:292355/1‑149 (MQ=255)                                                                                                                                       aGCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTATGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTTCTTCAACCGTGGTTCTTTCCGAGCAGCGCGTGGTAAATCTCGGTATCCCAGAGTATc      >  2:224227/1‑149 (MQ=38)                                                                                                                                        gCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATcc     <  4:116136/149‑1 (MQ=255)                                                                                                                                        gCCCTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTCAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                               <  5:273074/123‑1 (MQ=255)                                                                                                                                           cTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCcca  >  4:314402/1‑149 (MQ=255)                                                                                                                                           cTTCGCGGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCcca  <  7:342387/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                TGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCG‑AGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCA  >  NC_002947/352924‑353209

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