Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,741,246	Δ1 bp	intergenic (+483/+270)	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,741,245	0	A	.	89.6%	167.3 / 18.9	48	intergenic (+482/+271)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (25/18);  minor base C (2/3);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-03

TTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCCGCCAACTGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCT  >  NC_002947/4741185‑4741340                                                               |                                                                                                  ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcgg                                  >  1:4667579/1‑125 (MQ=255)  tGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTCAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggc                                 >  1:4667515/1‑125 (MQ=255)    aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGgcgc                               <  1:4667320/125‑1 (MQ=255)     aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc                              >  1:4667718/1‑125 (MQ=255)      cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt                             <  2:4667322/125‑1 (MQ=255)      cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt                             <  1:4667324/125‑1 (MQ=255)       tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                            >  2:4667648/1‑125 (MQ=255)       tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                            >  2:4667693/1‑125 (MQ=255)         ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGa                          <  2:4667326/125‑1 (MQ=255)          gTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt                             >  2:4667584/1‑121 (MQ=255)           tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                        <  1:4667328/125‑1 (MQ=255)           tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                        <  2:4667327/125‑1 (MQ=255)            ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                       <  1:4667329/125‑1 (MQ=255)             taGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc                      >  2:4667678/1‑125 (MQ=255)             taGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc                      >  2:4667368/1‑125 (MQ=255)               ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                     <  2:4667330/124‑1 (MQ=255)                gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACgg                    >  2:4667864/1‑124 (MQ=255)                gTCAAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                     >  2:4667554/1‑123 (MQ=255)                 tCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                   <  2:4667331/124‑1 (MQ=255)                 tCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                  >  2:4667652/1‑125 (MQ=255)                  ccAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                 <  1:4667355/125‑1 (MQ=255)                   cAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                 >  1:4667472/1‑124 (MQ=255)                   cAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcac                <  2:4667375/125‑1 (MQ=255)                    aGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcac                >  2:4667491/1‑124 (MQ=255)                    aGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACAAGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcac                >  2:4667716/1‑124 (MQ=255)                    aGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcac                >  1:4667604/1‑124 (MQ=255)                     gCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac              >  2:4667814/1‑125 (MQ=255)                     gCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac              >  2:4667653/1‑125 (MQ=255)                     gCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac              <  1:4667429/125‑1 (MQ=255)                      ctctGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac              <  2:4667414/124‑1 (MQ=255)                      ctctGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAcc             >  1:4667592/1‑125 (MQ=255)                      ctctGGTCACCTCCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac              <  1:4667381/124‑1 (MQ=255)                       tctGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAcc             >  1:4667534/1‑124 (MQ=255)                        ctGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc           <  1:4667428/125‑1 (MQ=255)                        ctGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc           >  2:4667612/1‑125 (MQ=255)                         tGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc           >  2:4667611/1‑124 (MQ=255)                          ggTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGa         <  1:4667420/125‑1 (MQ=255)                           gTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGAt        <  1:4667379/125‑1 (MQ=255)                           gTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGAt        >  1:4667469/1‑125 (MQ=255)                           gTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGAt        >  1:4667757/1‑125 (MQ=255)                           gTCACCTTCGTGCTTGAG‑CA‑GCTGCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGa         <  2:4667380/124‑1 (MQ=255)                             cACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg      <  1:4667423/125‑1 (MQ=255)                             cACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg      <  1:4667394/125‑1 (MQ=255)                             cACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg      >  2:4667713/1‑125 (MQ=255)                              aCCTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAgg     >  2:4667528/1‑125 (MQ=255)                               ccTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAATCCAGGCAC‑AGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc    <  1:4667389/125‑1 (MQ=255)                               ccTTCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc    >  1:4667545/1‑125 (MQ=255)                                 ttCGTGCTTGAGCAG‑CTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCt  >  1:4667573/1‑125 (MQ=255)                                                               |                                                                                                  TTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCCGCCAACTGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCT  >  NC_002947/4741185‑4741340

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑