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.9%	149.5 / NA	47	intergenic (+56/+697)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (24/22);  minor base G (0/1);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_002947	4,740,820	0	G	T	97.8%	139.3 / NA	47	intergenic (+57/+696)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base T (23/22);  minor base A (1/0);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACC  >  NC_002947/4740677‑4740953                                                                                                                                                 ||                                                                                                                                      tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCg                                                                                                                                 >  4:346805/1‑149 (MQ=255) tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCg                                                                                                                                 >  3:273280/1‑149 (MQ=255) tcCAGCCCGTACAAGATCGGCCGTCCGCCCCAGCTGTACACCGGCGAGGAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCCGCC‑TTTTTTATTTGTCCAGTGCCg                                                                                                                                 >  2:48195/1‑149 (MQ=255) tcCAGCCCGTACAAGATCGGACGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCg                                                                                                                                 >  5:273465/1‑149 (MQ=255)   cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCg                                                                                                                                 <  8:324684/147‑1 (MQ=255)   cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCg                                                                                                                                 <  1:184012/147‑1 (MQ=255)   cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGc                                                                                                                               <  7:124710/149‑1 (MQ=255)         gTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTc                                                                                                                         >  1:290833/1‑149 (MQ=255)         gTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTc                                                                                                                         <  3:180036/149‑1 (MQ=255)         gTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACCCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctct                                                                                                                           >  1:295019/1‑147 (MQ=255)             aaGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCCGCC‑GTGTTTATTTGTCCTGTGCCGGCCGCTTCGCgg                                                                                                                     >  7:140929/1‑149 (MQ=255)               gATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGcc                                                                                                                   <  2:77026/149‑1 (MQ=255)                   ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGc                                                                                                               <  4:273280/149‑1 (MQ=255)                    gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCTGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGcc                                                                                                              <  6:273465/149‑1 (MQ=255)                               atgCTTTATACCGGCTAGAAGAAGCTCGACATCTTTGCCCTGAAATACCGCAAGTAATCTGCAGTAGCCGAACGCAAGAA‑‑GTCTGGCCTCGGGCAG‑CGTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAagag                                                                                                     <  4:298407/145‑1 (MQ=255)                                aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGcc                                                                                                  <  3:346805/149‑1 (MQ=255)                                  cTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCACTACCCGAACGCAAGAC‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCgg                                                                                                <  1:88095/149‑1 (MQ=255)                                       cacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAACTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCCGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACTAAACGGCAATGCACGCCAAGCTGCAGGcaca                                                                                            >  6:123513/1‑148 (MQ=255)                                                 cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaa                                                                                  >  8:33672/1‑148 (MQ=255)                                                   gacgcgcGACAGCGTTGCCCTGAAAGCCCGCAAGTAAGCTGCAGTAGCCTAACGCAAGAC‑‑GGCGGCCCTCGGGCTGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTATGAGGCCGGCACAGACACCACAAAAc                                                                               <  2:204937/146‑1 (MQ=255)                                                      gcgcGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCt                                                                             >  5:100023/1‑148 (MQ=255)                                                       cgcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCa                                                                           <  7:112242/149‑1 (MQ=255)                                                       cgccACTTCGTTGCCCTCAACCCCCCCACGTACGCGCCTGTACCCCACCCCAACAC‑‑GGCTCCCCGGGGGCCG‑CTTTTTTTTTTTGCCCGGTGCCGCCCCCTCCGCGGCCTTGCCCCCTACGCGGCCGGCACAGACACCACAAAACCTCa                                                                           <  2:290833/149‑1 (MQ=255)                                                          gACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                        >  4:226435/1‑117 (MQ=255)                                                          gACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                        <  3:226435/117‑1 (MQ=255)                                                             aTCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGc                                                                                                   <  2:294202/119‑1 (MQ=255)                                                                       tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccgc                                                           >  2:52917/1‑149 (MQ=255)                                                                        gAAAGACCGCAAGTAAGCTGCAGTAGCCGAATGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGCCACCACAAAACCTCAGTTTTTCTCCGCCCGCt                                                          >  1:37761/1‑149 (MQ=255)                                                                         aaaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTACGCCCGCTc                                                         >  7:182309/1‑149 (MQ=255)                                                                               cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGAGCCCGGCACCGCCACCACAACACCGCACTTTTTGACCCCCCCCCCACGCAg                                                   >  7:196566/1‑149 (MQ=255)                                                                                       aGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCttt                                             >  1:44697/1‑147 (MQ=255)                                                                                       aGCGGCAGTAGCCGAACGCTTGAA‑‑GGCTGCCCTGGGGCTG‑CTTTTTTTATTTGCCCTTGGCCGGCCTCTGCGGGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                           <  5:123513/149‑1 (MQ=255)                                                                                           gCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                        >  5:14066/1‑84 (MQ=38)                                                                                           gCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                        <  6:14066/84‑1 (MQ=38)                                                                                             aGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                                     >  7:242609/1‑149 (MQ=255)                                                                                                  ccGAACGCAAGTACCGCCTGCCCCCGGGCAG‑CGTTTTTTATTTGTCCTGGGCCGGCCTGTCCGCGGCTTTCCCCGCCATGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGa                                  <  1:303745/149‑1 (MQ=255)                                                                                                  ccGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAAc                                >  8:97205/1‑149 (MQ=255)                                                                                                     aaCGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                             <  7:97205/149‑1 (MQ=255)                                                                                                      aCGCAAGAA‑‑GGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCGCACGCAGCCCTGTCAGGGTATTGAACGgcg                            >  3:286848/1‑149 (MQ=255)                                                                                                          aagaa‑‑gGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc                        <  6:100023/149‑1 (MQ=255)                                                                                                                  gCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCcaccac                  <  3:118388/149‑1 (MQ=255)                                                                                                                    tGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCGGGCACAGCCACCCAAAAACCGCAGTTTTTTTCCGCCCGCGCACGCAGCCCTTTAAGGGTTTTGAACCGCGCATCCCCCCCg                 >  3:204024/1‑148 (MQ=255)                                                                                                                       ccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAG‑‑TTTCTCAGCCCGCTCACGCAGCCCGTTAAGGGTATTGAAAGGCGCATCCACCACGAACttg           >  3:285664/1‑149 (MQ=255)                                                                                                                        cgagggCAGCC‑TTTTTTATTTGACGTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACtt            <  7:39156/146‑1 (MQ=255)                                                                                                                          cGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCCTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgt          >  5:210491/1‑149 (MQ=255)                                                                                                                                   cttttttATTTGTCCTGTGCCGGCCTCCTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGccacc  <  6:210491/148‑1 (MQ=255)                                                                                                                                    ttttttATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTATCCGCCCCCTCACCCAGCCCATTCAGGGTATTGAAAGGCGCATCCACCACCAACTTGTTGGccacc  >  2:2305/1‑148 (MQ=255)                                                                                                                                                 ||                                                                                                                                      TCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAA‑‑GGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACC  >  NC_002947/4740677‑4740953

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