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%	173.5 / 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 (23/29);  total (23/29)

TCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCA  >  NC_002947/352922‑353209                                                                                                                                                  |                                                                                                                                                 tCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCggg                                                                                                                                                >  6:169267/1‑148 (MQ=255)    ggAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACg                                                                                                                                            >  4:114357/1‑149 (MQ=255)     gAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGc                                                                                                                                           <  4:32640/149‑1 (MQ=255)      aaGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGcc                                                                                                                                          >  5:18832/1‑149 (MQ=255)           cATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTc                                                                                                                                     >  1:38150/1‑149 (MQ=255)           cATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTc                                                                                                                                     >  5:85330/1‑149 (MQ=255)            aTGCCAGCGCTGTGCTGGCGGAGAGCATCTTCGACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATTGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTTCTTTGCATCCCTTCGCGGGCACGACCGCTcc                                                                                                                                    <  6:149153/149‑1 (MQ=255)            aTGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTcc                                                                                                                                    <  3:360716/149‑1 (MQ=255)             tGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTcc                                                                                                                                    >  3:224949/1‑148 (MQ=255)                   cgcgGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGt                                                                                                                             <  4:224949/149‑1 (MQ=255)                   ccaggTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGTTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                               <  5:34989/144‑1 (MQ=255)                        tGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑Ct                                                                                                                       <  5:169267/149‑1 (MQ=255)                               ggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑Atct                                                                                                               <  6:36597/149‑1 (MQ=255)                               ggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑Atct                                                                                                               >  5:36597/1‑149 (MQ=255)                                 ccAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGTG‑ATctct                                                                                                             <  2:38150/149‑1 (MQ=255)                                 ccAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATctct                                                                                                             <  6:85330/149‑1 (MQ=255)                                  cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCtg                                                                                                            >  6:157647/1‑149 (MQ=255)                                         ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCGCACCCCCAGGTTTG‑CTGAG‑ATCACAGTTGGAGc                                                                                                     >  2:282454/1‑149 (MQ=39)                                          tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCg                                                                                                    >  2:109437/1‑149 (MQ=255)                                           ccACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGa                                                                                                       <  8:64218/145‑1 (MQ=255)                                           ccACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCgg                                                                                                   >  7:129215/1‑149 (MQ=255)                                            cACTTTGGCGAGTACACAGTGCCTGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACTCCCGCTCCCACTGTTTTT‑CTGAG‑ATATCTGTTGGAGCgg                                                                                                   <  2:317269/148‑1 (MQ=38)                                             aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGc                                                                                                 <  3:114357/149‑1 (MQ=255)                                                    cGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTgg                                                                                                         >  4:15350/1‑134 (MQ=255)                                                    cGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTgg                                                                                                         <  3:15350/134‑1 (MQ=255)                                                         acacAGTGCCGTTGGCAATGGCCTACTTGGCTTCGCGCGGGAGGGGCGTTCGCTTATCAATTGACGGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCGCCCACAGTTTTG‑CTGAT‑ATCTCTGTGGTAGGGGGGGTGCCCGCGaa                                                                                      <  2:144552/148‑1 (MQ=21)                                                               tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCtgt                                                                                                           <  6:71480/121‑1 (MQ=255)                                                               tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCtgt                                                                                                           >  5:71480/1‑121 (MQ=255)                                                                      ggCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGGTCGCTGATCCAATGACCGGGGCTGCTTTGCAGCCCGTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑AGCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCgg                                                                        <  5:360454/149‑1 (MQ=38)                                                                             gCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGttt                                                                                                                           >  4:141841/1‑93 (MQ=255)                                                                             gCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGttt                                                                                                                           <  3:141841/93‑1 (MQ=255)                                                                                  cATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGACTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTg                                                               <  1:109437/147‑1 (MQ=255)                                                                                       cTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAAcc                                                       <  8:70093/149‑1 (MQ=255)                                                                                           gcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTg                                                    <  7:357875/148‑1 (MQ=255)                                                                                             gcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGtt                                                 <  5:157647/149‑1 (MQ=255)                                                                                                 gATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGc                                                                                             >  5:225922/1‑101 (MQ=255)                                                                                                 gATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGc                                                                                             <  6:225922/101‑1 (MQ=255)                                                                                                      tcgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTcc                                                                                                                                    <  5:25134/59‑1 (MQ=255)                                                                                                      tcgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTcc                                                                                                                                    >  6:25134/1‑59 (MQ=255)                                                                                                            gCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAgcgc                                  >  5:216198/1‑149 (MQ=255)                                                                                                             gtGATCCATTGAGCGGGGCTGCTTTGCAGCCGTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGcgcg                                 <  1:282454/148‑1 (MQ=255)                                                                                                             cTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGcgcg                                 <  8:129215/149‑1 (MQ=255)                                                                                                              tGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGt                                >  7:129352/1‑149 (MQ=255)                                                                                                                aTCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTgg                              >  7:187687/1‑149 (MQ=255)                                                                                                                aTCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCACCCACAGGTTTG‑CGGAG‑ATATCAGTGGGAGCGGGCCTGCCCGCCAAGGGCCGCGAAGCGGCCCCAAATGCATAAAACGAGGTTCTTGCCGAGCAACGCCTgg                              >  8:243048/1‑149 (MQ=25)                                                                                                                    aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAg                                                                                     <  6:83192/90‑1 (MQ=255)                                                                                                                    aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAg                                                                                     >  5:83192/1‑90 (MQ=255)                                                                                                                         ccGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCg                     <  6:18832/149‑1 (MQ=255)                                                                                                                                   tttGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcccc             >  2:124658/1‑147 (MQ=255)                                                                                                                                     cgcggCCCTTCGCGGGCACGCCCGCTCCCACAGAGAT‑‑CTCAGCAAACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa           <  6:397490/143‑1 (MQ=14)                                                                                                                                          cccTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTg                                             <  1:331732/108‑1 (MQ=255)                                                                                                                                            cTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCcca  >  6:170848/1‑149 (MQ=255)                                                                                                                                            cTTCGCGGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCcca  >  5:318631/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                 TCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTG‑CTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCA  >  NC_002947/352922‑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: ‑