Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,379,816	C→T	intergenic (‑58/+37)	ubiE ← / ← rmuC	bifunctional 2‑octaprenyl‑6‑methoxy‑1,4‑benzoquinone methylase and S‑adenosylmethionine:2‑DMK methyltransferase/predicted recombination limiting protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,379,816	0	C	T	100.0%	74.8 / NA	24	intergenic (‑58/+37)	ubiE/rmuC	bifunctional 2‑octaprenyl‑6‑methoxy‑1,4‑benzoquinone methylase and S‑adenosylmethionine:2‑DMK methyltransferase/predicted recombination limiting protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/12);  total (12/12)

CATCAAAAAATTGTTCCAGAAGTGTAACAGATTGGGCGTCGATGCCCTAGATTTCTACCCGGCTTAACTACTCCCAATGGGCTAGCGCGACTGCTGATTATATTCAT  >  minE/2379772‑2379878                                             |                                                               catcaAAAAATTGTTCCAGAAGTGTAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAactac                                      >  1:2984527/1‑71 (MQ=255)     aaaaaaTTGTTCCAGAAGTGTAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTcc                                   >  1:616917/1‑70 (MQ=255)     aaaaaaTTGTTCCAGAAGTGTAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTcc                                   >  1:3080452/1‑70 (MQ=255)     aaaaaaTTGGTCCAGAAGTGTAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTcc                                   >  1:4740125/1‑70 (MQ=255)        aaaTTGTTCCAGAAGTGTAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTc                                    >  1:3738079/1‑66 (MQ=255)        aaaTTGTTCCAGAAGTGTAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCAAt                               >  1:1446956/1‑71 (MQ=255)        aaaTTGTTCCAGAAGTGTAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCAAt                               >  1:3735214/1‑71 (MQ=255)                ccAGAAGTGTAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCGATGGGCTAg                        <  1:709948/70‑1 (MQ=255)                ccAGAAGTGTAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCAATGGGCTAg                        <  1:1891777/70‑1 (MQ=255)                ccAGAAGTGTAACAGATAGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCAATGGGCTAg                        <  1:4019153/70‑1 (MQ=255)                  aGAAGTGTAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCAATGGGCTAgcgc                     <  1:239944/71‑1 (MQ=255)                  aGAAGTGTAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCAATGGGCTAgcgc                     <  1:219010/71‑1 (MQ=255)                  aGAAGTGTAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCAATGGGCTAgcgc                     <  1:571568/71‑1 (MQ=255)                  aGAAGTGTAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCAATGGGCTAgcgc                     <  1:2064635/71‑1 (MQ=255)                  aGAAGTGTAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCAATGGGCTAgcgc                     <  1:1931532/71‑1 (MQ=255)                         tAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCAATGGGCTAGCGCGActgct              <  1:3015743/71‑1 (MQ=255)                         tAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCAATGGGCTAGCGCGActgct              <  1:3449953/71‑1 (MQ=255)                         tAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCAATGGGCTAGCGCGActgct              <  1:2520962/71‑1 (MQ=255)                         tAACAGATTGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCAATGGGCTAGCGCGActgct              <  1:4144635/71‑1 (MQ=255)                                 tGGGCGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCAATGGGCTAGCGCGACTGCTGATTatat      >  1:427967/1‑71 (MQ=255)                                 tGGGAGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCAATGGGCTAGCGCGACTGCTGATTatat      >  1:58821/1‑71 (MQ=255)                                     cGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTGCCAATGGGCTAGCGCGACTGCTGATTATATtcat  >  1:3454655/1‑71 (MQ=255)                                     cGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCAATGGGCTAGCGCGACTGCTGATTATATtcat  >  1:3195739/1‑71 (MQ=255)                                     cGTCGATGTCCTAGATTTCTACCCGGCTTAACTACTCCCAATGGGCTAGCGCGACTGCTGATTATATtcat  >  1:3332966/1‑71 (MQ=255)                                             |                                                               CATCAAAAAATTGTTCCAGAAGTGTAACAGATTGGGCGTCGATGCCCTAGATTTCTACCCGGCTTAACTACTCCCAATGGGCTAGCGCGACTGCTGATTATATTCAT  >  minE/2379772‑2379878

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