Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,677,185	T→G	6.0%	R249R (AGA→CGA) ‡	yhgA ←	predicted transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,677,185	0	T	G	6.0%	229.6 / 6.8	83	R249R (AGA→CGA) ‡	yhgA	predicted transposase
Reads supporting (aligned to +/- strand):  ref base T (46/32);  new base G (2/3);  total (48/35)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

CCGGATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGT  >  minE/2677121‑2677253                                                                 |                                                                     ccGGATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCg                                                                 >  1:2227312/1‑70 (MQ=255) ccGGATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCg                                                                 >  1:2233742/1‑70 (MQ=255) ccGGATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCg                                                                 >  1:4432514/1‑70 (MQ=255)  cGGATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGc                                                                >  1:3978426/1‑70 (MQ=255)   ggATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCaa                                                              <  1:2413551/71‑1 (MQ=255)       tCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTc                                                                   >  1:477769/1‑62 (MQ=255)       tCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGt                                                          >  1:4851246/1‑71 (MQ=255)           gCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCAt                                                       >  1:1679226/1‑70 (MQ=255)            cATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTa                                                     <  1:3722993/71‑1 (MQ=255)             attattTTGGCTATATGCAGGGC‑TTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTaaa                                                   >  1:1993270/1‑71 (MQ=255)              ttattTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTaaa                                                   <  1:3674607/71‑1 (MQ=255)              ttattTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTaaa                                                   <  1:2999804/71‑1 (MQ=255)                attTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAAct                                                 >  1:1442851/1‑71 (MQ=255)                attTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAAc                                                  >  1:1614431/1‑70 (MQ=255)                   ttGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAActct                                               <  1:4030237/70‑1 (MQ=255)                      gCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCtt                                           <  1:2238824/71‑1 (MQ=255)                         atatGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCaa                                                              >  1:2286552/1‑49 (MQ=255)                         atatGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTtgtg                                        >  1:3018182/1‑71 (MQ=255)                          tatGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCtt                                                                     <  1:865956/41‑1 (MQ=255)                          tatGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTgtgt                                       >  1:2640550/1‑71 (MQ=255)                            tGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTa                                                     >  1:663721/1‑55 (MQ=255)                            tGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAActc                                                <  1:2322354/60‑1 (MQ=255)                            tGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGttt                                     >  1:458701/1‑71 (MQ=255)                            tGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGttt                                     >  1:5674497/1‑71 (MQ=255)                            tGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGttt                                     >  1:2039009/1‑71 (MQ=255)                             gCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCttt                                                                    <  1:4568646/39‑1 (MQ=255)                             gCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTc                                    >  1:2855571/1‑71 (MQ=255)                              cAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCg                                   >  1:2114348/1‑71 (MQ=255)                               aGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTcc                                                                  >  1:5476994/1‑39 (MQ=255)                               aGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCgg                                  <  1:3123443/71‑1 (MQ=255)                                gggCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTc                                                         <  1:1161130/47‑1 (MQ=255)                                gggCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGttt                                     <  1:4605188/67‑1 (MQ=255)                                 ggCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGt                                 >  1:2284974/1‑70 (MQ=255)                                   cTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGTAATAGTCATTAAACTCTCCTTGTGTTTCGGTGaa                              >  1:665471/1‑71 (MQ=255)                                   cTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTaaa                                                   >  1:3730373/1‑50 (MQ=255)                                   cTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTg                                <  1:2654664/69‑1 (MQ=255)                                   cTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGaa                              >  1:5566757/1‑71 (MQ=255)                                    tttGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTa                                                     <  1:1186435/47‑1 (MQ=255)                                     ttGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGaa                              >  1:1695325/1‑69 (MQ=255)                                      tGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGt                           >  1:2339278/1‑71 (MQ=255)                                      tGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGt                           >  1:3608937/1‑71 (MQ=255)                                      tGGATTGTTCCACCTCCTGCCGCAATATTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGt                           >  1:3721111/1‑71 (MQ=255)                                        gATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAg                                                           >  1:1901807/1‑37 (MQ=255)                                        gATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGtt                          <  1:2158191/70‑1 (MQ=255)                                        gATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTc                         >  1:1466136/1‑71 (MQ=255)                                          ttGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCa                                                        <  1:1809776/38‑1 (MQ=255)                                          ttGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCgg                       <  1:3575475/71‑1 (MQ=255)                                          ttGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCgg                       <  1:2933328/71‑1 (MQ=255)                                           tGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTaaa                                                   >  1:1461287/1‑42 (MQ=255)                                            gTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCAtt                                                      >  1:4177365/1‑38 (MQ=255)                                            gTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTa                                                     <  1:4398553/39‑1 (MQ=255)                                            gTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGc                      <  1:1473989/70‑1 (MQ=255)                                            gTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGc                      <  1:2483895/70‑1 (MQ=255)                                             ttCCCACTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCa                     <  1:1604583/70‑1 (MQ=255)                                              tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTaaa                                                   >  1:2550433/1‑39 (MQ=255)                                              tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCg                        >  1:4082537/1‑66 (MQ=255)                                              tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCg                        >  1:783954/1‑66 (MQ=255)                                              tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCaa                    >  1:5549790/1‑70 (MQ=255)                                              tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAAc                   >  1:3991758/1‑71 (MQ=255)                                              tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAAc                   >  1:2568676/1‑71 (MQ=255)                                                ccccTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACg                  <  1:1969183/70‑1 (MQ=255)                                                 cccTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTACTTGTGTTTCGGTg                                <  1:4651551/55‑1 (MQ=255)                                                   ttcctGCCGTAGCCGTTCAGCAATAGTCATTAATCTCTCCTtttg                                        <  1:4111483/44‑1 (MQ=12)                                                   ttcctGCCGTAGCCGTTCAGCAATAGTCATTAATCTCTCCTTttgta                                      >  1:3008479/2‑46 (MQ=12)                                                   ctcctGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCaa                    <  1:3778954/65‑1 (MQ=255)                                                   ctcctGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCaa                    <  1:2492474/65‑1 (MQ=255)                                                   ctcctGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCaa                    <  1:1389454/65‑1 (MQ=255)                                                   ctcctGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCaa                    <  1:4730658/65‑1 (MQ=255)                                                   ctcctGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGt              >  1:1520338/1‑71 (MQ=255)                                                    tcctGCCGTAGCCGTTCAGCAATAGTCATTAATCTCTCCTTttgt                                       <  1:1771849/45‑1 (MQ=12)                                                    tcctGCCGTAGCCGTTCAGCAATAGTCATTAATCTCTCCTTttgt                                       <  1:5020024/45‑1 (MQ=12)                                                    tcctGCCGCAGCCGTTCAGCAATAGTCATTAATCTCTCCt                                            >  1:2142348/1‑40 (MQ=16)                                                        gCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGAt          >  1:5014097/1‑70 (MQ=255)                                                        gCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGAt          >  1:3202811/1‑70 (MQ=255)                                                         ccGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTgtgt                                       <  1:1110434/40‑1 (MQ=255)                                                         ccGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATaa        <  1:5446177/71‑1 (MQ=255)                                                           gCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATaaaa      >  1:1918759/1‑71 (MQ=255)                                                             aaTCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATc    <  1:2978835/71‑1 (MQ=255)                                                             aaTCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATc    <  1:1136212/71‑1 (MQ=255)                                                              aTCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCg   >  1:4838314/1‑71 (MQ=255)                                                              aTCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCg   >  1:1451415/1‑71 (MQ=255)                                                               tCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTg                                <  1:1744585/41‑1 (MQ=255)                                                               tCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATaaaa      <  1:4637351/67‑1 (MQ=255)                                                               tCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGt  >  1:799715/1‑71 (MQ=255)                                                                cTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTc                         >  1:208756/1‑47 (MQ=255)                                                                 |                                                                     CCGGATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGT  >  minE/2677121‑2677253

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