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%	149.6 / NA	45	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 (23/22);  total (23/22)

CTGGAAGGGCATGCCAGCGCGGTGCTGGC‑GGCCAGCATCTTCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑AGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGC‑GAAGCGGC‑CCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCC  >  NC_002947/352923‑353208                                                                                                                                                    |                                                                                                                                                cTGGAAGGGCATGCCAGCGCGGTGCTGGC‑GGCCAGCATCTTCCACTTTGGCGAGTAGCACAGTGCCGGAGGCCAAGGCATACATGGATTCGCGCGGTAGCGTCGTTCGGTGATCCATTGCCCGGGGCTGCTTTGCAGCCATTCCC‑GGGc                                                                                                                                              >  1:90052/1‑149 (MQ=255)   ggAAGGGCATGCCAGCGCGGTGCTGGC‑GGCCAGCATCTTCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACg                                                                                                                                           >  2:124457/1‑149 (MQ=255)       gggCATGCCAGCGCGGTGCTGGC‑GGCCAGCATCTTCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCGTGCCCg                                                                                                                                       <  7:224296/149‑1 (MQ=255)       ctggccgCCAGCACCGCGCTGGCATGCCAGCATCTTCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGgcccgccc                                                                                                                                        >  4:77370/7‑149 (MQ=255)          cATGCCAGCGCGGTGCTGGC‑GGCCAGCATCTTCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGgcccgcccccac                                                                                                                                    >  6:205910/1‑145 (MQ=255)           ccgCCAGCACCGCGCTGGC‑GGCCAGCATCTTCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTcc                                                                                                                                   >  6:145916/3‑149 (MQ=255)            tGCCAGCGCGGTGCTGGC‑GGCCAGCATCTTCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTccc                                                                                                                                  <  2:90052/149‑1 (MQ=255)                aGCGCGGTGCTGGC‑GGCCAGCATCTTCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAg                                                                                                                              >  8:158844/1‑149 (MQ=255)                  gaagatGCTGGC‑CGCCAGCATCTTCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGACACGCCCGCTcc                                                                                                                                   >  1:73014/6‑142 (MQ=255)                  gaagatGCTGGC‑CGCCAGCATCTTCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGACACGCCCGCTcc                                                                                                                                   <  2:73014/137‑1 (MQ=255)                   gcgGTGCTGGC‑GGCCAGCATCTTCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGtt                                                                                                                           <  5:148006/149‑1 (MQ=255)                     ggTGCTGGC‑GGCCAGCATCTTCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTg                                                                                                                         <  7:101683/149‑1 (MQ=255)                            gc‑ggcCAGCATCTTCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGAt                                                                                                                  <  5:145916/149‑1 (MQ=255)                                    gCATCTTCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtgtg                                                                                                           <  7:169128/149‑1 (MQ=255)                                     cATCTTCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                                          <  1:242675/149‑1 (MQ=255)                                     cATCTTCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                                          <  6:185347/149‑1 (MQ=255)                                         ttCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGc                                                                                                      >  3:226944/1‑149 (MQ=255)                                           ccACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTTCTGATATATGTGTGGGAGCgg                                                                                                    >  6:66529/1‑149 (MQ=38)                                                 tGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGt                                                                                                <  4:226944/147‑1 (MQ=255)                                                   gCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGGTGAGATCTCTGTGGGAGCGGGCGTGccc                                                                                            >  3:302592/1‑149 (MQ=255)                                                   gCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGccc                                                                                            >  3:297886/1‑149 (MQ=255)                                                              aGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGcc                                                                                  >  5:101043/1‑149 (MQ=255)                                                                tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGAtct                                                                                                                >  3:228382/1‑117 (MQ=255)                                                                tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGAtct                                                                                                                <  4:228382/117‑1 (MQ=255)                                                                          cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCCCGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGC‑GAAGCGGc‑cc                                                                    >  6:123284/1‑149 (MQ=255)                                                                          cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAgg                                                                                                                             <  6:84115/94‑1 (MQ=255)                                                                          cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAgg                                                                                                                             >  5:84115/1‑94 (MQ=255)                                                                               ccTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGGCACCCCCCCACCCCCACGTTTTGGGAGAAAGATGTGGGAGCGGGGGGGGCCCCGACGGGCCCC‑GCAGCGGCGCCCaa                                                                 >  2:311921/1‑149 (MQ=255)                                                                               ccTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGC‑GAAGCGGC‑CCCAAtt                                                               >  4:58355/1‑149 (MQ=255)                                                                                         ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGcc                                                                                  <  6:276404/122‑1 (MQ=255)                                                                                            gcgcgGGAGCGTCGTGCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGCCACGCCCGCTCCCACAGGTTTCCTTAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGC‑GAAGCGGC‑CCCAATTGCTTCAACCGTgg                                                  <  5:66529/149‑1 (MQ=38)                                                                                                       tcgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGC‑GAAGCGGC‑CCCAATTGCTTCAACCGTGGTTCTTGCCGag                                       <  7:242943/149‑1 (MQ=255)                                                                                                                 aTCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTTCCACAGGTTTGGTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGC‑GAAGCGGC‑CCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTgg                             >  3:74217/1‑149 (MQ=255)                                                                                                                   ccATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGCGCGTGCAACCGAAGGGTATCGGAAGCGGC‑CCCAATTGCTTCAACAGTGGTTGTTGCCGAGCAGCGCGTgg                             >  1:105417/1‑148 (MQ=25)                                                                                                                          cgggggCGGCTTTCCAGCCCTGCGC‑GGGCACGCCCGCCCCCACCGGTTTGCTTACATCTCTGGGGGCGGGGGGGTGCCCGCGAAGGGCCGC‑GAAGCGGC‑CCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                    <  5:205910/147‑1 (MQ=255)                                                                                                                           cGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCGGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGC‑GAAGCGGC‑CCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCt                      <  5:3737/146‑1 (MQ=255)                                                                                                                           cGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACCCCCCCACCCACAGGGTTGCTGAGATCTCAGTGGGAGCCGGCGTGCCCGCGAAGGGCCCC‑GAACCGGC‑CCCAAATGGATAAACCCGGGTTGATTCCGAGCAGCGCGGGGGTAAGCt                      >  6:3737/1‑146 (MQ=19)                                                                                                                                  gctTTGCAGCCCTTCGC‑GGGCACGCCCGCACCCACAGGTTTGCGGAGATCTCAGTGGGAGCGGGCGGGCCCGCGAAGGGCCGC‑GAAGCGGC‑CCCAAATGGTTCAACCCGGGGTATTGCCGACCCGCGCGTGGGAAACCTCCGTATcccc            >  6:292718/1‑149 (MQ=25)                                                                                                                                      tGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGC‑GAAGCGGC‑CCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa          <  6:101043/147‑1 (MQ=255)                                                                                                                                       gCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGC‑GAAGCGGC‑CCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGt        <  3:58355/148‑1 (MQ=255)                                                                                                                                        cAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGC‑GAAGCGGC‑CCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTAt      >  1:42115/1‑149 (MQ=255)                                                                                                                                         aGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGC‑GAAGCGGC‑CCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                    >  6:265681/1‑134 (MQ=255)                                                                                                                                         aGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGC‑GAAGCGGC‑CCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                    <  5:265681/134‑1 (MQ=255)                                                                                                                                         aGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGC‑GAAGCGGC‑CCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATc     <  2:12153/149‑1 (MQ=255)                                                                                                                                            ccTTCGC‑GGGCACGCCCGCTTCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGC‑GAAGCGGC‑CCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATcccc  <  4:74217/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                CTGGAAGGGCATGCCAGCGCGGTGCTGGC‑GGCCAGCATCTTCCACTTTGGCGAGTA‑CACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑AGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGC‑GAAGCGGC‑CCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCC  >  NC_002947/352923‑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: ‑