Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,740,819	2 bp→CT	intergenic (+56/+696)	gltA → / ← yeiW	citrate synthase/putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,740,819	0	T	C	97.8%	142.6 / ‑4.8	47	intergenic (+56/+697)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (25/20);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01
*	NC_002947	4,740,820	0	G	T	97.9%	145.3 / ‑3.2	47	intergenic (+57/+696)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (26/20);  minor base . (0/1);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

AATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCT‑CCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCA  >  NC_002947/4740670‑4740954                                                                                                                                                      ||                                                                                                                                         aaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCt                                                                                                                                           >  2:369512/1‑149 (MQ=255)        tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGGCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCg                                                                                                                                    >  8:132875/1‑149 (MQ=255)          cAGGCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGGCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCg                                                                                                                                    <  7:39104/147‑1 (MQ=255)             cccGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctc                                                                                                                               >  8:237655/1‑149 (MQ=255)                 tACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcg                                                                                                                           <  7:18234/149‑1 (MQ=255)                     aGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑c                                                                                                                      >  1:68124/1‑149 (MQ=255)                          ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGc                                                                                                                 <  8:348041/149‑1 (MQ=255)                            ccGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGccc                                                                                                               >  1:255611/1‑149 (MQ=255)                            ccGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGccc                                                                                                               <  2:68124/149‑1 (MQ=255)                              gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTGGCGG‑CCGTGGCcgc                                                                                                             >  1:158189/1‑149 (MQ=255)                              gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCTTAGCGGCCCTTGccgc                                                                                                              >  7:134826/1‑147 (MQ=255)                                 cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑GCAAGGCCGCGaa                                                                                                          <  3:379382/149‑1 (MQ=255)                                 cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGaa                                                                                                          <  8:351690/149‑1 (MQ=255)                                  gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑GCAAGGCCGCGAag                                                                                                         <  6:334056/149‑1 (MQ=255)                                  gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCcgcg                                                                                                            >  6:16504/1‑146 (MQ=255)                                   cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCT‑GTGCCGGCCTCTTCGCGG‑GCAAGGCCGCGAagag                                                                                                       <  1:369512/149‑1 (MQ=38)                                    gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAagag                                                                                                       <  5:16504/149‑1 (MQ=255)                                        gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGCCg                                                                                                   >  8:198187/1‑149 (MQ=255)                                           gTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTTCCCGCGAAGAGGCCGGca                                                                                                >  5:355188/1‑149 (MQ=255)                                            tACAGCGGCGAGCAGAAGCGCGACATCGTTGGCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑cc                                                                                                                     <  3:54296/127‑1 (MQ=255)                                            tACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTATTATTTGTCATTTGCCGGCCTCTTCGCGG‑cc                                                                                                                     >  4:54296/1‑127 (MQ=255)                                                 cGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGACGAACGCAAGAAGGCTGCCCTCGGGCAGAC‑TTTTTTATTTGTCCTGTGCCGGCCTCTACGCGG‑GCAAGGCCGCGCAGCGGCCGGCCCAggcc                                                                                          >  8:184450/1‑146 (MQ=255)                                                         aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTCGCCGAACGCAACAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCTAAGAGGCCGGCACAGACACCACaaaa                                                                                  <  6:178742/149‑1 (MQ=255)                                                           aaGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGC‑CCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAAcc                                                                                <  7:237655/149‑1 (MQ=255)                                                                               gAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGG‑CCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTAT‑CCGCCCGCt                                                           >  3:134369/1‑149 (MQ=255)                                                                                  aGACGGCAAGTAAGCTGCAGTAGCCGAACGCAAGACGGCTGCCCGCGGGCAGCC‑TTTTTTATTTGTCCTGGGCCGGCCTCTGCGCGG‑CCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCT‑CCGCCCGCTCAc                                                        <  6:355188/149‑1 (MQ=255)                                                                                     ccGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGGCTATTAGCGG‑CCTTGCCCGCGAAGAATCCGGCACAGAAACCGCAAAAACTCAGTTTTTCT‑CCCGCCCCTCACGcc                                                     >  8:271598/1‑148 (MQ=255)                                                                                            tAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGCCGGCACAGACCCCACAAAACCTCAGTTTTTCT‑ACGCCCGCTCACGCAGCCCtt                                               >  1:266173/1‑148 (MQ=255)                                                                                            tAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCATGCCCGCGAAGAGGCCGGCACAGACAACACAAAACCGCAGTTTTTAT‑CCGCCCCCACACGCAGcccat                                               >  8:23741/1‑146 (MQ=255)                                                                                              aGCTGCAGTACCCCACCGCAACAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGCCCTGTCCCGGCCTCTGCGCGG‑CCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCT‑CCGCCCGCTCACGCAGCCCttt                                              <  4:249609/147‑1 (MQ=255)                                                                                                cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCT‑CCGCCCGCTCACGCAGCCCTTTCAgg                                          <  7:198187/149‑1 (MQ=255)                                                                                                   cAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGCCGGCCCCGACACCACAAAACCTCAGGTTTTCT‑ACGCCCGCGCACGCAGCCCTGTCAGGGTa                                       >  7:187844/1‑149 (MQ=255)                                                                                                       aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGCCCTGTGCCGGCCTCTGCGCGG‑CCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCT‑CCGCCCGCTCACGCAGCCCTTTCAGGGTATTg                                    <  7:23741/148‑1 (MQ=255)                                                                                                         ccGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGTCTCTTCGCGG‑CCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCT‑CCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAAc                                 >  4:378655/1‑149 (MQ=255)                                                                                                         ccGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGCCGGCACAGACACCCCAAAACCTCAGTTTTTCT‑CCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAAc                                 >  6:365690/1‑149 (MQ=255)                                                                                                                cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGACGGCACAGCCACCACAAAACCACAGTTTTTAG‑CCGCCCGCTCACGAAGCCCTTTCAGGGTAGTGAACGGCGCAt                          >  8:270779/1‑149 (MQ=255)                                                                                                                   gaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCT‑CCGCCCGCTCACGCAGCCCTTTCAGGGTATTg                                    >  2:50289/1‑136 (MQ=255)                                                                                                                   gaagGCGGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCCGCCCGCTCACGCAGCCCTTTCAGGGTATTg                                    <  1:50289/137‑1 (MQ=38)                                                                                                                         tGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGTCGGCACAGACACCACAAAACCTCAGTTTTTCT‑CCGCCCGCTCACGCAGc                                                   >  5:150802/1‑115 (MQ=255)                                                                                                                         tGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGTCGGCACAGACACCACAAAACCTCAGTTTTTCT‑CCGCCCGCTCACGCAGc                                                   <  6:150802/115‑1 (MQ=255)                                                                                                                         tGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCT‑CCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACg                  >  7:356585/1‑148 (MQ=255)                                                                                                                         tGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCT‑ACGCCCGCTCACGCAGACCTTTCAGGGTATTGAACGGCGCATACACCACg                  >  3:257798/1‑148 (MQ=255)                                                                                                                           cccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCT‑CCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc               >  6:180667/1‑149 (MQ=255)                                                                                                                            ccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGCCGGCACAGACACCCCAAAACCTCAGTTTTTAT‑GCGCCCGCTGACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCCACACGacag              >  8:357754/1‑146 (MQ=255)                                                                                                                                gggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCT‑CCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt          <  5:180667/149‑1 (MQ=255)                                                                                                                                gggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCT‑CCGCCCGCTCACGCAGCCCTTTCAGGGTAt                                      <  5:377426/121‑1 (MQ=38)                                                                                                                                   cAGCC‑TTTTTTATTTGTCCTGTGCCGGTCTCTTCGCGG‑CCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCT‑CCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGc       <  3:378655/149‑1 (MQ=255)                                                                                                                                         ttttttATTTGTCCTGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCTCAGTTTTTCT‑CCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGccgcca  >  7:84054/1‑149 (MQ=255)                                                                                                                                                      ||                                                                                                                                         AATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGG‑CCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCT‑CCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCA  >  NC_002947/4740670‑4740954

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