Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,951,377	G→C	50.0%	intergenic (‑113/+48)	oguA ← / ← PP_2585	8‑oxoguanine deaminase/alpha‑ketoglutaric semialdehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,951,377	0	G	C	50.0%	‑6.1 / 19.3	14	intergenic (‑113/+48)	oguA/PP_2585	8‑oxoguanine deaminase/alpha‑ketoglutaric semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base C (4/3);  total (7/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 = 5.65e-01

GAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTC  >  NC_002947/2951235‑2951524                                                                                                                                               |                                                                                                                                                    gagCGTTGTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGACGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCtc                                                                                                                                              >  2:442325/1‑149 (MQ=255)      tttttGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGcc                                                                                                                                         <  1:321921/150‑1 (MQ=255)           gCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGaa                                                                                                                                    <  2:400722/150‑1 (MQ=255)           gCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGaa                                                                                                                                    <  2:478190/150‑1 (MQ=255)              gggTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGAGTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAagag                                                                                                                                 <  1:442325/150‑1 (MQ=39)                                                                  gATGTCATAGGCGCTCGGTATTGGAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTcaca                                                                             <  1:88115/150‑1 (MQ=255)                                                                                                           aTTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGGACGGGTTGGAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGttttt                                    >  1:319523/1‑150 (MQ=255)                                                                                                                                          ctctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCgg                                                              >  1:343601/1‑93 (MQ=25)                                                                                                                                          ctctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCgg                                                              <  2:343601/93‑1 (MQ=25)                                                                                                                                           tctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGCTCGAAACGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGt    >  2:366049/1‑150 (MQ=255)                                                                                                                                           tctTCGCGGCTCTAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGc                                                  <  1:167226/104‑1 (MQ=255)                                                                                                                                           tctTCGCGGCTCTAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGc                                                  >  2:167226/1‑104 (MQ=38)                                                                                                                                           tctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGt    >  2:185213/1‑150 (MQ=39)                                                                                                                                             ttGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTc  >  1:354101/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                    GAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTC  >  NC_002947/2951235‑2951524

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: ‑