Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,655,626:1	+A	10.6%	intergenic (‑81/+71)	epd ← / ← tktA	D‑erythrose 4‑phosphate dehydrogenase/transketolase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,655,626	1	.	A	10.6%	170.9 / 14.1	47	intergenic (‑81/+71)	epd/tktA	D‑erythrose 4‑phosphate dehydrogenase/transketolase A
Reads supporting (aligned to +/- strand):  ref base . (17/25);  new base A (3/2);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCC  >  NC_002947/5655480‑5655769                                                                                                                                                     |                                                                                                                                                 gcgcAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGCTGTGGCATTCGCGGGCAAGCCCGCTCCCACAg‑g                                                                                                                                                 <  1:206124/149‑1 (MQ=255)   gcAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCTTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGtt                                                                                                                                              <  6:88706/149‑1 (MQ=255)   gcAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGtt                                                                                                                                              >  5:56913/1‑149 (MQ=255)      ggACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCa                                                                                                                                           <  1:275797/149‑1 (MQ=255)        acacAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATg                                                                                                                                         <  4:26887/149‑1 (MQ=255)            aGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAAc                                                                                                                                     >  5:74084/1‑149 (MQ=255)                   gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCtgtg                                                                                                                              <  4:3344/149‑1 (MQ=255)                    aTACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTgg                                                                                                                             >  6:140039/1‑149 (MQ=255)                     tACGGCCATAACCGTTGAGTGCAACTTTGTAGGTACGGTGGGTAATTTGGGGCTCTCGTAAGCTAACTATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACC‑Ctgtg                                                                                                                              <  1:265016/147‑1 (MQ=18)                          ccATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCgg                                                                                                                       <  5:72496/149‑1 (MQ=255)                            aTAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGc                                                                                                                     <  5:140039/149‑1 (MQ=255)                                ccGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGc                                                                                                                     >  8:136174/1‑145 (MQ=255)                                ccGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGc                                                                                                                     <  7:136174/145‑1 (MQ=255)                                     gAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCcgc                                                                                                            >  7:181132/1‑148 (MQ=255)                                      aGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCCGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCcgcg                                                                                                           <  2:184941/148‑1 (MQ=255)                                      aGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCcgcg                                                                                                           <  1:209841/148‑1 (MQ=255)                                            aCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAg                                                                                                    <  7:7402/149‑1 (MQ=255)                                                 gTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGa                                                                                                >  2:219882/1‑148 (MQ=255)                                                  tAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAC‑GGTTTACTGCAACCTGTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGcca                                                                                              >  2:3498/1‑149 (MQ=39)                                                    gggACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTGCCACAC‑GGGTTGCTGCAACCTCTGTGGGATCCGGCTTG‑CCCGCGTATTAAGGGGGCAcc                                                                                            >  5:24305/1‑149 (MQ=255)                                                     ggACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGTCAccg                                                                                           >  5:90446/1‑149 (MQ=255)                                                     ggACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCATCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACAccg                                                                                           >  5:107645/1‑149 (MQ=255)                                                       aCGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGaa                                                                                                     >  6:153496/1‑137 (MQ=255)                                                       aCGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGaa                                                                                                     <  5:153496/137‑1 (MQ=255)                                                             gggCATTCGGGGATCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑GGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCCGCGAATGAAGGCGACACCGCCCGCAg                                                                                    <  1:3498/149‑1 (MQ=255)                                                               gCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACAccgcccg                                                                                       >  2:79399/1‑143 (MQ=255)                                                               gCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATc                                                                                 <  1:79399/149‑1 (MQ=255)                                                                  ttCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCtt                                                                               >  4:209444/1‑148 (MQ=255)                                                                        gCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCa                                                                        <  6:291088/149‑1 (MQ=255)                                                                         cacACGACACAGTTCAATTGTCGTCAGCCTGCGTAATTTGCGG‑TGGGGCATTCCCGGGCAAGCCCGCTCCCACAG‑GGTTTCATCCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAg                                                                       <  1:166626/147‑1 (MQ=38)                                                                           cACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCCACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGCCACCGCCCGCAGATATTTTAATCCGTc                                                                     >  8:270294/1‑149 (MQ=38)                                                                            aCGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGt                                                                      <  6:74084/147‑1 (MQ=255)                                                                                gaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCca                                                                <  6:56913/149‑1 (MQ=255)                                                                                   aGTTCAATTGTCTTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCcagca                                                             <  6:107645/149‑1 (MQ=255)                                                                                   aGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGTCACCGCCCGCAGATCTTTAAATCAGTCTTCcagca                                                             <  6:90446/149‑1 (MQ=255)                                                                                                         ggAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGTCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGAt                                       <  5:155792/149‑1 (MQ=255)                                                                                                           attttGCGG‑GGGGCCATTCGCGGCAACCCCCGCTCCCACTG‑GGTTTCATCCAACCTGTGTGGTAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATg                                      <  5:48455/146‑1 (MQ=255)                                                                                                             tttGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTtc                                   >  6:82947/1‑149 (MQ=255)                                                                                                                         gCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGaa                        <  1:219882/149‑1 (MQ=255)                                                                                                                               gcgGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATg                                      <  7:170179/129‑1 (MQ=255)                                                                                                                               gcgGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATg                                      >  8:170179/1‑129 (MQ=255)                                                                                                                                  ggCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAActc               <  5:82947/149‑1 (MQ=255)                                                                                                                                   gCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAActct              >  8:133778/1‑149 (MQ=255)                                                                                                                                     aaGCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTc            >  3:139279/1‑149 (MQ=37)                                                                                                                                               ccACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAActct              >  5:114264/1‑137 (MQ=37)                                                                                                                                               ccACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAActct              <  6:114264/137‑1 (MQ=37)                                                                                                                                               ccACAGAGGTTGCATGAAACC‑ATGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGcc  >  7:56360/1‑149 (MQ=38)                                                                                                                                               ccACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGcc  <  8:181132/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GCGCAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGG‑TGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTG‑CCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCC  >  NC_002947/5655480‑5655769

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