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

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

GCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTG  >  minE/1695573‑1695696                                                              |                                                               gCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCtt                                                        <  1:2935748/70‑1 (MQ=255) gCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCtt                                                        <  1:2377829/70‑1 (MQ=255) gCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCtt                                                        <  1:1953106/70‑1 (MQ=255) gCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCtt                                                        <  1:258037/70‑1 (MQ=255) gCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCtt                                                        <  1:63615/70‑1 (MQ=255) gCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCtt                                                        <  1:1151674/70‑1 (MQ=255)  ccAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCTTCa                                                      <  1:2522891/71‑1 (MQ=255)  ccAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCTTCa                                                      <  1:650191/71‑1 (MQ=255)   cAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCTTCa                                                      <  1:304623/70‑1 (MQ=255)             gAGCTGTATTCCGGGTTAGCCGCAATCGCGAGATGATAATGGCATTTCTACAACTCTTCa                                                      <  1:2858006/60‑1 (MQ=255)                       ccGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTcg                                 >  1:2819441/1‑71 (MQ=255)                       ccGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTcg                                 >  1:2182774/1‑71 (MQ=255)                              aGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGa                                         >  1:4716054/1‑56 (MQ=255)                              aGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGa                                         >  1:4653393/1‑56 (MQ=255)                              aGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGa                                         >  1:407351/1‑56 (MQ=255)                              aGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGa                                         >  1:3656892/1‑56 (MQ=255)                              aGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGa                                         >  1:759741/1‑56 (MQ=255)                              aGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGa                                         >  1:2895433/1‑56 (MQ=255)                              aGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGa                                         >  1:1665404/1‑56 (MQ=255)                              aGCCACAATCGCGACATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGa                                         >  1:148510/1‑56 (MQ=255)                                        gcgACATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAgcgc                >  1:1128348/1‑71 (MQ=255)                                        gcgACATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCAGACGCGCTCGCGACCAAACACCAgcgc                >  1:650376/1‑71 (MQ=39)                                            cATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACAcc                     >  1:1684507/1‑62 (MQ=255)                                            cATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCAcc            >  1:371165/1‑71 (MQ=255)                                            cATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCAcc            >  1:1009067/1‑71 (MQ=255)                                            cATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCAcc            >  1:4863785/1‑71 (MQ=255)                                            cATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCAcc            >  1:1167241/1‑71 (MQ=255)                                            cATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCAcc            >  1:3487221/1‑71 (MQ=255)                                            cATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCAcc            >  1:3166431/1‑71 (MQ=255)                                            cATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCAcc            >  1:2190760/1‑71 (MQ=255)                                             aTGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCg           >  1:1292407/1‑71 (MQ=255)                                             aTGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCg           >  1:507135/1‑71 (MQ=255)                                               gATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCggg         >  1:2332686/1‑71 (MQ=255)                                               gATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCggg         >  1:1782149/1‑71 (MQ=255)                                               gATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCggg         >  1:1166037/1‑71 (MQ=255)                                                   aTGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTa      <  1:1387880/70‑1 (MQ=255)                                                   aTGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTa      <  1:19136/70‑1 (MQ=255)                                                      gaaTTTCTACAACTCTTCATTGGTCAGGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGGATTTg  >  1:4414669/3‑71 (MQ=38)                                                      gCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTg  >  1:3274310/1‑71 (MQ=255)                                                      gCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTg  >  1:4310799/1‑71 (MQ=255)                                                      gCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTg  >  1:1563063/1‑71 (MQ=255)                                                      gCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTg  >  1:1914977/1‑71 (MQ=255)                                                              |                                                               GCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTG  >  minE/1695573‑1695696

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