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	98.0%	158.7 / ‑4.2	49	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);  major base G (23/25);  minor base T (0/1);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

CTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCA  >  NC_002947/352923‑353209                                                                                                                                                 |                                                                                                                                                cTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCa                                                                                                                                             >  5:75168/1‑149 (MQ=255)  tGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCAc                                                                                                                                            >  5:324244/1‑149 (MQ=255)           ccgCCAGCACCGCGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTcc                                                                                                                                   >  1:262922/3‑149 (MQ=255)              ccAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCcac                                                                                                                                <  2:209622/149‑1 (MQ=255)              ccAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCca                                                                                                                                 <  7:75671/148‑1 (MQ=255)                aGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGcc                                                                                                                                         >  4:82512/1‑138 (MQ=255)                aGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGcc                                                                                                                                         <  3:82512/138‑1 (MQ=255)                aGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                              >  6:96772/1‑149 (MQ=255)                aGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                              <  2:346371/149‑1 (MQ=255)                      tggCTGTCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATAGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGc                                                                                                                        <  3:15670/147‑1 (MQ=255)                         cTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGttt                                                                                                                          >  4:358152/1‑144 (MQ=255)                         cTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGttt                                                                                                                          <  3:358152/144‑1 (MQ=255)                         cTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                                     <  8:5109/149‑1 (MQ=255)                          tGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                                       >  1:74976/1‑146 (MQ=255)                              ggcCAGCATCTGCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑Atct                                                                                                               <  3:226632/149‑1 (MQ=255)                                 cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCtg                                                                                                            >  7:177798/1‑149 (MQ=255)                                  aGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCGCGCCCACCGGTTTGGTGAG‑GTCTgtgt                                                                                                           >  1:156237/1‑149 (MQ=39)                                        ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGc                                                                                                     >  4:245399/1‑149 (MQ=255)                                          ccACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCgg                                                                                                   >  7:328298/1‑149 (MQ=255)                                          ccACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTGCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTggta                                                                                                       <  2:238311/145‑3 (MQ=255)                                           cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑Atct                                                                                                               >  6:3031/1‑136 (MQ=255)                                           cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑Atct                                                                                                               <  5:3031/136‑1 (MQ=255)                                           cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCgg                                                                                                   <  2:74976/148‑1 (MQ=255)                                            aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGc                                                                                                 <  8:328298/149‑1 (MQ=255)                                                        acacAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCTGGCACGCCCGCTCCCACAGAGAT‑CTCAGCAAACCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                      <  6:359416/148‑1 (MQ=21)                                                                 cGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTgg                                                                                                         >  6:270388/1‑120 (MQ=255)                                                                 cGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTgg                                                                                                         <  5:270388/120‑1 (MQ=255)                                                                 cGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCTGCGAAGGGCCGCGaa                                                                            <  5:96772/149‑1 (MQ=255)                                                                        cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCt                                                                                                                                     <  8:1503/86‑1 (MQ=255)                                                                        cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCt                                                                                                                                     >  7:1503/1‑86 (MQ=255)                                                                              cTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTg                                                               >  7:76315/1‑149 (MQ=255)                                                                                      cTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGTAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAAcc                                                       <  1:85398/149‑1 (MQ=255)                                                                                       ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCCCAACCACAGGGTTTCTGAG‑ATATATGTGGGAGCGGGCGGGCCCGCCAAGGGCCCCGCAGCGGCCCCCAATGCTTCAAccc                                                      >  2:123982/1‑148 (MQ=255)                                                                                          gcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCCCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTgg                                                   <  2:297134/149‑1 (MQ=255)                                                                                               ggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCtt                                              <  2:262922/149‑1 (MQ=255)                                                                                                    gtggtTCGCTTATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGGGGGAGCGGGCGGGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGa                                         <  1:123982/146‑1 (MQ=38)                                                                                                     tcgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCCCAACCCCAGGGTTGGTGAG‑ATATCTTTTGGTCCCGGGGTGCCCCCGCAGGGCCCCGCAGCGGGCCCAAATGCGTAAAACCTGGTTATTTCCGag                                        >  8:308967/1‑149 (MQ=255)                                                                                                          cGCTGATCCATTGTCCGGGGGTGCTTTGCAGCCGTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGCTCTTGCCGagcagc                                    <  7:308967/148‑1 (MQ=38)                                                                                                               aTCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTgg                              <  3:245399/149‑1 (MQ=255)                                                                                                               aTCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                               <  6:75168/148‑1 (MQ=255)                                                                                                                    ttGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCAAAGGGTTGGTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAg                         >  1:80575/1‑149 (MQ=255)                                                                                                                          ggggCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCGACCCCAGGGTTGCGGAG‑ATCTAAGTGGGAGCGGGCGGGCCCGCCCAGGGCCCCGAACCGGCCCCCATTTGATGAACCGAGGTTCGTGCCGAGCCGCCCGAGGGTAAGcagcgg                   >  6:195761/1‑144 (MQ=255)                                                                                                                             aggatgttGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATc                <  1:80997/143‑1 (MQ=255)                                                                                                                              ctgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAAATGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcc               >  6:141947/1‑149 (MQ=255)                                                                                                                               tgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagca                                      <  5:178765/125‑1 (MQ=255)                                                                                                                               tgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagca                                      >  6:178765/1‑125 (MQ=255)                                                                                                                                 ctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCg            <  8:76315/149‑1 (MQ=255)                                                                                                                                  tttGCAGCCCTTCGCGGGCACCCCCGCAACCACAGGTTTGCTGAG‑ATCTCTGTGGGAACGGGCCTGCCCCCGAAGGGCCGCGAAGCGGCCCCAATAGCTTAAACCCTGGTTCATGCCGAGCAGCGCGTGGTAAAGCTAGGTATcccc             >  8:281359/1‑147 (MQ=37)                                                                                                                                           cTTCGCGGGCACCCCCCCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGGGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGGGGTTGTTGGCGAGGAGCGCGTGGTAAAGCTCGGTATACCCGAGTATccccc  >  4:17460/1‑148 (MQ=38)                                                                                                                                                 |                                                                                                                                                CTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCA  >  NC_002947/352923‑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: ‑