Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	454,994	G→A	10.4%	L31L (CTG→TTG) ‡	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	454,994	0	G	A	10.4%	226.9 / 13.5	77	L31L (CTG→TTG) ‡	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (33/36);  new base A (4/4);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-01

CTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCG  >  minE/454930‑455058                                                                 |                                                                 cTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:61552/71‑1 (MQ=255) cTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:2228850/71‑1 (MQ=255)      tAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGc                                                       <  1:1913923/71‑1 (MQ=255)      tAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGc                                                       <  1:4097912/71‑1 (MQ=255)       aaCTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGcc                                                      >  1:3240464/1‑71 (MQ=255)          tCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCg                                                     >  1:3214852/1‑69 (MQ=255)          tCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGt                                                    >  1:3580465/1‑70 (MQ=255)           cGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTa                                                   >  1:2942984/1‑70 (MQ=255)            gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaa                                                  <  1:2377919/70‑1 (MQ=255)            gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTa                                                   <  1:756171/69‑1 (MQ=255)            gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTa                                                   <  1:136861/69‑1 (MQ=255)            gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTa                                                   <  1:2656889/69‑1 (MQ=255)               ggCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGc                                              >  1:3212056/1‑71 (MQ=255)                     aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:4046149/51‑1 (MQ=255)                     aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:1336543/51‑1 (MQ=255)                     aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:3381634/51‑1 (MQ=255)                     aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                         >  1:3014661/1‑70 (MQ=255)                     aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                         >  1:272347/1‑70 (MQ=255)                     aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                         <  1:700930/70‑1 (MQ=255)                     aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                         <  1:1583992/70‑1 (MQ=255)                     aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAat                                        <  1:704828/71‑1 (MQ=255)                      gTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                       >  1:1330484/1‑71 (MQ=255)                      gTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                       >  1:4131068/1‑71 (MQ=255)                       tATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:1607978/49‑1 (MQ=255)                         tCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:1134709/47‑1 (MQ=255)                         tCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTc                                    >  1:1368659/1‑71 (MQ=255)                           gTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTcc                                   >  1:2573458/1‑70 (MQ=255)                            ttGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:55879/44‑1 (MQ=255)                               tATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGc                                                               <  1:3343843/38‑1 (MQ=255)                                aTTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                       >  1:3959543/1‑61 (MQ=255)                                aTTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                       >  1:2968897/1‑61 (MQ=255)                                 ttttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaa                                          >  1:2994322/1‑57 (MQ=255)                                  tttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                           >  1:3699829/1‑71 (MQ=255)                                   ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCa                                                          >  1:1472529/1‑39 (MQ=255)                                   ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGc                                              >  1:1468493/1‑51 (MQ=255)                                   ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                           >  1:1661899/1‑70 (MQ=255)                                   ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                           <  1:3326154/70‑1 (MQ=255)                                   ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                           <  1:4121510/70‑1 (MQ=255)                                    tAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                           >  1:1307684/1‑69 (MQ=255)                                      gTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTATAAGCTAAAata                                       <  1:3876936/55‑1 (MQ=255)                                         gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                       >  1:960131/1‑52 (MQ=255)                                         gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcgc                     >  1:2511377/1‑70 (MQ=255)                                           ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTa                                                   <  1:2313457/38‑1 (MQ=255)                                           ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAAc                  <  1:1794726/71‑1 (MQ=255)                                             ggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTaa                                                  >  1:3891473/1‑37 (MQ=37)                                              gTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGa                <  1:1915790/70‑1 (MQ=255)                                               tGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGa                               <  1:2368093/54‑1 (MQ=255)                                               tGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGAc               >  1:1052665/1‑70 (MQ=255)                                               tGAGATAGAGAATCAG‑AGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGa                               <  1:206121/53‑1 (MQ=255)                                                 agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGata              >  1:3795297/1‑67 (MQ=255)                                                 agaTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCATCGACa              <  1:819345/69‑1 (MQ=255)                                                  gaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata            <  1:3559243/70‑1 (MQ=255)                                                   aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACtg                                            <  1:1679394/37‑2 (MQ=19)                                                   aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACtg                                            <  1:3218011/37‑2 (MQ=19)                                                   aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACtg                                            <  1:775266/37‑2 (MQ=19)                                                   aTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAg           >  1:1310295/1‑70 (MQ=255)                                                   aTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGt          <  1:1541908/71‑1 (MQ=255)                                                    tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTcc                                   <  1:2200551/45‑1 (MQ=255)                                                    tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa         <  1:3613702/71‑1 (MQ=255)                                                    tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa         <  1:9308/71‑1 (MQ=255)                                                    tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa         <  1:83469/71‑1 (MQ=255)                                                     agagAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGAtcgcca                   >  1:352824/1‑58 (MQ=255)                                                     agagAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAat                                        >  1:2821213/1‑39 (MQ=255)                                                     agagAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgc                       >  1:1399834/1‑56 (MQ=255)                                                     agagAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa         >  1:2988942/1‑70 (MQ=255)                                                      gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGc                          >  1:3460408/1‑52 (MQ=255)                                                       agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa       >  1:708108/1‑70 (MQ=255)                                                       agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa       >  1:581079/1‑70 (MQ=255)                                                        gaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatatc     >  1:3956993/1‑70 (MQ=255)                                                        gaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACAt             >  1:822525/1‑63 (MQ=255)                                                         aaTCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGAc               >  1:1784888/1‑60 (MQ=255)                                                          aTCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAg    <  1:1540939/70‑1 (MQ=255)                                                           tCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTc                                    <  1:315200/37‑1 (MQ=255)                                                           tCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTc                                    <  1:3549167/37‑1 (MQ=255)                                                           tCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTGCCAGATTTGCAGCGCAACGACATAGTAAATAgcg  >  1:1085761/1‑71 (MQ=255)                                                           tCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAAt      <  1:3237565/67‑1 (MQ=255)                                                           tCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcg  >  1:72838/1‑71 (MQ=255)                                                             agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCa                         >  1:3692836/1‑46 (MQ=255)                                                             agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACAt             >  1:1588407/1‑58 (MQ=255)                                                             agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACAt             >  1:1355340/1‑58 (MQ=255)                                                             agcagGGCACGCATGCCGTAAAAGCTAAAATATTCACAGATTTGCa                         >  1:4076241/1‑46 (MQ=255)                                                               cagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATa            >  1:1383647/1‑57 (MQ=255)                                                                 gggCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa              >  1:3761063/1‑53 (MQ=255)                                                                 |                                                                 CTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCG  >  minE/454930‑455058

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