Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,451,978	Δ1 bp	coding (318/651 nt)	atoA →	acetyl‑CoA:acetoacetyl‑CoA transferase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,451,977	0	A	.	100.0%	60.3 / NA	15	coding (317/651 nt)	atoA	acetyl‑CoA:acetoacetyl‑CoA transferase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (0/15);  total (0/15)

GCTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACGAAGAAGCAAACCTCGCGAACTGGGTAGTGCCTGGGAAAATGGTGCCCGGTATGGGTGGC  >  minE/1451918‑1452032                                                            |                                                        gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:367059/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1647081/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1735944/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1776283/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1793405/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1838332/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1841850/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:229979/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:245162/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:322212/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1736407/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:42955/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:437297/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:518047/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:522094/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:641388/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:810237/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:811153/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:860491/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1590562/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1226955/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1231896/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:154786/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1241306/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1499585/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1469185/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1459196/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:133860/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1387974/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1436803/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1406263/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaagaa                                                        >  1:1607830/1‑61 (MQ=255) gcTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACgaaga                                                         >  1:354972/1‑60 (MQ=255)                                                     gACGAAG‑AGCAAACCTCGCGAACTGGGTAGTGCCTGGGAAAATGGTGCCCGGTATGGGTGgc  <  1:1392199/62‑1 (MQ=255)                                                     gACGAAG‑AGCAAACCTCGCGAACTGGGTAGTGCCTGGGAAAATGGTGCCCGGTATGGGTGgc  <  1:358735/62‑1 (MQ=255)                                                     gACGAAG‑AGCAAACCTCGCGAACTGGGTAGTGCCTGGGAAAATGGTGCCCGGTATGGGTGgc  <  1:1868504/62‑1 (MQ=255)                                                     gACGAAG‑AGCAAACCTCGCGAACTGGGTAGTGCCTGGGAAAATGGTGCCCGGTATGGGTGgc  <  1:1387264/62‑1 (MQ=255)                                                     gACGAAG‑AGCAAACCTCGCGAACTGGGTAGTGCCTGGGAAAATGGTGCCCGGTATGGGTGgc  <  1:1378072/62‑1 (MQ=255)                                                     gACGAAG‑AGCAAACCTCGCGAACTGGGTAGTGCCTGGGAAAATGGTGCCCGGTATGGGTGgc  <  1:1445848/62‑1 (MQ=255)                                                     gACGAAG‑AGCAAACCTCGCGAACTGGGTAGTGCCTGGGAAAATGGTGCCCGGTATGGGTGgc  <  1:1016325/62‑1 (MQ=255)                                                     gACGAAG‑AGCAAACCTCGCGAACTGGGTAGTGCCTGGGAAAATGGTGCCCGGTATGGGTGgc  <  1:1730279/62‑1 (MQ=255)                                                     gACGAAG‑AGCAAACCTCGCGAACTGGGTAGTGCCTGGGAAAATGGTGCCCGGTATGGGTGgc  <  1:726458/62‑1 (MQ=255)                                                     gACGAAG‑AGCAAACCTCGCGAACTGGGTAGTGCCTGGGAAAATGGTGCCCGGTATGGGTGgc  <  1:78588/62‑1 (MQ=255)                                                     gACGAAG‑AGCAAACCTCGCGAACTGGGTAGTGCCTGGGAAAATGGTGCCCGGTATGGGTGgc  <  1:1656667/62‑1 (MQ=255)                                                     gACGAAG‑AGCAAACCTCGCGAACTGGGTAGTGCCTGGGAAAATGGTGCCCGGTATGGGTGgc  <  1:1191430/62‑1 (MQ=255)                                                     gACGAAG‑AGCAAACCTCGCGAACTGGGTAGTGCCTGGGAAAATGGTGCCCGGTATGGGTGgc  <  1:811780/62‑1 (MQ=255)                                                     gACGAAG‑AGCAAACCTCGCGAACTGGGTAGTGCCTGGGAAAATGGTGCCCGGTATGGGTGgc  <  1:856056/62‑1 (MQ=255)                                                     gACGAAG‑AGCAAACCTCGCGAACTGGGTAGTGCCTGGGAAAATGGTGCCCGGTATGGGTGgc  <  1:1185071/62‑1 (MQ=255)                                                            |                                                        GCTAATCCGTGGCGGTCATATTGATGCCTGCGTGCTCGGCGGTTTGCAAGTAGACGAAGAAGCAAACCTCGCGAACTGGGTAGTGCCTGGGAAAATGGTGCCCGGTATGGGTGGC  >  minE/1451918‑1452032

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