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%	165.8 / NA	52	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 (28/24);  total (28/24)

TCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCC  >  NC_002947/352922‑353207                                                                                                                                                  |                                                                                                                                             tCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTAGATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCggg                                                                                                                                            <  4:134937/148‑1 (MQ=255)   tGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCAc                                                                                                                                         >  5:405999/1‑149 (MQ=255)            aTGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACACTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTcc                                                                                                                                <  5:70124/149‑1 (MQ=255)                cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCcaca                                                                                                                            >  4:260611/1‑149 (MQ=255)                    gcgGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGtt                                                                                                                        <  4:83967/149‑1 (MQ=255)                      ggTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTg                                                                                                                      >  2:23738/1‑149 (MQ=255)                        tGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                                    <  7:175178/149‑1 (MQ=255)                         gCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACAGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                                    >  8:313647/1‑148 (MQ=255)                              cggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAtc                                                                                                              <  6:34049/149‑1 (MQ=255)                                   aGCATCTTCCACTTTGGCGAGTACACAGTGGCGGATGGCAAGGCCTACATGGCTTCGCGCGGGTTCGTGGTTCGCTGATCCATTGACCGGGGCTGGTTTGCAGCCGTTCGCGGGCTCGCCCGCTCCCACAGGTTTGCTGAGATCTCtgt                                                                                                         <  2:94964/149‑1 (MQ=255)                                     cATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTgg                                                                                                       >  4:25037/1‑149 (MQ=255)                                     cATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTCCCCCAGGTTTGCTGAGAtatcgttgg                                                                                                       >  3:354527/1‑144 (MQ=39)                                     cATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCGCACCCACAGCGTTGCTGAGATCTCTGTgg                                                                                                       >  8:233374/1‑149 (MQ=39)                                            cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCgg                                                                                                 <  3:260611/148‑1 (MQ=255)                                             aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCGCCCACAGGGTTGGTGAGATCTCTGTGGGAGCGGGc                                                                                               >  5:40405/1‑149 (MQ=39)                                                 tGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGt                                                                                             <  3:25037/147‑1 (MQ=255)                                                   gCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGccc                                                                                         >  8:9420/1‑149 (MQ=255)                                                          caccGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTGCGCGGGCACGCCCGCTCCCACAGGTTTGCTGACAGCTCTGTGGGAGCGGGCGTGCCCGCGAAg                                                                                   <  5:316379/148‑1 (MQ=39)                                                           acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTgaga                                                                                                                >  2:9700/1‑118 (MQ=255)                                                           acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTgaga                                                                                                                <  1:9700/118‑1 (MQ=255)                                                           acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                 >  6:182409/1‑149 (MQ=255)                                                               tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATATCTGTGGGAGCGGGCGTGGCCGCGAAGGGCcgc                                                                             >  3:185429/1‑149 (MQ=255)                                                                   ggAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAg                                                                         >  8:385012/1‑149 (MQ=255)                                                                       gCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGTTGCTTTCCAGCCCTGCGCGGGCACGCCCGCTCCCACAGCTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGc                                                                     <  1:266268/149‑1 (MQ=39)                                                                        ccATGGCCTACTTGGCTTCGCGGGGGAGCGTCTTTCGCTGATCCTTTGACCGGGGCTGTTTTGCTGCCGTGCGCGGGCACGCCCGCTCCCACTGGTTTGCTTAGATCTCTGGGGTAGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGcc                                                                    <  4:185429/149‑1 (MQ=255)                                                                         cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCACGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGACGGGCCGCGAAGAGGgcc                                                                   >  8:171140/1‑149 (MQ=38)                                                                         cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTGCCACAGGTTTGCTGAGATGTCTGTGGGAGCGGGCGTGCCCGCGCAGGGCCGCGAAGCggg                                                                     >  7:257993/1‑146 (MQ=39)                                                                                   aTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCtt                                                          >  5:156289/1‑148 (MQ=255)                                                                                        ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCTCCCACAAGTTTTCTGAGAACTCTTTGGGAGCGGGCGTGGCCGCGCAGGGCCCCGCAGCCGCCCCAAATGCTTCAAccc                                                    >  4:395201/1‑148 (MQ=37)                                                                                          cgcgcgGGTTCGTCGTGCGCTGATCCATGGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTg                                                  <  4:354527/149‑1 (MQ=255)                                                                                          cgcgcgGGCTCGTCGGGCGCTGATCCATTGACCGGTGCTGCTTTGCAGCCCGTCGCGGGCCCGCCCGCTCCCACCGGTTTGCTGAGATCTCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTg                                                  <  7:357766/149‑1 (MQ=38)                                                                                          cgcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGTAGCGGCCCCAATTGCTTCAACCg                                                    <  6:405999/147‑1 (MQ=255)                                                                                                            gCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCCCTCCCCCAGGGTTGGTGAGATATCTGTGGGAGCGGGCGTGCCCGCGAACGGCCGCGAAGCGGCCCCAATTGCTTCAACCGGGGTTCTTGCCGAGCCgcgc                                >  5:16331/1‑149 (MQ=38)                                                                                                            gCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCCCGAAGCGGCCCCAATTGCTTCAACCGTGGTTATTTCCGAGCAtcgc                                >  8:382336/1‑149 (MQ=38)                                                                                                             cTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGcgcg                               <  7:233374/149‑1 (MQ=255)                                                                                                              tGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCACCCACAAGTTTGCTGAGATCTAAGTGGGAGCGGGCGTGCCCGCGAACGGCCCCCAAGCGGCCCCAATAGCATCAAACGTGGTTCTTGCCGAGCAGCGCGt                              >  4:94241/1‑149 (MQ=37)                                                                                                               gATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagcagc                                  <  4:318543/144‑1 (MQ=255)                                                                                                               gATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATATCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCCACCCTGGTTCTTGCCGAGCAGcgcggt                             >  3:318543/1‑147 (MQ=255)                                                                                                                 tCCTTTGCCCGGGGCGGTTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACTGGTTTGCTGATATCTCTGGGGGAGGGGGGGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTgg                            <  6:16331/148‑1 (MQ=38)                                                                                                                    aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaaa                        <  7:9420/149‑1 (MQ=255)                                                                                                                    aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATATCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTa                          >  4:346014/1‑147 (MQ=255)                                                                                                                    aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGCGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGACCGGCCGCGCAGCGGCCCCAATGGCTTCAACCGGGGTTCTTGCCGCGCAGCGCGTGGTaaa                        >  3:274425/1‑149 (MQ=38)                                                                                                                     ttGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCTCCCACAGGTTTGCTGAAAACTAAGTGGGAACGCGCGAGCCCGCCAAGGGCCGCGAACCGGCCCCAATTGCTTCAACCCTGGTTCTTTCCGAGCAGCGCCTGGTa                          >  6:116449/1‑146 (MQ=25)                                                                                                                      tgcccGGGGCTGCTTTGCAGCCCTGCGCGGGCACGCCGGCTCCCACTGGTTTGGTTCTAGCTGTGTGGGTGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAg                       <  6:40405/145‑1 (MQ=37)                                                                                                                        aCCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTc                    >  2:132468/1‑149 (MQ=255)                                                                                                                         ccGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                   <  6:156289/149‑1 (MQ=255)                                                                                                                         ccGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                   <  7:385012/149‑1 (MQ=255)                                                                                                                         ccGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCAGAGATCTCTGTGGGAGCGGGCGTCCCCCCGAAGGGCCGCGAAGCGGCCCCAATTGCTTAAACCGTGGTTCATGACGAGAAGCGCGGGGTAAAtctcc                   >  6:422714/1‑148 (MQ=37)                                                                                                                         ccGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCGCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                   <  7:171140/149‑1 (MQ=255)                                                                                                                          cGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTTCGATCTCTGTGGTAGCGGGCTTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                   <  5:422714/148‑1 (MQ=38)                                                                                                                               ctgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcc             >  1:44285/1‑149 (MQ=255)                                                                                                                                tgctTTGCAGCCCTTCGCGGGCACGCCCGCAACCACAGGTTTGGAGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCCCCAAGCGGCCCCAAATTCATAAACCCTGGTTCATGCCGAGCACCGCCAGGTAAACCACGGTATacc            >  2:127962/1‑149 (MQ=21)                                                                                                                                          cccTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATccc  <  5:182409/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                             TCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCC  >  NC_002947/352922‑353207

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑