Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	455,013	G→A	10.2%	S24S (AGC→AGT)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	455,013	0	G	A	10.2%	295.2 / 17.5	98	S24S (AGC→AGT)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (50/38);  new base A (6/4);  total (56/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

GGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCA  >  minE/454944‑455083                                                                      |                                                                       ggCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGc                                                                       >  1:3212056/1‑71 (MQ=255)       aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                                                  <  1:1583992/70‑1 (MQ=255)       aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                                                  >  1:272347/1‑70 (MQ=255)       aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                                                  >  1:3014661/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)           tCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTc                                                             >  1:1368659/1‑71 (MQ=255)             gTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTcc                                                            >  1:2573458/1‑70 (MQ=255)                  aTTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                                                >  1:2968897/1‑61 (MQ=255)                  aTTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                                                >  1:3959543/1‑61 (MQ=255)                   ttttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaa                                                                   >  1:2994322/1‑57 (MQ=255)                    tttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                                                    >  1:3699829/1‑71 (MQ=255)                     ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGc                                                                       >  1:1468493/1‑51 (MQ=255)                     ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                                                    <  1:4121510/70‑1 (MQ=255)                     ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                                                    <  1:3326154/70‑1 (MQ=255)                     ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                                                    >  1:1661899/1‑70 (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)                             ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAAc                                           <  1:1794726/71‑1 (MQ=255)                                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:83469/71‑1 (MQ=255)                                      tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                                  <  1:9308/71‑1 (MQ=255)                                      tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                                  <  1:3613702/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: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: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)                                                   aggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGc                     >  1:1171949/2‑71 (MQ=255)                                                    ggCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCt                    <  1:667631/71‑1 (MQ=255)                                                     gCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAAt                               <  1:1128546/59‑1 (MQ=255)                                                     gCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCtg                   >  1:3684837/1‑71 (MQ=255)                                                      cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                                  <  1:1366525/55‑1 (MQ=255)                                                      cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                  <  1:2515024/71‑1 (MQ=255)                                                      cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                  >  1:3738395/1‑71 (MQ=255)                                                      cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                  <  1:465468/71‑1 (MQ=255)                                                      cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                  <  1:2733018/71‑1 (MQ=255)                                                       aCGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGa                                         <  1:2021355/47‑1 (MQ=255)                                                       aCGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                  <  1:2510045/70‑1 (MQ=255)                                                        cGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                >  1:2848709/1‑71 (MQ=255)                                                           aTGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGc             >  1:2283777/1‑71 (MQ=255)                                                            tGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACg                                          <  1:1244630/41‑1 (MQ=255)                                                            tGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                  <  1:1522903/65‑1 (MQ=255)                                                            tGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                  <  1:1527764/65‑1 (MQ=255)                                                              ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                >  1:4071431/1‑49 (MQ=255)                                                              ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                >  1:3060851/1‑49 (MQ=255)                                                              ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                  <  1:3222207/63‑1 (MQ=255)                                                              ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                  <  1:2692561/63‑1 (MQ=255)                                                              ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                >  1:2610202/1‑65 (MQ=255)                                                              ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgtg          >  1:3096739/1‑71 (MQ=255)                                                              ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgtg          >  1:1851137/1‑71 (MQ=255)                                                               cGTAAAAGCTAAACTATTCCCAGATTTGCAGCGCAAAGACATAGTAAATAGCGCGCGGCTGTGATGCgtgt         <  1:2854933/71‑1 (MQ=255)                                                               cGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                                  <  1:4016739/46‑1 (MQ=255)                                                               cGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgtgt         <  1:3482689/71‑1 (MQ=255)                                                               cGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgtgt         >  1:849753/1‑71 (MQ=255)                                                                gTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAGATAGCGCGCGGCTGTGATGCGTGtt        >  1:3803075/1‑71 (MQ=255)                                                                gTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGtt        >  1:1078163/1‑71 (MQ=255)                                                                gTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGtt        <  1:1966428/71‑1 (MQ=255)                                                                 tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa                                       <  1:1175312/39‑1 (MQ=255)                                                                 tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                >  1:1808109/1‑62 (MQ=255)                                                                 tAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGtttt      >  1:4093893/1‑70 (MQ=255)                                                                   aaaGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAt     >  1:1366680/1‑71 (MQ=255)                                                                   aaaGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAt     >  1:3352262/1‑71 (MQ=255)                                                                     aGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                >  1:3966265/1‑42 (MQ=255)                                                                     aGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                      >  1:2812259/1‑52 (MQ=255)                                                                      gCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCa  <  1:799288/71‑1 (MQ=255)                                                                      |                                                                       GGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCA  >  minE/454944‑455083

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