Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,695,295	G→A	Q239* (CAG→TAG)	yfhQ ←	predicted methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,695,295	0	G	A	100.0%	242.6 / NA	73	Q239* (CAG→TAG)	yfhQ	predicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (64/9);  total (64/9)

CGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGC  >  minE/1695239‑1695351                                                         |                                                         cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:1849614/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:2595288/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:2526491/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:2303425/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:2252500/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:2673935/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:2220864/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:272956/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:2143164/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:3024001/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:1992734/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:3033536/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:1942307/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:3092312/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:3107017/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:436530/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:454115/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:3193843/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:553209/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:559544/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:784255/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:1322485/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:1279139/1‑71 (MQ=255) cGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAg                                            >  1:569562/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:959163/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:809868/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:2612197/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:702882/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:700066/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:2706663/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:2830506/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:2924598/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:2953214/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:48865/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:48626/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:2596698/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:1711428/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:1051961/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:1129553/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:1214636/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:1673582/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:1761161/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                      >  1:2079121/1‑71 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGa                                       >  1:2896888/1‑70 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGa                                       >  1:199944/1‑70 (MQ=255)       aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGa                                       >  1:902465/1‑70 (MQ=255)                       aaCGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGatat                       >  1:3151581/1‑70 (MQ=255)                       aaCGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATAtt                      >  1:2704424/1‑71 (MQ=255)                           gCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAAc                  <  1:2122288/71‑1 (MQ=255)                           gCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAAc                  <  1:2477105/71‑1 (MQ=255)                           gCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAAc                  <  1:930425/71‑1 (MQ=255)                           gCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAAc                  <  1:83852/71‑1 (MQ=255)                           gCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAAc                  <  1:581715/71‑1 (MQ=255)                           gCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAAc                  <  1:1256017/71‑1 (MQ=255)                           gCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAAc                  <  1:1395362/71‑1 (MQ=255)                           gCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAAc                  <  1:1824577/71‑1 (MQ=255)                           gCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAAc                  <  1:2148797/71‑1 (MQ=255)                                    cGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACt                 >  1:459946/1‑63 (MQ=255)                                    cGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACt                 >  1:2033098/1‑63 (MQ=255)                                    cGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACt                 >  1:3011549/1‑63 (MQ=255)                                    cGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACt                 >  1:1975470/1‑63 (MQ=255)                                    cGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACt                 >  1:1085143/1‑63 (MQ=255)                                    cGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACt                 >  1:1097744/1‑63 (MQ=255)                                    cGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACt                 >  1:392348/1‑63 (MQ=255)                                    cGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACt                 >  1:505404/1‑63 (MQ=255)                                    cGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACt                 >  1:1559603/1‑63 (MQ=255)                                    cGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACt                 >  1:1623167/1‑63 (MQ=255)                                    cGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACt                 >  1:2180794/1‑63 (MQ=255)                                           gTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCgg    >  1:1041270/1‑69 (MQ=255)                                           gTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgc  >  1:1808159/1‑71 (MQ=255)                                           gTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgc  >  1:2243220/1‑71 (MQ=255)                                           gTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgc  >  1:1502727/1‑71 (MQ=255)                                           gTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgc  >  1:815282/1‑71 (MQ=255)                                                         |                                                         CGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGC  >  minE/1695239‑1695351

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