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

CCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGC  >  minE/1695229‑1695355                                                                   |                                                             ccTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:1167511/71‑1 (MQ=255) ccTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:2934092/71‑1 (MQ=255) ccTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:2864032/71‑1 (MQ=255) ccTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:2355038/71‑1 (MQ=255) ccTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:1553074/71‑1 (MQ=255)  cTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:1237486/70‑1 (MQ=255)  cTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:524796/70‑1 (MQ=255)  cTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:2076397/70‑1 (MQ=255)                 aCGCGTTGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                          >  1:179802/1‑71 (MQ=255)                 aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                          >  1:179323/1‑71 (MQ=255)                 aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                          >  1:2094377/1‑71 (MQ=255)                 aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGaa                                          >  1:3338682/1‑71 (MQ=255)                   gcgTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:3320862/53‑1 (MQ=255)                    cgTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:1415911/52‑1 (MQ=255)                           tCAGGCAACGGCTGTTATTCAGCCTTGTTACCTTTATTCtactgc                                                          <  1:162582/45‑1 (MQ=38)                             aGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:4048747/43‑1 (MQ=255)                             aGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:2219561/43‑1 (MQ=255)                               gCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:3494190/41‑1 (MQ=255)                               gCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:894494/41‑1 (MQ=255)                                cAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGatat                           >  1:4037748/1‑71 (MQ=255)                                cAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGatat                           >  1:972673/1‑71 (MQ=255)                                cAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGatat                           >  1:3890131/1‑71 (MQ=255)                                cAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGatat                           >  1:3562439/1‑71 (MQ=255)                                cAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGatat                           >  1:938262/1‑71 (MQ=255)                                cAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGatat                           >  1:3309348/1‑71 (MQ=255)                                cAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGatat                           >  1:1211951/1‑71 (MQ=255)                                  aCGGCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:1645802/38‑1 (MQ=255)                                   cGGCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:3504872/37‑1 (MQ=255)                                   cGGCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:2976491/37‑1 (MQ=255)                                     gCTGTTATTCGGCCTTGTTACCTTTATTCtactgc                                                          <  1:587898/35‑1 (MQ=255)                                         ttATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGa                              <  1:851354/59‑1 (MQ=255)                                         ttATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGa                              <  1:3476080/59‑1 (MQ=255)                                         ttATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGa                              <  1:2765277/59‑1 (MQ=255)                                         ttATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGa                              <  1:3533227/59‑1 (MQ=255)                                         ttATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGa                              <  1:4225419/59‑1 (MQ=255)                                         ttATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGa                              <  1:418376/59‑1 (MQ=255)                                         ttATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGa                              <  1:4140065/59‑1 (MQ=255)                                         ttATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGa                              <  1:3659062/59‑1 (MQ=255)                                         ttATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGa                              <  1:917084/59‑1 (MQ=255)                                         ttATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGa                              <  1:3897605/59‑1 (MQ=255)                                         ttATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGa                              <  1:4024780/59‑1 (MQ=255)                                           aTTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCttg                >  1:1723566/1‑69 (MQ=255)                                           aTTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTgg                >  1:804830/1‑71 (MQ=255)                                           aTTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTgg                >  1:3253810/1‑71 (MQ=255)                                           aTTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTgg                >  1:3106268/1‑71 (MQ=255)                                           aTTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTgg                >  1:3058672/1‑71 (MQ=255)                                           aTTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTg                 >  1:3648919/1‑70 (MQ=255)                                               ggCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCttt            <  1:3140452/71‑1 (MQ=255)                                               ggCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCttt            <  1:1248034/71‑1 (MQ=255)                                               ggCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCttt            <  1:660218/71‑1 (MQ=255)                                               ggCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCttt            <  1:1388704/71‑1 (MQ=255)                                               ggCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCttt            <  1:4329054/71‑1 (MQ=255)                                               ggCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCttt            <  1:517762/71‑1 (MQ=255)                                                gCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCTCAGGATATTCAACTCCTGGCttt            <  1:1174700/70‑1 (MQ=255)                                                  cTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCg         >  1:3335520/1‑71 (MQ=255)                                                  cTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCg         >  1:3318384/1‑71 (MQ=255)                                                  cTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCg         >  1:3897726/1‑71 (MQ=255)                                                        aCCTTTATTCTACTGCTAAATAGAAGCCAGAATCCCGCTCATGa                              <  1:1880650/44‑1 (MQ=25)                                                         ccTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgcgcgc  <  1:3982546/71‑1 (MQ=255)                                                         ccTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgcgcg   <  1:4297715/70‑1 (MQ=255)                                                         ccTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgcgcg   <  1:4194741/70‑1 (MQ=255)                                                         ccTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgcgcg   <  1:2537571/70‑1 (MQ=255)                                                         ccTTTATTCTACTGCTCAATAGAAGACAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgcgcgc  <  1:3561557/71‑1 (MQ=255)                                                                   |                                                             CCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGC  >  minE/1695229‑1695355

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