Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,000,186	G→A	G179D (GGT→GAT)	ydgD →	predicted peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,000,186	0	G	A	96.0%	156.4 / ‑1.5	50	G179D (GGT→GAT)	ydgD	predicted peptidase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (29/19);  total (31/19)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

ATTACGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTG  >  minE/1000122‑1000215                                                                    |                              aTTACGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTa                          >  1:1239801/1‑70 (MQ=255)    aCGCCGTTGCCGTTATTTTTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTa                          >  1:12160/1‑70 (MQ=255)    aCGCCGTTGCCGTTATTTTTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTa                          >  1:469835/1‑70 (MQ=255)    aCGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAGGCGGCAGATCGTAAAGt                      >  1:236041/1‑71 (MQ=255)    aCGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTaaa                        >  1:252190/1‑69 (MQ=255)    aCGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAg                       >  1:1645773/1‑70 (MQ=255)    aCGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                      >  1:236772/1‑71 (MQ=255)    aCGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                      >  1:948982/1‑71 (MQ=255)    aCGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                      >  1:839354/1‑71 (MQ=255)    aCGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                      >  1:80331/1‑71 (MQ=255)    aCGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                      >  1:3045534/1‑71 (MQ=255)    aCGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                      >  1:2960881/1‑71 (MQ=255)    aCGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                      >  1:2867280/1‑71 (MQ=255)    aCGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                      >  1:2759334/1‑71 (MQ=255)    aCGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                      >  1:2259202/1‑71 (MQ=255)    aCGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                      >  1:1938671/1‑71 (MQ=255)    aCGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                      >  1:1221183/1‑71 (MQ=255)    aCGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                      >  1:1250721/1‑71 (MQ=255)               gttA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTcaggcag           >  1:1122563/1‑71 (MQ=255)               gttA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTcaggcag           >  1:1523885/1‑71 (MQ=255)               gttA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTcaggcag           >  1:1840818/1‑71 (MQ=255)               gttA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTcaggcag           >  1:2875773/1‑71 (MQ=255)               gttA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTcaggcag           >  1:2874249/1‑71 (MQ=255)               gttA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTcaggcag           >  1:438888/1‑71 (MQ=255)               gttA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTcaggcag           >  1:2763701/1‑71 (MQ=255)               gttA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTcaggcag           >  1:2607876/1‑71 (MQ=255)               gttA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTcaggcag           >  1:498895/1‑71 (MQ=255)               gttA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTcaggcag           >  1:2276696/1‑71 (MQ=255)               gttA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTcaggcag           >  1:2044237/1‑71 (MQ=255)               gttA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTcaggcag           >  1:2021304/1‑71 (MQ=255)               gttA‑‑‑TTAGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTcaggcag           >  1:418206/1‑71 (MQ=255)                           gggAGATAAAGCCGTGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:1553516/71‑1 (MQ=255)                           gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:685806/71‑1 (MQ=255)                           gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:695142/71‑1 (MQ=255)                           gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:930388/71‑1 (MQ=255)                           gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:2938361/71‑1 (MQ=255)                           gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:3212664/71‑1 (MQ=255)                           gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:137818/71‑1 (MQ=255)                           gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:167757/71‑1 (MQ=255)                           gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:1689132/71‑1 (MQ=255)                           gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:2100773/71‑1 (MQ=255)                           gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:1722554/71‑1 (MQ=255)                           gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:1931087/71‑1 (MQ=255)                            ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGTCAGGCTACCCTg  <  1:155167/70‑1 (MQ=255)                            ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:1347426/70‑1 (MQ=255)                            ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:1555060/70‑1 (MQ=255)                            ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:1030931/70‑1 (MQ=255)                                                       ttAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:531203/43‑1 (MQ=255)                                                             gCGGCTGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:810977/37‑1 (MQ=38)                                                             gCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg  <  1:574447/37‑1 (MQ=255)                                                                    |                              ATTACGCCGTTGCCGTTA‑‑‑TTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTG  >  minE/1000122‑1000215

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