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%	148.2 / 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 (21/24);  total (21/24)

GAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑AGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTAT  >  NC_002947/352926‑353204                                                                                                                                               |                                                                                                                                          gAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGc                                                                                                                                    <  8:92381/149‑1 (MQ=255)     tggccgCCAGCACCGCGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGc                                                                                                                                >  7:101163/6‑149 (MQ=39)     ggCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCg                                                                                                                                 >  4:142444/1‑148 (MQ=255)        aTGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGcccccacc                                                                                                                             >  7:148870/1‑144 (MQ=255)                  ggTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTg                                                                                                                   <  2:257227/149‑1 (MQ=255)                    tGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                                 <  3:142444/149‑1 (MQ=255)                       tGGCGGCCAGCATCTTCCACTTTGGCTATGACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTgag                                                                                                              <  8:354917/149‑1 (MQ=255)                         gcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTgag                                                                                                              <  7:332597/147‑1 (MQ=255)                              cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtg                                                                                                       >  6:311815/1‑149 (MQ=255)                                     ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGc                                                                                                >  7:372266/1‑149 (MQ=255)                                      tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTcc                                                                                                                             >  5:40697/1‑119 (MQ=255)                                      tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTcc                                                                                                                             <  6:40697/119‑1 (MQ=255)                                       ccACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGACCGCTCCCAAAGGTTTGCATCGATCTCTGTTGGCGCgg                                                                                              >  4:219607/1‑149 (MQ=38)                                             tGGCGAGTACTCAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCAGTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGGTTGCTGAGATCTCTGTGGGAGCGGGCGt                                                                                          <  3:295026/147‑1 (MQ=255)                                             tGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGt                                                                                          <  8:101163/147‑1 (MQ=255)                                             tGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGt                                                                                          <  3:159362/147‑1 (MQ=255)                                               gCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCCCGCCCGCTCCCACAGGt                                                                                                                      >  3:384710/1‑117 (MQ=255)                                               gCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGt                                                                                                                      <  6:247896/117‑1 (MQ=255)                                               gCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGt                                                                                                                      >  5:247896/1‑117 (MQ=255)                                               gCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGt                                                                                                                      <  4:384710/117‑1 (MQ=255)                                                 gAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCCCGCCCGCGCCCACAGGGTTGCTGAGATCGCTGTGGGAGCGGGCGTGCCcgc                                                                                    >  2:333177/1‑149 (MQ=39)                                                         aGTGCCGTAGGCCAAGGCCTACAGGGCTTCGCGCGGCATCGGCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGcc                                                                            <  8:24717/149‑1 (MQ=38)                                                                gAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGc                                                                     <  5:311815/149‑1 (MQ=255)                                                                gAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGc                                                                     >  5:344052/1‑149 (MQ=255)                                                                   gCCACGGCCTACATGGCTGCGCGGGGGATCGTCGTTCGCTGAGCCATTCACGGGGGCTGCTTTCCTGCCCTGCGC‑GGGCACGCCCGCTCCCCCCGGTTTGGTTACATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGc                                                                  <  2:15659/149‑1 (MQ=25)                                                                     cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGGCCCCCCCCCACCCACAGGTTTGGGTAGAGATGTGGGGGGGCGGGGGGGCCCCCCAAGGGCCCCGACCCGGcc                                                                 >  2:377392/1‑149 (MQ=255)                                                                     cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGccc                                                                >  8:310373/1‑149 (MQ=255)                                                                         gCCTACATGGCTTCGCGCGGGAGCGTCGTGCGCTGATCCATTGCCCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAAt                                                            <  1:333177/149‑1 (MQ=255)                                                                            tACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGc                                                         >  1:219289/1‑149 (MQ=255)                                                                             aCATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCCCTCCCCCAGGTTTGCGTAGATGTCTGTGGGCGCGGGCGTGCCCGCCAAGGGCCGCGAAACGGCCCCCATTTCt                                                        >  7:340411/1‑149 (MQ=37)                                                                              cATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGc                                                         <  8:372266/147‑1 (MQ=255)                                                                               aTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCcaca                                                                                                                         >  8:169791/1‑82 (MQ=255)                                                                               aTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCcaca                                                                                                                         <  7:169791/82‑1 (MQ=255)                                                                                        cgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGt                                             <  7:310373/149‑1 (MQ=255)                                                                                                    gtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGag                                   >  5:89133/1‑147 (MQ=255)                                                                                                        ggggATCCTTTTACCGGGGCGGCTTGGCAGCCCTTGGC‑GGGCCCGCCCGCTCCCACAGGTTGGCGGCGTTCTCGGTGGGAGCGGGGGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAgcgc                             <  3:219607/146‑1 (MQ=37)                                                                                                               cATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                                 >  3:156126/1‑58 (MQ=255)                                                                                                               cATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                                 <  4:156126/58‑1 (MQ=255)                                                                                                                 ttGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGCTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCCATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAg                    >  3:346157/1‑149 (MQ=255)                                                                                                                 ttGACCGGGGCTGCTTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATATCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTTCTTAAACCGTGGTTCTTGACGAGCATCGCGTGGTAAAg                    >  3:240389/1‑149 (MQ=37)                                                                                                                              ctTTGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGTTTTGCTGAGTTCTCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcccc        <  8:340411/148‑1 (MQ=255)                                                                                                                               tgtGCCGCCCTGCGC‑GGGCCCGCCCGCCCCCACTGGTGTGCTCACTTCGGTGGGGGCGGGGGGGGCCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa      <  7:319864/147‑1 (MQ=255)                                                                                                                                 tGCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa      <  5:72167/147‑1 (MQ=255)                                                                                                                                  gCAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGt    <  6:89133/148‑1 (MQ=255)                                                                                                                                   cAGCCCTTCGC‑GGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTAt  >  1:110147/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                          GAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGC‑AGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTAT  >  NC_002947/352926‑353204

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