Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,996,170	G→A	intergenic (‑21/‑377)	speC ← / → yqgA	ornithine decarboxylase, constitutive/predicted inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,996,170	0	G	A	100.0%	325.4 / NA	93	intergenic (‑21/‑377)	speC/yqgA	ornithine decarboxylase, constitutive/predicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (47/46);  total (47/46)

GGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATACCGGACAAGGTGGAAAACCCTGCCCC  >  minE/1996126‑1996207                                             |                                      ggcggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCgg                        >  1:477083/1‑60 (MQ=255) ggcggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCgg                        >  1:1648108/1‑60 (MQ=255) ggcggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCgg                        >  1:3386990/1‑60 (MQ=255) ggcggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCgg                        >  1:3297296/1‑60 (MQ=255) ggcggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCgg                        >  1:1524259/1‑60 (MQ=255) ggcggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCgg                        >  1:1187864/1‑60 (MQ=255) ggcggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCgg                        >  1:1289588/1‑60 (MQ=255) ggcggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCgg                        >  1:2016030/1‑60 (MQ=255) ggcggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCgg                        >  1:265673/1‑60 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGaa                              >  1:1841950/1‑51 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACaa                    >  1:2344892/1‑61 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAg                   >  1:1641539/1‑62 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAg                   >  1:1027275/1‑62 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAg                   >  1:1860696/1‑62 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAg                   >  1:18730/1‑62 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAg                   >  1:2228624/1‑62 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAg                   >  1:2418103/1‑62 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAg                   >  1:2842516/1‑62 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAg                   >  1:301983/1‑62 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAg                   >  1:3048296/1‑62 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAg                   >  1:3151892/1‑62 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAg                   >  1:317397/1‑62 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAg                   >  1:3437450/1‑62 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAg                   >  1:3547896/1‑62 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAg                   >  1:1474579/1‑62 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAg                   >  1:1344484/1‑62 (MQ=255)    ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAg                   >  1:1119516/1‑62 (MQ=255)    ggcAATATTCATTGATTTCATAAACGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAg                   >  1:2041557/1‑62 (MQ=255)                  aTTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCt       >  1:1708835/1‑60 (MQ=255)                  aTTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCTGc     >  1:1568663/1‑62 (MQ=255)                  aTTTCATAAGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCTGc     >  1:2540701/1‑62 (MQ=255)                         aaGCGCAAACCCGTTTCGGGATAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:3485223/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:3322618/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:2804303/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:3300960/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:1103911/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:1227459/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:3130653/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:1236335/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:2979434/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:1314382/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:2815161/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:979927/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:1083690/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:3422238/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:343606/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:1073488/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:3607289/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:3630127/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:402600/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:434443/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:1019118/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:502292/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:578360/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:665948/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:1608408/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:2281515/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:1567074/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:1509829/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:1452547/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:2006945/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:1430392/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:2022993/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:2197255/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:2206340/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:2697985/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:2319025/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:2618421/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:2577128/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:2393237/51‑1 (MQ=255)                         aaGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:2560304/51‑1 (MQ=255)                          aGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:1376263/50‑1 (MQ=255)                          aGCGCAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:1057547/50‑1 (MQ=255)                           gcgcAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:954620/49‑1 (MQ=255)                             gcAAACCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:1405855/47‑1 (MQ=255)                                aaCCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCTGcccc  >  1:2059936/1‑51 (MQ=255)                                aaCCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCTGcccc  >  1:2664276/1‑51 (MQ=255)                                aaCCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCTGcccc  >  1:652622/1‑51 (MQ=255)                                aaCCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCTGcccc  >  1:607109/1‑51 (MQ=255)                                aaCCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCTGcccc  >  1:593772/1‑51 (MQ=255)                                aaCCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCTGcccc  >  1:2904379/1‑51 (MQ=255)                                aaCCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCTGcccc  >  1:1999146/1‑51 (MQ=255)                                aaCCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCTGcccc  >  1:3537317/1‑51 (MQ=255)                                aaCCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCTGcccc  >  1:2535991/1‑51 (MQ=255)                                aaCCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCTGcccc  >  1:2521445/1‑51 (MQ=255)                                aaCCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCTGcccc  >  1:343198/1‑51 (MQ=255)                                aaCCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCTGcccc  >  1:2286061/1‑51 (MQ=255)                                aaCCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCTGcccc  >  1:338915/1‑51 (MQ=255)                                aaCCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCTGcccc  >  1:1383306/1‑51 (MQ=255)                                aaCCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAACCCTGccc   >  1:2291450/1‑50 (MQ=255)                                aaCCCGTTTCGGGAAAGTAAGAATACCGGACAAGGTGG‑‑AACCCTGcccc  >  1:590857/1‑49 (MQ=255)                                  cccGTTTCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:3063320/42‑1 (MQ=255)                                      tttCGGGAAAGTAAGAATACCGGACAAGGTGGAAAAcc         <  1:1294736/38‑1 (MQ=255)                                             |                                      GGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATACCGGACAAGGTGGAAAACCCTGCCCC  >  minE/1996126‑1996207

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