Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,586,033:1	+G	100%	intergenic (+116/+99)	PP_4061 → / ← PP_4063	hypothetical protein/long‑chain fatty acid‑‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,586,033	1	.	G	94.1%	48.5 / ‑2.9	17	intergenic (+116/+99)	PP_4061/PP_4063	hypothetical protein/long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (4/12);  total (5/12)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GG‑TAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCC  >  NC_002947/4585927‑4586116                                                                                                             |                                                                                       aCGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATg                                             >  7:274807/1‑149 (MQ=39)  cGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgc                                            <  6:2032/149‑1 (MQ=39)   ggCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                           <  5:385823/149‑1 (MQ=39)    gCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcgg                                          <  8:274807/149‑1 (MQ=39)    gCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcgg                                          <  5:244508/149‑1 (MQ=39)                aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTgg                              <  1:89105/149‑1 (MQ=39)                   gCTGATCTGCGCTGCTCCATTGACAGGTTGAGGATGACCGGGGGCTT‑GGCAAGGTGAAGATCTTGCATTGCTCGGGCCGGCCGCTTCGCGGG‑TACACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAg                           <  3:249431/149‑1 (MQ=255)                   gCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAg                           <  4:147157/149‑1 (MQ=39)                       aTCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGa                         <  1:236007/147‑1 (MQ=39)                            cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGtt                  >  2:139983/1‑149 (MQ=38)                              ctgctCCATTGAAAGGTTGAGGAAGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCTTTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTa                <  7:251884/149‑1 (MQ=38)                                 ctCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAc               <  1:139983/147‑1 (MQ=38)                                    cATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCcgc           >  7:28983/1‑148 (MQ=38)                                         aaaGGTGGAGTATCACCGGG‑GCTTGGGCAAGGTTAAGATGTTGAATTGATCGGGCCGGCCTCTTCGCGGG‑AAAACCCGCTCCTACCCGGGACC‑GGGCCGGCTTTTAATGCGGCGCTGCCCCTGGTTGGAGCGGTTTTACCCGCGAaga      <  8:28983/148‑1 (MQ=255)                                            ggTTGAGGATGCCCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACCCGGGA‑CCGGGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggcc  >  2:47261‑M2/1‑142 (MQ=255)                                                                tt‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCt                                <  2:184932/99‑1 (MQ=25)                                                                tt‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCt                                >  1:184932/1‑99 (MQ=25)                                                                                                             |                                                                                       ACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GG‑TAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCC  >  NC_002947/4585927‑4586116

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