Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,555,181	C→T	50.0%	R120R (CGC→CGT)	dsdA →	D‑serine ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,555,181	0	C	T	50.0%	‑5.2 / 38.8	44	R120R (CGC→CGT)	dsdA	D‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (16/6);  new base T (15/7);  total (31/13)
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

AACTGTTACTGAAAAAAGATAGCCATTTGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTGGCTCTGGAAGCGGGGTTGCTGACG  >  minE/1555126‑1555244                                                        |                                                                aaCTGTTACTGAAAAAAGATAGCCATTTGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGATTTa                                                      >  1:3296450/1‑67 (MQ=255) aaCTGTTACTGAAAAAAGATAGCCATTTGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGATTTa                                                      >  1:1629463/1‑67 (MQ=255)             aaaaaTACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGt                                                >  1:757434/1‑61 (MQ=255)             aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        >  1:1146509/1‑67 (MQ=255)             aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        >  1:846293/1‑67 (MQ=255)             aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        >  1:3312805/1‑67 (MQ=255)             aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        >  1:2680598/1‑67 (MQ=255)             aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        >  1:1856553/1‑67 (MQ=255)             aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        >  1:1700188/1‑67 (MQ=255)             aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        >  1:1131745/1‑67 (MQ=255)             aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        >  1:1526967/1‑67 (MQ=255)             aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGcc                                         >  1:2861561/1‑67 (MQ=255)             aaaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGcc                                         >  1:2705092/1‑67 (MQ=255)              aaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        <  1:999877/68‑3 (MQ=255)              aaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        <  1:1277331/68‑3 (MQ=255)              aaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        <  1:2445274/68‑3 (MQ=255)              aaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        <  1:346316/68‑3 (MQ=255)              aaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        <  1:3353819/68‑3 (MQ=255)              aaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        <  1:1606249/68‑3 (MQ=255)              aaaaGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGccc                                        <  1:3098251/68‑3 (MQ=255)                 aGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTg                                            >  1:2240034/1‑61 (MQ=255)                 aGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGcccatg                                     >  1:1939818/1‑66 (MQ=255)                 aGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGcccatg                                     >  1:1901940/1‑66 (MQ=255)                 aGACAGCCATTTGCCGATTTCCGGCTCCATAAAAGCACGTGGCGGGATTTATGAAGTCCTGGcccatg                                     >  1:3225326/1‑66 (MQ=255)                            tGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGAtt                                                        >  1:1925736/1‑38 (MQ=255)                            tGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACt                           >  1:1973960/1‑67 (MQ=255)                            tGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTgg                         >  1:293777/1‑69 (MQ=255)                            tGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTgg                         >  1:813951/1‑69 (MQ=255)                            tGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTgg                         >  1:2421764/1‑69 (MQ=255)                                    tCCGGCTCCATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTGGCTCTGGaa                 >  1:262062/1‑69 (MQ=255)                                           ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCAGACGCAGaaaa                              >  1:2143537/1‑49 (MQ=255)                                           ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              <  1:2529459/49‑1 (MQ=255)                                           ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              <  1:18646/49‑1 (MQ=255)                                           ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              >  1:662595/1‑49 (MQ=255)                                           ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              >  1:366951/1‑49 (MQ=255)                                           ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              <  1:2548866/49‑1 (MQ=255)                                           ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              <  1:2719207/49‑1 (MQ=255)                                           ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              >  1:2826272/1‑49 (MQ=255)                                           ccATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              >  1:2925413/1‑49 (MQ=255)                                              tAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              <  1:1373960/46‑1 (MQ=255)                                                aaaGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTGGCTCTGGAAGCGGGGTTGCTg     >  1:552345/1‑69 (MQ=255)                                                aaaGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTGGCTCTGGAAGCGGGGTTGCTg     >  1:3016701/1‑69 (MQ=255)                                                  aGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGaaaa                              <  1:91914/42‑1 (MQ=255)                                                   ggaCGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTGGCTCTGGAAGCGGGGTTGCTGACg  >  1:1840782/3‑69 (MQ=255)                                                        |                                                                AACTGTTACTGAAAAAAGATAGCCATTTGCCCATTTCCGGCTCCATAAAAGCACGCGGCGGGATTTATGAAGTCCTGGCACACGCAGAAAAACTGGCTCTGGAAGCGGGGTTGCTGACG  >  minE/1555126‑1555244

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