Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,741,260	2 bp→GC	intergenic (+497/+255)	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,260	0	C	G	92.9%	29.1 / 0.6	28	intergenic (+497/+256)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base G (15/11);  total (16/12)
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.69e-01
*	NC_002947	4,741,261	0	A	C	96.3%	21.1 / ‑4.3	28	intergenic (+498/+255)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (15/11);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

AACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCTCCGCCAACTGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAG  >  NC_002947/4741188‑4741368                                                                           ||                                                                                                             aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGgcgc                                                           <  1:4023928/125‑1 (MQ=255)  aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc                                                          <  1:4023930/125‑1 (MQ=255)   cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt                                                         <  1:4023931/125‑1 (MQ=255)     gggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc                                                          <  1:4023933/122‑1 (MQ=255)      ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTg                                                       >  1:4024103/1‑124 (MQ=255)      ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGa                                                      >  1:4024073/1‑125 (MQ=255)      ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGa                                                      >  1:4023972/1‑125 (MQ=255)        tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                                                    <  1:4023934/125‑1 (MQ=255)         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                                                    >  1:4024074/1‑124 (MQ=255)         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                                                    >  1:4024097/1‑124 (MQ=255)         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                                                   <  1:4023935/125‑1 (MQ=255)         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACATCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                                                   <  1:4023936/125‑1 (MQ=255)           aGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                                 >  1:4024216/1‑125 (MQ=255)              tCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                                              >  1:4024160/1‑125 (MQ=255)                cAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcac                                            >  1:4024033/1‑125 (MQ=255)                 aGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcac                                            >  1:4024110/1‑124 (MQ=255)                  gCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                                          >  1:4024376/1‑125 (MQ=255)                  gCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                                          >  1:4024189/1‑125 (MQ=255)                      tGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgcg                                      <  1:4023953/125‑1 (MQ=255)                        gTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGATACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGAt                                    <  1:4023969/125‑1 (MQ=255)                            ccTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAATCCAGGCAC‑AGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc                                <  1:4023976/125‑1 (MQ=255)                             cTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc                               >  1:4024134/1‑125 (MQ=255)                             cTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc                               <  1:4023951/125‑1 (MQ=255)                                                      acttgcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTc        >  1:4024159/7‑124 (MQ=255)                                                      acttgcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTc        >  1:4024194/7‑124 (MQ=255)                                                       cttgcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCAc      >  1:4024172/6‑125 (MQ=255)                                                         tgcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACt     <  1:4023954/121‑1 (MQ=255)                                                           caTTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAg  >  1:4024251/2‑125 (MQ=255)                                                                           ||                                                                                                             AACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCTCCGCCAACTGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAG  >  NC_002947/4741188‑4741368

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