Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,988,027	C→G	65.7%	G197G (GGG→GGC)	yggN ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,988,027	0	C	G	65.7%	30.1 / 24.6	35	G197G (GGG→GGC)	yggN	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (6/6);  new base G (12/11);  total (18/17)
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

TGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGC  >  minE/1987966‑1988078                                                              |                                                    tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCa                                              >  1:3389359/1‑69 (MQ=255)   ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCa                                              <  1:603638/67‑1 (MQ=255)   ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCa                                              <  1:2886503/67‑1 (MQ=255)   ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCa                                              <  1:2843180/67‑1 (MQ=255)   ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCa                                              <  1:3015982/67‑1 (MQ=255)   ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCa                                              <  1:2994024/67‑1 (MQ=255)   ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCa                                              <  1:809466/67‑1 (MQ=255)   ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCa                                              <  1:708541/67‑1 (MQ=255)   ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCa                                              <  1:2564967/67‑1 (MQ=255)   ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCa                                              <  1:2191520/67‑1 (MQ=255)    cAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCa                                              <  1:2325624/66‑1 (MQ=255)     aaaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCa                                              <  1:2425083/65‑1 (MQ=255)       aCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGcc                                                    >  1:1886428/1‑57 (MQ=255)            tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCAGGCTTCCCAg                                    >  1:344852/1‑68 (MQ=255)            tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCAGGCTTCCCAGc                                   >  1:344048/1‑69 (MQ=255)            tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCAGGCTTCCCAGc                                   >  1:3046796/1‑69 (MQ=255)            tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCAGGCTTCCCAGc                                   >  1:3386884/1‑69 (MQ=255)            tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCAGGCTTCCCAGc                                   >  1:515411/1‑69 (MQ=255)            tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCAGGCTTCCCAGc                                   >  1:2897481/1‑69 (MQ=255)            tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCAGGCTTCCCAGc                                   >  1:956486/1‑69 (MQ=255)            tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCAGGCTTCCCAGc                                   >  1:1345836/1‑69 (MQ=255)            tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCAGGCTTCCCAGc                                   >  1:1288411/1‑69 (MQ=255)            tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCAGGCTTCCCAGc                                   >  1:1284245/1‑69 (MQ=255)            tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCAGGCTTCCCAGc                                   >  1:1198952/1‑69 (MQ=255)            tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGGCCGCCCAGGCTTCCCAGc                                   >  1:1044333/1‑69 (MQ=255)                  tCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCa                                              >  1:3056880/1‑52 (MQ=255)                  tCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTc                             >  1:485668/1‑69 (MQ=255)                  tCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTc                             >  1:1171556/1‑69 (MQ=255)                   cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTc                             >  1:3307528/1‑68 (MQ=255)                        tCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGt                          <  1:3096975/66‑1 (MQ=255)                        tCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGt                          <  1:2921472/66‑1 (MQ=255)                        tCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGt                          <  1:590867/66‑1 (MQ=255)                        tCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGt                          <  1:1550748/66‑1 (MQ=255)                        tCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCACCAGGCTTCCCAGCACATTCTGt                          <  1:1497147/66‑1 (MQ=255)                         ccTGCTTTTTCCACTCGGTTTGGTTTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGt                          <  1:2887876/65‑1 (MQ=255)                                             tGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGc  >  1:1357932/1‑69 (MQ=255)                                                              |                                                    TGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGC  >  minE/1987966‑1988078

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