Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,741,275	A→G	100%	intergenic (+466/+286)	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,275	0	A	G	86.8%	27.2 / ‑3.1	14	intergenic (+466/+286)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (5/7);  total (6/8)
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
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCCGCCAACTG‑CATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGG  >  NC_002947/4741169‑4741382                                                                                                            |                                                                                                               cGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTG‑CATTCGTGAATCCATGC‑ACAGGCCCCGGa                                                                   <  2:305486/150‑1 (MQ=255)  gCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGATTTAGCAG‑CTTCTGCGACTTG‑CATTCGTGAATCCAGGC‑ACAGGCCCCGGac                                                                  <  1:419558/150‑1 (MQ=255)  gCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGGATTGAACTGGGGATAGGTCCAGCTCTGGTCACCTTCGTGCTTGCGCGG‑CTTCTGCGATTTG‑CATTCGCGAC‑TCCAGGCACAGGCCCCGGac                                                                  >  2:194565/1‑150 (MQ=255)                      gtgCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTG‑CATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTc                                              >  1:191434/1‑150 (MQ=255)                         catcaCGGAAGCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTG‑CATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGAc                                           >  2:565923/4‑150 (MQ=255)                           gcaCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTG‑CATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCt                                         >  2:465214/1‑150 (MQ=255)                                           cgTCACCGGCCATGGCGTATTGAACTGCGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTG‑CATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                         <  1:239731/150‑1 (MQ=255)                                               aCCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTG‑CATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                     <  1:565923/150‑1 (MQ=255)                                               aCCGGCCATGGCGTATTGAAATGGGTATAGGTCCAGCTATGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTG‑CATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGACAGCGGCGCCTTGATc                     >  1:176853/1‑150 (MQ=255)                                                  ggCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTCATTCCG‑CA‑GATTTTGCGACTGTCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                  <  1:140991/150‑1 (MQ=255)                                                     cATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTG‑CATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc               >  2:425161/1‑150 (MQ=255)                                                         gtgTATTGAACTGCGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTCAGCAG‑CTTCTGCGACTTG‑CATTCGTGAA‑TGCAGGCACAGGCCCCGGACACGTCCTCCTCCAGCGCCTCTACCTTGGCCAGCCGCGCCTTGATCACAGTCacac           <  2:176853/148‑1 (MQ=255)                                                              ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTG‑CATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGa      <  1:425161/150‑1 (MQ=255)                                                                  aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACTTG‑CATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAgg  <  1:465214/150‑1 (MQ=255)                                                                                                            |                                                                                                               CGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCCGCCAACTG‑CATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGG  >  NC_002947/4741169‑4741382

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