Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	668,428	C→A	T589K (ACG→AAG)	uup →	fused predicted transporter subunits and ATP‑binding components of ABC superfamily

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	668,428	0	C	A	100.0%	112.0 / NA	38	T589K (ACG→AAG)	uup	fused predicted transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (22/16);  total (22/16)

CAGCGCGAACTGGAGCAGCTACCGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGACGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGTGCTTGCTGATA  >  minE/668361‑668493                                                                    |                                                                  cAGCGCGAACTGGAGCAGCTACCGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGc                                                                 <  1:2113193/70‑1 (MQ=255) cAGCGCGAACTGGAGCAGCTACCGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGc                                                                 <  1:2038486/70‑1 (MQ=255)  aGCGCGAACTGGAGCAGCTACCGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCaa                                                               <  1:110991/71‑1 (MQ=255)   gcgcGAACTGGAGCAGCTACCGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAATCCCTACAGAAGCAAg                                                              <  1:1074870/71‑1 (MQ=255)    cgcgAACTGGAGCAGCTACCGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGt                                                             >  1:2454192/1‑71 (MQ=255)    agcgAACTGGAGCAGCTACCGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGt                                                             >  1:1072690/2‑71 (MQ=255)     gcgAACTGGAGCAGCTACCGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTg                                                            <  1:995853/71‑1 (MQ=255)     gcgAACTGGAGCAGCTACCGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTg                                                            <  1:2529374/71‑1 (MQ=255)             gAGCAGCTACCGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGc                                                    <  1:2056196/71‑1 (MQ=255)               gcagcTACCGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGct                                                   >  1:1649156/1‑70 (MQ=255)                      ccGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGCTTCcttctt                                           <  1:1037864/71‑1 (MQ=255)                      ccGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGCTTCcttctt                                           >  1:881389/1‑71 (MQ=255)                      ccGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGCTTCcttct                                            >  1:1956934/1‑70 (MQ=255)                      ccGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGCTTCcttct                                            <  1:376697/70‑1 (MQ=255)                      ccGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGCTTCcttct                                            <  1:2185933/70‑1 (MQ=255)                              gCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGcttcc                                                >  1:256711/1‑58 (MQ=255)                              gCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGcttcc                                                >  1:2347050/1‑58 (MQ=255)                              gCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGc                                   >  1:2910328/1‑71 (MQ=255)                                  gAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGcttc                                                 >  1:1120005/1‑53 (MQ=255)                                    aGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCAt                              >  1:1189476/1‑70 (MQ=255)                                       tCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGc                                   >  1:2773456/1‑62 (MQ=255)                                           gAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAg                        <  1:611141/69‑1 (MQ=255)                                           gAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAg                        <  1:2084607/69‑1 (MQ=255)                                                gAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCaaa                 <  1:2589485/71‑1 (MQ=255)                                                  aGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCaaa                 >  1:1008245/1‑69 (MQ=255)                                                      ggAAGCCCTACAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGtgc           >  1:351165/1‑71 (MQ=255)                                                          gCCCTACAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACTCAAAAAGtgcttgc       <  1:2073463/71‑1 (MQ=255)                                                              tACAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGTGCTTGCTGa    >  1:1481588/1‑70 (MQ=255)                                                              tACAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGTGCTTGCTGa    >  1:1870227/1‑70 (MQ=255)                                                              tACAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGTGCTTGCTGa    >  1:277765/1‑70 (MQ=255)                                                               aCAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCa                               >  1:2676054/1‑42 (MQ=255)                                                               aCAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCAt                              >  1:2172822/1‑43 (MQ=255)                                                               aCAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAg              >  1:2117841/1‑59 (MQ=255)                                                                cAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGtgct          >  1:1707107/1‑62 (MQ=255)                                                                cAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGTGCTTGCTGata  <  1:156143/70‑1 (MQ=255)                                                                 aGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCg                                 <  1:1858113/38‑1 (MQ=255)                                                                 aGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCaaaaa               >  1:2051165/1‑56 (MQ=255)                                                                 aGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGtgcttg        >  1:998300/1‑63 (MQ=255)                                                                    |                                                                  CAGCGCGAACTGGAGCAGCTACCGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGACGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGTGCTTGCTGATA  >  minE/668361‑668493

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