Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	454,992	C→T	10.4%	L31L (CTG→CTA) ‡	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	454,992	0	C	T	10.4%	229.6 / 13.5	77	L31L (CTG→CTA) ‡	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (33/35);  new base T (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 = 9.97e-01

CGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCG  >  minE/454928‑455058                                                                 |                                                                   cgcTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                   >  1:2224203/1‑66 (MQ=255)   cTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:61552/71‑1 (MQ=255)   cTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:2228850/71‑1 (MQ=255)        tAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGc                                                       <  1:4097912/71‑1 (MQ=255)        tAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGc                                                       <  1:1913923/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)              gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                   >  1:640421/1‑53 (MQ=255)              gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaa                                                  <  1:2377919/70‑1 (MQ=255)              gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTa                                                   <  1:136861/69‑1 (MQ=255)              gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTa                                                   <  1:2656889/69‑1 (MQ=255)              gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTa                                                   <  1:756171/69‑1 (MQ=255)                 ggCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGc                                              >  1:3212056/1‑71 (MQ=255)                    gTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                   >  1:3908257/1‑47 (MQ=255)                    gTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                   >  1:1883470/1‑47 (MQ=255)                       aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                   >  1:3398756/1‑44 (MQ=255)                       aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:3381634/51‑1 (MQ=255)                       aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:4046149/51‑1 (MQ=255)                       aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:1336543/51‑1 (MQ=255)                       aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                         <  1:1583992/70‑1 (MQ=255)                       aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                         >  1:3014661/1‑70 (MQ=255)                       aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                         >  1:272347/1‑70 (MQ=255)                       aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                         <  1:700930/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)                           tCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                   >  1:1014040/1‑40 (MQ=255)                           tCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:1134709/47‑1 (MQ=255)                           tCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTc                                    >  1:1368659/1‑71 (MQ=255)                             gTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                   >  1:3745011/1‑38 (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:775266/37‑2 (MQ=19)                                                     aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACtg                                            <  1:3218011/37‑2 (MQ=19)                                                     aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACtg                                            <  1:1679394/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:581079/1‑70 (MQ=255)                                                         agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa       >  1:708108/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:3549167/37‑1 (MQ=255)                                                             tCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTc                                    <  1:315200/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:1355340/1‑58 (MQ=255)                                                               agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACAt             >  1:1588407/1‑58 (MQ=255)                                                               agcagGGCACGCATGCCGTAAAAGCTAAAATATTCACAGATTTGCa                         >  1:4076241/1‑46 (MQ=255)                                                                 cagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATa            >  1:1383647/1‑57 (MQ=255)                                                                 |                                                                   CGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCG  >  minE/454928‑455058

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑