Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,655,641	Δ1 bp	16.7%	intergenic (‑96/+56)	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,641	0	T	.	16.7%	113.3 / 20.1	42	intergenic (‑96/+56)	epd/tktA	D‑erythrose 4‑phosphate dehydrogenase/transketolase A
Reads supporting (aligned to +/- strand):  ref base T (17/18);  new base . (5/2);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

CGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACC  >  NC_002947/5655497‑5655789                                                                                                                                                  |                                                                                                                                                      cGATACGGCCATAACCGTGGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCtgt                                                                                                                                                   <  2:253936/149‑1 (MQ=255)  gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGGCGGTGTCGCCT‑TCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTgg                                                                                                                                                 <  4:154663/149‑1 (MQ=17)  gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCtgtg                                                                                                                                                  <  6:304851/149‑1 (MQ=255)        gCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCCCTCCCAAAC‑GGGTTGATGCAACCTCTGTGGGGGCg                                                                                                                                            >  4:5806/1‑149 (MQ=255)         ccATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCgg                                                                                                                                           <  6:133461/149‑1 (MQ=255)         ccATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCgg                                                                                                                                           >  7:219243/1‑149 (MQ=255)          cATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCggg                                                                                                                                          <  1:310060/149‑1 (MQ=255)           aTAACCGTTGATTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGc                                                                                                                                         <  8:209697/149‑1 (MQ=255)              aCCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCtt                                                                                                                                       <  4:129978/148‑1 (MQ=255)                 gTTGAGTGCATCTTTGTAGGGTCGGTGGGGCATTGGGGGCTCACGAGACAGTTCAATTGTCGTCTGCCTGCGGATTTTGCTGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGTGCGGGCTTGccc                                                                                                                                   <  6:55377/149‑1 (MQ=38)                      gTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTgg                                                                                                                                                 <  8:94809/130‑1 (MQ=255)                      gTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTgg                                                                                                                                                 >  7:94809/1‑130 (MQ=255)                       tggaaCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GCTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAAt                                                                                                                             >  7:230015/4‑149 (MQ=255)                                                     tggCTCACTAGACAGTTCAATTGTCGTCAGCCTGCTGAATTTGCTGTGTGGCATTCGGGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAAt                                                                                              <  8:230015/148‑1 (MQ=255)                                                       gCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACAccgc                                                                                                              >  4:110996/1‑131 (MQ=255)                                                       gCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACAccgc                                                                                                              <  3:110996/131‑1 (MQ=255)                                                                     tCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCACCAAAAG‑TATTTAATGAAACCTCTGATCTAGCGGACTTGCCAGCAAATT‑AAGCAGCCACCGCCCCCAGAGATTTAAATCAGTCTTccaccacca                                                                              >  2:75483/1‑146 (MQ=255)                                                                           gtcgtcAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCg                                                                         <  8:219243/148‑1 (MQ=255)                                                                                               gCGGTGTGCCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTtctc                                                     <  3:5806/148‑1 (MQ=255)                                                                                                cGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCAGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTcc                                                    <  7:229236/148‑1 (MQ=255)                                                                                                             gcgGGCAAGCCGTCTACCTCTG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCTGTACCGAGGATGTTCTCCAggg                                                <  2:52423/139‑1 (MQ=255)                                                                                                             gcgGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAggg                                                >  1:52423/1‑139 (MQ=255)                                                                                                                 gCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCaa                                       >  6:324470/1‑144 (MQ=255)                                                                                                                 gCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCaa                                       <  5:324470/144‑1 (MQ=255)                                                                                                                   aaGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTc                                >  5:153366/1‑149 (MQ=255)                                                                                                                       ccGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCa                                                                      <  5:43277/107‑1 (MQ=255)                                                                                                                       ccGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCa                                                                      >  6:43277/1‑107 (MQ=255)                                                                                                                           tCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATg                        >  3:168111/1‑149 (MQ=38)                                                                                                                            cccACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTc                                                                                      <  4:280331/86‑1 (MQ=255)                                                                                                                            cccACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTc                                                                                      >  3:280331/1‑86 (MQ=255)                                                                                                                             ccACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCCACACCGCAAATTCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAAc                            <  2:118828/144‑1 (MQ=17)                                                                                                                             ccACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCCACACCGCAAATTCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAAc                            >  1:118828/1‑144 (MQ=17)                                                                                                                             ccACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGcc                      >  7:250837/1‑149 (MQ=38)                                                                                                                             caacaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGCTCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGACGATGTTCTCCAGGGTGAAGCCAACCACTTCGAACAATgac                      >  1:168890/3‑147 (MQ=255)                                                                                                                               acagaggTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg                     >  1:355917/6‑148 (MQ=38)                                                                                                                               acaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg                     >  8:32505/1‑148 (MQ=255)                                                                                                                                  g‑ggTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGc                 >  6:360957/1‑149 (MQ=255)                                                                                                                                      atgcATGAAACC‑CTGTCGGAGCCGGCTTCCCCGCGAATG‑AAGTCGACACCGCCCCCAGATCATTAAATAACTCATCCAGCACCTCATCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTAGAACAATGACGAACACGgc             >  2:185892/4‑149 (MQ=255)                                                                                                                                           tGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGcc          >  3:228439/1‑148 (MQ=255)                                                                                                                                              aaCCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTtctc                                                     <  3:207489/102‑1 (MQ=255)                                                                                                                                              aaCCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTtctc                                                     >  4:207489/1‑102 (MQ=255)                                                                                                                                              aaCCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACt      >  7:217565/1‑149 (MQ=255)                                                                                                                                                 ctctGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCAc   >  6:237878/1‑149 (MQ=255)                                                                                                                                                  tctGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCAcc  >  5:242564/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                      CGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACC  >  NC_002947/5655497‑5655789

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