Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,741,243	3 bp→TGC	intergenic (+480/+271)	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,243	0	G	T	100.0%	32.1 / NA	11	intergenic (+480/+273)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (4/7);  total (4/7)
*	NC_002947	4,741,244	0	C	G	90.9%	30.1 / NA	11	intergenic (+481/+272)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (4/6);  minor base T (0/1);  total (4/7)
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
*	NC_002947	4,741,245	0	A	C	90.9%	29.7 / ‑6.0	11	intergenic (+482/+271)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (4/6);  total (4/7)
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

CTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGAACCTCCGCCAACTGCAATTCGTGAACCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTC  >  NC_002947/4741127‑4741362                                                                                                                     |||                                                                                                                       ctcctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCggac                                                                                           <  7:23773/148‑4 (MQ=255)     ttGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCggacacgt                                                                                       >  1:143142/1‑141 (MQ=255)                aaTGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGc                                                                            <  8:93753/148‑1 (MQ=255)                          aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGcgcctcgacct                                                                  >  1:27986/1‑141 (MQ=255)                                         gTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt                                                   <  2:158540/148‑1 (MQ=255)                                                gCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc                                            <  6:178836/148‑1 (MQ=255)                                                                    gTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCt                        <  2:81578/148‑1 (MQ=255)                                                                            gcAGCGCTGGTCAGCTTCGTGCTTGTGAAGCTTCTGCGACTTTAATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCTCCTGGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg                <  8:103749/147‑1 (MQ=255)                                                                              aGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGa              <  2:143142/148‑1 (MQ=255)                                                                                      tCACCTTCGTGCTTGAGCAGCTTCTGCGCCTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCa      >  5:92406/1‑148 (MQ=255)                                                                                          cTTCGTGCTTGAGCAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTc  >  2:105469/1‑148 (MQ=255)                                                                                                                     |||                                                                                                                       CTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGAACCTCCGCCAACTGCAATTCGTGAACCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTC  >  NC_002947/4741127‑4741362

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