Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,068,195	A→T	26.7%	intergenic (‑97/+24)	plsC ← / ← PP_1845	1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase/enoyl‑CoA hydratase/isomerase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,068,195	0	A	T	26.7%	19.7 / 7.0	15	intergenic (‑97/+24)	plsC/PP_1845	1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase/enoyl‑CoA hydratase/isomerase family protein
Reads supporting (aligned to +/- strand):  ref base A (5/6);  new base T (2/2);  total (7/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 = 6.46e-01

AAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGC  >  NC_002947/2068058‑2068338                                                                                                                                           |                                                                                                                                                aaaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAg                                                                                                                                     <  2:444045/150‑1 (MQ=255)   aTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTTCGCGGGCTTGCCCGCGAAAggg                                                                                                                                    <  1:599259/150‑1 (MQ=255)                       tACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAgg                                                                                                                                    >  1:14123/1‑129 (MQ=255)                       tACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAgg                                                                                                                                    <  2:14123/129‑1 (MQ=255)                                                  tCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTTGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGt                                                                                    <  2:506886/150‑1 (MQ=255)                                                                                             tggGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGGCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGATTGGTGtt                                         <  1:447845/150‑1 (MQ=255)                                                                                                gtcAAACCGTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCGCCAGCATCCCTTTGCGGACGTTCGGCGTGGTGTTCtt                                      >  1:244174/3‑150 (MQ=255)                                                                                                ggCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAg                                                                                                                      >  1:593373/1‑70 (MQ=255)                                                                                                ggCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAg                                                                                                                      <  2:593373/70‑1 (MQ=255)                                                                                                                   cTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGc                   <  1:520940/150‑1 (MQ=255)                                                                                                                               ttttCGCGGGCTTGCCCGCG‑AAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTtgt                                    >  1:312000/2‑121 (MQ=255)                                                                                                                               ttttCGCGGGCTTGCCCGCG‑AAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTtgt                                    <  2:312000/120‑1 (MQ=255)                                                                                                                                tttCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTcc       >  1:459032/1‑150 (MQ=255)                                                                                                                                  tCGCGGGCTTGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCcg     >  1:172400/1‑150 (MQ=255)                                                                                                                                     cgGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCgcgc  >  2:360875/1‑150 (MQ=255)                                                                                                                                           |                                                                                                                                                AAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGC  >  NC_002947/2068058‑2068338

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