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%	146.9 / NA	44	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 (25/19);  total (25/19)

AGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCC  >  NC_002947/352938‑353208                                                                                                                                  |                                                                                                                                                 aGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCcacac                                                                                                                               >  3:212170/1‑148 (MQ=255)      ggTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCATTCGCGGGCACGCCCGCTCCCACAGGTTTg                                                                                                                          >  1:439198/1‑149 (MQ=255)          cTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑ga                                                                                                                     <  6:286766/149‑1 (MQ=255)                 ccAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATctct                                                                                                              >  4:146587/1‑149 (MQ=255)                   aGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCtgt                                                                                                            <  2:1487/149‑1 (MQ=255)                       tCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCggg                                                                                                                                                >  1:119193/1‑110 (MQ=255)                       tCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCggg                                                                                                                                                <  2:119193/110‑1 (MQ=255)                        cTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAg                                                                                                       >  3:144958/1‑149 (MQ=255)                           ccACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGa                                                                                                        <  3:430776/145‑1 (MQ=255)                           ccACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCgg                                                                                                    >  2:146350/1‑149 (MQ=255)                            cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCggg                                                                                                   >  3:409307/1‑149 (MQ=255)                            cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCgg                                                                                                    <  4:353046/148‑1 (MQ=255)                             aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGc                                                                                                  <  3:146587/149‑1 (MQ=255)                                 tGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGt                                                                                                <  4:144958/147‑1 (MQ=255)                                   gCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGccc                                                                                            <  6:424867/149‑1 (MQ=255)                                     gAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCACCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCcg‑c                                                                                         >  2:226580/1‑149 (MQ=255)                                         acacAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGGAAGCCCCCGCCCCCCGGGTTTGT‑GAGATGTGGGGGGGTGCGGGCGTGCCcg‑cgcac                                                                                     >  3:209627/1‑146 (MQ=255)                                          cacaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAg                                                                                     <  1:100718/148‑1 (MQ=255)                                               tGCCGGCGGCCAAGGCCTACCGGGCTGCGCGCGGGATCGGCGTTCGCTGATCCATTGCCCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCCGGTTTGCT‑GAGAGCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCcgc                                                                               <  4:409307/149‑1 (MQ=37)                                                gCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGcccc                                                                               >  6:389719/1‑146 (MQ=255)                                                 ccGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCcgcg                                                                              >  7:299650/1‑148 (MQ=255)                                                         cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGt                                                                                                <  8:319248/123‑1 (MQ=255)                                                         cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGt                                                                                                >  7:319248/1‑123 (MQ=255)                                                         cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑ccc                                                                    >  7:292405/1‑149 (MQ=255)                                                             gCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTgg                                                                                                          >  5:67562/1‑109 (MQ=255)                                                             gCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTgg                                                                                                          <  6:67562/109‑1 (MQ=255)                                                                tACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGGCCCCCCCCCCCCCCACAGTTGGT‑GAGAAAAATGTTGGGGCCGGCGGGGCCG‑CCAAGCGCCCCCAAGAGC‑CCCCacaat                                                              >  4:2523/1‑144 (MQ=255)                                                                     ggCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCACCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCCTGCCCG‑CGAAGGGCCGCGAAGCGG‑CCCCCAATGCTTCaa                                                        >  2:128448/1‑149 (MQ=39)                                                                             gcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑CCCCAATTGCtt                                                           >  8:224188/1‑138 (MQ=255)                                                                             gcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑CCCCAATTGCtt                                                           <  7:224188/138‑1 (MQ=255)                                                                                        gtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCCCCGCCCCCGGGTTGCTGGGGAAAAAGGTGGGCCCGGGC‑TGCCCGCCGACGGGCCGC‑AAGCGGCCCCCAATGGGTGCAACCGGGGTGATTGccaccc                                      >  2:381227/1‑145 (MQ=255)                                                                                           cGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTGAACCGTGGTTCTTGGCGAGCcgcg                                  >  3:345661/1‑149 (MQ=255)                                                                                             cTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTCCCACAGGGTTTGT‑GAGATATATGTGGGCGCGGGGGTGGCCG‑CGCAGGGCCGCGAAGCGG‑CCCCAATTGCTTGAACCGGGGGTATTGCCGAGCAGcgcg                                >  4:47035/1‑149 (MQ=25)                                                                                                  ccATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCCCACCCACAGGGTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑GCCCCAATGCTTAAACCGGGGGTATTGGCGAGCCGCCCGGGGt                            >  6:5280/1‑148 (MQ=25)                                                                                                    aTTGACCGGGGCTGCTTTGCAGCCCTTCGGGGGCTCGCCCGCTCCCACTGTTTTGCT‑GATATCTCTGTGGGAGGGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaaa                         <  8:292405/149‑1 (MQ=38)                                                                                                         ccGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                    <  8:299650/149‑1 (MQ=255)                                                                                                            gggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTGTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTa                 >  3:111127/1‑149 (MQ=255)                                                                                                            gggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGGT‑GGGATCTCGGTGGGGGCGGGCGTGCCCG‑CGCAGGGCCGCGCAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGGGCAGCGCGGGGTAAAGCTCGGTa                 >  3:102138/1‑149 (MQ=37)                                                                                                                      gCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTGGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGt        <  5:5280/148‑1 (MQ=255)                                                                                                                         gCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATcc    <  4:102138/149‑1 (MQ=255)                                                                                                                          cccTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATccc   >  4:60151/1‑149 (MQ=255)                                                                                                                           ccTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcccc            <  6:335801/139‑1 (MQ=255)                                                                                                                           ccTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcccc            >  5:335801/1‑139 (MQ=255)                                                                                                                           ccTTCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑CAGATCTCTGTGGGAGGGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGCAAAGCTCGGTATCCCCGAGTATcccc  <  4:212170/149‑1 (MQ=255)                                                                                                                             ttCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTc                     <  7:263117/128‑1 (MQ=255)                                                                                                                             ttCGCGGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTc                     >  8:263117/1‑128 (MQ=255)                                                                                                                                  |                                                                                                                                                 AGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCT‑GAGATCTCTGTGGGAGCGGGCGTGCCCG‑CGAAGGGCCGCGAAGCGG‑CCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCC  >  NC_002947/352938‑353208

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