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	93.5%	151.6 / ‑5.9	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);  major base G (17/26);  minor base C (2/1);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 5.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCC  >  NC_002947/352923‑353207                                                                                                                                                   |                                                                                                                                              cTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCggg                                                                                                                                             >  2:174967/1‑147 (MQ=255)     aaGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACg                                                                                                                                         <  5:357783/147‑1 (MQ=255)      aGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGc                                                                                                                                        >  4:143482/1‑147 (MQ=255)          cATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACCGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGc                                                                                                                                    <  7:131154/147‑1 (MQ=255)              ccAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTccc                                                                                                                                >  4:33061/1‑147 (MQ=255)                      gTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGttt                                                                                                                        <  3:143482/147‑1 (MQ=255)                      gTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGttt                                                                                                                        <  5:87813/147‑1 (MQ=255)                        gCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGc                                                                                                                      <  1:25011/147‑1 (MQ=255)                          ttgcggcCAGCATCTTCCACTTTGGCGATTACAGAGTGCCGCAGGCCAAGGCCTACATGGCTGCGCGCGGGATCGTCGTTCGCGGATCCATTGAGCGGGTCTGCTTTGCAGCCCTTC‑‑GCGTGCACGCCCGCTCCCACAGTTTTGCt                                                                                                                     <  6:93811/144‑1 (MQ=255)                          tGGCGGCCAGCATCTTCCACTTTGGCGCGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTg                                                                                                                    <  8:333418/147‑1 (MQ=255)                          tGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGGTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGACCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTg                                                                                                                    <  6:381213/147‑1 (MQ=255)                           ggcggcCAGCATTTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                                   >  7:59836/1‑147 (MQ=255)                             cggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTgaga                                                                                                                 <  1:174967/147‑1 (MQ=255)                                         tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGTATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCTGGCACGACCGCTCACACAGGTGTGCTGAGATCTCTGT‑GGGAg                                                                                                    >  5:83575/1‑147 (MQ=39)                                          ccACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGc                                                                                                   >  2:68546/1‑147 (MQ=255)                                           cACTTTGGCGAGTACACAGTGCCGGAGGCCAATGACTACATGGCTTCGCGCGGGAGCGTCGTTCTCTGATCCATTGACCGGGTCTGCTTTGCAGCCCGTC‑‑GCGGGAACGCCCGCTCCCACAGGTTTGCTTAGATCTCTTT‑gggtg                                                                                                    <  6:379867/145‑3 (MQ=255)                                             cTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTGCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCggg                                                                                                <  6:212998/147‑1 (MQ=255)                                                 ggCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCCCGCCCACAGGTTTGCTGAGATCTGTGT‑GGGAGCGGGCggg                                                                                            >  8:279536/1‑145 (MQ=255)                                                     aGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTGCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCcg                                                                                        <  6:64793/147‑1 (MQ=255)                                                             gTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGc                                                                                <  6:382673/147‑1 (MQ=255)                                                             gTGACGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCACCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGc                                                                                >  2:99503/1‑147 (MQ=255)                                                                    aGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAg                                                                         <  2:250327/147‑1 (MQ=255)                                                                         aacgCCTACATGCCTTCGCGCGGGATCGTCGTTCGCTGATCAATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCAACCCCCCTCCCACAGGTTTGGTGAGATCTCTGT‑GGGAGCGGGCGGGCCCGCGACTGGACGGCAAGCggac                                                                    >  4:270042/4‑145 (MQ=25)                                                                                  aTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTGTGT‑GGGAGCGGGCGTGGCCGCGAAGGGCCGCGCAGCGGCCCCAATTGCt                                                           >  1:132374/1‑147 (MQ=255)                                                                                  aTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGc                                                            >  1:127783/1‑146 (MQ=255)                                                                                   tGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCCGGCACCCCCCCTCCCCCACGTTTTCAGAGATAACTTT‑GGGAGCCGGCGTTCCCCCCACGGGCCCCCCAGCCGCCCCAAGTg                                                             >  6:23847/1‑144 (MQ=255)                                                                                     gCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGGGCCGGCACGCCCGCTACCACAGGTTTGCTGAGATATATGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCt                                                           >  1:111623/1‑146 (MQ=38)                                                                                       ttCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAAc                                                      >  1:70944/1‑147 (MQ=255)                                                                                        tCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAAGTTTTCTGAGATCTCTGT‑TGGAGCGGGCGTTCCCGCGAAGGGCCCCGAAGCGGCCCCAATTGATTaaaaa                                                     >  1:332397/1‑145 (MQ=37)                                                                                         cgcgcgGGATCGTCGTTCGCTGATCCATTGCCCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCCCAACCACAGGTTTTCTGCGATATCTGT‑GGGAGCCGGCGAGCCCGCGCAGGGCCCCGCACCCGCCCCCAGTGCTTCAACCg                                                    >  8:176571/1‑147 (MQ=25)                                                                                            gcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTgg                                                 <  3:33061/147‑1 (MQ=255)                                                                                                  tcgtcgTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTg                                           <  7:279536/147‑1 (MQ=255)                                                                                                      cgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGt                                                   >  2:93157/1‑135 (MQ=255)                                                                                                      cgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGt                                                   <  1:93157/135‑1 (MQ=255)                                                                                                      cgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGa                                       <  2:127783/147‑1 (MQ=255)                                                                                                             tGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAgcgc                                >  2:259723/1‑147 (MQ=255)                                                                                                                tCCTTTCACCGGGGCTGCTTTGCTGCCCTTC‑‑GCGGGCACGCCCGCTCCCACCGGTTTGCTGAGATCTCTGT‑GGGTGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                             <  3:194538/147‑1 (MQ=38)                                                                                                                         gggggCTGCTTTGCCGCCCTGC‑‑GCGGGCACGCCCGCGCCCCCAGGTTTGCGGATATCTCTGG‑GGGCGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTc                    <  5:227294/146‑1 (MQ=38)                                                                                                                         cGGGGCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCCCAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCGTCAACCGTGGTTCTTGCCGCGCAGCGCGGGGTAAAGCTc                    >  4:258749/1‑147 (MQ=38)                                                                                                                            ggCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCCCACGTTTGCTGAGAACTCTGT‑GGGAGCGGGCGTGCCCGCGAAAGGCCGCGAAGCGGCCCCAAATGCTTCAAACGTGGTTCATGCCGAGCACCGCGGGGTAAAGCACgg                  >  6:348589/1‑146 (MQ=37)                                                                                                                            ggCTGCTTTGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGt                 <  1:99503/147‑1 (MQ=255)                                                                                                                                  tttGCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATcccc           <  1:68546/147‑1 (MQ=255)                                                                                                                                     gCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCg          <  2:111623/145‑1 (MQ=255)                                                                                                                                     gCAGCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCg          <  3:258749/145‑1 (MQ=255)                                                                                                                                        gCCCTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGGTTCAACCGGGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTAt     >  5:207232/1‑147 (MQ=255)                                                                                                                                          ccTTC‑‑GCGGGCACGCCCTGTCCCACAGGTTGCCTGAGATCTCTGG‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATcc   <  4:224822/147‑1 (MQ=38)                                                                                                                                           cTTC‑‑GCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATccc  <  1:259723/147‑1 (MQ=255)                                                                                                                                           cTTC‑‑GCCGGCACGCCCGCGCCCAACGTTTTCCTGACATC‑CAGTGGGGCGGGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCATCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATccc  <  6:223616/147‑1 (MQ=25)                                                                                                                                                   |                                                                                                                                              CTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTC‑‑GCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGT‑GGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCC  >  NC_002947/352923‑353207

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