Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,741,302	A→G	intergenic (+493/+259)	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,302	0	A	G	77.8%	30.9 / 3.2	19	intergenic (+493/+259)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (8/6);  total (10/9)
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

ATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCTCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCC  >  NC_002947/4741166‑4741438                                                                                                                                             |                                                                                                                                         aTCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGcccc                                                                                                                              <  2:697757/150‑1 (MQ=255)                               caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCtt                                                                                                >  1:650771/1‑150 (MQ=255)                                aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGCA‑TAGCGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                                                               <  1:691642/150‑1 (MQ=255)                                          gcggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggcg                                                                                     <  2:650771/150‑1 (MQ=255)                                              cgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                                                                                 <  1:307065/150‑1 (MQ=255)                                                                         ggTATAGGTCCAGCTCTGCTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCtt                                                      >  2:287692/1‑150 (MQ=255)                                                                          gTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTg                                                     >  1:46244/1‑150 (MQ=255)                                                                             taGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                                                               <  1:642137/105‑1 (MQ=255)                                                                             taGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                                                               >  2:642137/1‑105 (MQ=255)                                                                              aGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggc                                                                                      <  2:16883/113‑1 (MQ=255)                                                                              aGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggc                                                                                      >  1:16883/1‑113 (MQ=255)                                                                                 tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTg                                                     >  1:630322/1‑143 (MQ=255)                                                                                 tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTg                                                     <  2:630322/143‑1 (MQ=255)                                                                                         tGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg                                                           >  1:10824/1‑129 (MQ=255)                                                                                         tGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg                                                           <  2:10824/129‑1 (MQ=255)                                                                                           gTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAc                                    >  1:591391/1‑150 (MQ=255)                                                                                                                          acttGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCa                                                                                           <  1:472592/63‑1 (MQ=255)                                                                                                                          acttGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCa                                                                                           >  2:472592/4‑66 (MQ=255)                                                                                                                               cATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcgtcc  >  2:32088/1‑150 (MQ=255)                                                                                                                                             |                                                                                                                                         ATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCTCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCC  >  NC_002947/4741166‑4741438

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