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%	106.3 / NA	36	Q239* (CAG→TAG)	yfhQ	predicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/21);  total (15/21)

GGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGT  >  minE/1695227‑1695361                                                                     |                                                                   ggCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctact                                                                  >  1:2392975/1‑71 (MQ=255)       aTAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCaa                                                            >  1:1490966/1‑71 (MQ=255)       aTAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCa                                                             <  1:2005172/70‑1 (MQ=255)           gACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAg                                                         >  1:991446/1‑70 (MQ=255)           gACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAg                                                         >  1:2027488/1‑70 (MQ=255)           gACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGa                                                        <  1:78373/71‑1 (MQ=255)              gCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGc                                                     <  1:1763265/71‑1 (MQ=255)              gCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGc                                                     <  1:330064/71‑1 (MQ=255)                 ttACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCa                                                   <  1:284234/70‑1 (MQ=255)                 ttACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAg                                                  <  1:3020052/71‑1 (MQ=255)                  tACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAg                                                  >  1:2925050/1‑70 (MQ=255)                  tACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAg                                                  >  1:2888929/1‑70 (MQ=255)                   aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATa                                                          >  1:1483112/1‑61 (MQ=255)                   aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                                >  1:2167469/1‑71 (MQ=255)                   aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGa                                                 <  1:2372249/70‑1 (MQ=255)                       gTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGcc                                                    >  1:2252470/1‑63 (MQ=255)                       gTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGcc                                                    >  1:2386551/1‑63 (MQ=255)                                 gCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAg                                                         <  1:2484373/48‑1 (MQ=255)                                   aaCGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAgg                                     <  1:1008848/66‑1 (MQ=255)                                    aCGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCa                                                   >  1:77124/1‑51 (MQ=255)                                           ttATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCt                        <  1:168646/71‑1 (MQ=255)                                           ttATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGAATATTCAACTCCt                        >  1:2121485/1‑71 (MQ=255)                                             aTTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATTTTCAACTCCTgg                      >  1:2223534/1‑71 (MQ=255)                                                cGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCt                    <  1:1393582/70‑1 (MQ=255)                                                 ggCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCttt                  <  1:1807110/71‑1 (MQ=255)                                                 ggCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCttt                  <  1:429389/71‑1 (MQ=255)                                                  gCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTc                 <  1:2849424/71‑1 (MQ=255)                                                   ccTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTcc                <  1:2971369/71‑1 (MQ=255)                                                   ccTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTcc                <  1:2999542/71‑1 (MQ=255)                                                   ccTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTcc                >  1:387119/1‑71 (MQ=255)                                                      tGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCtt                   <  1:1782123/65‑1 (MQ=255)                                                       gTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCgcgc                                        >  1:1923412/1‑43 (MQ=255)                                                       gTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCgg              <  1:140912/69‑1 (MQ=255)                                                             tttATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgcgc          <  1:2312185/67‑1 (MQ=255)                                                               tATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCAc      <  1:2316169/69‑1 (MQ=255)                                                               tATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCAc      <  1:1674644/69‑1 (MQ=255)                                                                  tCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGt  >  1:438888/1‑70 (MQ=255)                                                                     |                                                                   GGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGT  >  minE/1695227‑1695361

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