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%	157.9 / NA	47	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 (24/23);  total (24/23)

TTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCAGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACC  >  NC_002947/352921‑353211                                                                                                                                                    |                                                                                                                                                   ttCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCggg                                                                                                                                                  <  4:137460/149‑1 (MQ=255)      gAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGc                                                                                                                                             >  6:12068/1‑149 (MQ=255)                ccAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCca                                                                                                                                   >  6:371950/1‑148 (MQ=255)                  aGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑g                                                                                                                               <  7:160765/149‑1 (MQ=255)                      cgGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGttt                                                                                                                           <  3:294958/149‑1 (MQ=255)                      cgGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGttt                                                                                                                           <  4:366764/149‑1 (MQ=255)                         tGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCt                                                                                                                        >  8:77684/1‑149 (MQ=255)                           cTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCGTGCCCGcgaag                                                                                                                                     >  8:82/1‑132 (MQ=255)                           cTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCGTGCCCGcgaag                                                                                                                                     <  7:82/136‑5 (MQ=255)                           cTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTga                                                                                                                      <  7:77684/149‑1 (MQ=255)                            tGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCCAGCCCCC‑ACCCCCACGGTTTgggc                                                                                                                       >  4:189541/1‑145 (MQ=255)                              gcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTgaga                                                                                                                    <  3:86617/148‑1 (MQ=255)                                ggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGATTCGGGCGGGTTCGTCGTTCGCTGATCCATTGACGGGGGCTGCTTTGCAGC‑CCTTCGCGTGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGAtct                                                                                                                 <  1:233167/149‑1 (MQ=38)                                    aGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCtgt                                                                                                             <  7:130213/149‑1 (MQ=255)                                    aGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCtgt                                                                                                             <  1:229708/149‑1 (MQ=255)                                    aGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGGTTGCTGAGATCTCtgt                                                                                                             >  2:241562/1‑149 (MQ=255)                                     gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGTCTTCTCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCtgtg                                                                                                            <  3:140632/149‑1 (MQ=255)                                       aTCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTggg                                                                                                          <  8:84434/149‑1 (MQ=255)                                          ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTgg                                                                                                           <  7:61982/145‑1 (MQ=255)                                          ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTgg                                                                                                           >  8:61982/1‑145 (MQ=255)                                            ccACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAg                                                                                                        <  3:311369/146‑1 (MQ=255)                                               cTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCg                                                                                                  >  3:268074/1‑149 (MQ=255)                                                      gAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCCCGCCCCCGACCCCAA‑GGTTTG‑TGAGATATCTGTTGGCGCCGGCGGGccccc                                                                                           >  6:220440/1‑147 (MQ=255)                                                         tACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCTCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCcgcg                                                                                          >  6:368251/1‑147 (MQ=255)                                                          acacAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                        <  8:3367/148‑1 (MQ=255)                                                           cacaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTAGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAgg                                                                                      >  1:85396/1‑149 (MQ=255)                                                            acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                     >  6:150475/1‑149 (MQ=255)                                                                  ccGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGa                                                                                                         <  1:79030/123‑1 (MQ=255)                                                                  ccGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑CGTTTGCTGAGATCTCTGTGGGa                                                                                                         >  2:79030/1‑123 (MQ=255)                                                                    ggAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAg                                                                             >  6:258838/1‑149 (MQ=255)                                                                          cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGccc                                                                       >  4:283035/1‑149 (MQ=255)                                                                           aaGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGccc                                                                       >  3:207683/1‑148 (MQ=255)                                                                                            gcgcgGGATCGTCGTTCGCTGATCCATTCACCGGGGCTGCTTTGCAGC‑CCTGCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTgg                                                     <  5:220440/149‑1 (MQ=255)                                                                                                             gCTGATCCATTGACCGGGGCTGCTTTGCAGCGCCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                                                   >  1:39612/1‑87 (MQ=25)                                                                                                             gCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGc                                                                                                   <  2:39612/86‑1 (MQ=255)                                                                                                                 aTCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTgg                                <  3:283035/149‑1 (MQ=255)                                                                                                                 aTCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCAACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTgg                                >  6:355633/1‑149 (MQ=255)                                                                                                                   ccATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGt                               >  2:13101/1‑148 (MQ=255)                                                                                                                     aTTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTgg                                                     <  2:221829/124‑1 (MQ=255)                                                                                                                     aTTGACCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTgg                                                     >  1:221829/1‑124 (MQ=255)                                                                                                                         aCCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCt                                                               <  3:216119/110‑1 (MQ=255)                                                                                                                         aCCGGGGCTGCTTTGCAGC‑CCTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCt                                                               >  4:216119/1‑110 (MQ=255)                                                                                                                                  gctTTGCTGC‑CTTCCGCGGGCACGCCCGC‑TCCCACT‑GGTTTGCTGAGAGCTCTGTGGGTGGGGGGGGGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcccc               <  6:42600/149‑1 (MQ=255)                                                                                                                                      tGCAGC‑CCTTCGCGGGCACGCCCGC‑GCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa             <  5:12068/147‑1 (MQ=255)                                                                                                                                           c‑ccTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa             >  6:279180/1‑142 (MQ=255)                                                                                                                                           c‑ccTTCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa             <  5:279180/142‑1 (MQ=255)                                                                                                                                               ttCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCccac   >  6:322813/1‑149 (MQ=255)                                                                                                                                                tCGCGGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGCAGCGGCCCGAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTTAAGCTCGGTATCCCCGAGTATCcccccc  >  7:176860/1‑149 (MQ=38)                                                                                                                                                    |                                                                                                                                                   TTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGC‑CCTTCGCAGGCACGCCCGC‑TCCCACA‑GGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACC  >  NC_002947/352921‑353211

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